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.

GCGCAAAGGCCAGTGGCAACCAGCCGGTTTGTGGGCGCAGTTTTTCCTGGCCGACATAGTGTGCCCAGCCGCCACCGCTTTGCCCGACACAACCGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAA  >  NC_000913/1540805‑1541348
                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                                
        gCCATTGGCAACCAGCCGGTTTGTGGGCGCAGTTTTTCCTGGCCGAACGAGTGTGACCAGCCGCCACCGCTGTGCCCGACCCAACCGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTa                                                                                                                                                                                                                                                                        <  4:90901/275‑1 (MQ=255)
                       ccGGTTTGTGGGCGCAGTTTTTCCTGGCCGACATAGTGTGCCCAGCCGCCACCGCTATGCCCGACACAACCGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATCTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATACTCATCCAGCGCCAATGCGTTTTATGGGCGGTATAGGCCAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTACTTTGCTCACCCCCGTCTGGGGTTTTCGGTTGAAtttgtgg                                                                                                                                                                                                                                                         >  4:466085/1‑271 (MQ=255)
                       ccGGTTTGTGGGCGCAGTTTTTCCTGGCCGACATAGTGTGCCCAGCAGCCCCCGCATTACCCGACCCAACCGCAGAAGATCGGCATATTGATCATACCGCGGTAGTTCATGTCCATTGGATACACGTGGTTAACACCTGCGCAGTGGATAATCATCGAGCGCTCATGCGTTTTCTGGGCGGTATCGGCAAATTCAAGGGCAAGGGTTTCGATTTAACTGCCGGGGACGCCAGGTTATTGCCCAACCCACGCCGGGGTGTGCGGGTTGACTTctg                                                                                                                                                                                                                                                          >  4:508/1‑272 (MQ=255)
                                                  ccGACATAGTGTGCCCAGCCGCCACCGCTTTGCCCGACACAACCGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATTCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTc                                                                                                                                                                                                                             >  2:405268/1‑275 (MQ=255)
                                                       aTAGTGTGCCCAGCCGCCACCGCTTTGCCCGACACAACCGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTNTTATGGGCGGTAACGGCAAATTCACGGGCGATGGTTTCACTATACTGGCGCGGAACACCGGTAATTTGCTCCCCNCCGGATGGGGGGTTCGGTTTGATTTTAGCGNNATCTTTTACCTGTTTTTCATTTTCCac                                                                                                                                                                                                                          >  3:75338/1‑272 (MQ=255)
                                                               cccAGCCGCCACCGCTTTGCCCGACACAACCGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATc                                                                                                                                                                                                                 >  2:78572/1‑274 (MQ=255)
                                                                   gccgccACCGCTTTGCCCGACACAACCGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGAGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGAACGCCGGTAATTTGCTCAACCCAGTCTGGGGTGTACGGTTTGATTTCAGCGTCATCTTTCTCACTGTTTTCATCTTCACGCCcccgttgcaggc                                                                                                                                                                                                            >  4:384209/1‑267 (MQ=255)
                                                                   gccacgacggtgttGCCCGACAAACCCGGAGAGGCTCAGCAAGTGGATCAACCTGGGGTTGGTCATGTGATTGTGATACCAGTGGTTACCCGCTGGCGCGAGGATAAGCACCCGGGCACAAGGCGTTTTATAGGGGGTATCGGCAAATTCACGGGCGATGGTTTCAAGATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGCGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGc                                                                                                                                                                                                            <  4:283967/263‑1 (MQ=255)
                                                                                            ccGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTACTTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCAGTACTTCGCCAGTACCAAATCATa                                                                                                                                                                                   >  1:115687/1‑275 (MQ=255)
                                                                                            ccGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCACGTGATACCAGTGGTTAACACCTGCCCCGAGGATAATCATCGAGCGCCCATGCGTTTTAAGGGCGGTATGCGCAAATTCCCGGGCGATGGTTTCAATCTACTGGCCCGGCACGCCGGTAATTTGCTCACCCCAGTCTGGGGTGGACGGTTTGAGTTCAGCGTAATATTTCGTCCGGTTTTCATCTTCCAGCCCGCCATCGAGGCCGTAATTCACACGTCCCAAATCATa                                                                                                                                                                                   >  4:452965/1‑275 (MQ=255)
                                                                                                 gaagaTCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGGACCAAATCATAAACGc                                                                                                                                                                              >  1:114923/1‑275 (MQ=255)
                                                                                                            tatTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTc                                                                                                                                                                                                                                                            >  4:251949/1‑187 (MQ=255)
                                                                                                                gATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCg                                                                                                                                                                                                                                                        <  2:251950/187‑1 (MQ=255)
                                                                                                                 atcatcCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGgt                                                                                                                                                              >  2:319269/1‑275 (MQ=255)
                                                                                                                         gcgGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCa                                                                                                                                                       <  2:374532/274‑1 (MQ=255)
                                                                                                                                ttCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTg                                                                                                                                                                                                                                                                               <  8:249620/148‑1 (MQ=255)
                                                                                                                                ttCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGAACACAGGACAGGTGTTGCCATCGACTa                                                                                                                                                >  2:79879/1‑274 (MQ=255)
                                                                                                                                    tGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGCGTTGCCATCGACTAACGTc                                                                                                                                           <  2:231781/275‑1 (MQ=255)
                                                                                                                                    tGTCCATGGGATTCCGGTGTTTAACACCTGCGCCGAAGATGATCATCGAGGGCCCATGCGTTTTATGGGCGGGATCGTCAAATTGACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGTCTGGGGTGTACGGTTTTATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTc                                                                                                                                           <  4:79878/275‑1 (MQ=255)
                                                                                                                                     gTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTc                                                                                                                                           >  2:371728/1‑274 (MQ=255)
                                                                                                                                     gTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTc                                                                                                                                           >  4:319852/1‑274 (MQ=255)
                                                                                                                                       ccATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCCGTACCAAATCATACACGCTGACCACCGGACAGGTGTTGCCATCGCCTACCGTCAgg                                                                                                                                        >  1:69033/1‑275 (MQ=255)
                                                                                                                                         aTGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGgtgt                                                                                                                                                                                                                                                                         >  6:1075670/1‑145 (MQ=255)
                                                                                                                                         aTGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGCCTGGGgtg                                                                                                                                                                                                                                                                          <  8:1090200/144‑1 (MQ=255)
                                                                                                                                              gtACCAGTGGTTAAAACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGt                                                                                                                                            <  2:388838/263‑1 (MQ=255)
                                                                                                                                                   aGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTc                                                                                                                                                                                                                                                            <  7:918649/148‑1 (MQ=255)
                                                                                                                                                     tGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTc                                                                                                                                                                                                                                                            <  5:379923/146‑1 (MQ=255)
                                                                                                                                                        ttAACACCTGCGCCGAGGATAATCATCGAGCGCCCACGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAgg                                                                                                                                                                                                                                                                                  >  6:999267/1‑121 (MQ=255)
                                                                                                                                                         tAACACCTGCGCCTAGGATAATCAGCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGTCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTgg                                                                                                                      <  1:34506/275‑1 (MQ=255)
                                                                                                                                                            cacCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAACTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATc                                                                                                                                                                                                                                                   >  5:795807/1‑148 (MQ=255)
                                                                                                                                                             aacTGCGCCGAGGATAAGAATCGTGCGCACATACGTTTTATGGGCGGTATCGGCAAATTCACGTGCGAGGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTAACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGgcgc                                                                                                                   <  2:538287/272‑1 (MQ=255)
                                                                                                                                                               cTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTACTTTGCTCACCCCAGGCCGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCCTCTTCCAGCTCGCGATCGAGCCCGTAATTCGCCAGTACCAAATCATAAACCCTTACCACCGGACAGGTTTCGCCATCGCCGAACGTCAGGTTTTTAACGGGCCATTGGCCCAcc                                                                                                                >  4:300767/1‑275 (MQ=255)
                                                                                                                                                                 gcgcCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGaccac                                                                                                                                                                       >  4:144935/1‑218 (MQ=255)
                                                                                                                                                                        ggATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGtt                                                                                                                                                                                                                                      >  5:447586/1‑148 (MQ=255)
                                                                                                                                                                        ggATAATCATCGAGAGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGa                                                                                                       <  2:480417/275‑1 (MQ=255)
                                                                                                                                                                          aTAATCCTCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCATGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGTTTTTTAACGGGCAATTGGCGCACCAGCACCGGAtt                                                                                                     <  2:529983/275‑1 (MQ=255)
                                                                                                                                                                          aTAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTgg                                                                                                                                                                                                                                                                              >  7:1070921/1‑107 (MQ=255)
                                                                                                                                                                          aTAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGtttt                                                                                                                                                                                                                                    <  6:447586/148‑1 (MQ=255)
                                                                                                                                                                          aTAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGt                                                                                                                                                                                                                                       <  8:278444/145‑1 (MQ=255)
                                                                                                                                                                          aTAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGt                                                                                                                                                                                                                                       >  7:278444/1‑145 (MQ=255)
                                                                                                                                                                             atcatcGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTgg                                                                                                                                                                                                                                                                              >  8:215975/1‑104 (MQ=255)
                                                                                                                                                                              tcatcGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGCCTGGGgtg                                                                                                                                                                                                                                                                          <  7:88405/107‑1 (MQ=255)
                                                                                                                                                                              tcatcGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGCCTGGGgtg                                                                                                                                                                                                                                                                          >  8:88405/1‑107 (MQ=255)
                                                                                                                                                                                 tcGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGt                                                                                                                                                                                                                                                                    >  6:157295/1‑110 (MQ=255)
                                                                                                                                                                                 tcGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGt                                                                                                                                                                                                                                                                    <  5:157295/110‑1 (MQ=255)
                                                                                                                                                                                 tcGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTc                                                                                                                                                                                                                             >  7:738684/1‑148 (MQ=255)
                                                                                                                                                                                 tcGAGCGCCCATGCGTTTCATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTc                                                                                                                                                                                                                             >  7:671894/1‑148 (MQ=255)
                                                                                                                                                                                                 tttATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTc                                                                                                                                                                                                                                                            <  7:314569/102‑1 (MQ=255)
                                                                                                                                                                                                   tATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGCCAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCCCCAGCCCCGGGTTGTGTTTTACGCTGCGAACATGCgg                                                                             >  2:18339/1‑274 (MQ=255)
                                                                                                                                                                                                   tATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCgg                                                                             >  2:264025/1‑274 (MQ=255)
                                                                                                                                                                                                           gTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGAATAACGTCAGGTTTTTAACGTGCAATTGGTGCACCCGCCCCGGCTTTGATTTTACGCGGCGAAAATGCGGAATTTAAAt                                                                    >  1:44423/1‑275 (MQ=255)
                                                                                                                                                                                                             aTCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGcc                                                                                                                                                                                                 >  5:49115/1‑148 (MQ=255)
                                                                                                                                                                                                             aTCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGcc                                                                                                                                                                                                 >  5:630886/1‑148 (MQ=255)
                                                                                                                                                                                                              tCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCa                                                                                                                                                                                                <  8:671894/148‑1 (MQ=255)
                                                                                                                                                                                                                 gCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTa                                                                                                                                                                                                        >  5:482798/1‑137 (MQ=255)
                                                                                                                                                                                                                       tCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTc                                                                                                                                                                                                                             <  8:1008881/110‑1 (MQ=255)
                                                                                                                                                                                                                            ggCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTgg                                                                                                                                                                                                                                                                              <  7:222346/57‑1 (MQ=255)
                                                                                                                                                                                                                            ggCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATaa                                                                                                                                                                                  >  8:677472/1‑148 (MQ=255)
                                                                                                                                                                                                                            ggCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGa                                                   <  2:136484/275‑1 (MQ=255)
                                                                                                                                                                                                                             gCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTgg                                                                                                                                                                                                                                                                              >  8:626430/1‑56 (MQ=255)
                                                                                                                                                                                                                             gCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTgg                                                                                                                                                                                                                                                                              <  7:626430/56‑1 (MQ=255)
                                                                                                                                                                                                                              cGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGttt                                                                                                                                                                                                                                                                  <  8:1012943/67‑1 (MQ=255)
                                                                                                                                                                                                                              cGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGttt                                                                                                                                                                                                                                                                  >  7:1012943/1‑67 (MQ=255)
                                                                                                                                                                                                                                aTGGTTTCATTATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGt                                                                                                                                                                                                                                       >  7:68765/1‑91 (MQ=255)
                                                                                                                                                                                                                                aTGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGt                                                                                                                                                                                                                                       <  8:68765/91‑1 (MQ=255)
                                                                                                                                                                                                                                 tGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGcc                                              <  2:319853/275‑1 (MQ=255)
                                                                                                                                                                                                                                  ggTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTg                                                                                                                                                                            <  7:677472/148‑1 (MQ=255)
                                                                                                                                                                                                                                    tttCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGgtgt                                                                                                                                                                                                                                                                         >  6:30699/1‑54 (MQ=255)
                                                                                                                                                                                                                                    tttCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGgtgt                                                                                                                                                                                                                                                                         <  5:30699/54‑1 (MQ=255)
                                                                                                                                                                                                                                     ttCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGAt                                                                                                                                                                                                                  >  6:260506/1‑107 (MQ=255)
                                                                                                                                                                                                                                      tCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGacca                                                                                                                                                                        <  6:384275/148‑1 (MQ=255)
                                                                                                                                                                                                                                          attaCCGGCGTGCCGCGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTg                                                                                                                                                                                                                                                                 >  6:902604/3‑56 (MQ=21)
                                                                                                                                                                                                                                          attaCCGGCGTGCCGCGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTg                                                                                                                                                                                                                                                                 <  5:902604/54‑1 (MQ=21)
                                                                                                                                                                                                                                             gaCTGGCGCGGTCCGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCTCGATCGAGGCCGTAATTCGCCAGTACCAAATCAGAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCaa                                    <  4:264024/273‑1 (MQ=255)
                                                                                                                                                                                                                                              aCTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACa                                                                                                                                                                 <  6:49115/148‑1 (MQ=255)
                                                                                                                                                                                                                                                tGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACCATGCGAAAATGCGGAATTTCAATGCCGCCAAACTAGGGTAAGGCCACGCATGCAACAGc                                >  4:139818/1‑275 (MQ=255)
                                                                                                                                                                                                                                                 ggCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTa                                                                                                                                                                                                        >  7:174301/1‑106 (MQ=255)
                                                                                                                                                                                                                                                   cgcgGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGCTGGTGTTGTACGCTGCGAAAATGCGGATTTTCAATGCCGCAAAAGTAGGGGAAGGCCACGCCTGAAAGAGcagc                             >  4:133247/1‑272 (MQ=255)
                                                                                                                                                                                                                                                   cgcgGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGCCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTTTGTTTTACGTTGCGAAACTCAGGATTTTAAATGCCGCCAATGTAGGGGAAGGCCACGCCTCCAACAGCGTc                             >  4:123600/1‑275 (MQ=255)
                                                                                                                                                                                                                                                     cgGCACGCCGGTAATTTGCCCACACCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGACAGTACCAAATCAGAAACGCTGACCACCGGCCAGGTGTTGCCATCGACTAACGTCAGGTTCTAAACGGGCAATCGGCGCCACAGCACCGGCTAGTGTTTTACGCTGCCAAAATGCGGCCTGTCAACGCCGCAAAAGGAGGGGAAGGCCCCGCCTGCAACCGCGGCAt                           >  4:280166/1‑275 (MQ=255)
                                                                                                                                                                                                                                                       gCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGt                         >  4:375368/1‑275 (MQ=255)
                                                                                                                                                                                                                                                       gCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGCACCAAATCATAAACGCTGCCCACCGGACAGGTGTTGCCATCGACTAACTTCAGGTTTTTAACGGACACTTGACGCACAAGCACCGGATTGTGTTTTACGCTCCGAAAATGCGGATTTTCAATGCAGCCAAAGTAGGGGAACGCCAACCCTGCAACAGCGTCATGt                         >  1:99936/1‑275 (MQ=255)
                                                                                                                                                                                                                                                         aggCCGGTAATTGGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTg                       <  2:299263/273‑1 (MQ=255)
                                                                                                                                                                                                                                                          cGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTAATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCGGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGa                      <  2:123601/275‑1 (MQ=255)
                                                                                                                                                                                                                                                            ccGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGc                                                                                                                                                                                                                                             >  5:218252/1‑59 (MQ=255)
                                                                                                                                                                                                                                                                 tAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCAGCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAg                 <  2:508/274‑1 (MQ=255)
                                                                                                                                                                                                                                                                        cTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCt                                                                                                                                                                             <  8:1137079/111‑1 (MQ=255)
                                                                                                                                                                                                                                                                         acACCCCAGCCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCAt                                                                                                                                                      <  7:297686/132‑1 (MQ=255)
                                                                                                                                                                                                                                                                         acACCCCAGCCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCAt                                                                                                                                                      >  8:297686/2‑133 (MQ=255)
                                                                                                                                                                                                                                                                          cACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAgg                                                                                                                                                                                                             >  8:388468/1‑77 (MQ=255)
                                                                                                                                                                                                                                                                          cACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTa                                                                                                                                                                                                        >  5:68368/1‑82 (MQ=255)
                                                                                                                                                                                                                                                                          cACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTa                                                                                                                                                                                                        <  6:68368/82‑1 (MQ=255)
                                                                                                                                                                                                                                                                          cACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATa                                                                                                                                                                                   <  6:487571/103‑1 (MQ=255)
                                                                                                                                                                                                                                                                           aCCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGttt                                                                                                                                     <  8:738684/148‑1 (MQ=255)
                                                                                                                                                                                                                                                                              ccAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCa                                                                                                                                                                                     >  7:1133701/1‑97 (MQ=255)
                                                                                                                                                                                                                                                                              ccAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGgtgt                                                                                                                                                            <  8:1128638/122‑1 (MQ=255)
                                                                                                                                                                                                                                                                              ccAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGgtgt                                                                                                                                                            >  7:1128638/1‑122 (MQ=255)
                                                                                                                                                                                                                                                                              ccAGCCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATaaa                                                                                                                                                                                 <  8:164123/101‑1 (MQ=255)
                                                                                                                                                                                                                                                                              ccAGCCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATaaa                                                                                                                                                                                 >  7:164123/1‑101 (MQ=255)
                                                                                                                                                                                                                                                                                aGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAAc                                                                                                                                >  8:309403/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                                
GCGCAAAGGCCAGTGGCAACCAGCCGGTTTGTGGGCGCAGTTTTTCCTGGCCGACATAGTGTGCCCAGCCGCCACCGCTTTGCCCGACACAACCGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAA  >  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: