Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,948,180 T→A 11.2% *151C (TGA→TGT nudB ← dihydroneopterin triphosphate pyrophosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,948,1800TA11.2% 187.3 / 13.1 71*151C (TGA→TGTnudBdihydroneopterin triphosphate pyrophosphatase
Reads supporting (aligned to +/- strand):  ref base T (30/32);  new base A (4/4);  total (34/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.93e-01

CCGTAACCTTCGTAGATGATGGTTTCCATGTTTGCATCATCATCACCGCCCACACCACGTGCAATTGCGCGGTTCAGTGTGTCACGGGTCATGTTGTTAGACAGTGCTTTATCCACCGCCGCACGCAGACGCGGGTTAGCGTCCGGATCGCCACCGCCCAGCTTAGCCGCGGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTCAATTAAGGTCAGTTGTTCAGCGACAACATCAATGGTGACCTCTTCCTTTACTTCGCGCA  >  NC_000913/1947880‑1948472
                                                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                                    
ccGTAACCTTCGTAGATGATGGTTTCCATGTTTGCATCATCATCACCGCCCACACCACGTGCAATTGCGCGGTTCAGTGTGTCACGGGTCATGTTGTTAGACAGTGCTTTATCCACCGCCGCACGCAGACGCGGGTTAGCGTCCGGATCGCCACCGCCCAGCTTAGCCGCGGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATACCTCCAGATAGCCCTGCCCGtt                                                                                                                                                                                                                                                                                                      >  1:149693/1‑301 (MQ=255)
 ccgggttgtctttgttgatgGTTTCCATGTTTGCATCATCATCACCGCCCACGCCACGTGCAGTTGCGCGGTTCAGTGTGTCACGGGTCATGTTGTTAGACAGTGCTTTATCCACCGCCGCACGCAGACGCGGGTTAGCGTCCGGATCGCCACCGCCCAGCTTAGCCGCGGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTc                                                                                                                                                                                                                                                                                                     <  1:102908/286‑1 (MQ=255)
 cGTAACCTTCGTAGATGATGGTTTCCATGTTTGCATCATCATCACCGCCCACACCCCGTGCAATTGCGCGGTTCAGTGTGTCACGGGTCATGTTGTTAGACAGTGCTTTATCCACCGCCGCACGCAGACGCGGGTTAGCGTCCGGATCGCCACCGCCCAGCTTAGCCGCGGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTc                                                                                                                                                                                                                                                                                                     <  1:133199/301‑1 (MQ=255)
       tttCTTAGATGTTGGTTTCCATGTTTGCATCATCATCACCGCCCACACCACGTGCAATTGCGCGGTTCAGTGTGTCACGGGTCATGTTGTTAGACAGTGCTTTATCCACCGCCGCACGCAGACGCGGGTTAGCGTCCGGATCGCCACCGCCCAGCTTAGCCGCGGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAg                                                                                                                                                                                                                                                                                               <  2:268007/300‑1 (MQ=255)
        ttCGTAGATGATGGTTTCCATGTTTGCATCATCATCACCGCCCACACCACGTGCAATTGCGCGGTTCAGTGTGTCACGGGTCATGTTGTTAGACAGTGCTTTATCCACCGCCGCACGCAGACGCGGGTTAGCGTCCGGATCGCCACCGCCCAGCTTAGCCGCGGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTAAgaacgg                                                                                                                                                                                                                                                                                              >  2:103032/1‑296 (MQ=255)
        ttCGTAGATGATGGTTTCCATGTTTGCATCATCATCACCGCCCACACCACGTGCAATTGCGCGGTTCAGTGTGTCACGGGTCATGTTGTTAGACAGTGCTTTATCCACCGCCGCACGCAGACGCGGGTTAGCGTCCGGATCGCCACCGCCCAGCTTAGCCGCGGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGAACAGGCAgg                                                                                                                                                                                                                                                                                              <  2:488566/301‑2 (MQ=255)
         tCGTAGATGTTGGTTTCCATGTTTGCATCATCATCACCGCCCACACCACGTGCAATGGCGCGGTTCAGTGTGTCACGGGTCATGTTGTTAGACAGTGCTTTATCCACCGCCGCACGCAGACGCGGGTTAGCGTCCGGATCGCCACCGCCCAGCTTAGCCGCGGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGAACAGGCAggg                                                                                                                                                                                                                                                                                             <  2:7738/301‑3 (MQ=255)
         tCGTAGATGATGGTTTCCATGTTTGCATCATCATCACCGCCCACACCACGTGCAATTGCGCGGTTCAGTGTGTCACGGGTCATGTTGTTAGACAGTGCTTTATCCACCGCCGCACGCAGACGCGGGTTAGCGTCCGGATCGCCACCGCCCAGCTTAGCCGCGGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCg                                                                                                                                                                                                                                                                                             <  2:461854/301‑1 (MQ=255)
             gggtgattGTGTCCATTTTTGCATCATCATCTGCGGCCAACGAACGGGCAAGTGCTCGGTTTAGTGGGTCACGGGTTATGTTTTTAGATAGGGCTGTATCCACAGCCGAACGCAGACGAGGGGTAGCGGCCGGATCGCCACGGTGCAGCTTAGCCGCGGTTACCAGCTCACGAATGATGTTAGTGAAGAGTTTACCGCGCTTAGCATCCTGCGCAGCGTTACGATGTCTGGTGTTGGCTCATTTACTATGACCTGACATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTaa                                                                                                                                                                                                                                                                                         <  1:497832/298‑1 (MQ=255)
              gttgttgGTTTGCATTTTTGCATCCTCATCCCCGCCCACACCCAGCGCAATTGCGCGGTTCAGTGGGTCACGGGTCATGTTGTTAGACAGTGCTTTATCCACCGCCGCACGCAGACGCGGGTTAGCGTCCGGATCGCCACCGCCCAGCTTAGCCGCGGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAAt                                                                                                                                                                                                                                                                                        <  1:437626/296‑1 (MQ=255)
               ttgatgGTTTCCATGTTTGCATCATCATCACCGCCCACACCACGTGCAATTGCGCGGTTCAGTGTGTCACGGGTCATGTTGTTAGACAGTGCTTGATCCACCGCCGCACGCAGACGCGGGTTAGCGTCCGGATCGCCACCGCCCAGCTTAGCCGCGGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGAACAGGCAgggctatct                                                                                                                                                                                                                                                                                       <  2:413885/300‑9 (MQ=255)
                          cATTTTTGCTTCATTATCACCGCCCACACCACGTGCAATTGCGCGGTTCAGTGTGTCACGTGTCATGTTGTTAGACAGTGCTTTATCCACCGCCGCACGCAGACGCGGGTTAGCGTCCGGATCGCCACCGCCCAGCTTAGCCGCGGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTc                                                                                                                                                                                                                                                                            <  2:569106/301‑1 (MQ=255)
                               ttGCATCATCATCACCGCCCACACCACGTGCAATTGCGCGGTTCAGTGTGTCACGGGTCATGTTGTTAGACAGTGCTTTATCCACCGCCGCACGCAGACGCGGGTTAGCGTCCGGATCGCCACCGCCCAGCTTAGCCGCGGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATcc                                                                                                                                                                                                                                                                       >  2:409386/1‑300 (MQ=255)
                                tGCATCATCATCACCGCCCACACCACGTGCAATTGCGCGGTTCAGTGTGTCACGGGTCATGTTGTTAGACAGTGCTTTATCCACCGCCGCACGCAGACGCGGGTTAGCGTCCGGATCGCCACCGCCCAGCTTAGCCGCGGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGc                                                                                                                                                                                                                                                                      >  2:280437/1‑301 (MQ=255)
                                  catcatcatCACCGCCCACACCACGGGCAATTGCGCGGTTCAGTGTGTCACGGGTCATGTTGTTAGTCAGTGCTTTATCCGCCGCCGCACGCAGACGCGGGTTAGCGTCCGGATCGCCACCGCCCAGCTTAGCCGCGGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCt                                                                                                                                                                                                                                                                    <  1:594331/301‑1 (MQ=255)
                                       tcatcaCCGCCCACACCACGTGCAATTGCGCGGTTCAGTGTGTCACGGGTCATGTTGTTAGACAGTGCTTTATCCACCGCCGCACGCAGACGCGGGTTAGCGTCCGGATCGCCACCGCCCAGCTTAGCCGCGGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCCTAAAAATATCTCAAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCgg                                                                                                                                                                                                                                                               >  1:284462/1‑301 (MQ=255)
                                           caCCGCCCACACCACGTGCAATTGCGCGGTTCAGTGTGTCACGGGTCATGTTGTTAGACAGTGCTTTATCCACCGCCGCACGCAGACGCGGGTTAGCGTCCGGATCGCCACCGCCCAGCTTAGCCGCGGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGc                                                                                                                                                                                                                                                           <  2:210827/301‑1 (MQ=255)
                                           caCCGCCCACACCACGTGCAATTGCGCGGTTCAGTGTGTCACGGGTCATGTTGTTAGACAGTGCTTTATCCACCGCCGCACGCAGACGCGGGTTAGCGTCCGGATCGCCACCGCCCAGCTTAGCCGCGGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGc                                                                                                                                                                                                                                                           >  2:222045/1‑301 (MQ=255)
                                                 ccACACCACGTGCAATTGCGCGGTTCAGTGTGTCACGGGTCATGTTGTTAGACAGTGCTTTATCCACCGCCGCACGCAGACGCGGGTTAGCGTCCGGATCGCCACCGCCCAGCTTAGCCGCGGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAgg                                                                                                                                                                                                                                                     <  2:195612/301‑1 (MQ=255)
                                                  cacaCCACGTGCAATTGCGCGGTTCAGTGTGTCACGGGTCATGTTGTTAGACAGTGCTTTATCCACCGCCGCACGCAGACGCGGGTTAGCGTCCGGATCGCCACCGCCCAGCTTAGCCGCGGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACaaagaacagaagagcaaacgtctgaactccaga                                                                                                                                                                                                                                                    >  2:337354/1‑271 (MQ=255)
                                                                     cgGTTCAGTGTGTCACGGGTCATGTTGTTAGACAGTGCTTTATCCACCGCCGCACGCAGACGCGGGTTAGCGTCCGGATCGCCACCGCCCAGCTTAGCCGCGGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCTATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTCCAAACTGTTCAATCGCCTACCGGTTGCCCACGGCCTTACTTACcgccgccgcc                                                                                                                                                                                                                                 >  1:51436/1‑300 (MQ=255)
                                                                            gtgtgtCACGGGTCATGTTGTTAGACAGTGCTTTATCCACCGCCGCACGCAGACGCGGGTTAGCGTCCGGATCGCCACCGCCCAGCTTAGCCGCGGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGgcg                                                                                                                                                                                                                          >  2:71691/1‑301 (MQ=255)
                                                                                      ggTCATGTTGTTAGACAGTGCTTTATCCACCGCCGCACGCAGACGCGGGTTAGCGTCCGGATCGCCACCGCCCAGCTTAGCCGCGGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCCGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAACCTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCCCAGCCGGCCCATCAAGccc                                                                                                                                                                                                                >  1:436965/1‑300 (MQ=255)
                                                                                           tgttgtTAGACAGTGCTTTATCCACCGCCGCACGCAGACGCGGGTTAGCGTCCGGATCGCCACCGCCCAGCTTAGCCGCGGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGt                                                                                                                                                                                                           >  2:259763/1‑301 (MQ=255)
                                                                                                                 cACCGCCGCACGCAGACGCGGGTTAGCGTCCGTATCGCCACCGCCCAGCTTAGCCGCGGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAAc                                                                                                                                                                                     <  2:51504/301‑1 (MQ=255)
                                                                                                                                 ggcgGGTTGGCGTCCGGATCGCCACTGCCTATCTTATCCTCGGTTTCCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCGGCGCTGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTGGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGAGAAAACGATCTGCCGCTCGTGc                                                                                                                                                                     <  1:48996/300‑1 (MQ=255)
                                                                                                                                 cgcgGGTTAGCGTCCGGATCGCCACCGCCCAGCTTAGCCGCGGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGc                                                                                                                                                                     >  2:620784/1‑301 (MQ=255)
                                                                                                                                  gctggTTAGCGTCCGGATCGCCACCGCCCAGCTTAGCCGCGGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACCGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCg                                                                                                                                                                    <  2:249511/298‑1 (MQ=255)
                                                                                                                                           cGTCCGGATCGCCACCGCCCAGCTTAGCCGCGGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAACCGCaa                                                                                                                                                           >  1:171685/1‑301 (MQ=255)
                                                                                                                                                ggATCGCCACCGCCCAGCTTAGCCGCGGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCCGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAAAGAACTCACGCCCGTGCGGCAGCGCacaccca                                                                                                                                                      >  1:223087/1‑295 (MQ=255)
                                                                                                                                                      ccACCGCCCAGCTTAGCCGCGGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCaa                                                                                                                                                >  2:176264/1‑300 (MQ=255)
                                                                                                                                                       caacgccgtggttttgcgcgGTTACCATGTGACGAATGATGTTAGTGAAGATTGTACCGCGCTTAGCATCCTGCGGAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTGCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGGCGGTTGCTCCAGGACTGAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACGTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCATg                                                                                                                                               <  1:365970/285‑1 (MQ=255)
                                                                                                                                                       cACCGCCCAGCTTAGCCGCGGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCAGATCggaagagcacacgtctgaactccagt                                                                                                                                               >  2:361916/1‑276 (MQ=255)
                                                                                                                                                          cGCCCAGCTTAGCCGCGGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCACTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACAATCTGCCGCACGTGCGGAAACGCAAGACAGAACCATGAtt                                                                                                                                            >  2:422358/1‑301 (MQ=255)
                                                                                                                                                                           tttACCAGGTCTCGAATTATTTTAGTGAAGTTTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACg                                                                                                                           <  1:221802/300‑1 (MQ=255)
                                                                                                                                                                            ttACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCATGTTCAGGAAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTCATGAGCGCCGCCGCAGCAGGCGCATACAGCCACTTGTAAGCCAGATGTTCATTAAAAACGCACTGCACCCCGTTCAGAACCACAACAAAGAACCATCAATCCGTATTCCGCATccccc                                                                                                                          >  1:588246/1‑297 (MQ=255)
                                                                                                                                                                                 agttCATGATTGGTTTTGGTGAATATTTTCCCGCGCTTGGCATCCTGGGCAGATTTACGATGTCTGGTGTTGGCCCATTTACTATGACCGGACATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCACTCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCGgc                                                                                                                     <  1:421795/298‑1 (MQ=255)
                                                                                                                                                                                        gAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCCGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCATGAATCCGTATTACGCGTCACCCCCGGCCCACAGc                                                                                                              >  1:593335/1‑301 (MQ=255)
                                                                                                                                                                                           tGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGaa                                                                                                           >  2:548957/1‑300 (MQ=255)
                                                                                                                                                                                                          ttttACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCCAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTCCGGAAGCGCAAGACAGAACCATGATTCCGTATTACCCCTCACCCCCGgccacaagcgaccacccaaacctgaaa                                                                                            >  1:632463/1‑276 (MQ=255)
                                                                                                                                                                                                            ttACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTTCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCCCATAGCGATAACGTAAATGTGaaaaa                                                                                          >  1:504801/1‑301 (MQ=255)
                                                                                                                                                                                                              aCCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCACAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGaaaaaaa                                                                                        >  2:423726/1‑300 (MQ=255)
                                                                                                                                                                                                                 gcgcTTAGCATCCTGCGCAGCTTTACGATGTTTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTc                                                                                     <  2:149844/301‑1 (MQ=255)
                                                                                                                                                                                                                   tcTTAGCATCCTGCGCAGCTTTACTATGTCTGGTGTTGTCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACCAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCaa                                                                                   <  2:437529/300‑1 (MQ=255)
                                                                                                                                                                                                                     ttAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAAc                                                                                 <  2:223330/301‑1 (MQ=255)
                                                                                                                                                                                                                           tCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCATGATTCCATATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTAcc                                                                           >  2:244345/1‑301 (MQ=255)
                                                                                                                                                                                                                                      tttGCGATGTTGGGTGTTTGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACGGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACa                                                                <  1:176067/301‑1 (MQ=255)
                                                                                                                                                                                                                                         aCGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATGACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCACAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTc                                                             <  1:423163/301‑1 (MQ=255)
                                                                                                                                                                                                                                          cGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTCa                                                            <  2:594096/301‑1 (MQ=255)
                                                                                                                                                                                                                                                cTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGAAAGTCAATTaaa                                                      >  1:77862/1‑300 (MQ=255)
                                                                                                                                                                                                                                                 tGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTCAATTAAgg                                                     <  2:505473/301‑1 (MQ=255)
                                                                                                                                                                                                                                                 ggtgtttGGCGCGTTTGTTATGTCCTGCCCTAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCAGCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTCAATTAAgg                                                     <  1:408823/296‑1 (MQ=255)
                                                                                                                                                                                                                                                   ttgtTTGCCCATTTACTATTACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTCAATTAAGGTc                                                   <  2:284765/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                 tcTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTCAATTAAGGTCAGTTGTTCAGCGac                                     <  2:254983/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                       cccTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTCAATTAAGGTCAGTTGTTCAGCGACAACATc                               <  2:77952/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                        ggagttcagacatttgctcttccgatctctgcctgTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTCAATTAAGGTCAGTTGTTCAGCGACAACATCa                              <  2:65942/273‑1 (MQ=255)
                                                                                                                                                                                                                                                                         cccgatctaaaagttttcCAGTTAGCCCTGCCTGTTCAGGCAGCGTTAATTACAACCTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTCAATTAAGGTCAGTTGTTCAGCGACAACATCaa                             <  1:341552/285‑1 (MQ=255)
                                                                                                                                                                                                                                                                               tAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCCTAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCCTACCCTGACAGTCAATTAAGATCAGTTGTTCAGCGCCAAAAACACAGGTGa                       >  1:564848/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                               ataaaaTGTGTCCTGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCAGCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACGCTACCGTGCGCTGACAGTCAATTAAGGTCAGTTGTTCAGCGACAACATCAATGGTGa                       <  1:71601/299‑1 (MQ=255)
                                                                                                                                                                                                                                                                               aaaaaaTATTTCCGTGTAGCCCTGCCTGTTCAGGCAGCGTTAATTACAATCTGTTCAAGCGCCTGCCGGGTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTCAATTAAGGTCAGTTGTTCAGCGACAACATCAATGGTGa                       <  1:619926/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                 aaaaTATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTCAATTAAGGTCAGTTGTTCAGCGACAACATCaaaaatcgga                     >  2:342015/1‑293 (MQ=255)
                                                                                                                                                                                                                                                                                  aacagtttgtaattaacgctgcctgAACAGGCAGGGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTCAATTAAGGTCAGTTGTTCAGCGACAACATCAATGGTAAccc                    >  1:230302/19‑300 (MQ=255)
                                                                                                                                                                                                                                                                                      ttctcCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCCGACAGTCAATTAAGGTCAGTTGTTCAGCGACAACATCAATGGTGACCTCTTc                <  2:171882/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                        ctcccgttggCCTTGCCTGTTCAGGCAGCGTTAATTACAAACGGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGGGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTCAATTAAGGTCAGTTGTTCAGCGACAACATCAATGGTGACCTCTTCCt              <  1:548284/292‑1 (MQ=255)
                                                                                                                                                                                                                                                                                        ctcCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAACGATGACGTAAATGAGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTCAATTAAGGTCAGTTGTTCAGCGACAACAACAAAGGAGAACTCTTACt              >  2:403318/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                            aattaacgctgcctgAACAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTCAATTAAGGTCAGTTGTTCAGCGACAACATCAATGGTGACATCTTCCTTTaa          >  2:161259/9‑300 (MQ=255)
                                                                                                                                                                                                                                                                                            aGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTCAATTAAGGACAGTTGTTCAGCGCAAACACCAATGGTAACCTCTTCCTTTat          >  1:642867/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                             gATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTCAATTAAGGTCAGTTGTTCAGCGACAACATCAATGGTGACCTCTTCCTTTACt         >  2:381853/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                             attaacgctgcctgAACAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTCAATTAAGGTCAGTTGTTCAGCAACAACATCAATGGTGACCTCttactttact         >  2:374671/8‑301 (MQ=255)
                                                                                                                                                                                                                                                                                                aacgctgcctgAACAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTCACTTAAGGTAAGTTGTTCAGCGACAACATCAATGGAGACCCCATCCTTAACTTcc      >  1:562497/5‑300 (MQ=255)
                                                                                                                                                                                                                                                                                                  cgctgcctgAACAGGTAGCGTTAATTACAAACGGTTCAGTCGCCTGCCGGTTGCTCCAGGACTTAGTGTGCGCCGCCGAAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCAGGATGCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTCAATTAAGGTCAGTTGTTCAGCGACAACATCAATGGTGACCTCTTCCTTTACTTcgcg    <  1:374175/299‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                    ctgcctgTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTCAATTAAGGTCAGTTGTTCAACGAcaacaacaagaacggaaaagcgacgcataaggaaaa  >  1:65852/1‑274 (MQ=255)
                                                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                                    
CCGTAACCTTCGTAGATGATGGTTTCCATGTTTGCATCATCATCACCGCCCACACCACGTGCAATTGCGCGGTTCAGTGTGTCACGGGTCATGTTGTTAGACAGTGCTTTATCCACCGCCGCACGCAGACGCGGGTTAGCGTCCGGATCGCCACCGCCCAGCTTAGCCGCGGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTCAATTAAGGTCAGTTGTTCAGCGACAACATCAATGGTGACCTCTTCCTTTACTTCGCGCA  >  NC_000913/1947880‑1948472

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: