Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,568,770 G→T 17.0% N186K (AAC→AAA)  glgC ← glucose‑1‑phosphate adenylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,568,7700GT17.0% 170.7 / 6.6 89N186K (AAC→AAA) glgCglucose‑1‑phosphate adenylyltransferase
Reads supporting (aligned to +/- strand):  ref base G (54/17);  new base T (2/14);  total (58/31)
Fisher's exact test for biased strand distribution p-value = 3.59e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.85e-01

TACGCAAGAGAGCGGGAACGGGTGCGCATAGGCCAGACCGGCTTCGGTGATCTTGGGAATCAAATCTTTGCCAAAGTCGTGGCTGGAGTTCTCATCGCGATCGTCTTCTTCCAGCAGTTCATACAGATAGTCGGCGTCAAAGACGTAGATACCCATA‑CTCGCCAGAGATTTGCTCGGATCGTTCGGCATTGACGGCGGGTTAGCAGGTTTTTC‑AA‑CGAA‑TTCGATAATTTTATCGTTCTCATCAACCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCATACAAGCAACGGTGCAACGTGCGCCTTTTTCGACGTGATCGATAAGCATACGCGAGTAGTCTTGCTTGTAGATATGGTCGCCCGCCAGGATCACCACGTATTCCGCTTTATAACGGCGGATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCCGCGATACCAGTTTTCCCCTTTCATT  >  NC_000913/3568536‑3569007
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                             
tACGCAAGAGAGCGGGAACGGGTGCGCATAGGCCAGACCGGCTTCGGTGATCTTGGGAATCAAATCTTTGCCAAAGTCGTGGCTGGAGTTCTCATCGCGATCGTCTTCTTCCAGCAGTTCATACAGATAGTCGGCGTCAAAGACGTAGATACCCATA‑CTCGCCAGAGATTTGCTCGGATCGTTCGGCATTGACGGCGGGTTAGCAGGGTTTTC‑AA‑CGAA‑TTCGATAATTTTATCGttt                                                                                                                                                                                                                                            >  2:1114953/1‑237 (MQ=255)
 aCGCAAGAGAGCGGGAACGGGTGCGCATAGGCCAGACCGGCTTCGGTGATCTTGGGAATCAAATCTTTGCCAAAGTCGTGGCTGGAGTTCTCATCGCGATCGTCTTCTTCCAGCAGTTCATACAGATAGTCGGCGTCAAAGACGTAGATACCCATA‑CTCGCCAGAGATTTGCTCGGATCGTTCGGCATTGACGGCGGGTTAGCAGGTTTTTT‑AA‑CGAA‑TTCGATAATTTTATCGTtc                                                                                                                                                                                                                                            >  1:367435/1‑237 (MQ=255)
  cGCAAGAGAGCGGGAACGGGTGCGCATAGGCCAGACCGGCTTCGGTGATCTTGGGAATCAAATCTTTGCCAAAGTCGTGGCTGGAGTTCTCATCGCGATCGTCTTCTTCCAGCAGTTCATACAGATAGTCGGCGTCAAAGACGTAGATACCCATA‑CTCGCCAGAGATTTGCTCGGATCGTTCGGCATTGACGGCGGGTTAGCAGGTTTTTC‑AA‑CGAA‑TTCGATAATTTTATCGTTCtca                                                                                                                                                                                                                                         >  2:1049388/1‑239 (MQ=255)
          agCGGGAACGGGTGCGCATAGGCCAGACCGGCTTCGGTGATCTTGGGAATCAAATCTTTGCCAAAGTCGTGGCTGGAGTTCTCATCGCGATCGTCTTCTTCCAGCAGTTCATACAGATAGTCGGCGTCAAAGACGTAGATACCCATA‑CTCGCCAGAGATTTGCTCGGATCGTTCGGCATTGACGGCGGGTTAGCAGGGTTTTC‑AA‑CGAA‑TTCGATAATTTTATCGTTCTCATCAAcca                                                                                                                                                                                                                                  >  2:485034/1‑237 (MQ=255)
          agCGGGAACGGGTGCGCATAGGCCAGACCGGCTTCGGTGATCTTGGGAATCAAATCTTTGCCAAAGTCGTGGCTGGAGTTCTCATCGCGATCGTCTTCTTCCAGCAGTTCATACAGATAGTCGGCGTCAAAGACGTAGATACCCATA‑CTCGCCAGAGATTTGCTCGGATCGTTCGGCATTGACGGCGGGGTAGCAGGTTTTTC‑AA‑CGGA‑TTCGATAATTTTATCGTTCTCATCAACCg                                                                                                                                                                                                                                  >  2:995100/1‑238 (MQ=255)
            cGGGAACGGGTGCGCATAGGCCAGACCGGCTTCGGTGATCTTGGGAATCAAATCTTTGCCAAAGTCGTGGCTGGAGTTCTCATCGCGATCGTCTTCTTCCAGCAGTTCATACAGATAGTCGGCGTCAAAGACGTAGATACCCATA‑CTCGCCAGAGATTTGCTCGGATCGTTCGGCATTGACGGCGGGTTAGCAGGTTTTTC‑AA‑CGAA‑TTCGATAATTTTATCGTTCTCATCAACCGCCa                                                                                                                                                                                                                               >  2:233004/1‑239 (MQ=255)
               gAACGGGTGCGCATCGGCCAGACCGGCTTCGGTGATCTTGGGAATCAAATCTTTGCCAAAGTCGTGGCTGGAGTTCTCATCGCGATCGTCTTCTTCCAGCAGTTCATACAGATAGTCGGCGTCAAAGACGTAGATACCCATA‑CTCGCCAGAGATTTGCTCGGATCGTTCGGCATTGACGGCGGGTTAGCAGGTTTTTT‑AA‑GGCA‑TTCGATAATTTTATCGTTCTTATCaaacacccat                                                                                                                                                                                                                             >  1:91866/1‑230 (MQ=255)
               gAACGGGTGCGCATAGGCCAGACCGGCTTCGGTGATCTTGGGAATCAAATCTTTGCCAAAGTCGTGGCTGGAGTTCTCATCGCGATCGTCTTCTTCCAGCAGTTCATACAGATAGTCGGCGTCAAAGACGTAGATACCCATA‑CTCGCCAGAGATTTGTTCGGATCGTTCGGCATTGACGGCGGGTTGACAGGGTTTTC‑AA‑CCGA‑TTCGAAAATTTTTTCCTTCTCATCAAACGcccat                                                                                                                                                                                                                             >  1:476697/1‑235 (MQ=255)
               gAACGGGTGCGCATAGGCCAGACCGGCTTCGGTGATCTTGGGAATCAAATCTTTGCCAAAGTCGTGGCTGGAGTTCTCATCGCGATCGTCTTCTTCCAGCAGTTCATACAGATAGTCGGCGTCAAAGACGTAGATACCCATA‑CTCGCCAGAGATTTGCTCGGATCGTTCGGCATTGACGGCGGGTTAGCAGGGTTTTC‑AA‑CGGA‑TTCGATAATTTTATCGTTCTCATCAACCGCCATa                                                                                                                                                                                                                             >  2:309691/1‑238 (MQ=255)
               gAACGGGTGCGCATAGGCCAGACCGGCTTCGGTGATCTTGGGAATCAAATCATTGCCACAGTCGTGGCTGGCGTTCTCATCGCGATCGTCTTCTTCCCGCAGTTCATACAGATCGCCGGCGTCAAAGCCGTAGATACCCATA‑CTCGCCCGATATTTGCTCGGATCGTTCGGCCTTGACGGCGGTTTAGCCAGTTTTTC‑AT‑CCTG‑TTCGATACTTTTTATTTTCTCATCCACCCCCCTa                                                                                                                                                                                                                             >  1:290104/1‑238 (MQ=255)
                    ggTGCGCATAGGCCAGACCGGCTTCGGTGATCTTGGGAATCAAATCTTTGCCAAAGTCGTGGCTGGAGTTCTCATCGCGATCGTCTTCTTCCAGCAGTTCATACAGATAGTCGGCGTCAAAGACGTAGATACCCATA‑CTCGCCAGAGATTTGCTCGGATCGTTCGGCATTGACGGCGGGTTAGCAGGTTTTTC‑AA‑CGAA‑TTCGGTAATTTTATCCTTCTTATCAACCGCCATAAAGCCa                                                                                                                                                                                                                       >  1:700883/1‑239 (MQ=255)
                    ggTGCGCATAGGCCAGACCGGCTTCGGCGATCTTGGGAATCAAATCTTTGCCAAAGTCGTGGCTGGAGTTCTCATCGCGATCGTCTTCTTCCAGCAGTTCATACAGATAGTCGGCGTCAAAGACGTAGATACCCATA‑CTCGCCAGAGATTTGCTCGGATCGTTCGGCATTGACGGCGGGTTGGCGGGGTTTTC‑AA‑CGAC‑TTTGCTCATTTTTTCGTTCTCGTTCAACGCCATacaccca                                                                                                                                                                                                                       >  1:262286/1‑233 (MQ=255)
                           aTAGGCCAGCCGGGTTTGGGTGATTTTGGGAATCAAATCTTTGCCAAAGTGGTGGCTGGAGTTCTCATCGCGATCTTCTTCTTCCAGCATTTCATACAGATAGTCGGCGTCAAAGACGTAGATCCCCATA‑CCCGCCAGAGATTTGCTCGGATCGTTCGGCATTGACGGCGGGTTAGCAGGTTTTTC‑AA‑CAAA‑TTCGATAATTTTATCTTTCTCATCAACCGCCATAACGCCAAATGCg                                                                                                                                                                                                                 <  2:570219/238‑1 (MQ=255)
                             aGGCCAGACCGGCTTCGGTGATCTTGGGAATCAAATCTTTGCCAAAGTCGTGGCTGGAGTTCTCATCGCGATCGTCTTCTTCCAGCAGTTCATACAGATAGTCGGCGTCAAAGACGTAGATACCCATA‑CTCGCCAGAGATTTGCTCGGATCGTTCGGCATTGACGGCGGGTTAGCAGGTTTTTC‑AA‑CGAA‑TTCGATAATTTTATCGGTCTCATCAACCGCCATAACGCCAAATGCGGAg                                                                                                                                                                                                              >  2:337468/1‑239 (MQ=255)
                               gCCAGACCGGCTTCGGTGATCTTGGGAATCAAATCTTTGCCAAAGTCGTGGCTGGAGTTCTCATCGCGATCGTCTTCTTCCAGCAGTTCATACAGATAGTCGGCGTCAAAGACGTAGATACCCATA‑CTCGCCAGAGATTTGCTCGGATCGTTCGGCATTGACGGCGGGTTAGCAGGTTTTTC‑AA‑CGAA‑TTCGATAATTTTATCGTTCTCATCAACCGCCATAACGCCAAATGCGGAg                                                                                                                                                                                                              >  2:66472/1‑237 (MQ=255)
                               gCCAGACCGGCTTCGGTGATCTTGGGAATCAAATCTTTGCCAAAGTCGTGGCTGGAGTTCTCATCGCGATCGTCTTCTTCCAGCAGTTCATACAGATAGTCGGCGTCAAAGACGTAGATACCCATA‑CTCGCCAGAGATTTGCTCGGATCGTTCGGCATTGACGGCGGGTTAGCAGGGTTTTC‑AA‑CGAG‑TTCGATAATTTTATCGTTCTCATCAAACGCCATAAAGCCAACTGCGGAggg                                                                                                                                                                                                            >  1:798668/1‑238 (MQ=255)
                                    aCCGGCTTGGGTGCTCTTGGGAATCAAATCTTTGACAAAGGCGCGGCTGGAGATCCCATCCGGATCGTCATCTTCCAGCAGTTCATACAGATAGTCGGCGGCAAAGACGTATATACCCATA‑CTCGCCAAAGATCTGCTAGGATCGTTCGGCATTTACGGAGGACTACCGGGATTTTC‑TA‑CGCA‑TATGATGATTTTAACTTTCTTATCAACCGCATTAAAACCCAATGcaggggatgtc                                                                                                                                                                                                        >  1:122031/1‑228 (MQ=255)
                                      cgtgtccGGTGATCTGGGGAATCAATTCTTTGCCAAGGTGGGGGCGGGAGTTCCCATCGCGAGCGTCTTCTCCAAGCAGTTCATACAGAGGGTCGGAGTAAAACACGTAGACACCCATC‑CTCGCTAAAATTTTCCGCGGATCTTTCGGCATGGGCGG‑GGGTTTGCCAGTTTTTCAAA‑CAAA‑TTGAATATTTTTTTCGTTCGCTCCACCCGCCATAACGCCAAATGCGGAGGCTTCTTCa                                                                                                                                                                                                     <  1:939480/233‑1 (MQ=255)
                                         tctcgtgtttgttGGAAATCAATTTTTTCCAAAAGTCGGGGCGGGAGTCCTCATCGCGACCGTCTCCTTCCAGCGGTTAATAAAGATGGGCGGCTTCAAAGACGTAGATACCCATA‑CTCGCCAGAGTTTTGCTCGTATCGTCCGGCATTGACGGCGGGTTAGCAGGTTTTAC‑AA‑CGAA‑TCCGATATTTTTATCTTTCTCCACAACCGCCATAACGCCAAATGCGGAGGCTTCTTCAAt                                                                                                                                                                                                   <  1:1049388/227‑1 (MQ=255)
                                         cttgcggtttttgggcacacAAATCTTGGCCAAAGTCTGGGCTGGGGTTCCCTGCGAGATCGTCTTTTCCCACGTGTTCATACAGTTCGTGGGCGTCAAGGAGGTAGATCCCCATT‑CTCGCCAGAGATTGGCTGGGATCTGTCGGATTTGCCGGCGGTTTAGCAGTTTTTCC‑AA‑CAAA‑TGCGATAATTTTATCTTTCTCACCAACCGCCATAACGCCAAATGCGGAGGCTTCTTCAAt                                                                                                                                                                                                   <  1:410144/219‑1 (MQ=255)
                                         ctgcGGTGTTCTTGGAAAACAAATCTTTGCAAAAGTCGTGGCTGGAGTTCTCATCGCGATCGTCTTCTCCCAGCATTTCATACAGTTAGTGGGGTTAAAAGACGTAGATCCCCATC‑CCCGCAAGAGTTTTGCTCGGTTCGTTCGGCTTTGAGGGGGGTTTAGCAGTTTTTTC‑AA‑CGAA‑TTCGATATTTTTATCTTTCTCACCAACCGCCATAACGCCAAATGCGGAGGCTTCTTCAAt                                                                                                                                                                                                   <  1:269853/235‑1 (MQ=255)
                                         cTTCGGTGATCTTGGGAATCAAATCTTTGCCAAAGTCGTGGCTGGAGTTCTCATCGCGATCGTCTTCTTCCAGCAGTTCATACAGATAGTCGGCGTCAAAGACGTAGATACCCATA‑CTCGCCAGAGATTTGCTCGGATCGTTCGGCATTGACGGCGGGTTAGCAGGTTTTTC‑AA‑CGAA‑TTCGATAATTTTATCGTTCTCATCAACCGCCATAACGCCAAATGCGGAGGCTTCTTCAATc                                                                                                                                                                                                  >  2:772587/1‑239 (MQ=255)
                                         cTTCGGTGATCTTGGGAATAAAATTTTTGCAAAAGTCGTGGCTGGAGTTCTCACCGCGATCGTCTTCTCCCAGCAGTTCATACAGATAGTCGGCGTAAAAGACGTAGATCCCCATC‑CTCGCCAGAGATTTGCTGGGATCGTCCGGCATTGACGGGGGGTTAGCAGGTTTTTC‑AA‑CGAA‑TTCGATATTTTTATCGTTCTCATCAACCGCCATAACGCCAAATGCGGAGGCTTCTTCAAt                                                                                                                                                                                                   <  2:367435/238‑1 (MQ=255)
                                             tgTGATCTTGGTAATCAAGTCTTTGCAAAAGCGGGGGGTGTAGTCCTCATCTCGGTCGTCTTCTCCCATCTGTTCACACAGTTGGTCGGGTGCAAAGGGGTAGATCCCATTC‑CGCCCCAGAGTTTGGCCCGCATCGTTCGGCATTTCCGGGGGGTTAGCATTTTTTTC‑AA‑CGAA‑TGCTATATTTTTATCGTCCTCCCCAACCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGtt                                                                                                                                                                                              <  1:866029/238‑2 (MQ=255)
                                                   cTTGGGAATCAAATCTTTGCCAAAGTCGTGGCTGGAGTTCTCATCGCGATCGTCTTCTTCCAGCAGTTCATACAGATAGTCGGCGTCAAAGACGTAGATACCCATA‑CTCGCCAGAGATTTGCTCGGATCGTTCGGCATTGACGGCGGGTTAGCAGGTTTTTC‑AA‑CGAA‑TTCGATAATTTTATCGTTCTCATCAACCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCa                                                                                                                                                                                        >  2:639805/1‑239 (MQ=255)
                                                      gggAATCAAATCTTTGCCAAAGTCGTGGCTGGAGTTCTCATCGCGATCGTCTTCTTCCAGCAGTTCATACAGATAGTCGGCGTCAAAGACGTAGATACCCATA‑CTCGCCAGAGATTTGCTCGGATCGTTCGGCATTGACGGCGGGTTAGCAGGTTTTTC‑AA‑CGGA‑TTCGATAATTTTATCGTTTTCATCAACCGCCATAACGCCAAATGCGGAGGCTTTTTCAATCGGTATTGGccttc                                                                                                                                                                                     >  2:581048/1‑235 (MQ=255)
                                                         aaTCAAATCTTTGCCAAAGTCGTGGCTGGAGTTCTCATCGCGATCGTCTTCTTCCAGCAGTTCATACAGATAGTCGGCGTCAAAGACGTAGATACCCATA‑CTCGCCAGAGATTTGCTCGGATCGTTCGGCATTGACGGCGGGTTAGCAGGTTTTTC‑AA‑CGAA‑TTCGATAATTTTATCGTTCTCATCAACCGCCATAACGCCAAATGCGGAGGCTTTTTCACTCGGTACTGGCATACaa                                                                                                                                                                                   >  2:223845/1‑238 (MQ=255)
                                                         aaTCAAATCTTTGCCAAAGTCGTGGCTGGAGTTCTCATCGCGATCGTCTTCTTCCAGCAGTTCATACAGATAGTCGGCGTCAAAGACGTAGATACCCATA‑CTCGCCAGAGATTTGCTCGGATCGTTCGGCATTGACGGCGGGTTAGAAGGGTTTTC‑AA‑CGAA‑TTCGGTAATTTTTTCGTTCTCATCAAACGACATAAAGCCACATTCGGAGGCTTCTTTGATCGGGACTGGCAAAAAAg                                                                                                                                                                                  >  1:1013877/1‑239 (MQ=255)
                                                           tCAAATCTTTGCCAAAGTCGTGGCTGGAGTTCTCATCGCGATCGTCTTCTTCCAGCAGTTCATACAGATAGTCGGCGTCAAAGACGTAGATACCCATA‑CTCGCCAGAGATTTGCTCGGATCGTTCGGCATTGACGGCGGGGTAGCCGGTTTTTC‑CA‑CGAA‑TTCGATAATTTTATCGTTCTCATCAAACGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGGACTGGCATACAAGCa                                                                                                                                                                                >  1:503991/1‑239 (MQ=255)
                                                            aaattCTTTGCAAAAGCCGGGGCTGGTGTTCTCCTCGCGCTCGTCTTCTTCCAGCAGTTCATAAATTGAGGCGGCTACAAAGCCGTAGATCCCCATC‑CTCCCAAGATATTTGTGCGGATCGTTCGGCATTGAGGGGGGGTTAGCAGGTTTTAC‑AA‑CAAA‑TTCGATATTTTTATCTTTCTCAACAACCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCATACAAGCaa                                                                                                                                                                               <  1:444645/235‑1 (MQ=255)
                                                                 cttgcccaagttcggggtatgtgTTCTCATCGCGATCTTCTTCTCCCAACAGTTCATACAGTTAGTCGTAGTAAAAGGCGTAGTTCCCCTTA‑CTCGCCAGAGATTTGCCCGGTTCTTTCGGTATGGACGGGGGGTTAGCGGTTTTTAC‑AA‑CAAA‑TTCGATAATTTTATCTTTCTCATCAACCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCATACAAGCAACGGTg                                                                                                                                                                          <  1:772587/217‑1 (MQ=255)
                                                                    tGCCAAAGTCGTGGCTGGAGTTCTCATCGCGATCGTCTTCTTCCAGCAGTTCATACAGATAGTCGGCGTCAAAGACGTAGATACCCATA‑CTCGCCAGAGATTTGCTCGGATCGTTCGGCATTGACGGCGGGTTAGCAGGTTTTTC‑AA‑CGAA‑TTCGATAATTTTATCGTTCTCATCAACCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGGACTGGCATACAAGCAACGGTGCa                                                                                                                                                                        >  2:236360/1‑238 (MQ=255)
                                                                    tGCCAAAGTCGTGGCTGGAGTTCTCATCGCGATCGTCTTCTTCCAGCAGTTCATACAGATAGTCGGCGTCAAAGACGTAGATACCCATA‑CTCGCCAGAGATTTGCTCGGATCGTTCGGCATTGACGGCGGGTTAGCAGGGTTTTC‑AA‑CGAA‑TTCGATAATTTTATCGTTCTCATCAACCGCCATAACGCCAAATGCGGAGGCTTCTTCTATCGGGACTGGCATACAAGAAACGGTGCa                                                                                                                                                                        >  2:701206/1‑238 (MQ=255)
                                                                               tGGCTGGAGTTCTCATCGCGATCGTCTTCTTCCAGCAGTTCATACAGATAGTCGGCGTCAAAGACGTAGATACCCATA‑CTCGCCAGAGATTTGCTCGGATCGTTCGGCATTGACGGCGGGTTAGCAGGGTTTTC‑AA‑CGAA‑TTCGATAATTTTATCGTTCTCATCAACCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCATACAAGCAACGGTGCAACGTGCGCCttt                                                                                                                                                            >  2:1066062/1‑239 (MQ=255)
                                                                               tGGCTGGAGTTCTCATCGCGATCGTCTTCTTCCAGCAGTTCATACAGATAGTCGGCGTCAAAGACGTAGATACCCATA‑CTCGCCAGAGATTTGCTCGGATCGTTCGGCATTGACGGCGGGTTAGCAGGGTTTTC‑AA‑CGAA‑TTCGATAATTTTATCGTTCTCATCAACCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCATACAAGCAACGGTGCAACGTGCGCCttt                                                                                                                                                            >  2:536438/1‑239 (MQ=255)
                                                                                     gAGTTCTCATCGCGATCGTCTTCTTCCAGCAGTTCATACAGATAGTCGGCGTCAAAGACGTAGATACCCATA‑CTCGCCAGAGATTTGCTCGGATCGTTCGGCATTGACGGCGGGGTAGGAGGGTTTTC‑AA‑CGAA‑TTCGATAATTTTATCGTTCTCATCAACCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGGACTGGCATACAAGCAACGGGGCAACGTGCGCCCTTTTCGAc                                                                                                                                                      >  1:235280/1‑239 (MQ=255)
                                                                                                     cGTCTTCTTCCAGCAGTTCATACAGATAGTCGGCGTCAAAGACGTAGATACCCATA‑CTCGCCAGAGATTTGCTCGGCTCGTTCGGCATTGACGGCGGGTTAGCAGGTTTTTC‑AA‑CGAA‑TTCGATAATTTTATCGTTCGCAACAACCGCCATAACGCCACATGCGGAGGCTTCTTCAATCCGTATTGGc                                                                                                                                                                                         >  2:413853/1‑188 (MQ=255)
                                                                                                      tgctgcttcCAGTAGTTCATAAGTTTGGTCGGCTTAAAAGCCGTAAACCCCCATC‑CGCCCCAGAATGTTGCGGGGTTCGTCTGGTATTGGCGGGGGGTTAGCAGTTTTTAC‑AG‑AAAA‑TTGGATATTTTTAGGGTTCGCACGACCCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCATACAAGCAACGGTGCAACGTGCGCCTTTTTCGACGTGATCGATAAGCATa                                                                                                                                      <  1:254873/233‑1 (MQ=255)
                                                                                                        cttctcCCAGCGGTTCATACAGATAGTCGGCGTCAAAGACTTAGATCCCCATC‑CTCGCCAGAGATTTGCTCGGATCGTTCGGCATGGACGGGGGTTTAGCGGTTTTTTC‑AC‑CAAA‑TTCTATAATTTTATCTTTCCCATCAACCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCATACAAGCAACGGTGCAACGTGCGCCTTTTTCGACGTGATCGATAAGCATAcgc                                                                                                                                   <  2:1065736/239‑1 (MQ=255)
                                                                                                         ttcttcCAGCAGTTCATACAGATAGTCGGCGTCAAAGACGTAGATACCCATA‑CTCGCCAGAGATTTGCTCGGATCGTTCGGCATTGACGGCGGGTTAGCAGGTTTTTC‑AA‑CGAA‑TTCGATAATTTTATCGTTCTCATCAACCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCATACAAGCAACGGTGCAACGTGCGCCTTTTTCGACGTGATCGATAAGCATAcgcg                                                                                                                                  >  2:700256/1‑239 (MQ=255)
                                                                                                           cttcCAGCATTTCATACAGATAGTCGGCGTAAAAGAGTTAGACACCCTTC‑CTCGCAAGAGATTTGGTCGGATCGTTCGGCATGGCCGGCGGTTTAGCAGTTTTTTC‑AA‑CAAA‑TTCGATAATTTTATCGTTCTCACCAACCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCATACAAGCAACGGTGCAACGTGCGCCTTTTTCGACGTGATCGATAAGCATACGCGag                                                                                                                                <  1:337468/239‑1 (MQ=255)
                                                                                                           ctcccAGCAGTAAATACAGATGGGCGGCGTCAAAGAGTGAGACACCCACA‑CCCGCCAGAGTTTGGCCGGGATCGTTCGGCATTGACGGCGGTTAAGAAGGTTTTCC‑AA‑CGAA‑TTCTATATTTTTATCGTTCTCACCAACCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCATACAAGCAACGGTGCAACGTGCGCCTTTTTCGACGTGATCGATAAGCATACGCGag                                                                                                                                <  1:223845/236‑1 (MQ=255)
                                                                                                              ccAGCAGTTCATACAGATAGTCGGCGTCAAAGACGTAGATACCCATA‑CTCGCCAGAGATTTGCTCGGATCGTTCGGCATTGACGGCGGGTTAGCAGGTTTTTC‑AA‑CGAA‑TTCGATAATTTTATCGTTCTCATCAACCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCATACAAGCAACGGGGCAACGGGCGCCTTTTTCGACGTGATCGATAAGCATACGCGagtag                                                                                                                             >  1:390983/1‑239 (MQ=255)
                                                                                                               caggagTACATAAAGGTGGTGGGCGTAAAAGAGTTAGACACCCTTA‑CTCGCAAGATATTTGCTGGGAGGTTTGGGCATTGAGGGGGGTTTAGGAGTTTTTTC‑AA‑CGAA‑TTCGATATTTTTATCTTTCTCACCAACCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCATACAAGCAACGGTGCAACGTGCGCCTTTTTCGACGTGATCGATAAGCATACGCGagtagt                                                                                                                            <  1:123841/239‑1 (MQ=255)
                                                                                                                agcagTTCATACAGATAGTCGGCGTCAAAGACGTAGATACCCATA‑CTCGCCAGAGATTTGCTCGGATCGTTCGGCATTGACGGCGGGTTAACAGGTTTTTC‑AA‑CGAA‑TTCGATAATTTTATCGTTCTCTTCAACGGACATAACGCCAAATACGGAGGCTTCTTCAATCGGTACTGGCATACAAGCAACGGAGCAACGTGCGCCTTTTTCACCGTGATCGGTAAGCATAcgcggggag                                                                                                                             >  2:291762/1‑232 (MQ=255)
                                                                                                                 acaattCATAAAGATAGGCGGCGAAAAAGTGGTAGATCCCCATC‑CTCGCGGGAAATTTGGCAGGATCGTTCGGCTTTTACGGGGGGTAAGGAGTTTTTAA‑AA‑CAAA‑TGCGAAAATTTTATCGTCCTCACCAACCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCATACAAGCAACGGTGCAACGTGCGCCTTTTTCGACGTGATCGATAAGCATACGCGAGTAGTct                                                                                                                          <  1:335537/235‑1 (MQ=255)
                                                                                                                   agttaaaaaagagatgcGGCGTCAAAGACGTAGATCCCCATC‑CTCGCAAGAGATTGGCTCGGTTCGTTCGGCATTGAGGGGGGTTTAGCGGTTTTTTC‑AA‑CGAA‑TTCGATAATTTTACCTTTCTCACAAACCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCATACAAGCAACGGTGCAACGTGCGCCTTTTTCGACGTGATCGATAAGCATACGCGAGTAGTc                                                                                                                           <  2:664823/220‑1 (MQ=255)
                                                                                                                      taaaaCAGATGTTCGGAGTAAAAGCCGTAAATACCCATC‑CTCGCCAGAGTTTTGCCCGGATCGTCCGGCATTGACGGGGGGTAAGCAGTTTTTCA‑AA‑CAAA‑TTCGATAATTTTATCGTTCTCATCAACCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCATACAAGCAACGGTGCAACGTGCGCCTTTTTCGACGTGATCGATAAGCATACGCGAGTAGTcttgctt                                                                                                                     <  2:692985/235‑1 (MQ=255)
                                                                                                                      tCATAAAGATGGCCGGCGAAAAAGACGTAAATCCCCATC‑CCCGCCAGAGTTTTGCTGGGATCGTTCGGTATTGCGGGCGGTTTAGCATTTTTTTC‑AA‑CGAA‑TTCGATAATTTTTTCGTTCTCATCAACCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCATACAAGCAACGGTGCAACGTGCGCCTTTTTCGACGTGATCGATAAGCATACGCGAGTAGTcttgctt                                                                                                                     <  2:700883/239‑1 (MQ=255)
                                                                                                                       cATATAGATGGGCGGCATAAAAAACTTAGA‑CCCCATGCCTCGCCAGAGTTTGACTCGGTTGGCTGGTCATTGACGGCGTATACACAGTTTTTCC‑AA‑CAAT‑TTTTATATTTTTATCTTTCCCACCACCCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTTGCATACAAGCATCGTTGCATCGTGCGCCTTTGTCGATGTGATCGATAAGCATACGCGAGTAGTcttgctt                                                                                                                     <  1:309691/238‑1 (MQ=255)
                                                                                                                            agttgggcGGCGTAAAAGCCGTAAATCCCCATC‑CTCGCCAGAGTTTTGCGCGGATCGTTCGGTATCGACGGGGGGTTAGCAGGTTTTTC‑AA‑AAAC‑TTGGATATTTTTATCGTTCCCAACAACCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCATACAAGCAACGGTGCAACGTGCGCCTTTTTCGAc                                                                                                                                                      <  1:594262/193‑1 (MQ=255)
                                                                                                                               tAGTCGGCGTCAAAGACGTAGATACCCATA‑CTCGCCAGAGATTTGCTCGGATCGTTCGGCATTGACGGCGGGTTAGCAGGTTTTTC‑AA‑CGAA‑TTCGATAATTTTATCGTTCTCATCAACCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCGTACAAGCAACGGTGCAACGTGCGCCTTTTTCGACGTGATCGATAAGCATACGCGAGTAGTCTTGCTTGTAGATATg                                                                                                            >  2:697813/1‑239 (MQ=255)
                                                                                                                               tAGTCGGCGTCAAAGACGTAGATACCCATA‑CTCGCCAGAGATTTGCTCGGATCGTTCGGCATTGACGGCGGGTTAGCAGGTTTTTC‑AA‑CGAA‑TTCGATAATTTTATCGTTCTCATCAACCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCCTACAAGCAACGGTGCAACGTGCGCCTTTTTCGACGTGATCGATAAGCATACGCGAGTAGTCTTGCTTGTAGGTATg                                                                                                            >  2:148753/1‑239 (MQ=255)
                                                                                                                               tAGTCGGCGTCAAAGACGTAGATACCCATA‑CTCGCCAGAGATTTGCTCGGATCGTTCGGCATTGACGGCGGGTTAGCAGGTTTTTC‑AA‑CGAA‑TTCGATAATTTTATCGTTCTCATCAACCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCATACAAGCAACGGTGCAACGTGCGCCTTTTTCGACGTGATCGATAAGCATACGCGAGTAGTCTTGCTTGTAGATATg                                                                                                            >  2:120230/1‑239 (MQ=255)
                                                                                                                               aagacGGCGTCAAAGACGTAGATACCCATA‑CGCGCCAGAGATTTGCTCGGATCGTTCGGCATTGACGGCGGGTTAGCAGGTTTTTC‑AA‑CGAA‑TTCGATAATATAATCGGTCCCACCAACCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCATACAAGCAACGGTGCAACGGGCGCCTTTTTCGACGTGATCGATAAGCATACGCGAGTAGTCTTGCTTGTAGATATg                                                                                                            >  1:379018/5‑239 (MQ=255)
                                                                                                                                    ggCGTCAAAGACGTAGATACCCATA‑CTCGCCAGAGATTTGCTCGGATCGTTCGGCATTGACGGCGGGTTAGCAGGTTTTTC‑AA‑CGAA‑TTCGATAATTTTTTCGTTCTCATCAACCGCCCTAACGCCAAATGCGGGGGCTTCTTCAATCGGTACTGGCCTACAAGCAACGGCGCAACGTGCGCCTTTTTCTACGTGCTCGGTTAGCCTACGCGGGTAGTCTTCCTTGTAGGTATGGTcg                                                                                                        >  1:533006/1‑238 (MQ=255)
                                                                                                                                    ggCGTCAAAGACGTAGATACCCATA‑CTCGCCAGAGATTTGCTCGGATCGTTCGGCATTGACGGCGGGTTAGCAGGTTTTTC‑AA‑CGAA‑TTCGATAATTTTATCGTTCTCATCAACCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCATACAAGCAACGGTGCAACGTGCGCCTTTTTCGACGTGATCGATAAGCATACGCGAGTAGTCTTGCTTGTAGATATGGTcg                                                                                                        >  2:218671/1‑238 (MQ=255)
                                                                                                                                    ggCGTCAAAGACGTAGATACCCATA‑CTCGCCAGAGATTTGCTCGGATCGTTCGGCATTGACGGCGGGTTAGCAGGTTTTTC‑AA‑CGAA‑TTCGATAATTTTATCGTTCTCATCAACCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCATACAAGCAACGGTGCAACGCGCGCCTTTTTCGACGTGATCGATAAGCATACGCGAGTAGTCTTGCTTGTAGATATGGGcg                                                                                                        >  2:914370/1‑238 (MQ=255)
                                                                                                                                     ggtaaaaaGCCGTAGATCCCCACA‑CCCCCAAGAGTTTTGCTCGGATCGTTCGGCATTGAGGGGGGTTTAGCAGTTTTTTA‑AA‑CAAA‑TTATATATTTTTATCGTTCTAAACACCCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCATACAAGCAACGGTGCAACGTGCGCCTTTTTCGACGTGATCGATAAGCATACGCGAGTAGTCTTGCTTGTAGATATGGTcgc                                                                                                       <  2:364424/233‑1 (MQ=255)
                                                                                                                                      cGTCAAAGACGTAGATACCCATA‑CTCGCCAGAGATTTGCTCGGATCGTTCGGCATTGACGGCGGGTTAGCAGGTTTTTC‑AA‑CGAA‑TTCGATAATTTTATCGTTCTCATCAACCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCATACAAGCAACGGTGCAACGTGCGCCTTTTTCGACGTGATCGATa                                                                                                                                            >  2:594262/1‑200 (MQ=255)
                                                                                                                                      cGTCAAAGACGTAGATACCCATA‑CTCGCCAGAGATTTGCTCGGATCGTTCGGCATTGACGGCGGGTTAGCAGGTTTTTC‑AA‑CGAA‑TTCGATAATTTTATCGTTCTCATCAACCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCATACAAGCAACGGTGCAACGTGCGCCTTTTTCGACGTGATCGATAAGCATACGCGAGTAGTCTTGCTTGTAGATATGGTcgccc                                                                                                     >  2:1091881/1‑239 (MQ=255)
                                                                                                                                      cGTCAAAGACGTAGATACCCATA‑CTCGCCAGAGATTTGCTCGGATCGTTCGGCATTGACGGCGGGTTAGCAGGTTTTTC‑AA‑CGAA‑TTCGATAATTTTATCGTTCTCATCAACCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCATACAAGCAACGGTGCAACGTGCGCCTTTTTCGACGTGATCGATAAGCATACGCGAGTAGTCTTGCTTGTAGATATGGTcgccc                                                                                                     >  1:940980/1‑239 (MQ=255)
                                                                                                                                      cGTCAAAGACGTAGATACCCATA‑CTCGCCAGAGATTTGCTCGGATCGTTCGGCATTGACGGCGGGTTAGCAGGTTTTTC‑AA‑CGAA‑TTCGATAATTTTATCGTTCTCATCAACCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCATACAAGCAACGGTGCAACGTGCGCCTTTTTCGACGTGATCGATAAGCATACGCGAGTAGTCTTGCTTGTAGATATGGTcgccc                                                                                                     >  2:878661/1‑239 (MQ=255)
                                                                                                                                      cGTCAAAGACGTAGATACCCATA‑CTCGCCAGAGATTTGCTCGGATCGTTCGGCATTGACGGCGGGTTAGCAGGGTTTTC‑AA‑CGAA‑TTCGATAATTTTATCGTTCTCATCAACCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCATACAAGCAACGGTGCAACGTGCGCCTTTTTCGACGTGATCGATAAGCATACGCGAGTAGTCTTGCTTGTAGATATGGTcgccc                                                                                                     >  1:233524/1‑239 (MQ=255)
                                                                                                                                        tCAAAGACGTAGATACCCATA‑CTCGCCAGAGATTTGCTCGGATCGTTCGGCATTGACGGCGGGTTAGCAGGTTTTTC‑CA‑CGAA‑TTCGATAATTTTATCGTTCTCATCAACCCCCATAACGCCAAATGCGGAGGGTTTTTCAATCGGTACTGGCATACAAGCCAAGGTGCAACGCGCGCCTTTTTCGACGGGATCGATAAGCATACGCGGGTAGTCTTGCTTGTTGATATGGTagcccg                                                                                                    >  1:747335/1‑238 (MQ=255)
                                                                                                                                                      aCCCATA‑CTCGCCAGAGATTTGCTCGGATCGTTCGGCATTGACGGCGGGTTAGCAGGTTTTTC‑AA‑CGAA‑TTCGATAATTTTATCGTTCTCATCAACCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCATACAAGCAACGGTGCAACg                                                                                                                                                                     >  2:884817/1‑159 (MQ=255)
                                                                                                                                                        ccATA‑CTCGCCAGAGATTTGCTCGGATCGTTCGGCATTGACGGCGGGTTAGCAGGTTTTTC‑AA‑CGAA‑TTCGATAATTTTATCGTTCTCATCAACCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCATACAAGCAACGGTGCAACGTGCGCCTTTTTCGACGTGATCGATAAGCATACGCGGGTAGTCTTGCTTGTAGATATGGTCGCCCGCCAGGATCACCACGTAt                                                                                   >  2:433363/1‑239 (MQ=255)
                                                                                                                                                        ccATA‑CTCGCCAGAGATTTGCTCGGATCGTTCGGCATTGACGGCGGGTTAGCAGGTTTTTC‑AA‑CGAA‑TTCGATAATTTTATCGTTCTCATCAACCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCATACAAGCAACGGTGCAACGTGCGCCTTTTTCGACGTGATCGATAAGCATACGCGAGTAGTCTTGCTTGTAGCTATGGGCGCCCGCCAGGATCACCCCGGAt                                                                                   >  1:707409/1‑239 (MQ=255)
                                                                                                                                                         cATA‑CCCCCCAGAGATTTGCTCGGATCGTTCGGCTTGGACGGGGGGTTAGCAGGTTTTTC‑AA‑CAAA‑TTCGATATTTTTATCTTCCTAACAAACCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCATACAAGCAACGGTGCAACGTGCGCCTTTTTCGACGTGATCGATAAGCATACGCGAGTAGTCTTGCTTGTAGATATGGTCGCCCGCCAGGATCACCACGTAtt                                                                                  <  2:1203076/239‑1 (MQ=255)
                                                                                                                                                           ttaccccccaggtttttGCTGGGTTTGTTCGGCTTTGCCGGGGTGTAAGCTTTTTTTAC‑AA‑CTAT‑TCCGTAATTTTTGCCGTTCTCAACAACCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCATACAAGCAACGGTGCAACGTGCGCCTTTTTCGACGTGATCGATAAGCATACGCGAGTAGTCTTGCTTGTAGATATGGTCGCCCGCCAGGATCACCACGTATTc                                                                                 <  1:558684/225‑1 (MQ=255)
                                                                                                                                                           tA‑CTCGCCAGAGATTTGCTCGGATCGTTCGGCATTGACGGCGGGTTAGCAGGGTTTTC‑AA‑CGAA‑TTCGATAATTTTATCGTTCTCATCAACCGCCATAAAGCCAAATGCGGAGGCATCTTCAATCGGCACTGGCATACAAGCAAAGGTGCAACGTGCGGCTTTTTTGACGTGATCGATAAGGAATAGCGGGTAGACTTGGTTGAAGACACGGTCGCCCGTCAGGAACACCACGGATTc                                                                                 >  1:387460/1‑238 (MQ=255)
                                                                                                                                                                  ccAGAGATTTGCTCGGATCGTTCGGCATTGACGGCGGGTTAGCAGGTTTTTC‑AA‑CGAA‑TTCGATAATTTTATCGTTCTCATCAACCGCCATAACGCCAAATGCGGAGGCTTCTTCCATCGGTACTGGCATACAAGCAACGGTGCAACGTGCGCCTTTTTCGACGTGATCGATAAGCATACGCGAGTAGTCTTTCTTGTAGATATGGTCGCCCCCCAGGGTCACCCCGTATTCCGCTTta                                                                          >  1:978615/1‑239 (MQ=255)
                                                                                                                                                                   cAGAGATTTGCTCGGATCGTTCGGCATTGACGGCGGGTTAGCAGGTTTTTC‑AA‑CGAA‑TTCGATAATTTTATCGTTCTCATCCACCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCATACAAGCAACGGTGCAACGTGCGCCTTTTTCGGCGAGATCGATAAGCATACGCGAGTAGTCTTGCTTGTAGATATGGTCGCCCGCCCGGATCACCACGTATTCCGCTTtat                                                                         >  2:1002207/1‑239 (MQ=255)
                                                                                                                                                                   cAGAGATTTGCTCGGATCGTTCGGCATTGACCGCGGGCTAGCAGGTTTTTC‑AA‑CGAA‑TTCGATAATTTTATCGTTCTCATCAACCGCCATAACGCCAAATGCGGTGGCTTTTTCAATCGGAAATGGCATACAAGCAACGGTGCAACGTGCGCCTTTTTCGACTGGATCGATAATCCTACGCGTGTAGTCTTGCTTGTCGATAGGGTCGCCCACCAGACTCATCACGTATTCCGCTTtat                                                                         >  1:798273/1‑239 (MQ=255)
                                                                                                                                                                      agaTTTGCTCGGATCGTTCGGCATTGACGGCGGGTTAGCAGGTTTTTC‑AA‑CGAA‑TTCGATAATTTTATCGTTCTCATCAAACGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCATACAAGCAACGGTGCAACGTGCGCCTTTTTCGACGTGATCGATAAGCATACGCAAGTAGTCTTGCTTGTAGATATGGTCGCCCGCCAGGATCACCACGTATTCCGCTTTATAAc                                                                      >  2:171711/1‑239 (MQ=255)
                                                                                                                                                                             ctcgtttCGTCCGCCATTGCCGGAGGTTTAGCAGTTTTTAC‑AA‑CGAT‑TCAGATATTTTTATCTTTCTAACCACCCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCATACAAGCAACGGTGCAACGTGCGCCTTTTTCGACGTGATCGATAAGCATACGCGAGTAGTCTTGCTTGTAGATATGGTCGCCCGCCAGGATCACCACGTATTCCGCTTTATAACGGCGGa                                                                <  1:639805/232‑1 (MQ=255)
                                                                                                                                                                                    cggtcgccatggggggggggTGAGCAGG‑TTTTC‑AA‑CCAATTTCGATAATTTTATCTTCCTCATCAACCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCATACAAGCAACGGTGCAACGTGCGCCTTTTTCGACGTGATCGATAAGCATACGCGAGTAGTCTTGCTTGTAGATATGGTCGCCCGCCAGGATCACCACGTATTCCGCTTTATAACGGCGGATAATGTc                                                         <  2:533006/221‑1 (MQ=255)
                                                                                                                                                                                            aTTGACGGCGGGTTAGCAGGTTTTTC‑AA‑CGAA‑TTCGATAATTTTATCGTTCTCATCAACCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCATACAAGCAACGGTGCAACGTGCGCCTTTTTCGACGTGATCGATAAGCATACGCGAGTAGTCTTGCTTGTAGATATGGTCGCCCGCCAGGATCACCACGTATTCCGCTTTATAACGGCGGATAATGTCGAGGTTTTg                                                >  1:1053366/1‑239 (MQ=255)
                                                                                                                                                                                                aCGGCGGGTTAGCAGGTTTTTC‑AA‑CGAA‑TTCGATAATTTTATCGTTCTCATCAACCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCATACAAGCAACGGTGCAACGTGCGCCTTTTTCGACGTGATCGATAAGCATACGCGAGTAGTCTTGCTTGTAGATATGGTCGCCCGCCAGGATCACCACGTATTCCGCTTTATAACGGCGGATAATGTTGAGGGTTTGGGTg                                            >  2:1077997/1‑239 (MQ=255)
                                                                                                                                                                                                aCGGCGGGTTAGCAGGTTTTTC‑AA‑CGAA‑TTCGATAATTTTATCGTTCTCATCAACCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCATACAAGCAACGGTGCAACGTGCGCCTTTTTCGACGTGATCGATAAGCATACGCGAGTAGTCTTGCTTGTAGATATGGTCGCCCGCCAGGATCACCACGTATTCCGCTTTATAACGGCGGATAATGTCGAGGGTTTtgg                                              >  2:868255/1‑237 (MQ=255)
                                                                                                                                                                                                       gcttgtgtgTTTTTC‑AACCGAT‑GTCAAAATTTTTCTCGTTCCTATCACCCCCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCATACAAGCAACGGTGCAACGTGCGCCTTTTTCGACGTGATCGATAAGCATACGCGAGTAGTCTTGCTTGTAGATATGGTCGCCCGCCAGGATCACCACGTATTCCGCTTTATAACGGCGGATAATGTCGAGGTTTTGGGTGACCGc                                       <  1:578408/230‑1 (MQ=255)
                                                                                                                                                                                                          aGCAGGTTTTTC‑AA‑CGAA‑TTCGATAATTTTATCGTTCTCATCAACCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCATACAAGCAACGGTGCAACGTGCGCCTTTTTCGACGTGATCGATAAGCATACGCGAGTAGTCTTGCTTGTAGATATGGTCGCCCGCCAGGATCACCACGTATTCCGCTTTATAACGGCGGATAATGACGAGGTTTTGGGTCACCCCATcc                                   >  2:1098120/1‑237 (MQ=255)
                                                                                                                                                                                                                tttttC‑AA‑CAAA‑TTCAATATTTTTTTCTTTCTCAACACCCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCATACAAGCAACGGTGCAACGTGCGCCTTTTTCGACGTGATCGATAAGCATACGCGAGTAGTCTTGCTTGTAGATATGGTCGCCCGCCAGGATCACCACGTATTCCGCTTTATAACGGCGGATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGTGc                            <  2:235280/239‑1 (MQ=255)
                                                                                                                                                                                                                 ttttC‑AA‑CGAA‑TTCGATAATTTTATCGTTCTCATCAACCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCATACAAGCAACGGTGCAACGTGCGCCTTTTTCGACGTGATCGATAAGCATACGCGAGTAGTCTTGCTTGTAGATATGGTCGCCCGCCAGGATCACCACGTATTCCGCTTTATAACGGCGGATAATGGCGAGGTTTTGGGTGTACGCATTTGCgggggc                           >  1:345510/1‑235 (MQ=255)
                                                                                                                                                                                                                  ttac‑AA‑CGAA‑TTCGATAATTTTATTGTCCCCACAAACCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCATACAAGCAACGGTGCAACGTGCGCCTTTTTCGACGTGATCGATAAGCATACGCGAGTAGTCTTGCTTGTAGATATGGTCGCCCGCCAGGATCACCACGTATTCCGCTTTATAACGGCGGATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCcg                          <  2:926438/236‑1 (MQ=255)
                                                                                                                                                                                                                             a‑TTCGATAATTTTATCGTTCTCATCAACCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCATACAAGCAACGGTGCAACGTGCGCCTTTTTCGACGTGATCGATAAGCATACGCGAGTAGTCTTGCTTGTAGATATGGTCGCCCGCCAGGATCACCACGTATTCCGCTT‑‑TAACGGCGGATAATGTCGAGGTCTTGGGTGACCGCATCCGCGGTGCCGCGATaaccggt               >  1:609569/1‑233 (MQ=255)
                                                                                                                                                                                                                                    agtttttAGCGTTCCTACCATCCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCATACAAGCAACGGTGCAACGTGCGCCTTTTTCGACGTGATCGATAAGCATACGCGAGTAGTCTTGCTTGTAGATATGGTCGCCCGCCAGGATCACCACGTATTCCGCTTTATAACGGCGGATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCCGCGATACCAGTTTTc            <  1:29971/235‑1 (MQ=255)
                                                                                                                                                                                                                                      tttttCTCGTTCTCAACAGCCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCATACAAGCAACGGTGCAACGTGCGCCTTTTGCGACGTGATCGATAAGCATACGCGAGTAGTCTTGCTTGTAGATATGGTCGCCCGCCAGGATCACCACGTATTCCGCTT‑‑TAACGGCGGATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCCGCGATACCAGTTTTCCCCtt       <  2:609569/238‑1 (MQ=255)
                                                                                                                                                                                                                                        tatactgTTCTCATAAACCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCATACAAGCAACGGTGCAACGTGCGCCTTTTTCGACGTGATCGATAAGCATACGCGAGTAGTCTTGCTTGTAGATATGGTCGCCCGCCAGGATCACCACGTATTCCGCTTTATAACGGCGGATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCCGCGATACCAGTTTTCCCCtt       <  2:993333/233‑1 (MQ=255)
                                                                                                                                                                                                                                             cGTCCTCACCAGCCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCATACAAGCAACGGTGCAACGTGCGCCTTTTTCGACGTGATCGATAAGCATACGCGAGTAGTCTTGCTTGTAGATATGGTCGCCCGCCAGGATCACCACGTATTCCGCTTTATAACGGCGGATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCCGCGATACCAGTTTTCCCCTTTCAtt  <  2:210703/239‑1 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                             
TACGCAAGAGAGCGGGAACGGGTGCGCATAGGCCAGACCGGCTTCGGTGATCTTGGGAATCAAATCTTTGCCAAAGTCGTGGCTGGAGTTCTCATCGCGATCGTCTTCTTCCAGCAGTTCATACAGATAGTCGGCGTCAAAGACGTAGATACCCATA‑CTCGCCAGAGATTTGCTCGGATCGTTCGGCATTGACGGCGGGTTAGCAGGTTTTTC‑AA‑CGAA‑TTCGATAATTTTATCGTTCTCATCAACCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCATACAAGCAACGGTGCAACGTGCGCCTTTTTCGACGTGATCGATAAGCATACGCGAGTAGTCTTGCTTGTAGATATGGTCGCCCGCCAGGATCACCACGTATTCCGCTTTATAACGGCGGATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCCGCGATACCAGTTTTCCCCTTTCATT  >  NC_000913/3568536‑3569007

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