Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,542,110 A→G 26.9% S161P (TCC→CCC)  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,542,1100AG26.9% 144.6 / 13.2 88S161P (TCC→CCC) narZnitrate reductase 2 (NRZ), alpha subunit
Reads supporting (aligned to +/- strand):  ref base A (38/24);  new base G (2/23);  total (40/48)
Fisher's exact test for biased strand distribution p-value = 3.51e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCGGTACGTCGGTTTGCTCGCCCCAGGTCATCGGCGAGGCGGGCGGCAGGTCGCAATACCAGTC‑GTAGAAACTTAAACAGGTG‑CCGCCAAGCAGCGACAGATAACGCGTTCCGGCGGCGTAAGAAACCATCGACATCGCCGGGATCGGCGAGAAACCGGCAACGCGATCCGGGCCGTAG‑GTTTTGATGGTCC‑AGA‑CGTTAGCGGCGGCAATCAGCTGGTTTAGTTCC‑TGCCAGTTGGAGCGGATAAA‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCTTCGCGCCACAGTTCAATCAGTCGTTTACGAATGAGCGGGTATTTCAGGCGGTTAGCGCTGTAAAGATACCAGGAGTAACTTGCGCCACGCGGGCAGCCGCGAGGTTCATGATTGGGCAGGTCAGGGC  >  NC_000913/1541872‑1542346
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                             
tCCGGTACGTCGGTTTGCTCGCCCCAGGTCATCGGCGAGGCGGGCGGCAGGTCGCAATACCAGTC‑GTAGAAACTTAAACAGGTG‑CCGCCAAGCAGCGACAGATAACGCGTTCCGGCGGCGTAAGAAACCATCGACATCGCCGGGATCGGCGAGAAACCGGCAACGCGATCCGGGCCGTAG‑GTTTTGATGGTCC‑AGA‑CGTTAGCGGCGGCAATCAGCTGGTTTAGTTCC‑TGCCAGTTGGa                                                                                                                                                                                                                                               >  1:145468/1‑239 (MQ=255)
tCCGGTACGTCGGTTTGCTCGCCCCAGGTCACCGGCGAGGCGGGCGGCAGGTCGCAATACCAGTC‑GTAGAAACTTAAACAGGTG‑CCGCCAAGCAGCGACAGATAACGCGTTCCGGCGGCGTAAGAAACCATCGACATCGCCGGGATCGGCGAGAAACCGGCAACGCGATCCGGGCCGTAG‑GTTTTGATGGTCC‑AGA‑CGTTAGCGGCGGCAATCAGCTGGTTTAGTTCC‑TGCCAGTTGGa                                                                                                                                                                                                                                               <  1:1088706/239‑1 (MQ=255)
         tCGGTTTGCTCGCCCCAGGTCATCGGCGAGGCGGGCGGCAGGTCGCAATACCAGTC‑GTAGAAACTTAAACAGGTG‑CCGCCAAGCAGCGACAGATAACGCGTTCCGGCGGCGTAAGAAACCATCGACATCGCCGGGATCGGCGAGAAACCGGCAACGCGATCCGGGCCGTAG‑GGTTTGATGGTCC‑AGA‑CGTTAGCGGCGGCAATCAGCTGGTTTAGTTCC‑TGCCAGTTGGAGCGGGTaa                                                                                                                                                                                                                                       >  1:648790/1‑238 (MQ=255)
          cGGTTTGCTCGCCCCAGGTCATCGGCGAGGCGGGCGGCAGGTCGCAATACCAGTC‑GTAGAAACTTAAACAGGTG‑CCGCCAAGCAGCGACAGATAACGCGTTCCGGCGGCGTAAGAAACCATCGACATCGCCGGGATCGGCGAGAAACCGGCAACGCGATCCGGGCCGTAG‑GTTTTGATGGTCC‑AGA‑CGTTAGCGGCGGCAATCAGCTGGTTTAGTTCC‑TGCCAGTTGGAGCGGATaaa                                                                                                                                                                                                                                      >  1:190364/1‑238 (MQ=255)
               tGCTCGCCCCAGGTCATCGGCGAGGCGGGCGGCAGGTCGCCATACCAGTC‑GTAGAAACTTAAACAGGTG‑CCGCCAAGCAGCGACAGATAACGCGTTCCGGCGGCGTAAGAAACCATCGACATCGCCGGGATCGGCGAGGAACCGGCAACGCGATCCCGGCCGTAG‑GTTTTGATGGTCA‑CGA‑CGTTAGCGGCGGCAAACAAATGGGTTAGTTCCTTCCCATTTGGGGCGGAGAAA‑ccgc                                                                                                                                                                                                                                 >  2:459760/1‑236 (MQ=255)
                   cGCCCGAGGTGGTCGGGGAGGCGGGCGGCGGGTCGAAATCCAAGCG‑GTAGAAACTAAAACAGGGC‑CCGCCAAGCGGCGACAAATAACGGTTTCCGGGGGCGAAAGAACCCACCGACTTCGCGGGGATGGGCGAGAAACCGGCAGCGCAACCCGGGCCGTAG‑GTTTGGATGGTCC‑AGC‑CGTTAGCGGCGGCAATCAGCTGGTTTAGTTCC‑TCCCAGTTGGAGCGGATAAA‑CCCGCCgcgc                                                                                                                                                                                                                           <  2:396825/239‑1 (MQ=255)
                   cGCCCCGGGCCATCGGGGGGGCGGGCGGCAGGTCAAAATACCAGTC‑GTAGAACCTAAAACAGGTG‑CCGCCAAGCGGCGACAGATACCGCGTTCCGGCGGCGAAAGAAACCACCGACACCGCCGGGATCGGCGAGAAACCGGCAACGCGATCCGGGCCGTAG‑GTTTTGATGGTCC‑AGA‑CTTTAGCGGCGGCAATCAGCTGGTTTAGTTCC‑TGCCAGTGGGAGCGGATAAA‑CCCGCCgcgc                                                                                                                                                                                                                           <  2:1013135/239‑1 (MQ=255)
                           ggtcacgggggagggggcggGAGGTCGAAATACCAGCG‑GTAGAAACTAAAACAGGTG‑CCCCCAAGCAGCGACAGAAGACGGGTTCCGGCGGGTTAAAAACCCATCGGAATCGCGGGGATCGGCGAGAAACCGGCAACGCGATCCGGGCGGTAG‑GTTTTGATGGTCC‑AGC‑CGTTAGCGGCGGCAATCAGCTGGTTTATTTCC‑TGCCAGTTGGAGCGGATAAA‑CCCGCCGCGCCCACGCa                                                                                                                                                                                                                    <  2:477377/224‑1 (MQ=255)
                             gagggggggggcgggggGGAGGCCCCAAT‑CCAGGCGGTAGAAACTAAAACAGGGG‑CCGCCAAGCGGGGCAAGATAACGCGTTCGGGGGGCGTAAGAACCCATCGAAACCGCCGGGATCGGCGAGAAACCGGAAACGCGATCCGGGCCGTAG‑GTTTTGATGGTCC‑AGA‑CGTTAGCGGCGGCAATCAGCTGTTTTAGTTCC‑TGCCAGTTGGAGCGAAAAAA‑CCCGCCGCGCCCACGCACtt                                                                                                                                                                                                                 <  2:287276/230‑1 (MQ=255)
                                cGGTGAGGCGGGCGGCAGGTCGCAATACCAGTC‑GTAGAAACTTAAACAGGTG‑CCGCCAAGCAGCGACAGATAACGCGTTCCGGCGGCGTAAGAAACCATCGACATCGCCGGGATCGGCGAGAAACCGGCAACGCGATCCGGGCCGTAG‑GGTTTGATGGTCC‑AGA‑CGTTAGCGGCGGCAATCAGCTGGTTTAGTTCC‑TGCCAGTTGGAGCGGATAAA‑CCCGCCGCGCCCCCGCACttgtt                                                                                                                                                                                                              >  1:646218/1‑239 (MQ=255)
                                cGGCGAGGCGGGCGGCAGGTCGCAATACCAGTC‑GTAGAAACTTAAACAGGTG‑CCGCCAAGCAGCGACAGATAACGCATTCCGGCGGCGTAAGAAACCATCGACATCGCCGGGATCGGCGAGAAACCGGCAACGCGATCCGGGCCGTAG‑GTTTTGATGGTCC‑AGA‑CGTTAGCGGCGGCAATCAGCTGGTTTAGTTCC‑TGCCCGTTGGAGCGGATAAA‑CCCGCCGCGCCCACGCACttgtt                                                                                                                                                                                                              >  1:203740/1‑239 (MQ=255)
                                                  gTCGCAATACCAGTC‑GTAGAAACTTAAACAGGTG‑CCGCCAAGCAGCGACAGATAACGCGTTCCGGCGGCGTAAGAAACCATCGACATCGCCGGGATCGGCGAGAAACCGGCAACGCGATCCGGGCCGTAG‑GTTTTGATGGTCC‑AGA‑CGTTAGCGGCGGCAATCAGCTGGGTTAGTTCC‑TGCCCGGTTGAGCGGGTAAA‑CCCGCCGCGCCCACACCCCTGTTTTTCGCTCCAGCAGtttt                                                                                                                                                                                            >  2:19044/1‑239 (MQ=255)
                                                  gTCGCAATACCAGTC‑GTAGAAACTTAAACAGGTG‑CCGCCAAGCAGCGACAGATAACGCGTTCCGGCGGCGTAAGAAACCATCGACATCGCCGGGATCGGCGAGAAACCGGCAACGCGATCCGGGCCGTAG‑GTTTTGATGGGCC‑AGA‑CGTTAGCGGCGGCAATCAGCTGGTTTAGTTTC‑TGCCAGTTGGAGCGGGTAAA‑CCCGCCCCGCCCACCCACTTGTTTGTAGCTCAGGCACtttt                                                                                                                                                                                            >  2:322363/1‑239 (MQ=255)
                                                       aaaaCCAGTC‑GTAAAAACTAAAACAGGGG‑CCGCCAAGCAGCGACAGATAAGGCGTCCCGGCGGCGGAAGAACCCAGCGACATCGCCGGGATCGGGGAAAAACCGGCAACGCGATCCGGGCCGTAG‑GTTTTGATGGTCC‑AGA‑CGTTAGCGGCGGCAATCAGCTGGTTTAGTTCC‑TGCCAGTTGGAGCGGATAAC‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGa                                                                                                                                                                                       <  2:1311234/236‑1 (MQ=255)
                                                            cAGTC‑GTAGAAACTTAAACAGGTG‑CCGCCAAGCAGCGACAGATAACGCGTTCCGGCGGCGTAAGAAACCATCGACATCGCCGGGATCGGCGAGAAACCGGCAACGCGATCCGGGCCGTAG‑GTTTTGATGGTCC‑AGA‑CGTTAGCGGCGGCAATCCGCCGGTTTAGTTCC‑TGCCAGTTGGAGCGGATAAA‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCCCTTTTGCGGATCGt                                                                                                                                                                                   >  1:738141/1‑238 (MQ=255)
                                                              gTC‑GTAGAAACTTAAACAGGTG‑CCGCCAAGCAGCGACAGATAACGCGTTCCGGCGGCGTAAGAAAACATCGACATCGACGGGATCGGCGAGAAACCGGCAACGCGATCCGGGCCGTAG‑GTTTTGATGGTCC‑AGA‑CGTTAGCGGCGCCAATCAACTGGTTTAGTTCC‑TGCCAGTTGGGGCGGGTAAA‑CCCGCCCCGCCCACGCACTTGTTTGTGGGTCAGTGACCTTTGCGCATCGTTct                                                                                                                                                                                >  2:1296748/1‑238 (MQ=255)
                                                                  gTAGAAACTTAAACAGGTG‑CCGCCAAGCAGCGACAGATAACGCGTTCCGGCGGCGTAAGAAACCATCGACATCGCCGGGATCGGCGAGAAACCGGCAACGCGATCCGGGCCGTAG‑GTTTTGATGGTCC‑AGA‑CGTTAGCGGCGGCAATCAGCTGGTTTAGTTCC‑TGCCAGTTGGAGCGGATAAA‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATaa                                                                                                                                                                             >  2:343825/1‑239 (MQ=255)
                                                                          ttAAACAGGTG‑CCGCCAAGCAGCGACAGATAACGCGTTCCGGCGGCGTAAGAAACCATCGACATCGCCGGGATCGGCGAGAAACCGGCAACGCGATCCGGGCCGTAG‑GTTTTGATGGTCC‑AGA‑CGTTAGCGGCGGCAATCAGCTGGTTTAGGTCC‑TGCCAGTTGGAGCGGATAAA‑CCCGCCGCGCCCACGCCCTTGTTTTTAGATCAGGCACTTTTGCGGGTCGTTCATAAGCGATGc                                                                                                                                                                      >  2:1144667/1‑238 (MQ=255)
                                                                           taaaagaggcccccacAAGCGGCGACAGATACCGGGTGCCGGCGGCGAAAGAAACCATCGACACCGCCGGGATCGGGGAGAACCCGGCACGGCGATCGGGGCGGTGG‑GTTTTGATGGTCC‑AGA‑GGTTGGCGGCGGCAATCGGCTGGTTTGGTGCC‑TGCCGGTGGGGGCGGATAAC‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGcc                                                                                                                                                                     <  2:146311/224‑1 (MQ=255)
                                                                           tAAACAGGGGCCCCCCAAGCCGCCACAGAAACCCGGTTCGGGCGGCGAAAAAAACCAGCGACAGCGCGGGGAGGGGCGAAAAACCGGCAGCGCGACCGGGGCCGGAG‑GTTTGGAGGGCCC‑AGA‑CGTTAGCGGCGGCAATCAGTTGTTTTAGTTCC‑TGCCAGTGGGAGCGGAAAAC‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGcc                                                                                                                                                                     <  2:1322891/239‑1 (MQ=255)
                                                                           aaaaacggggcccgccAAGCAGCGACAGATACCGCGTTCCGGCGGCGAAAGAACCCATCAACATCGCGGGGATCGGGAAGAAGCCGGCAACGCGACCGGGGCCGTGG‑GTTTTGAGGGCCA‑AGC‑CGTAGGGGGGGGAAATCAGCGGTTTTAGTCCC‑TCCCGGTGGGAGGGGATAAC‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGcc                                                                                                                                                                     <  2:98263/228‑1 (MQ=255)
                                                                             aaCAGGTG‑CCGCCAAGCAGCGACAGATAACGCGTTCCGGCGGCGTAAGAAACCATCGACATCGCCGGGATCGGCGAGAAACCGGCAACGCGATCCGGGCCGTAG‑GTTTTGATGGTCC‑AGA‑CGTTAGCGGCGGCAATCAGCTGGTTTAGTTCC‑TGCCAGTTGGAGCGGATAAA‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCGCTTTTGCGGATCGTTCATAATCGATGCCCa                                                                                                                                                                   >  2:645497/1‑238 (MQ=255)
                                                                               cgggTC‑CCGCCAGGGGGCGAAAAAAAACGGGTGCGGGGGGGGTAAGAAACCACCGACCCCGCGGGGTTCGGGGAAAACCCGGAAACGGGATCCGGGCGGTAG‑GTTTGGTGGGCCC‑AGA‑CGTTGGCGGCGGCAACCGGCTGGTTTAGTTCC‑TGCCAGTTGGGGGGGATAAC‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGc                                                                                                                                                                <  2:1340241/237‑1 (MQ=255)
                                                                                cggggcccgcaAGGCGGCGAAAGAAAACGCGTTCCGGCGGCGAAAAAACCCACCGACATCGCGGGGAGGGGCGAAAACCCGGCAACGCGACCGGGGCGGTAG‑TTTTTGAGGGCCC‑AGA‑CGTTGGGGGCGGCAATCAGCTGTTTTAGTTCC‑TGCCAGTGGGGGCGGATAAA‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGc                                                                                                                                                                <  2:646218/233‑1 (MQ=255)
                                                                                  gTG‑CCGCCAAGCAGCGACAGATAACGCGTTCCGGCGGCGTAAGAAACCATCGACATCGCCGGGATCGGCGAGAAACCGGCAACGCGATCCGGGCCGTAG‑GTTTTGATGGTCC‑AGA‑CGTTAGCGGCGGCAATCAGCTGGTTTAGTTCC‑TGCCAGTTGGAGCGGATAAA‑Cccgcc                                                                                                                                                                                                                               >  1:188415/1‑173 (MQ=255)
                                                                                  gTG‑CCGCCAAGCAGCGACAGATAACGCGTTCCGGCGGCGTAAGAAACCATCGACATCGCCGGGATCGGCGAGAAACCGGCAACGCGATCCGGGCCGTAG‑GTTTTGATGGTCC‑AGA‑CGTTAGCGGCGGCAATCAGCTGGTTTAGTTCC‑TGCCAGTTGGAGCGGATAAA‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCCCTTTTGCGGATCGTTCCTAATCGATGCCCACGCCAg                                                                                                                                                             >  2:1306564/1‑239 (MQ=255)
                                                                                  gTG‑CCGCCAAGCAGCGACAGATAACGCGTTCCGGCGGCGTAAGAAACCATCGACATCGCCGGGATCGGCGAGAAACCGGCAACGCGATCCGGGCCGTAG‑GTTTTGATGGTCC‑AGA‑CGTTAGCGGCGGCAATCAGCTGGTTTAGTTCC‑TGCCAGTTGGAGCGGATAAA‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAg                                                                                                                                                             >  1:546914/1‑239 (MQ=255)
                                                                                    g‑ccgccAAGCAGCGACAGATAACGCGTTCCGGCGGCGTAAGAAACCATCGACATCGCCGGGATCGGCGAGAAACCGGCAACGCGATCCGGGCCGTAG‑GTTTTGATGGTCC‑AGA‑CGTTAGCGGCGGCAATCAGCTGGTTTAGTTCC‑TGCCAGTTGGAGCGGATAAA‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAgg                                                                                                                                                            >  2:55631/1‑237 (MQ=255)
                                                                                    g‑ccgccAAGCAGCGACAGATAACGCGTTCCGGCGGCGTAAGAAACCATCGACATCGCCGGGATCGGCGAGAAACCGGCAACGCGATCCGGGCCGTAG‑GTTTTGATGGTCC‑AGA‑CGTTAGCGGCGGCAATCAGCTGGTTTAGTTCC‑TGCCAGTTGGAGCGGATAAA‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGt                                                                                                                                                            >  1:31487/1‑238 (MQ=255)
                                                                                        gccAAGCAGCGACAGATACGGCGTTCCGGCGGCGAAAGAACCCATCGACACCGCCGGGATCGGCGAGAACCCGGCACCGCGATCCGGGCCGTAG‑GTTTGGATGGTCC‑AGA‑CGTTAGCGGCGGCAATCAGCGGGTTTAGTCCC‑TGCCAGTGGGAGCGGATAAC‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCg                                                                                                                                                        <  2:431097/239‑1 (MQ=255)
                                                                                           agccagcGACAGAGAGCGTGTTCGGGGGGCGAAAAAACCCCTCGAATTCGCCGAGACGGGCGAGAAACGGGCACCGCGTTCCGGGCTGAGG‑GTTTTAATGGCCC‑AGA‑CGTTAGCGGCGCCAACCAGCGGTTTTATTTCC‑TGCAATTGGGAGCGGATAAC‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGAt                                                                                                                                                     <  2:219054/236‑1 (MQ=255)
                                                                                             cggcggACAGATACCGGGTCCCGGCGCCGAAAAAAACCACCGCCATCGCGGGGACGGGGGAGAACCCGGCACCGCCACCCGGGCGTTAG‑GTTTTGATGGCCC‑AGA‑CGTTAGCGGCGCCAACCGGCTGTTTTAGTTCC‑TGCCAGTGGGAGCGGATAAA‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGAt                                                                                                                                                     <  2:137203/232‑1 (MQ=255)
                                                                                               agcGACAGATAACGCGTTCGGGCGGCGTAAGAACCCAGCGACAGCGCGGGGATCGGCGAGAACCCGGCAACGCGATCCGGGCGGTAG‑GTTTTGAGGGCCC‑AGA‑CGTTGGCGGCGGCAACCAGCGGGTTTAGTTCC‑TGCCAGTGGGAGCGGATAAA‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGc                                                                                                                                                  <  2:648790/238‑1 (MQ=255)
                                                                                                       aTAACGCGTTCCGGCGGCGTAAGAAACCATCGACATCGCCGGGATCGGCGAGAAACCGGCAACGCGATCCGGGCCGTAG‑GTTTTGATGGTCC‑AGA‑CGTTAGCGGCGGCAATCAGCTGGTTTAGTTCC‑TGCCAGTTGGAGCGGATAAA‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCt                                                                                                                                         >  1:396565/1‑239 (MQ=255)
                                                                                                        tAACGCTTTCCGGCGGGGAAAGAAACCACCGACACCGCCGGGACCGGGAAAAAACCGGCAACGCGTTCGGGGCCGTGG‑TTTTTGATGTCCC‑AGA‑CTTTAGGGGCGGCAATCAGCTGTTTTAGTTCC‑TGCCAGTTGGAGCGGATAAA‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCtt                                                                                                                                        <  2:1120386/239‑1 (MQ=255)
                                                                                                        tAACGCGTTCCGGCGGCGTAAGAAACCATCGACATCGCCGGGATCGGCGAGAAACCGGCAACGCGATCCGGGCCGTAG‑GTTTTGATGGTCC‑AGA‑CGTTAGCGGCGGCAATCAGCTGGTTTAGTTCC‑TGCCAGTTGGAGCGGATAAA‑CCCGCCGCGCCCACGCACTTGTTTGTAGCCCAGGCACTTTTGCGGATCGTTCATAATCGATCCCCACGCCAGTACTGGATCGCTGTGTTTCtt                                                                                                                                        >  2:479709/1‑239 (MQ=255)
                                                                                                           cgcgTTCCGGCGGCGTAAGAAACCATCGACATCGCCGGGATCGGCGAGAAACCGGCAACGCGATCCGGGCCGTAG‑GTTTTGATGGTCC‑AGA‑CGTTAGCGGCGGCAATCAGCTGGTTTAGTTCC‑TGCCAGTTGGAGCGGATAAA‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTg                                                                                                                                           >  1:137203/1‑233 (MQ=255)
                                                                                                           cgcgTTCCGGCGGCGTAAGAAACCATCGACATCGCCGGGATCGGCGAGAAACCGGCAACGCGATCCGGGCCGTAG‑GTTTTGATGGTCC‑AGA‑CGTTAGCGGCGGCAATCAGCTGGTTTAGTTCC‑TGCCAGTTGGAGCGGATAAA‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAg                                                                                                                                     >  1:84913/1‑239 (MQ=255)
                                                                                                                tCCGGCGGCGTAAGAAACCATCGACATCGCCGGGATCGGCGAGAAACCGGCAACGCGATCCGGGCCGTAG‑GTTTTGATGGTCC‑AGA‑CGTTAGCGGCGGCAATCAGCTGGTTTAGTTCC‑TGCCAGTTGGAGCGGATAAA‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGTGCt                                                                                                                                 >  1:1349664/1‑238 (MQ=255)
                                                                                                                 ccGGCGGCGTAAGAAACCATCGACATCGCCGGGATCGGCGAGAAACCGGCAACGCGATCCGGGCCGTAG‑GTTTTGATGGTCC‑AGA‑CGTTAGCGGCGGCAATCAGCTGGTTTAGTTCC‑TGCCAGTTGGAGCGGATAAA‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCttg                                                                                                                               >  1:984061/1‑238 (MQ=255)
                                                                                                                 ccGGCGGCGTAAGAAACCATCGACATCGCCGGGATCGGCGAGAAACCGGCAACGCGATCCGGGCCGTAG‑GTTTTGATGGTCC‑AGA‑CGTTAGCGGCGGCAATCAGCTGGTTTAGTTCC‑TGCCAGTTGGAGCGGATAAA‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCTTc                                                                                                                               >  2:142840/1‑239 (MQ=255)
                                                                                                                 ccGGCGGCGTAAGAAACCATCGACATCGCCGGGATCGGCGAGAAACCGGCAACGCGATCCGGGCCGTAG‑GTTTTGATGGTCC‑AGA‑CGTTAGCGGCGGCAATCAGCTGGTTTAGTTCC‑TGCCAGTTGGAGCGGATAAA‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGGGTTGCTTGAGGGCTTc                                                                                                                               >  1:1338453/1‑239 (MQ=255)
                                                                                                                 ccGGCGGCGTAAGAAACCATCGACATCGCCGGGATCGGCGAGAAACCGGCAACGCGATCCGGGCCGTAG‑GTTTTGATGGTCC‑AGA‑CGTTAGCGGCGGCAATCAGCTGGCTTAGTTTC‑TGCCAGTTGGAGCGGATAAA‑CCCGCCGCGCCCCCGCACCTGTTTGTAGCTCCGGCACTTTTGCGGATCGTTCATAATAGATGCCCACGCCAGTCCCGGATCGCTGTGTCGCTTgtagtcctc                                                                                                                               >  2:988427/1‑231 (MQ=255)
                                                                                                                  ggggggcgTAAGAACCCAGCGCCACCGCCGGGAGCGGCGAGAAACCGGAAACGCGACCCGGGCGGGAG‑GTTTGGAGGGTCC‑AGA‑CTTTAGCGGCGGCAATCAGCTGGTTTAGTTCC‑TGCAAGTGGGGGCGGATAAC‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCTTcg                                                                                                                              <  2:251052/235‑1 (MQ=255)
                                                                                                                      ggggaaaaaaaccaacgcactgcccGGGAGCGGCGAAAACCCGCCAACGCTATCGGGGCCGTAG‑TTTTTGTGGGCCC‑AGC‑GGTTGGCGGCGCCACTCGCCTGGTTTAGTTCG‑TGCAAGTTGGAGCGGAAAAA‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCTTCgcgc                                                                                                                           <  2:1092108/215‑1 (MQ=255)
                                                                                                                      ggcgtaggaaccagcggACTTCGCCGGGAGGGGGGAGAAACCGGCAAGCCGACCCGGGCCGTGG‑GTTTTGATGGCCC‑AGC‑CGTTGGCGGCGGCAATCAGCTGTTTTAGTCCC‑TGCAAGTGGGGGCGGATAAC‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCTTCgcgc                                                                                                                           <  2:203740/222‑1 (MQ=255)
                                                                                                                      ggcaaaaGAAACCAGGGACGGCGCCGGGACCGCGGAAAACCGGGAAACGCGA‑CCGGGCCGTTG‑GGTTTGGGGGGCCAAGC‑CGTTGGGGGGGCCAATCAGCGGTTTGATTCCC‑GGCCAGTGGGGGCGGAAAAC‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCTTCgcgc                                                                                                                           <  2:712052/233‑1 (MQ=255)
                                                                                                                        gcctaaaaacccctcgactccgcGGGGTGGGGGAAAAACCCGAAACCCGGTCCCGGGCGTTGTGTTTGGTTGGCCC‑AAA‑CGTTAGCGGCGGCAACCAGCTGTTTTAGTCCC‑TGCCATTGGGGGCGGATAAC‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCTTCGCGcc                                                                                                                          <  2:373085/216‑1 (MQ=255)
                                                                                                                          tAAGAAACCATCGACATCGCCGGGATCGGCGAGAAACCGGCAACGCGATCCGGGCCGTAG‑GTTTTGATGGTCC‑AGA‑CGTTAGCGGCGGCAATCAGCTGGTTTAGTTCC‑TGCCAGTTGGAGCGGATAAA‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCTTCGCGCCACAg                                                                                                                      >  1:1311333/1‑239 (MQ=255)
                                                                                                                           aaGAAACCACCGACATGGCGGGGACCGGCGGAAAACGGGCACCGGGATCGGGGCCGTAG‑GTTTGAATGTCCC‑AGC‑CGTTAGGGGGGCAAATCGGCTGGTTTAGTCCC‑TGCCAGTTGGGGGGGATAAC‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCTTCGCGCCACAGt                                                                                                                     <  2:1301997/239‑1 (MQ=255)
                                                                                                                                    tCGACATCGCCGGGATCGGCGAGAAACCGGCAACGCGATCCGGGCCGTAG‑GTTTTGATGGTCC‑AGA‑CGTTAGCGGCGGCAATCAGCTGGTTTAGTTCC‑TGCCAGTTGGAGCGGATAAA‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCTTCGCGCCCCAGTTCAATCAg                                                                                                             >  1:701024/1‑238 (MQ=255)
                                                                                                                                         aTCGCCGGGATCGGCGAGAAACCGGCAACGCGATCCGGGCCGTAG‑GTTTTGATGGTCC‑AGA‑CGTTAGCGGCGGCAATCAGCTGGTTTAGTTCC‑TGCCAGTTGGAGCGGATAAA‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCTTCGCGCCACAGTTCAATCAGTCGttt                                                                                                       >  2:1084244/1‑239 (MQ=255)
                                                                                                                                                  aggggggagaAACGGGCACGGGGACCGGGGCGGAGG‑GTTTTGAGGGCCA‑AGC‑GGTAAGCGGGGGAAACCGGCTGGTTAAGTTCC‑GGCCGGTGGGGGCGAATAAC‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCTTCGCGCCACAGTTCAATCAGTCGTTTACGAATGa                                                                                               <  2:84913/232‑1 (MQ=255)
                                                                                                                                                   gggggaaaaaaCCGGCACGGCGACCCGGGCGGTAG‑GTTTGGATGGTCA‑AGA‑GGTAAGGGGCGGAAAGCAGCGGGTTTATTTCC‑TGCAAGTTGGGGCGGATAAA‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCTTCGCGCCACAGTTCAATCAGTCGTTTACGAATGAg                                                                                              <  2:738141/230‑1 (MQ=255)
                                                                                                                                                      gggaaaaaCCGGCGCGGCGAGCGGGGCGGTGG‑TTTTTGTGGGTCA‑AGA‑CGTGAGCGGCGGCAATCGGTGGGTTTAGTTCC‑GGCCAGTGGGGGCGGATAAC‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCTTCGCGCCACAGTTCAATCAGTCGTTTACGAATGAGCgg                                                                                           <  2:31487/233‑1 (MQ=255)
                                                                                                                                                        gagaACCCGGCACGGGGATCGGGGCCGTAG‑GTTTTGATGG‑GC‑AGACCGTTGGCGGCGGCAATCGGCTGGTTTAGTTCC‑TGCCAGTTGGGGCGAATAAC‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCTTCGCGCCACAGTTCAATCAGTCGTTTACGAATGAGCGGGt                                                                                         <  2:701024/238‑1 (MQ=255)
                                                                                                                                                                   aaggcgATCCGGGCCGTGG‑TTTTTGATGGTCC‑AGA‑CGTGAGCGGCGGCAATCAGCTGGTTTAGTTCC‑GGCCAGTTGGAGCGGATAAA‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCTTCGCGCCACAGTTCAATCAGTCGTTTACGAATGAGCGGGTATTTCAGGCGGt                                                                             <  1:1378543/236‑1 (MQ=255)
                                                                                                                                                                      gcgATCCGGGCCGTAG‑GTTTTGATGGTCC‑GAA‑CGTTAGCGGCGGCAATAAGTTGTTTTAGTTCC‑TGCCAGTTGGAGCGGATAAC‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCTTCGCGCCACAGTTCAATCAGTCGTTTACGAATGAGCGGGTATTTCAGGCGGt                                                                             <  1:1262320/236‑1 (MQ=255)
                                                                                                                                                                      cgacgctgggCGGTGG‑GTTTGGAGGGCCC‑AGA‑CGTTGGCGCCGGCAACCAGC‑GGTTTTGTTCC‑CGGCCGTTGGGGGCGATAAACCCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCTTCGCGCCACAGTTCAATCAGTCGTTTACGAATGAGCGGGTATTTCAGGCGGTTAg                                                                          <  2:190364/232‑1 (MQ=255)
                                                                                                                                                                         gcccGGGGTGTGG‑GTTTGGGTGGTCC‑AGG‑CTTGAGCGGGGGCAATCAGGTGTTTTATTCCC‑TGCCAGTTGGAGGGAATAAC‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCTTCGCGCCACAGTTCAATCAGTCGTTTACGAATGAGCGGGTATTTCAGGCGGTTAgcg                                                                        <  2:3630/236‑1 (MQ=255)
                                                                                                                                                                         acccGGGCCGTAG‑TTTTTGAGGGCCC‑AGG‑CGTTAGCGGGGGCAATCAGCTGGTTTATTTCC‑TGCCAGTGGGGGCGGATAAA‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCTTCGCGCCACAGTTCAATCAGTCGTTTACGAATGAGCGGGTATTTCAGGCGGTTAgcg                                                                        <  1:263366/236‑1 (MQ=255)
                                                                                                                                                                          tCCGGGCCTTAG‑GTTTTGATGGTCC‑AGC‑CTTTAGCGGCGGCAATCGGTGGGTTTAGTTCC‑TGCCAGTTGGAGCGAATAAA‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCTTCGCGCCACAGTTCAATCAGTCGTTTACGAATGAGCGGGTATTTCAGGCGGTTAGCGCt                                                                      <  1:1115534/239‑1 (MQ=255)
                                                                                                                                                                          tCCGGGCCGTAG‑GTTTTGATGGTCC‑AGA‑CGTTAGCGGCGGCAATCAGCTGGTTTAGTTCC‑TGCCAGTTGGAGCGGATAAA‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCTTCGCGCCACAGTTCAATCAGTCGTTTACGAATGAGCGGGTATTTCAGGCGGTTAgc                                                                         >  1:180481/1‑236 (MQ=255)
                                                                                                                                                                                    cg‑gTTTTGAGGGTCC‑GGC‑CGTTACCGCCGCCAACGAGTGGTTTTAGTCCT‑TGCCAGTTGGAGCGGATAAA‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCTTCGCGCCACAGTTCAATCAGTCGTTTACGAATGAGCGGGTATTTCAGGCGGTTAGCGCTGTAAAGATAc                                                            <  1:370743/238‑1 (MQ=255)
                                                                                                                                                                                    aG‑GTTTTGATGGTCC‑AGA‑CGTTAGCGGCGGCAATCAGCTGGTTTAGTTCC‑TGCCAGTTGGAGCGGATAAA‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCTTCGCGCCACAGTTCAATCAGTCGTTTACGAATGAGCGGGTATTTCAGGCGGTTAGCGCTGTAAAGATa                                                             >  2:485192/1‑238 (MQ=255)
                                                                                                                                                                                      gttttgGATGGCCC‑AGG‑CGTTAGGGGCGGCAATCAGGTGTTTTAGTTCC‑GGCCAGTTGGGGCGAATAAC‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCTTCGCGCCACAGTTCAATCAGTCGTTTACGAATGAGCGGGTATTTCAGGCGGTTAGCGCTGTAAAGATAc                                                            <  2:1311333/236‑1 (MQ=255)
                                                                                                                                                                                          ttGATGGTCC‑AGA‑CGTTAGCGGCGGCAATCAGCTGGTTTAGTTCC‑TGCCAGTTGGAGCGGATAAA‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGa                                                                                                                                                      >  1:985473/1‑144 (MQ=255)
                                                                                                                                                                                          ttGATGGTCC‑AGA‑CGTTAGCGGCGGCAATCAGCTGGTTTAGTTCC‑TGCCAGTTGGAGCGGATAAA‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCTTCGCGCCACAGTTCAATCAGTCGTTTACGAATGAGCGGGTATTTCAGGCGGTTAGCGCTGTAAAGATACCAGGa                                                       >  2:1017524/1‑239 (MQ=255)
                                                                                                                                                                                          ttGATGGTCC‑AGA‑CGTTAGCGGCGGCAATCAGCTGGTTTAGTTCC‑TGCCAGTTGGAGCGGATAAA‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCTTCGCGCCACAGTTCAATCAGTCGTTTACGAATGAGCGGGTATTTCAGGCGGTTAGCGCTGTAAAGATACCAGGa                                                       >  1:70935/1‑239 (MQ=255)
                                                                                                                                                                                          ttGAGGTTCC‑AGA‑CGTTGGCGGCGGCAACCAGGTGTTTTAGTTCC‑GGCCAGTTGGGGCGGATAAA‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCTTCGCGCCACAGTGCAATCAGTCGTTTACGAATGAGCGGGTATTTCAGGCGGTTAGCGCTGTAAAGATACCAGGa                                                       <  1:479709/239‑1 (MQ=255)
                                                                                                                                                                                                gTCC‑AGA‑CGTTAGCGGCGGCAATCAGCTGGTTTAGTTCC‑TGCCAGTTGGAGCGGATAAA‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCTTCGCGCCACAGTTCAATCAGTCGTTTACGAATGAGCGGGTATTTCAGGCGGTTAGCGCTGTAAAGATACCAGGAGTAACt                                                 >  2:789161/1‑239 (MQ=255)
                                                                                                                                                                                                       a‑CGTTAGCGGCGGCAATCAGCTGGTTTAGTTCC‑TGCCAGTTGGAGCGGATAAA‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCTTCGCGCCACAGTTCAATCAGTCGTTTACGAATGAGCGGGTATTTCAGGCGGTTAGCGCTGTAAAGATACCAGGAGTAACTTGCGcc                                           >  1:1213369/1‑239 (MQ=255)
                                                                                                                                                                                                             aGCGGCGGCAATCAGCTGGTTTAGTTCC‑TGCCAGTTGGAGCGGATAAA‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTTCTTGAGGGCTTAGCGCCACCGTTCAATCAGTCGTTTCCGACTGAGCGGGTATTTCAGGCGGTTAGCGGTGTAAAGATACCAGGAGGCACTTGCGGCAcgcg                                      >  1:1397183/1‑239 (MQ=255)
                                                                                                                                                                                                             aGCGGCGGCAATCAGCCGGTTTAGTTCC‑TGCCAGCTGGAGCGGATAAA‑CCCGCCGCGCCCACGCACTTGCTTGTAGCTCAGGCACTTTTTCGGATCGTTCATAATCGATGCCCACGCCAGCACCGGGTCGCTGTGTTGCTTGAGGGCTTCGCGCCACAGCTCAATCCGCCGTTTACGAATGCGCGGGTATTTCAGGCGGTTAGCGACGTGAAGATACACGGAGGAACTTgctcgacccc                                      >  2:916516/1‑231 (MQ=255)
                                                                                                                                                                                                                    gcAATCAGCTGGTTTAGTTCC‑TGCCAGTTGGAGCGGATAAA‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCTTCGCGCCACAGTTCAATCAGTCGTTTACGAATGAGCGGGTATTTCAGGCGGTTAGCGCTGTAAAGATACCAGGAGTAACTTGCGCCACGCGGGCAGcc                               >  1:1047331/1‑239 (MQ=255)
                                                                                                                                                                                                                     aaaTCAGCGGGTTTTGTTCC‑GGCCAGTTGGAGCGGATAAC‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCTTCGCGCCACAGTTCAATCAGTCGTTTACGAATGAGCGGGTATTTCAGGCGGTTAGCGCTGTAAAGATACCAGGAGTAACTTGCGCCACGCGGGCAGCcg                              <  1:887595/238‑1 (MQ=255)
                                                                                                                                                                                                                       ataggttGGTTTAGTCCC‑GGCCGGTGGGGGCGGATAAC‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCTTCGCGCCACAGTTCAATCAGTCGTTTACGAATGAGCGGGTATTTCAGGCGGTTAGCGCTGTAAAGATACCAGGAGTAACTTGCGCCACGCGGGCAGCcgcg                            <  2:984061/233‑1 (MQ=255)
                                                                                                                                                                                                                       agcAGCTGGTTAAGTGCC‑TGCCATTTGGAGCGGATAAA‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCTTCGCGCCACAGTTCAATCAGTCGTTTACGAATGAGCGGGTATTTCAGGCGGTTAGCGCTGTAAAGATACCAGGAGTAACTTGCGCCACGCGGGCAGCcgcg                            <  1:825953/237‑1 (MQ=255)
                                                                                                                                                                                                                       accggtgggTTTATTTCG‑TGCCGTTGGGAGGGAATAAC‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCTTCGCGCCACAGTTCAATCAGTCGTTTACGAATGAGCGGGTATTTCAGGCGGTTAGCGCTGTAAAGATACCAGGAGTAACTTGCGCCACGCGGGCAGCcgcg                            <  2:180481/232‑1 (MQ=255)
                                                                                                                                                                                                                       aTCAGCTGGTTTGGTGCC‑GGCGAGTTGGGGCGGAAACA‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCTTCGCGCCACAGTTCAATCAGTCGTTTACGAATGAGCGGGTATTTCAGGCGGTTAGCGCTGTAAAGATACCAGGAGTAACTTGCGCCACGCGGGCAGCcgcg                            <  2:496662/239‑1 (MQ=255)
                                                                                                                                                                                                                       aTCAGCTGGGTTATTCCG‑TGCCAGTTGGTGGGGATAAT‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCTTCGCGCCACAGTGCAATCAGTCGTTTACGAATGAGCGGGTATTTCAGGCGGTTAGCGCTGTAAAGATACCAGGAGTAACTTGCGCCACGCGGGCAGCcgcg                            <  1:117745/239‑1 (MQ=255)
                                                                                                                                                                                                                       aTCAGCTGCTTTAGTTCG‑TGCCAGTTGGAGCGAATAAT‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCTTCGCGCCACAGTTCAATCAGTCGTTTACGAATGAGCGGGTATTTCAGGCGGTTAGCGCTGTAAAGATACCAGGAGTAACTTGCGCCACGCGGGCAGCcgcg                            <  1:197524/239‑1 (MQ=255)
                                                                                                                                                                                                                                tttAGTTCC‑TGCCTGTTGGAGCGGATAAA‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCTTCGCGCCACAGTTCAATCAGTCGTTTACGAATGAGCGGGTATTTCAGGCGGTTAGCGCTGTAAAGATACCAGGAGTAACTTGCGCCACGCGGGCAGCCGCGAGGTTCAt                    <  1:343825/238‑1 (MQ=255)
                                                                                                                                                                                                                                 ttAGTCCC‑TGCCAGGTGGGCGGGAAAAC‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCTTCGCGCCACAGTTCAATCAGTCGTTTACGAATGAGCGGGTATTTCAGGCGGTTAGCGCTGTAAAGATACCAGGAGTAACTTGCGCCACGCGGGCAGCCGCGAGGTTCATGa                  <  2:936071/239‑1 (MQ=255)
                                                                                                                                                                                                                                  ttggttctggcgggTTGGGGCGGATAAC‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCTTCGCGCCACAGTTCAATCAGTCGTTTACGAATGAGCGGGTATTTCAGGCGGTTAGCGCTGTAAAGATACCAGGAGTAACTTGCGCCACGCGGGCAGCCGCGAGGTTCATGa                  <  2:868654/226‑1 (MQ=255)
                                                                                                                                                                                                                                      tCC‑TGCGAGTTGGGGCGAATAAC‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCATTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCTTCGCGCCACAGTTCAATCAGTCGTTTACGAATGAGCGGGTATTTCAGGCGGTTAGCGCTGTAAAGATACCAGGAGTAACTTGCGCCACGCGGGCAGCCGCGAGGTTCATGATTggg             <  1:1306564/239‑1 (MQ=255)
                                                                                                                                                                                                                                             cAGTGGGGGCGGATAAC‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCTTCGCGCCACAGTTCAATCAGTCGTTTACGAATGAGCGGGTATTTCAGGCGGTTAGCGCTGTAAAGATACCAGGAGTAACTTGCGCCACGCGGGCAGCCGCGAGGTTCATGATTGGGCAGGt        <  1:322363/238‑1 (MQ=255)
                                                                                                                                                                                                                                             cAGTGGGGGCGAATAAA‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCTTCGCGCCACAGTTCAATCAGTCGTTTACGAATGAGCGGGTATTTCAGGCGGTTAGCGCTGTAAAGATACCAGGAGTAACTTGCGCCACGCGGGCAGCCGCGAGGTTCATGATTGGGCAGGt        <  1:823504/238‑1 (MQ=255)
                                                                                                                                                                                                                                                ttGGAGCGGATAAA‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCTTCGCGCCACAGTTCAATCAGTCGTTTACGAATGAGCGGGTATTTCAGGCGGTTAGCGCTGTAAAGATACCAGGAGTAACTTGCGCCACGCGGGCAGCCGCGAGGTTCATGATTGGGCAGGTCAgg    >  1:141076/1‑239 (MQ=255)
                                                                                                                                                                                                                                                   gAGCGGATAAC‑CCCGCCGCGCCCACTCATTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCTTCGCGCCACAGTTCAATCAGTCGTTTACGAATGAGCGGGTATTTCAGGCGGTTAGCGCTGTAAAGATACCAGGAGTAACTTGCGCCACGCGGGCAGCCGCGAGGTTCATGATTGGGCAGGTCAGGGc  <  1:909160/238‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                             
TCCGGTACGTCGGTTTGCTCGCCCCAGGTCATCGGCGAGGCGGGCGGCAGGTCGCAATACCAGTC‑GTAGAAACTTAAACAGGTG‑CCGCCAAGCAGCGACAGATAACGCGTTCCGGCGGCGTAAGAAACCATCGACATCGCCGGGATCGGCGAGAAACCGGCAACGCGATCCGGGCCGTAG‑GTTTTGATGGTCC‑AGA‑CGTTAGCGGCGGCAATCAGCTGGTTTAGTTCC‑TGCCAGTTGGAGCGGATAAA‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCTTCGCGCCACAGTTCAATCAGTCGTTTACGAATGAGCGGGTATTTCAGGCGGTTAGCGCTGTAAAGATACCAGGAGTAACTTGCGCCACGCGGGCAGCCGCGAGGTTCATGATTGGGCAGGTCAGGGC  >  NC_000913/1541872‑1542346

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