Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,541,077 G→C 6.7% A505G (GCC→GGC)  narZ ← nitrate reductase 2 (NRZ), alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,541,0770GC 6.7% 511.2 / 22.8 150A505G (GCC→GGC) narZnitrate reductase 2 (NRZ), alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (73/63);  new base C (4/6);  total (80/70)
Fisher's exact test for biased strand distribution p-value = 5.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.79e-03

Only 100 of 154 total aligned reads displayed.

GCGCAAAGGCCAGTGGCAACCAGCCGGTTTGTGGGCGCAGTTTTTCCTGGCCGACATAGTGTGCCCAGCCGCCACCGCTTTGCCCGACACAACCGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATAT‑ACTGGCGCGGCACGCC‑GGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAA  >  NC_000913/1540805‑1541348
                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                                
        gCCATTGGCAACCAGCCGGTTTGTGGGCGCAGTTTTTCCTGGCCGAACGAGTGTGACCAGCCGCCACCGCTGTGCCCGACCCAACCGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATAT‑ACTGGCGCGGCACGCC‑GGTAATTTGCTCACCCCAGGCTGGGGTGTa                                                                                                                                                                                                                                                                        <  4:90901/275‑1 (MQ=255)
                       ccGGTTTGTGGGCGCAGTTTTTCCTGGCCGACATAGTGTGCCCAGCCGCCACCGCTATGCCCGACACAACCGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATCTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATACTCATCCAGCGCCAATGCGTTTTATGGGCGGTATAGGCCAATTCACGGGCGATGGTTTCAATAT‑ACTGGCGCGGCACGCC‑GGTACTTTGCTCACCCCCGTCTGGGGTTTTCGGTTGAAtttgtgg                                                                                                                                                                                                                                                         >  4:466085/1‑271 (MQ=255)
                       ccGGTTTGTGGGCGCAGTTTTTCCTGGCCGACATAGTGTGCCCAGCAGCCCCCGCATTACCCGACCCAACCGCAGAAGATCGGCATATTGATCATACCGCGGTAGTTCATGTCCATTGGATACACGTGGTTAACACCTGCGCAGTGGATAATCATCGAGCGCTCATGCGTTTTCTGGGCGGTATCGGCAAATTCAAGGGCAAGGGTTTCGATTTAACTGCCG‑GGGACGCCAGGTTA‑TTGCCCAACCCACGCCGGGGTGTGCGGGTTGACTTctg                                                                                                                                                                                                                                                          >  4:508/1‑272 (MQ=255)
                                                  ccGACATAGTGTGCCCAGCCGCCACCGCTTTGCCCGACACAACCGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATTCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATAT‑ACTGGCGCGGCACGCC‑GGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTc                                                                                                                                                                                                                             >  2:405268/1‑275 (MQ=255)
                                                       aTAGTGTGCCCAGCCGCCACCGCTTTGCCCGACACAACCGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTNTTATGGGCGGTAACGGCAAATTCACGGGCGATGGTTTCACTAT‑ACTGGCGCGGAACACC‑GGTAATTTGCTCCCCNCCGGATGGGGGGTTCGGTTTGATTTTAGCGNNATCTTTTAC‑CTGTTTTTCATTTTCCac                                                                                                                                                                                                                          >  3:75338/1‑272 (MQ=255)
                                                               cccAGCCGCCACCGCTTTGCCCGACACAACCGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATAT‑ACTGGCGCGGCACGCC‑GGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATc                                                                                                                                                                                                                 >  2:78572/1‑274 (MQ=255)
                                                                   gccgccACCGCTTTGCCCGACACAACCGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGAGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATAT‑ACTGGCGCGGAACGCC‑GGTAATTTGCTCAACCCAGTCTGGGGTGTACGGTTTGATTTCAGCGTCATCTTTCTCACTG‑TTTTCATCTTCACGCCcccgttgcaggc                                                                                                                                                                                                            >  4:384209/1‑267 (MQ=255)
                                                                   gccacgacggtgttGCCCGACAAACCCGGAGAGGCTCAGCAAGTGGATCAACCTGGGGTTGGTCATGTGATTGTGATACCAGTGGTTACCCGCTGGCGCGAGGATAAGCACCCGGGCACAAGGCGTTTTATAGGGGGTATCGGCAAATTCACGGGCGATGGTTTCAAGAT‑ACTGGCGCGGCACGCC‑GGTAATTTGCTCACCCCAGGCTGGCGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGc                                                                                                                                                                                                            <  4:283967/263‑1 (MQ=255)
                                                                                            ccGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATAT‑ACTGGCGCGGCACGCC‑GGTACTTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCAGTACTTCGCCAGTACCAAATCATa                                                                                                                                                                                   >  1:115687/1‑275 (MQ=255)
                                                                                            ccGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCACGTGATACCAGTGGTTAACACCTGCCCCGAGGATAATCATCGAGCGCCCATGCGTTTTAAGGGCGGTATGCGCAAATTCCCGGGCGATGGTTTCAATCT‑ACTGGCCCGGCACGCC‑GGTAATTTGCTCACCCCAGTCTGGGGTGGACGGTTTGAGTTCAGCGTAATATTTCGTCCGG‑TTTTCATCTTCCAGCCCGCCATCGAGGCCGTAATTCACACGTCCCAAATCATa                                                                                                                                                                                   >  4:452965/1‑275 (MQ=255)
                                                                                                 gaagaTCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATAT‑ACTGGCGCGGCACGCC‑GGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGGACCAAATCATAAACGc                                                                                                                                                                              >  1:114923/1‑275 (MQ=255)
                                                                                                            tatTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATAT‑ACTGGCGCGGCACGCC‑GGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTc                                                                                                                                                                                                                                                            >  4:251949/1‑187 (MQ=255)
                                                                                                                gATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATAT‑ACTGGCGCGGCACGCC‑GGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCg                                                                                                                                                                                                                                                        <  2:251950/187‑1 (MQ=255)
                                                                                                                 atcatcCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATAT‑ACTGGCGCGGCACGCC‑GGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGgt                                                                                                                                                              >  2:319269/1‑275 (MQ=255)
                                                                                                                         gcgGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATAT‑ACTGGCGCGGCACGCC‑GGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCa                                                                                                                                                       <  2:374532/274‑1 (MQ=255)
                                                                                                                                ttCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATAT‑ACTGGCGCGGCACGCC‑GGTAATTTGCTCACCCCAGGCTg                                                                                                                                                                                                                                                                               <  8:249620/148‑1 (MQ=255)
                                                                                                                                ttCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATAT‑ACTGGCGCGGCACGCC‑GGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGAACACAGGACAGGTGTTGCCATCGACTa                                                                                                                                                >  2:79879/1‑274 (MQ=255)
                                                                                                                                    tGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATAT‑ACTGGCGCGGCACGCC‑GGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGCGTTGCCATCGACTAACGTc                                                                                                                                           <  2:231781/275‑1 (MQ=255)
                                                                                                                                    tGTCCATGGGATTCCGGTGTTTAACACCTGCGCCGAAGATGATCATCGAGGGCCCATGCGTTTTATGGGCGGGATCGTCAAATTGACGGGCGATGGTTTCAATAT‑ACTGGCGCGGCACGCC‑GGTAATTTGCTCACCCCAGTCTGGGGTGTACGGTTTTATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTc                                                                                                                                           <  4:79878/275‑1 (MQ=255)
                                                                                                                                     gTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATAT‑ACTGGCGCGGCACGCC‑GGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTc                                                                                                                                           >  2:371728/1‑274 (MQ=255)
                                                                                                                                     gTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATAT‑ACTGGCGCGGCACGCC‑GGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTc                                                                                                                                           >  4:319852/1‑274 (MQ=255)
                                                                                                                                       ccATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATAT‑ACTGGCGCGGCACGCC‑GGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCCGTACCAAATCATACACGCTGACCACCGGACAGGTGTTGCCATCGCCTACCGTCAgg                                                                                                                                        >  1:69033/1‑275 (MQ=255)
                                                                                                                                         aTGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATAT‑ACTGGCGCGGCACGCC‑GGTAATTTGCTCACCCCAGGCTGGGgtgt                                                                                                                                                                                                                                                                         >  6:1075670/1‑145 (MQ=255)
                                                                                                                                         aTGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATAT‑ACTGGCGCGGCACGCC‑GGTAATTTGCTCACCCCAGCCTGGGgtg                                                                                                                                                                                                                                                                          <  8:1090200/144‑1 (MQ=255)
                                                                                                                                              gtACCAGTGGTTAAAACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATAT‑ACTGGCGCGGCACGCC‑GGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGt                                                                                                                                            <  2:388838/263‑1 (MQ=255)
                                                                                                                                                   aGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATAT‑ACTGGCGCGGCACGCC‑GGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTc                                                                                                                                                                                                                                                            <  7:918649/148‑1 (MQ=255)
                                                                                                                                                     tGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATAT‑ACTGGCGCGGCACGCC‑GGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTc                                                                                                                                                                                                                                                            <  5:379923/146‑1 (MQ=255)
                                                                                                                                                        ttAACACCTGCGCCGAGGATAATCATCGAGCGCCCACGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATAT‑ACTGGCGCGGCACGCC‑GGTAATTTGCTCACCCCAgg                                                                                                                                                                                                                                                                                  >  6:999267/1‑121 (MQ=255)
                                                                                                                                                         tAACACCTGCGCCTAGGATAATCAGCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGTCGATGGTTTCAATAT‑ACTGGCGCGGCACGCC‑GGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTgg                                                                                                                      <  1:34506/275‑1 (MQ=255)
                                                                                                                                                            cacCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATAT‑ACTGGCGCGGCACGCC‑GGTAACTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATc                                                                                                                                                                                                                                                   >  5:795807/1‑148 (MQ=255)
                                                                                                                                                             aacTGCGCCGAGGATAAGAATCGTGCGCACATACGTTTTATGGGCGGTATCGGCAAATTCACGTGCGAGGGTTTCAATAT‑ACTGGCGCGGCACGCC‑GGTAATTTGCTAACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGgcgc                                                                                                                   <  2:538287/272‑1 (MQ=255)
                                                                                                                                                               cTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATAT‑ACTGGCGCGGCACGCC‑GGTACTTTGCTCACCCCAGGCCGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCCTCTTCCAGCTCGCGATCGAGCCCGTAATTCGCCAGTACCAAATCATAAACCCTTACCACCGGACAGGTTTCGCCATCGCCGAACGTCAGGTTTTTAACGGGCCATTGGCCCAcc                                                                                                                >  4:300767/1‑275 (MQ=255)
                                                                                                                                                                 gcgcCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATAT‑ACTGGCGCGGCACGCC‑GGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGaccac                                                                                                                                                                       >  4:144935/1‑218 (MQ=255)
                                                                                                                                                                        ggATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATAT‑ACTGGCGCGGCACGCC‑GGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑tt                                                                                                                                                                                                                                      >  5:447586/1‑148 (MQ=255)
                                                                                                                                                                        ggATAATCATCGAGAGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATAT‑ACTGGCGCGGCACGCC‑GGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGa                                                                                                       <  2:480417/275‑1 (MQ=255)
                                                                                                                                                                          aTAATCCTCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATAT‑ACTGGCGCGGCACGCC‑GGTAATTTGCTCACCCCATGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGTTTTTTAACGGGCAATTGGCGCACCAGCACCGGAtt                                                                                                     <  2:529983/275‑1 (MQ=255)
                                                                                                                                                                          aTAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATAT‑ACTGGCGCGGCACGCC‑GGTAATTTGCTCACCCCAGGCTgg                                                                                                                                                                                                                                                                              >  7:1070921/1‑107 (MQ=255)
                                                                                                                                                                          aTAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATAT‑ACTGGCGCGGCACGCC‑GGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑tttt                                                                                                                                                                                                                                    <  6:447586/148‑1 (MQ=255)
                                                                                                                                                                          aTAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATAT‑ACTGGCGCGGCACGCC‑GGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑t                                                                                                                                                                                                                                       <  8:278444/145‑1 (MQ=255)
                                                                                                                                                                          aTAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATAT‑ACTGGCGCGGCACGCC‑GGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑t                                                                                                                                                                                                                                       >  7:278444/1‑145 (MQ=255)
                                                                                                                                                                             atcatcGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATAT‑ACTGGCGCGGCACGCC‑GGTAATTTGCTCACCCCAGGCTgg                                                                                                                                                                                                                                                                              >  8:215975/1‑104 (MQ=255)
                                                                                                                                                                              tcatcGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATAT‑ACTGGCGCGGCACGCC‑GGTAATTTGCTCACCCCAGCCTGGGgtg                                                                                                                                                                                                                                                                          <  7:88405/107‑1 (MQ=255)
                                                                                                                                                                              tcatcGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATAT‑ACTGGCGCGGCACGCC‑GGTAATTTGCTCACCCCAGCCTGGGgtg                                                                                                                                                                                                                                                                          >  8:88405/1‑107 (MQ=255)
                                                                                                                                                                                 tcGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATAT‑ACTGGCGCGGCACGCC‑GGTAATTTGCTCACCCCAGGCTGGGGTGTACGGt                                                                                                                                                                                                                                                                    >  6:157295/1‑110 (MQ=255)
                                                                                                                                                                                 tcGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATAT‑ACTGGCGCGGCACGCC‑GGTAATTTGCTCACCCCAGGCTGGGGTGTACGGt                                                                                                                                                                                                                                                                    <  5:157295/110‑1 (MQ=255)
                                                                                                                                                                                 tcGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATAT‑ACTGGCGCGGCACGCC‑GGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTc                                                                                                                                                                                                                             >  7:738684/1‑148 (MQ=255)
                                                                                                                                                                                 tcGAGCGCCCATGCGTTTCATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATAT‑ACTGGCGCGGCACGCC‑GGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTc                                                                                                                                                                                                                             >  7:671894/1‑148 (MQ=255)
                                                                                                                                                                                                 tttATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATAT‑ACTGGCGCGGCACGCC‑GGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTc                                                                                                                                                                                                                                                            <  7:314569/102‑1 (MQ=255)
                                                                                                                                                                                                   tATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATAT‑ACTGGCGCGGCACGCC‑GGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGCCAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCCCCAGCCCCGGGTTGTGTTTTACGCTGCGAACATGCgg                                                                             >  2:18339/1‑274 (MQ=255)
                                                                                                                                                                                                   tATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATAT‑ACTGGCGCGGCACGCC‑GGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCgg                                                                             >  2:264025/1‑274 (MQ=255)
                                                                                                                                                                                                           gTATCGGCAAATTCACGGGCGATGGTTTCAATAT‑ACTGGCGCGGCACGCC‑GGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGAATAACGTCAGGTTTTTAACGTGCAATTGGTGCACCCGCCCCGGCTTTGATTTTACGCGGCGAAAATGCGGAATTTAAAt                                                                    >  1:44423/1‑275 (MQ=255)
                                                                                                                                                                                                             aTCGGCAAATTCACGGGCGATGGTTTCAATAT‑ACTGGCGCGGCACGCC‑GGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGcc                                                                                                                                                                                                 >  5:49115/1‑148 (MQ=255)
                                                                                                                                                                                                             aTCGGCAAATTCACGGGCGATGGTTTCAATAT‑ACTGGCGCGGCACGCC‑GGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGcc                                                                                                                                                                                                 >  5:630886/1‑148 (MQ=255)
                                                                                                                                                                                                              tCGGCAAATTCACGGGCGATGGTTTCAATAT‑ACTGGCGCGGCACGCC‑GGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCa                                                                                                                                                                                                <  8:671894/148‑1 (MQ=255)
                                                                                                                                                                                                                 gCAAATTCACGGGCGATGGTTTCAATAT‑ACTGGCGCGGCACGCC‑GGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTa                                                                                                                                                                                                        >  5:482798/1‑137 (MQ=255)
                                                                                                                                                                                                                       tCACGGGCGATGGTTTCAATAT‑ACTGGCGCGGCACGCC‑GGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTc                                                                                                                                                                                                                             <  8:1008881/110‑1 (MQ=255)
                                                                                                                                                                                                                            ggCGATGGTTTCAATAT‑ACTGGCGCGGCACGCC‑GGTAATTTGCTCACCCCAGGCTgg                                                                                                                                                                                                                                                                              <  7:222346/57‑1 (MQ=255)
                                                                                                                                                                                                                            ggCGATGGTTTCAATAT‑ACTGGCGCGGCACGCC‑GGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATaa                                                                                                                                                                                  >  8:677472/1‑148 (MQ=255)
                                                                                                                                                                                                                            ggCGATGGTTTCAATAT‑ACTGGCGCGGCACGCC‑GGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGa                                                   <  2:136484/275‑1 (MQ=255)
                                                                                                                                                                                                                             gCGATGGTTTCAATAT‑ACTGGCGCGGCACGCC‑GGTAATTTGCTCACCCCAGGCTgg                                                                                                                                                                                                                                                                              >  8:626430/1‑56 (MQ=255)
                                                                                                                                                                                                                             gCGATGGTTTCAATAT‑ACTGGCGCGGCACGCC‑GGTAATTTGCTCACCCCAGGCTgg                                                                                                                                                                                                                                                                              <  7:626430/56‑1 (MQ=255)
                                                                                                                                                                                                                              cGATGGTTTCAATAT‑ACTGGCGCGGCACGCC‑GGTAATTTGCTCACCCCAGGCTGGGGTGTACGGttt                                                                                                                                                                                                                                                                  <  8:1012943/67‑1 (MQ=255)
                                                                                                                                                                                                                              cGATGGTTTCAATAT‑ACTGGCGCGGCACGCC‑GGTAATTTGCTCACCCCAGGCTGGGGTGTACGGttt                                                                                                                                                                                                                                                                  >  7:1012943/1‑67 (MQ=255)
                                                                                                                                                                                                                                aTGGTTTCATTAT‑ACTGGCGCGGCACGCC‑GGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑t                                                                                                                                                                                                                                       >  7:68765/1‑91 (MQ=255)
                                                                                                                                                                                                                                aTGGTTTCAATAT‑ACTGGCGCGGCACGCC‑GGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑t                                                                                                                                                                                                                                       <  8:68765/91‑1 (MQ=255)
                                                                                                                                                                                                                                 tGGTTTCAATAT‑ACTGGCGCGGCACGCC‑GGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGcc                                              <  2:319853/275‑1 (MQ=255)
                                                                                                                                                                                                                                  ggTTTCAATAT‑ACTGGCGCGGCACGCC‑GGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTg                                                                                                                                                                            <  7:677472/148‑1 (MQ=255)
                                                                                                                                                                                                                                    tttCAATAT‑ACTGGCGCGGCACGCC‑GGTAATTTGCTCACCCCAGGCTGGGgtgt                                                                                                                                                                                                                                                                         >  6:30699/1‑54 (MQ=255)
                                                                                                                                                                                                                                    tttCAATAT‑ACTGGCGCGGCACGCC‑GGTAATTTGCTCACCCCAGGCTGGGgtgt                                                                                                                                                                                                                                                                         <  5:30699/54‑1 (MQ=255)
                                                                                                                                                                                                                                     ttCAATAT‑ACTGGCGCGGCACGCC‑GGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGAt                                                                                                                                                                                                                  >  6:260506/1‑107 (MQ=255)
                                                                                                                                                                                                                                      tCAATAT‑ACTGGCGCGGCACGCC‑GGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGacca                                                                                                                                                                        <  6:384275/148‑1 (MQ=255)
                                                                                                                                                                                                                                          att‑aCCGGCGTGCCGCGCC‑AGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTg                                                                                                                                                                                                                                                                 >  6:902604/3‑56 (MQ=21)
                                                                                                                                                                                                                                          att‑aCCGGCGTGCCGCGCC‑AGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTg                                                                                                                                                                                                                                                                 <  5:902604/54‑1 (MQ=21)
                                                                                                                                                                                                                                             gaCTGGCGCGGTCCGCC‑GGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCTCGATCGAGGCCGTAATTCGCCAGTACCAAATCAGAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCaa                                    <  4:264024/273‑1 (MQ=255)
                                                                                                                                                                                                                                              aCTGGCGCGGCACGCC‑GGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACa                                                                                                                                                                 <  6:49115/148‑1 (MQ=255)
                                                                                                                                                                                                                                                tGGCGCGGCACGCC‑GGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACCATGCGAAAATGCGGAATTTCAATGCCGCCAAACTAGGGTAAGGCCACGCATGCAACAGc                                >  4:139818/1‑275 (MQ=255)
                                                                                                                                                                                                                                                 ggCGCGGCACGCC‑GGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTa                                                                                                                                                                                                        >  7:174301/1‑106 (MQ=255)
                                                                                                                                                                                                                                                   cgcgGCACGCC‑GGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGCTGGTGTTGTACGCTGCGAAAATGCGGATTTTCAATGCCGCAAAAGTAGGGGAAGGCCACGCCTGAAAGAGcagc                             >  4:133247/1‑272 (MQ=255)
                                                                                                                                                                                                                                                   cgcgGCACGCC‑GGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGCCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTTTGTTTTACGTTGCGAAACTCAGGATTTTAAATGCCGCCAATGTAGGGGAAGGCCACGCCTCCAACAGCGTc                             >  4:123600/1‑275 (MQ=255)
                                                                                                                                                                                                                                                     cgGCACGCC‑GGTAATTTGCCCACACCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGACAGTACCAAATCAGAAACGCTGACCACCGGCCAGGTGTTGCCATCGACTAACGTCAGGTTCTAAACGGGCAATCGGCGCCACAGCACCGGCTAGTGTTTTACGCTGCCAAAATGCGGCCTGTCAACGCCGCAAAAGGAGGGGAAGGCCCCGCCTGCAACCGCGGCAt                           >  4:280166/1‑275 (MQ=255)
                                                                                                                                                                                                                                                       gCACGCC‑GGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGt                         >  4:375368/1‑275 (MQ=255)
                                                                                                                                                                                                                                                       gCACGCC‑GGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGCACCAAATCATAAACGCTGCCCACCGGACAGGTGTTGCCATCGACTAACTTCAGGTTTTTAACGGACACTTGACGCACAAGCACCGGATTGTGTTTTACGCTCCGAAAATGCGGATTTTCAATGCAGCCAAAGTAGGGGAACGCCAACCCTGCAACAGCGTCATGt                         >  1:99936/1‑275 (MQ=255)
                                                                                                                                                                                                                                                         aggCC‑GGTAATTGGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTg                       <  2:299263/273‑1 (MQ=255)
                                                                                                                                                                                                                                                          cGCC‑GGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTAATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCGGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGa                      <  2:123601/275‑1 (MQ=255)
                                                                                                                                                                                                                                                            cc‑GGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGc                                                                                                                                                                                                                                             >  5:218252/1‑59 (MQ=255)
                                                                                                                                                                                                                                                                 tAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCAGCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAg                 <  2:508/274‑1 (MQ=255)
                                                                                                                                                                                                                                                                        cTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCt                                                                                                                                                                             <  8:1137079/111‑1 (MQ=255)
                                                                                                                                                                                                                                                                         acACCCCAGCCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCAt                                                                                                                                                      <  7:297686/132‑1 (MQ=255)
                                                                                                                                                                                                                                                                         acACCCCAGCCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCAt                                                                                                                                                      >  8:297686/2‑133 (MQ=255)
                                                                                                                                                                                                                                                                          cACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAgg                                                                                                                                                                                                             >  8:388468/1‑77 (MQ=255)
                                                                                                                                                                                                                                                                          cACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTa                                                                                                                                                                                                        >  5:68368/1‑82 (MQ=255)
                                                                                                                                                                                                                                                                          cACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTa                                                                                                                                                                                                        <  6:68368/82‑1 (MQ=255)
                                                                                                                                                                                                                                                                          cACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATa                                                                                                                                                                                   <  6:487571/103‑1 (MQ=255)
                                                                                                                                                                                                                                                                           aCCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGttt                                                                                                                                     <  8:738684/148‑1 (MQ=255)
                                                                                                                                                                                                                                                                              ccAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCa                                                                                                                                                                                     >  7:1133701/1‑97 (MQ=255)
                                                                                                                                                                                                                                                                              ccAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGgtgt                                                                                                                                                            <  8:1128638/122‑1 (MQ=255)
                                                                                                                                                                                                                                                                              ccAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGgtgt                                                                                                                                                            >  7:1128638/1‑122 (MQ=255)
                                                                                                                                                                                                                                                                              ccAGCCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATaaa                                                                                                                                                                                 <  8:164123/101‑1 (MQ=255)
                                                                                                                                                                                                                                                                              ccAGCCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATaaa                                                                                                                                                                                 >  7:164123/1‑101 (MQ=255)
                                                                                                                                                                                                                                                                                aGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAAc                                                                                                                                >  8:309403/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                                
GCGCAAAGGCCAGTGGCAACCAGCCGGTTTGTGGGCGCAGTTTTTCCTGGCCGACATAGTGTGCCCAGCCGCCACCGCTTTGCCCGACACAACCGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATAT‑ACTGGCGCGGCACGCC‑GGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAA  >  NC_000913/1540805‑1541348

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