Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,823,288 C→T G161D (GGC→GAC)  recA ← DNA recombination/repair protein RecA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,823,2880CT97.9% 342.8 / ‑6.4 97G161D (GGC→GAC) recADNA recombination/repair protein RecA
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base T (50/44);  total (50/47)
Fisher's exact test for biased strand distribution p-value = 2.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

TTCGCTACCCACCACGTTTTCGCCCTCTTTCACCGCGCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGCCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATC  >  NC_000913/2823044‑2823530
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                    
ttCGCTACCCACCACGTTTTCGCCCTCTTTCACCGCGCCGATACGACGGATGTCGAGACGAACAGAGACGTAGAATTTCAGCGCGTTACCACCGGTAGTGGGTTCCGGGTTACCGCACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGTCCCATGTGAGAGTCGTCGAttt                                                                                                                                                                                                                                                >  2:503942/1‑251 (MQ=255)
    cTACCCACCACGTTTTCGCCCTCTTTCACCGCGCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCACGGCCCATGTGAGAGTCGTCGCTTTCGcc                                                                                                                                                                                                                                            >  1:527846/1‑251 (MQ=255)
       ccccccacGTTTTCGCCCTCTGTCACCGCGCCGATACGACGGATGTCGAGACGAACAGATGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGTTTTCCGGTTTACCGAACATCACATCTATTTTCGTACTGTTCTGGTTGATTAAGATCAGCAGCGTTTTTGACTGCTTCAGGTTACCCGCCTGCTTACGCATCGCGTGGCTCATCATTCTTGCTGCAGGGCCCATGTGAGAGTCGTCGATTTCGCCTTc                                                                                                                                                                                                                                         <  1:596541/247‑1 (MQ=255)
          cccacGTTTTCGCCCTCTTTCACCGCGCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGAt                                                                                                                                                                                                                                      <  1:214582/250‑1 (MQ=255)
              gttgttCGCCCTCTTTCTCCTGGCCGAGACGGCGGTTGTCGAGGCGAACGGAGGCTTGGTATTGCGTCGCGTGACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCTTGCGGATTTGGTTTATGAAGACCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGTCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTcc                                                                                                                                                                                                                                  <  2:237589/247‑1 (MQ=255)
               gTTTTCTCCCTCTTTCACCCGGCCGCTACGACGGATGGCGATACGAACAGTGGCGGAGAATTTCAGCGCGTGACCACCGTGAGTGGTTTCCGGGTTATCGAACATCACACCAATTTTCATACGGAGCTGGTTGATGAAGATCAGCAGCGTTTTGGACGGCTGCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTGCGCTTTCGAGTTCCg                                                                                                                                                                                                                                 <  2:346634/251‑1 (MQ=255)
                ttttCGCCCTCTTTCACCGCGCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGc                                                                                                                                                                                                                                >  1:523470/1‑251 (MQ=255)
                    cGCCCTCTTTCACCGCGCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTc                                                                                                                                                                                                                            <  1:596182/251‑1 (MQ=255)
                    cGCCCTCTATCACCGCTCCGATCCGACGGATGTCGAGAGGAACAGAGGCGTAGAATTTCAGCGCGTTACCCCCGGGAGTGTTTTCCGGGTTACCGAACACCCCACCAATTTTCATACGGATCTGGTTGATGAAGTTCAGCAGCGTGTTGGACTTCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTc                                                                                                                                                                                                                            <  2:444729/251‑1 (MQ=255)
                     ccccTCTTTCACCGCGCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTTATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGGGAGAGGCGTCGATTTCGCCTTCGATTTCCGCTTTCg                                                                                                                                                                                                                           <  1:63532/250‑1 (MQ=255)
                      cccTCTTTCACCGCGCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACAGCACACCAATTTTCATACGGATCTGGTTGAGGAAGATCAGCAGCGTGGTGGACTGCTTCAGGGTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATTGGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCgg                                                                                                                                                                                                                          <  1:544368/251‑1 (MQ=255)
                          ctTTCCCCGCGCCGAGACGACGGAGTTCGAGACGAACAGATGCGGAGAATTTCAGCGCGTTGCCACCGGTAGTGGTTTCCGGGTTTGCGAACATCACACCATTTTTCATACGTATCTTTTTGATGAAGATCAGCTGCGTGTTTGGCGGCTTCAGGTTTCCCGCCAGCTTACGCATCTCCTGGCTCAGCATACGTGCCGCAAGGCCCATGTGAGAGTCCCCTAGTTCGCCTTCGAGTGCCGCTTTCGGCg                                                                                                                                                                                                                        <  2:530074/249‑1 (MQ=255)
                            ttCACCGCGCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCACGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTc                                                                                                                                                                                                                      >  1:389923/1‑249 (MQ=255)
                              cACCGCGCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCATCAGTg                                                                                                                                                                                                                  <  1:163967/251‑1 (MQ=255)
                              cACCGCGACGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTg                                                                                                                                                                                                                  >  1:128224/1‑251 (MQ=255)
                                ccGCGCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTCCCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTgcc                                                                                                                                                                                                                <  1:421770/251‑1 (MQ=255)
                                        aTACGACGGATTGCGAGACGAACAGAGGCGTAGAATTTCAGCGTTTTACCACCGTTAGTTGTTTCGGGTTTACCGAACATCACCCCAATTTTCATACGTATCGGGTTGATGAATATCAGCAGCGTGTTGGACTGCTTCAGTTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGATAGTCGTCGATTTCGCCTTCGATTTCCGCTTACGGCGTCAGTGCCGCCACGGa                                                                                                                                                                                                        <  2:391034/251‑1 (MQ=255)
                                          acgactgatctggagaCGTACAGAGACGTAGAATTTCTGCTCGTTACCACCGGTAGTGGTGTCCGGTTTACCGCACATCACCCCAATTTTCCTCCGGATGTGGTTGAGGAAGATAATCCGCGTGTTGGACTGCTTCAGGTTCCCCGCCAGCTTACGCATCGCCTGGCTCATCATCCGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGCCAGTGTCGCTACGGAGt                                                                                                                                                                                                      <  2:532475/239‑1 (MQ=255)
                                            gacgGCTGTCGAGACGAGCAGAGGGGTAGAATTTCATCGAGTTACCACCGGAAGTGGTTTCCGGGATCCAGCACATCACACCAATTTTCATCCGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGCCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTc                                                                                                                                                                                                     >  2:288102/1‑250 (MQ=255)
                                             acgGATGTCGAGACGAACAGAGGCGTAGAATATCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATAAGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAggca                                                                                                                                                                                                    >  2:513722/1‑247 (MQ=255)
                                              cgGATGTCGTGACGAACAGAGGCGTAGAATTACAGCGCGTTACCTCCGGTAGTGGTTGCCGGGTTACCGAACATCACACCAATTTTCATACGTATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTTGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGGCAGTGCCGCCACGGAGTCAac                                                                                                                                                                                                  <  2:85057/251‑1 (MQ=255)
                                                          cGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTc                                                                                                                                                                                      >  2:483562/1‑251 (MQ=255)
                                                          cGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTc                                                                                                                                                                                      >  1:33079/1‑251 (MQ=255)
                                                             aCAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGCTAACGTCTAc                                                                                                                                                                                   >  2:181031/1‑251 (MQ=255)
                                                                 agGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTgcg                                                                                                                                                                               >  1:441491/1‑251 (MQ=255)
                                                                   gCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCTTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGcc                                                                                                                                                                             <  1:53860/251‑1 (MQ=255)
                                                                     gTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGcc                                                                                                                                                                             >  1:21274/1‑249 (MQ=255)
                                                                              cAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACCCCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACGGCTGCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCa                                                                                                                                                                  <  2:33079/251‑1 (MQ=255)
                                                                               aGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGcc                                                                                                                                                                   >  1:435374/1‑249 (MQ=255)
                                                                               aGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGcc                                                                                                                                                                   >  1:576645/1‑249 (MQ=255)
                                                                                   cgtttcacgcGGTAGTGGTGACCGTGTTCCCGAACATCACACCAATTTCCATACGGAGCTGGGTGAGGAAGAGCAGCAGCGTGTGGGACTGCTTCAGGTGCCCCGCCAGCTGACGCTTCGCCTGGCTCATCATTCGTGCCGGAAGGCCCATGTGAGAGTCGTCGATTTCGCTGTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGAGAACGGCTACGGCGCCAGAACGCGCCAggg                                                                                                                                                               <  2:23880/240‑1 (MQ=255)
                                                                                   cgTTACCACAGGTAGTGGTTCCCGGGTTACCGAACAGCACACCAATTTTCATACGGATTTGGTTGAGAAAGATCAGCAGCTTGTTGGACTGCTCCAGGTTCCCCGCCAGCTTACGCAGCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATACCGTCTACTGCGCCAGAACGCGCCAGGGCg                                                                                                                                                             <  2:43194/251‑1 (MQ=255)
                                                                                     ttACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTc                                                                                                                                                           >  1:322409/1‑251 (MQ=255)
                                                                                     ttACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTc                                                                                                                                                           >  2:839/1‑251 (MQ=255)
                                                                                       cccaccGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCTCCAGAACGCGCCAGGGCGTcac                                                                                                                                                         <  1:178392/250‑1 (MQ=255)
                                                                                          accGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTcac                                                                                                                                                         >  1:388076/1‑248 (MQ=255)
                                                                                            cGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGTAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGAtt                                                                                                                                                    >  1:234669/1‑251 (MQ=255)
                                                                                             ggTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCTTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGAttt                                                                                                                                                   <  1:839/251‑1 (MQ=255)
                                                                                              gTAGTGGATTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGGTTTCGCCTTCGATTTCCTCTTTCGTCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTcaccgctttc                                                                                                                                                  >  2:208078/1‑244 (MQ=255)
                                                                                               tAGTGCTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCCGCTTACCCCTCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGCGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGCTAACGTCTACTGCCCCCGAACGCGCCAGGGCGTCACAGATTTcc                                                                                                                                                 >  2:177312/1‑251 (MQ=255)
                                                                                                    gTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCAATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTCCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGc                                                                                                                                            >  2:406376/1‑251 (MQ=255)
                                                                                                    gTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGGGTTGGACGGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGc                                                                                                                                            <  2:523470/251‑1 (MQ=255)
                                                                                                      ttCCGGGTTACCTAACTTCAAACCAATTTTCAGACGGAGCTGGTTGAGGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGGGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGcc                                                                                                                                           <  2:454634/250‑1 (MQ=255)
                                                                                                      ttCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCt                                                                                                                                          >  1:205837/1‑251 (MQ=255)
                                                                                                         cGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCTTCACAGATTTCCAGTGCCTGCt                                                                                                                                       <  1:485704/251‑1 (MQ=255)
                                                                                                         cGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCt                                                                                                                                       >  1:130084/1‑251 (MQ=255)
                                                                                                              tACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCg                                                                                                                                  >  2:80080/1‑251 (MQ=255)
                                                                                                                ccGAACATCACCCCAATTTTCAGACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGTTTACCCGCCAGCTTACGCATCGCCTGTCTCATCATACGTGCCGCAAGGCCCAGGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGgt                                                                                                                                <  2:128224/251‑1 (MQ=255)
                                                                                                                ccGAACATCACACCAATTTTCATACGGATCTGGTTGCTGCAGATCAGAAGAGTGATGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCACCTGGATCCTCATACGTGCCGCAAGGCCACTGTCAGAGTCGTCGATTTCGCCTTCGCTTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCACCAGAACGAGCCAGGCCGTCACAGATTTCCAGAGCCTGATCGCCGgt                                                                                                                                >  1:362113/1‑251 (MQ=255)
                                                                                                                            ccAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGCTTTCCAGTGCCTGCTCGCCGGTGTCCGGCTgggg                                                                                                                    >  2:422103/1‑250 (MQ=255)
                                                                                                                            ccAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTggg                                                                                                                     >  1:12829/1‑250 (MQ=255)
                                                                                                                              aaTTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGc                                                                                                                  >  1:534913/1‑251 (MQ=255)
                                                                                                                                 tttCATACGGATCTGTTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTGCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGcac                                                                                                                <  1:73979/250‑1 (MQ=255)
                                                                                                                                            tCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAg                                                                                                                                                       >  2:529346/1‑200 (MQ=255)
                                                                                                                                            tCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAg                                                                                                                                                       <  1:529346/200‑1 (MQ=255)
                                                                                                                                               ggTTGATGAAGATCAGCAGCGTGTTGGACTGCTCCATGTTCCCCGCCAGCTGACGCATCGCCTGGCTCATCATACGTGCCTCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGGCAACGACGATAACGTCTTCTGCGCCAGAACGCGCCAGGGCGTCACATATTTCCAGGGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTc                                                                                                    <  1:533087/247‑1 (MQ=255)
                                                                                                                                                gTTGATGAAGCACACCAGCGTGTTGGACTGCATCAGGTTACCCGCAAGCTAACCCTTCGCCTGGCTCATCATACGTGCCGCAAGGACCATGTGCGAGACGTCGATTTCGCCTTCGATTTCCGCTCTCGGCGTCAGTGCCGCCACGGACTCAACGACGATAACGTCTACTGCGCCAGACCCCGCCAGGCCGTCACAGATTTCCAGTGCCTGCTCGCCGGGGTCCGGCTGGGAGCACAGCAGGGTGTCGatat                                                                                               >  2:72085/1‑251 (MQ=255)
                                                                                                                                                gTTGATGAAGATCAGCAGAGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGatat                                                                                               >  1:287247/1‑251 (MQ=255)
                                                                                                                                                attGATGAAGAGCAGCATCGGGTTGGACGGCGTTGGGTTCCCCGCCAGCTTACGCATCGGCTGGCGCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTGCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGGCGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCTCAGAAGGTTGTCGatat                                                                                               <  2:495102/250‑1 (MQ=255)
                                                                                                                                                  tgatgaAGATCAGCAGCGTGTTGGACTGCTTCAGGTTCCCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATc                                                                                              <  2:130084/250‑1 (MQ=255)
                                                                                                                                                        aGATCAGCAGCGTGTTGGCCTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCAACGccc                                                                                       <  1:153976/251‑1 (MQ=255)
                                                                                                                                                        aGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGccc                                                                                       >  1:583009/1‑251 (MQ=255)
                                                                                                                                                                  cGTGTGGGACTGCTTCAGGTTCCCCGCCAGCTTACGCTTCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTg                                                                             <  1:453251/251‑1 (MQ=255)
                                                                                                                                                                   gtgtTGGACTGCCTCAGGTTTCCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTTCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGc                                                                            <  1:208078/251‑1 (MQ=255)
                                                                                                                                                                     gtTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCgg                                                                          >  2:316211/1‑250 (MQ=255)
                                                                                                                                                                     gtTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGt                                                                          <  2:143971/251‑1 (MQ=255)
                                                                                                                                                                       tGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGt                                                                          >  1:567755/1‑249 (MQ=255)
                                                                                                                                                                       tGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGt                                                                          >  2:425979/1‑249 (MQ=255)
                                                                                                                                                                       tGGACTGCTTCAGGTTAACCGCAAGCTTACGCATCGCCTGGCTCATCATAAGTGCCGCCAGGCCCATGAGCGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCACCGACGATAAAGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACCGCAGGTTGTCGATATCGCCGCCAAGTTTACTTGCGtcg                                                                        >  2:559419/1‑249 (MQ=255)
                                                                                                                                                                        ggACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAg                                                                        <  1:483562/250‑1 (MQ=255)
                                                                                                                                                                             gCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCTTCGATTTCGCCTTCGATTTCCGCTTCCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTg                                                                    <  2:151701/249‑1 (MQ=255)
                                                                                                                                                                             gCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTggg                                                                  >  1:544604/1‑251 (MQ=255)
                                                                                                                                                                                  aGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGCTTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTcc                                                               >  2:109288/1‑249 (MQ=255)
                                                                                                                                                                                  aGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAg                                                             <  1:109288/251‑1 (MQ=255)
                                                                                                                                                                                    gTTACCCGCCAGCTTACGCTTCGCCTGGCTCATCATACGTGCCGCAAGACCCATGTGCGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAgcg                                                           >  2:408283/1‑251 (MQ=255)
                                                                                                                                                                                     ttACCCGCCAGCGTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAgcg                                                           <  2:388076/250‑1 (MQ=255)
                                                                                                                                                                                          cGCCAGCTTACGCCGCGCCTTGCTCATCATACGTGCCGCCAGGCCCCGTTTAGAGTCGCCGATGTCGCCTTCGATTTCCTCTTTCGGCGTCAGTGCCCCCACGGAGTCAACGACGATAACTTCTACTGCGCCAGAACGCGCCAGGGGCTCACAGAGTTCCAGCTCCTTCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGGTGTCGATTTCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGtt                                                     <  2:152855/251‑1 (MQ=255)
                                                                                                                                                                                           gCCAGCCTACGCATCGCCTGGCTCATAATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCATTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAAAGTCTACTGCGCCCGAACGCGCCAGGGCGTCACAGATTTCCCGTGCCTGCTCGCCGGTGGCCGGCTGGGAGCCCAGCCGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTc                                                    >  2:352036/1‑251 (MQ=255)
                                                                                                                                                                                            ccAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGTCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACATATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCTATAGCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCa                                                   <  1:513722/251‑1 (MQ=255)
                                                                                                                                                                                                    cGCATCGCCTGGCGCATCATACGTGCCGCAAGGCTCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGGCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGAt                                           <  2:556718/251‑1 (MQ=255)
                                                                                                                                                                                                      cATCGCCTGGCTCATCATACGTGCCGCAAGGCCCCTGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCCACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCCCAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGCTTGGGTCCAGCGCGTGTTCAGCATCGATa                                          >  1:54109/1‑250 (MQ=255)
                                                                                                                                                                                                              ggCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAg                                 <  1:519244/251‑1 (MQ=255)
                                                                                                                                                                                                                 tcatcaTACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCAACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGCCGCCCAGTTTACGTGCGTAGATTGGGTCCAGCCCGTGTTCAGCATCGATAACCGCACAGGtt                             >  2:400145/1‑251 (MQ=255)
                                                                                                                                                                                                                  catcatACGTGCCGCAAGGCCCATGTTAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCTTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTGCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTGACGTGCGTAGATTGGGGCCTGCGCGTGGTCATCATCGATGAACGCACAGGtt                             <  2:336024/250‑1 (MQ=255)
                                                                                                                                                                                                                  catcatACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGttt                            >  1:577976/1‑251 (MQ=255)
                                                                                                                                                                                                                     catACGTGCCGAAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGGCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGCTAAACGCACCGGTTTTAc                         >  2:330258/1‑251 (MQ=255)
                                                                                                                                                                                                                         cGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGACAGAACGCGCCAGGGCGACACAGATTTCCAGTGCCTGCTCGCCGGGGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCAAGCGCGTGATCAGCAGCGATAAACGCACAgggtatacatt                       >  2:334420/1‑241 (MQ=255)
                                                                                                                                                                                                                         cGTGCCGCAAGGCCCATGTGAGAGTCGGCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTc                     <  2:534913/251‑1 (MQ=255)
                                                                                                                                                                                                                          gTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCa                    <  1:587759/251‑1 (MQ=255)
                                                                                                                                                                                                                              cGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCAcgct                <  1:38516/251‑1 (MQ=255)
                                                                                                                                                                                                                              cGCAAGGCCCATGGGAGAGTCGTCGAGTTCTCCGTCGATTTCCGCTGTCGGGGGCAGTGCCGCCACGGAGTCAACGACGATAACGTCGACTGCGCCAGAACGCGCCAGGTCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGTCTGGGAGCGCAGCAGGATTGCGATATTGACACCCAGTTTACGTGTGTAGATTTGGGCCTGCGCTGGTTCAGCATCGATAAACGGACATTTTTTACCGTCAcgct                <  2:592355/251‑1 (MQ=255)
                                                                                                                                                                                                                                 aaGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCAcgctgcg             >  1:146790/1‑251 (MQ=255)
                                                                                                                                                                                                                                       cATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTgcgg       <  1:304915/251‑1 (MQ=255)
                                                                                                                                                                                                                                        aTGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTgcggc      >  1:192724/1‑251 (MQ=255)
                                                                                                                                                                                                                                          gtgAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGGTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTgaggcg     >  1:515270/1‑250 (MQ=255)
                                                                                                                                                                                                                                           tgAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGACAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGa  >  2:280213/1‑250 (MQ=255)
                                                                                                                                                                                                                                            gagggtcgtcGATTTCGCCTTCGATTTCCGCTGTAGTCGTCAGTACCGCCACGACGTCAACGACGATAACGTCTACTGAGCCAGAACGCGCCAGGGCGTCACAGATTTAAAGTGCCTGCTCGCAAGAGTCCGGCTGCGAGCAAAGCAGGTTGTAGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCCAGGGGTTTACCTTCACACTGCGCTgcgg       >  1:257883/5‑246 (MQ=255)
                                                                                                                                                                                                                                             agagTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGGGCCGCAACGGAGTCCACGACGATAACGTCTACTGCGACAGAACGCGCCAGGGCGACACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGACCACAGCAGATTGACGATATCGACGCCCAGTTTACGTGCATAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCGGCGCTGCGGCGATc  >  1:297043/1‑250 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                    
TTCGCTACCCACCACGTTTTCGCCCTCTTTCACCGCGCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGCCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATC  >  NC_000913/2823044‑2823530

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: