Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,285,261 Δ1 bp 5.2% coding (123/711 nt) narJ → molybdenum‑cofactor‑assembly chaperone delta subunit of nitrate reductase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,285,2610T. 5.2% 420.1 / 15.5 116coding (123/711 nt)narJmolybdenum‑cofactor‑assembly chaperone delta subunit of nitrate reductase 1
Reads supporting (aligned to +/- strand):  ref base T (53/56);  new base . (3/3);  total (56/60)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.79e-02

Only 100 of 116 total aligned reads displayed.

GTGAACTGGCACGGGAAGCCTTCCCGGAGAAAAATGGCTGCGGCTTTACCTTTGGTGATGGCTGCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTAAACAGCCGCGAATTGCCGGACCATCTGCCGCTGTATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAA  >  NC_000913/1284991‑1285522
                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                      
gTGAACTGGCACGGGAAGCCTTCCCGGAGAAAAATGGCTGCGGCTTTACCTTTGGTGATGGCTGCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCGCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATgcg                                                                                                                                                                                                                                                                     <  2:343657/274‑1 (MQ=255)
gTGAACTGGAACGGGAAGCCTTCCCGGAGAAAAATGGCTGCGGCTTTACCTTTGGTGATGGCTGCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATgcg                                                                                                                                                                                                                                                                     <  2:424541/274‑1 (MQ=255)
  aagatGGCACGGGAATCCTTCCCGGAGAAAAATTGAGGTGGCTTGACCTTTGGTGAGGGGTGCCAAGGTTCAGATACCAAATTCAAGCTGTTCAACGGCGGTCGTATCGATGCCATCGCTGTGACCAGCAAAACGGAGCCGCATACATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGAGGCTGCCTTATGGCAGCAGCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCt                                                                                                                                                                                                                                                                  <  4:156426/270‑1 (MQ=255)
                agacgtccctgcaaaaaaTGGCTGCGGCTTTACCTTTGGTGATGGCTGACACGGTTCAGATACCAAATTCATGCTGTTCAACATCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATTATCGAACTCGGGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTgcg                                                                                                                                                                                                                                                   <  1:105981/262‑1 (MQ=255)
                      cccGGAGAAAAATGGCAGCGGCTTTACCTTTGGTGATGGCTGCCACGGTTCAGCTACCAAATTCAATCTTTTCAACAGCCGTCGTAACGATGCCATCGATGTGACCAGCAAAACGGAGCCGAATCCATGATCGAACTCGTGCTTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTCTGGCAGCATCACCAACTGATGTTGGAGGCGATTGCCACGTACACATATCTGCCCAAAGAGGCTGCCCCTGCTCTGGGCATTTTCCTgacccatt                                                                                                                                                                                                                                              >  4:308165/1‑267 (MQ=255)
                        cGGAGAAAAATGGCTGCGGCTTTACCTTTGGTGATGGCTGCAACGGTTCAGATAACAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTaa                                                                                                                                                                                                                                           <  2:7739/275‑1 (MQ=255)
                         ggAGAAAAATGGCTGCGGCTTTACCTTTGGTGATGGATGCCACGGTTCAGATACCAAATTAAAACTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACAAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGACTTATGGCAGCATCAACAAGAGATGTTTAAGGTGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAac                                                                                                                                                                                                                                          <  2:95528/275‑1 (MQ=255)
                             aaaaaTGGCTGCGGCTTTACCTGCGGTGATGGCTGCGACGGGTCAGATACCAAATTCAATCTGTTCAACATCCGTCGTAGCGATGCCATAGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCTAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAacgacg                                                                                                                                                                                                                                      <  1:77654/275‑1 (MQ=255)
                                    gttGCGGCTTTACCTTTGGGGCTGCGTGCTCCGGTTCAGCTACTAAATTCAATCAGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGACCGAACTCGTGATTGTATCGCGTGTCCTGGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGGGATGTGTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTGCCTGCGCGATTTAACGACGATGGATc                                                                                                                                                                                                                               <  4:506909/273‑1 (MQ=255)
                                                tttttgggtgttgCTGGCACTGTGCAGATNCCAAGTCTCATCTGGTTCACAGCCGNCGTATCGATGCCATCGATTGGACCAGCAAAACGGGGCCGCATCCATGATCGAACTCGTGATTGTATAGCGTCTCCTTGAANATCCGGAGGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGccc                                                                                                                                                                                                                   <  3:97059/263‑1 (MQ=255)
                                                ccTTTGTGGATGGCTGCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGccc                                                                                                                                                                                                                   <  2:248167/275‑1 (MQ=255)
                                                acTTTGGTGAGGGCTGCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGccc                                                                                                                                                                                                                   <  2:165399/274‑1 (MQ=255)
                                                     ggTGATGGCTGCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGgcg                                                                                                                                                                                                              <  2:151924/275‑1 (MQ=255)
                                                      gTGATGGCTGCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGgcg                                                                                                                                                                                                              <  4:124674/274‑1 (MQ=255)
                                                          tGGCTGCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCCCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAACATGCCAAAAGAGGATGCCCATGCGCTGGGCCTTTCCCTCCGCGATTTAACGCCAATGGCACCGCTACATACACAGGCGCCGtc                                                                                                                                                                                                         >  4:91154/1‑274 (MQ=255)
                                                               gCCACGTTTCATATACCAAAGTCAACCTGTTCAACAGCCGTCGTATCGATGCCATCGATGGGACCAGAAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCg                                                                                                                                                                                                    <  4:259946/275‑1 (MQ=255)
                                                               gCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCCCGATGCCCAGGAGCAGTACAgg                                                                                                                                                                                                     >  4:236922/1‑273 (MQ=255)
                                                                ccACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGa                                                                                                                                                                                                   >  2:212400/1‑275 (MQ=255)
                                                                acACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGa                                                                                                                                                                                                   <  1:99431/274‑1 (MQ=255)
                                                                   cGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCg                                                                                                                                                                                                    <  2:236923/271‑1 (MQ=255)
                                                                                            ttCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCa                                                                                                                                                                       >  2:496204/1‑275 (MQ=255)
                                                                                                    ccGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGGGGCCCCGCCACGTCACGCTTGCt                                                                                                                                                               >  4:337857/1‑275 (MQ=255)
                                                                                                                 gCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGc                                                                                                                                                                                   >  4:161592/1‑242 (MQ=255)
                                                                                                                     tCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCcgcg                                                                                                                                                                               <  2:161593/242‑1 (MQ=255)
                                                                                                                                aGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCt                                                                                                                                                                                                                                                                  >  6:493365/1‑148 (MQ=255)
                                                                                                                                 gCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTg                                                                                                                                                                                                                                                                 <  5:493365/148‑1 (MQ=255)
                                                                                                                                         gAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCAtttt                                                                                                                                                                                                                                                         <  7:1036419/148‑1 (MQ=255)
                                                                                                                                          aGCCGCATCCATGATCGAACTCGTGATTGTATCGAGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGCGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGc                                                                                                                          <  2:318623/274‑1 (MQ=255)
                                                                                                                                             agCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGGCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAgg                                                                                                                      <  1:117838/274‑1 (MQ=255)
                                                                                                                                                            aCTTGTGTTTTTATCGCGTCTCCTTGAATACCCGGTGGCTGCCTTATGGCAGCATCAACAAGGGAAGTTGTAGGCGATTGCCGCGTCGAGAAATCTGCCAAAATAGGATGCCCATGCGATGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAGCTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGAGctgc                                                                                                       <  4:65579/275‑1 (MQ=255)
                                                                                                                                                                        aTCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACg                                                                                           <  2:337858/275‑1 (MQ=255)
                                                                                                                                                                            cgTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCCATTGCCGCGACCGCGGCCAGGCGTTGGTGGACCTGAGGGCGCCGTCAGagcg                                                                                       >  2:88658/1‑274 (MQ=255)
                                                                                                                                                                              tctcCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCa                                                                                                                                                                                                                  >  4:282594/1‑150 (MQ=255)
                                                                                                                                                                                  cTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGgcg                                                                                                                                                                                                              <  2:282595/150‑1 (MQ=255)
                                                                                                                                                                                    tGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGGGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCt                                                                               <  4:212399/275‑1 (MQ=255)
                                                                                                                                                                                    tGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCAtt                                                                                                                                                                                                                                                           <  5:425247/103‑1 (MQ=255)
                                                                                                                                                                                    tGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCAtt                                                                                                                                                                                                                                                           >  6:425247/1‑103 (MQ=255)
                                                                                                                                                                                       atatCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGa                                                                                                                                                                                                                                                <  8:761129/111‑1 (MQ=255)
                                                                                                                                                                                       atatCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGa                                                                                                                                                                                                                                                >  7:761129/1‑111 (MQ=255)
                                                                                                                                                                                        tatCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCAGCGCATGGGCAt                                                                                                                                                                                                                                                            >  7:1061/1‑98 (MQ=255)
                                                                                                                                                                                           ccGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGAACGGCTTGCAGtt                                                                        >  1:113602/1‑275 (MQ=255)
                                                                                                                                                                                            cGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGAACTTCTGGCGGCGTACGAGCAGCACGGCTTGCAGTTa                                                                       >  1:113158/1‑275 (MQ=255)
                                                                                                                                                                                              gATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGCTTGTGGACCTGCTGGCGCCGTACGAGCAGCACGGCTTGCAGTTaa                                                                      >  2:456977/1‑274 (MQ=255)
                                                                                                                                                                                               aTGCTGCCATAAGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGa                                                                                                                                                                                                   >  7:1015131/1‑148 (MQ=255)
                                                                                                                                                                                                    gcCTTATGGCAGCATCAAGAAGAGATGTTTGAGGCGATTGCCGCGTCGAAATATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATgg                                                                                                                                                                                                                                  >  8:717623/1‑112 (MQ=255)
                                                                                                                                                                                                    gcCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCCCGGCTGCTGTTTGAACATGTGCACAGCGAATCCCGAGACCGCGGCCAGGCTATGGTGGACCTGATGGCGCAGTACGAGACGCCCGGCTTGCAGGTAACCAGcgc                                                               >  4:188882/1‑273 (MQ=255)
                                                                                                                                                                                                    gcCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCGGGCGCAGTACGAGCCGCCCGGCTTGCCGTTAAAcgaccc                                                               >  4:171685/1‑270 (MQ=255)
                                                                                                                                                                                                    gcCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCAGCGCATGGGCAtc                                                                                                                                                                                                                                                           <  8:380005/87‑2 (MQ=255)
                                                                                                                                                                                                    gcCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCAGCGCATGGGCAtc                                                                                                                                                                                                                                                           >  7:380005/1‑86 (MQ=255)
                                                                                                                                                                                                     ccTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTACTGGCGCAGTACGAGAAGCCCGGCTTGCAGTTAAACAGCcgc                                                              >  2:518948/1‑275 (MQ=255)
                                                                                                                                                                                                      cTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCt                                                                                                                                                                                                                           <  8:108979/117‑1 (MQ=255)
                                                                                                                                                                                                       ttATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCAGCGCATGGGCAt                                                                                                                                                                                                                                                            >  5:917031/1‑83 (MQ=255)
                                                                                                                                                                                                       ttATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCAGCGCATGGGCAt                                                                                                                                                                                                                                                            <  6:917031/83‑1 (MQ=255)
                                                                                                                                                                                                           ggCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTaa                                                                                                                                                                                                                                           >  8:111793/1‑96 (MQ=255)
                                                                                                                                                                                                           ggCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTaa                                                                                                                                                                                                                                           <  7:111793/96‑1 (MQ=255)
                                                                                                                                                                                                             cagcaTCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGta                                                                                                                                                                                     <  8:1015131/148‑2 (MQ=255)
                                                                                                                                                                                                               gcaTCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCCGGCGCAGTACAGCGAACTGTTCGACCCTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGAGAATCTCGCGACCGCGTCCCGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGCCTGCCGGTCAAACAGCCGCGAATTg                                                        >  1:47132/1‑271 (MQ=255)
                                                                                                                                                                                                               gcaTCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGAcc                                                                                                                               >  4:63445/1‑200 (MQ=255)
                                                                                                                                                                                                               gcaTCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTAAACAGCCGCGAATTg                                                        >  2:54634/1‑271 (MQ=255)
                                                                                                                                                                                                               gcaTCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGCCCGCGCCACGTCCCTTGTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCCGCCAGGCGATGGTTGACCTGCCGGCGCAGTACGAGCAGCCCGGCAGTCCGTTAAACAGCCGCGAATTGCCgg                                                    >  4:158913/1‑275 (MQ=255)
                                                                                                                                                                                                                   caacaaGAGGTGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGgc                                                                                                   <  2:435602/224‑1 (MQ=255)
                                                                                                                                                                                                                   caacaaGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCgg                                                                                                                           <  2:63446/200‑1 (MQ=255)
                                                                                                                                                                                                                     acaaGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGcgcg                                                                                                                                                                                                                                                 >  5:1037423/1‑80 (MQ=255)
                                                                                                                                                                                                                     acaaGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGcgcg                                                                                                                                                                                                                                                 <  6:1037423/80‑1 (MQ=255)
                                                                                                                                                                                                                     acaaGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTAc                                                                                                                                                                                                        >  7:448459/1‑121 (MQ=255)
                                                                                                                                                                                                                       aaGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCTTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGgcgc                                                                                                                                                                                                             >  5:111013/1‑114 (MQ=255)
                                                                                                                                                                                                                       aaGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGAt                                                                                                                                                                                                                                <  8:934324/95‑1 (MQ=255)
                                                                                                                                                                                                                       aaGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGgcgc                                                                                                                                                                                                             <  6:111013/114‑1 (MQ=255)
                                                                                                                                                                                                                       aaGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCAAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGAt                                                                                                                                                                                                                                >  7:934324/1‑95 (MQ=255)
                                                                                                                                                                                                                          agaTGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCa                                                                                                                                                                                                                  <  8:1065852/106‑1 (MQ=255)
                                                                                                                                                                                                                          agaTGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCa                                                                                                                                                                                                                  >  5:311176/1‑106 (MQ=255)
                                                                                                                                                                                                                          agaTGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCa                                                                                                                                                                                                                  >  7:1065852/1‑106 (MQ=255)
                                                                                                                                                                                                                            aTGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTAAACAGCCGCGAATTGCCGGACCATCTg                                            >  4:226747/1‑270 (MQ=255)
                                                                                                                                                                                                                             tGTTTTAGGCGATGTCTGCGTGGAAAAACCTGCCAAAGAGGGGTGCCCAGTGGCTGGGCATTTGCCGGCGGGGTTTAACGTCGTTGGAGCCGCTCGATGCCCAGGGGAAGTAAAGCGAGCTGGTCGACCGGGGCCGGGCCACGTCACAGTGGCTGTTGGAACATGTGCACGGCGAATCCCGGGACCGCGGGCAGGGGATGGTGGACCTGCTGGCGCAGGACGGGAAGAACAGGTTGCAGTTAAACAGCCGCGAATTGCCGGACCATCTGCCGCt                                       <  4:456976/274‑1 (MQ=255)
                                                                                                                                                                                                                               tttGGGGTGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGATCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTAAACAGCCGCGAATTGCCGGACCATCTGCCGCTGt                                     <  4:509305/274‑1 (MQ=255)
                                                                                                                                                                                                                                    ggCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCcgcg                                                                                                                                                                               >  7:849117/1‑131 (MQ=255)
                                                                                                                                                                                                                                    ggCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCcgcg                                                                                                                                                                               <  8:849117/131‑1 (MQ=255)
                                                                                                                                                                                                                                      cGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGAc                                                                                                                                                                                         <  8:476322/119‑1 (MQ=255)
                                                                                                                                                                                                                                      cGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGAc                                                                                                                                                                                         >  7:476322/1‑119 (MQ=255)
                                                                                                                                                                                                                                        aTTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGa                                                                                                                                                                                                   >  8:333696/1‑107 (MQ=255)
                                                                                                                                                                                                                                        aTTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGa                                                                                                                                                                                                   <  7:333696/107‑1 (MQ=255)
                                                                                                                                                                                                                                        aTTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCTTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTAAACAGCCGCGAATTGCCGGACCATCTGCCGCTGTATCTGGAGt                            >  2:142697/1‑274 (MQ=255)
                                                                                                                                                                                                                                         ttGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCTCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGtt                                                                                                                                                                                             <  7:954925/112‑1 (MQ=255)
                                                                                                                                                                                                                                                     aaaaTCTGCCAAAAGGGGATGCCCATGCGCTGGGCAGTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGACCTGTTCGACCGTGGCAGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGTCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTAAACAGCCGCAAAGTGCCGGACCATCTGCCGCTGTATCTGGAGTACCTGGCGTAGCt               <  2:139722/274‑1 (MQ=255)
                                                                                                                                                                                                                                                      aaaTCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACAtg                                                                                                                                                    <  6:769884/140‑1 (MQ=255)
                                                                                                                                                                                                                                                       aaTCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGttt                                                                                                                                                           <  8:504594/132‑1 (MQ=255)
                                                                                                                                                                                                                                                       aaTCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGttt                                                                                                                                                           >  7:504594/1‑132 (MQ=255)
                                                                                                                                                                                                                                                          cTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGa                                                                                                                                                                                                   <  6:237957/89‑1 (MQ=255)
                                                                                                                                                                                                                                                          cTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGa                                                                                                                                                                                                   >  5:237957/1‑89 (MQ=255)
                                                                                                                                                                                                                                                            gCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACa                                                                                                                                                      >  8:109526/1‑132 (MQ=255)
                                                                                                                                                                                                                                                            gCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACa                                                                                                                                                      <  7:109526/132‑1 (MQ=255)
                                                                                                                                                                                                                                                               aaaaGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGaa                                                                                                                                         >  5:938560/1‑142 (MQ=255)
                                                                                                                                                                                                                                                               aaaaGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGaa                                                                                                                                         <  6:938560/142‑1 (MQ=255)
                                                                                                                                                                                                                                                                aaaGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCcgcg                                                                                                                                  >  7:834586/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                 aaGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCa                                                                                                                                                                                                                  <  5:858651/67‑1 (MQ=255)
                                                                                                                                                                                                                                                                 aaGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCa                                                                                                                                                                                                                  >  6:858651/1‑67 (MQ=255)
                                                                                                                                                                                                                                                                 aaGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTAAACAGCCGCGAATTGCCGGACCATCTGCCGCTGTATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGaa  <  2:98258/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                    agGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACcg                                                                                                                              >  6:865786/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                       aTGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGt                                                                                                                                                                                              >  6:632083/1‑81 (MQ=255)
                                                                                                                                                                                                                                                                       aTGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGt                                                                                                                                                                                              <  5:632083/81‑1 (MQ=255)
                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                      
GTGAACTGGCACGGGAAGCCTTCCCGGAGAAAAATGGCTGCGGCTTTACCTTTGGTGATGGCTGCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTAAACAGCCGCGAATTGCCGGACCATCTGCCGCTGTATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAA  >  NC_000913/1284991‑1285522

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