Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 354,812 C→T intergenic (+220/‑110) prpE → / → codB propionate‑‑CoA ligase/cytosine transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913354,8120CT99.1% 428.1 / NA 105intergenic (+220/‑110)prpE/codBpropionate‑‑CoA ligase/cytosine transporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (58/46);  minor base G (0/1);  total (58/47)
Fisher's exact test for biased strand distribution p-value = 4.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 105 total aligned reads displayed.

CGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCACTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGCAAGATAACAACTTTAGCCAGGGGCCAGTCCCGCAGTCGGCGCGGAAAGGGGTATTGGCATTGACGTTCGTCATGCTGGGATTAACCTTCTTTTCCGCCAGT  >  NC_000913/354579‑355029
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                         
cGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATtca                                                                                                                                                                                                                        <  2:127830/237‑1 (MQ=255)
      aaGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCatat                                                                                                                                                                                                                 <  1:425018/238‑1 (MQ=255)
            aGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATaaa                                                                                                                                                                                                              <  1:550452/235‑1 (MQ=255)
             ggTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAata                                                                                                                                                                                                         >  2:674588/1‑239 (MQ=255)
             ggTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAata                                                                                                                                                                                                         >  1:610369/1‑239 (MQ=255)
             ggTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACCGCTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATaaaa                                                                                                                                                                                                             >  2:1711337/1‑235 (MQ=255)
                         aTGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCTCCAcc                                                                                                                                                                                             >  2:2196418/1‑238 (MQ=255)
                                               cACTGTGTGCCGAAGCCGGATGCGTCGCTGACGCGTGTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCtt                                                                                                                                                                        <  2:330307/238‑1 (MQ=255)
                                                    tgtgCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAACCGATTGCTTTTTATCTTCAGATg                                                                                                                                                                  >  2:779162/1‑239 (MQ=255)
                                                    tgtgCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATg                                                                                                                                                                  >  2:1167323/1‑239 (MQ=255)
                                                    tgtgCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATg                                                                                                                                                                  >  2:423626/1‑239 (MQ=255)
                                                    tgtgCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCAACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAatatat                                                                                                                                                                                                      <  1:1195684/203‑1 (MQ=255)
                                                    tgtgCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCAACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAatatat                                                                                                                                                                                                      >  2:1195684/1‑203 (MQ=255)
                                                     gtgCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAg                                                                                                                                                                     <  2:645516/235‑1 (MQ=255)
                                                          ggtgtcggATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAACAATATATTTCCCCACGAAAACGATTGACTTTTATCTTCAGATGAATAgt                                                                                                                                                            >  2:1210152/8‑238 (MQ=255)
                                                           aTGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGa                                                                                                                                                            >  1:96230/1‑239 (MQ=255)
                                                            tGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGaa                                                                                                                                                           >  1:786263/1‑238 (MQ=255)
                                                                       aCGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGAttt                                                                                                                                               >  1:1469518/1‑239 (MQ=255)
                                                                        cGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCtt                                                                                                                                                                        <  1:1166303/213‑1 (MQ=255)
                                                                        cGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCtt                                                                                                                                                                        >  2:1166303/1‑213 (MQ=255)
                                                                        cGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGAtttt                                                                                                                                              >  2:1238683/1‑239 (MQ=255)
                                                                             aCGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTggg                                                                                                                                         >  1:2280908/1‑239 (MQ=255)
                                                                                   tttATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCaa                                                                                                                                   <  1:2188314/239‑1 (MQ=255)
                                                                                   tttATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCa                                                                                                                                    <  1:423626/238‑1 (MQ=255)
                                                                                   tttATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCa                                                                                                                                    <  2:786263/238‑1 (MQ=255)
                                                                                           gCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATc                                                                                                                                                                          >  2:1491872/1‑192 (MQ=255)
                                                                                           gCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATc                                                                                                                                                                          <  1:1491872/192‑1 (MQ=255)
                                                                                           gCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTATTAAACAGCaaa                                                                                                                           >  1:1875764/1‑239 (MQ=255)
                                                                                            ccTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCaaaa                                                                                                                          >  2:695209/1‑239 (MQ=255)
                                                                                                cGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAACCGATTCCTTTTTATCTTAAGATGAATAGAATGAGGCGGATTTTTTGGGTTTCAAACAGCACAAAtgg                                                                                                                      >  1:697554/1‑239 (MQ=255)
                                                                                                cGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATgcgg                                                                                                                                                      >  2:832465/1‑207 (MQ=255)
                                                                                                cGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATgcgg                                                                                                                                                      <  1:832465/207‑1 (MQ=255)
                                                                                                     cTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAAt                                                                                                                 >  1:1941287/1‑239 (MQ=255)
                                                                                                       gAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATgcgg                                                                                                                                                      <  2:491465/200‑1 (MQ=255)
                                                                                                       gAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATgcgg                                                                                                                                                      >  1:491465/1‑200 (MQ=255)
                                                                                                        aaCCGTAGGTCGGATAAGGCGCTCGCGTCGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTc                                                                                                              >  1:851395/1‑239 (MQ=255)
                                                                                                           cGTAGGTCGGATAAGGCGCTCGCGTCGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAAttt                                                                                                               >  1:692040/1‑235 (MQ=255)
                                                                                                                          gcgcTCGCGTCGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCttttt                                                                                                                                                                             >  1:1998194/1‑158 (MQ=255)
                                                                                                                          gcgcTCGCGTCGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCttttt                                                                                                                                                                             <  2:1998194/158‑1 (MQ=255)
                                                                                                                          gcgcTCGCGTCGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTTAAACAGCAAAAAGGGGGAATTTCGTGTCGCAAGATaaca                                                                                              >  2:1222758/1‑237 (MQ=255)
                                                                                                                          gcgcTCGCGTCGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCg                                                                                                             <  1:1798125/222‑1 (MQ=255)
                                                                                                                          gcgcTCGCGTCGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCg                                                                                                             >  2:1798125/1‑222 (MQ=255)
                                                                                                                          gcgcTCGCGTCGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGCAAGATaaca                                                                                              >  2:459953/1‑237 (MQ=255)
                                                                                                                           cgcTCGCGTCGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGCAAGATAACAACt                                                                                           >  1:1691153/1‑239 (MQ=255)
                                                                                                                             cTCGCGTCGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCttt                                                                                                                                                                               >  2:287374/1‑153 (MQ=255)
                                                                                                                             cTCGCGTCGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCttt                                                                                                                                                                               <  1:287374/153‑1 (MQ=255)
                                                                                                                              tCGCGTCGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGCAAGATAACAACttt                                                                                         <  2:2280908/238‑1 (MQ=255)
                                                                                                                               cgcgTCGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATaaa                                                                                                                                                                                                              <  1:1868387/120‑1 (MQ=255)
                                                                                                                               cgcgTCGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATaaa                                                                                                                                                                                                              >  2:1868387/1‑120 (MQ=255)
                                                                                                                               cgcgTCGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGCAAGATAACAACTTTa                                                                                        >  1:363988/1‑238 (MQ=255)
                                                                                                                                 cgTCGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGACAACGATTGCTTTTTATCTTAAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCCAACAGCAAAAAGGGGGAATTTCGTGTCGCAAGATAACACCTTAAGcc                                                                                     >  2:1649879/1‑239 (MQ=255)
                                                                                                                                   tCGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGCAAGATAACAACTTTAGCCa                                                                                    <  2:1774365/238‑1 (MQ=255)
                                                                                                                                    cGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCtt                                                                                                                                                                        >  1:2196102/1‑153 (MQ=255)
                                                                                                                                    cGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCtt                                                                                                                                                                        <  2:2196102/153‑1 (MQ=255)
                                                                                                                                    cGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGCAAGACAACAACTTTAACCAgg                                                                                  >  2:1271093/1‑239 (MQ=255)
                                                                                                                                       aTCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGCAAGATAACAACTTTAGCCAGGGGc                                                                               <  1:459953/239‑1 (MQ=255)
                                                                                                                                       aTCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGCAAGATAACAACTTTAGCCAGGGGc                                                                               <  2:1469518/239‑1 (MQ=255)
                                                                                                                                       aTCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGCAAGATAACAACTTTAGCCAGGGGc                                                                               <  1:674588/239‑1 (MQ=255)
                                                                                                                                         ccGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGCAAGATAACAACTTTAGCCAGGGGCCa                                                                             <  1:1238683/239‑1 (MQ=255)
                                                                                                                                                            ccTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGCAAGATa                                                                                                 <  1:661723/200‑1 (MQ=255)
                                                                                                                                                            ccTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGCAAGATa                                                                                                 >  2:661723/1‑200 (MQ=255)
                                                                                                                                                            ccTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGCAAGATAACAACTTTAGCCAGGGGCCAGTCCCGCAGTCGGCGCGGa                                                          >  2:175207/1‑239 (MQ=255)
                                                                                                                                                              tGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATtca                                                                                                                                                                                                                        >  1:602848/1‑79 (MQ=255)
                                                                                                                                                              tGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATtca                                                                                                                                                                                                                        <  2:602848/79‑1 (MQ=255)
                                                                                                                                                              tGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGCAAGATAACAACTTTAGCCAGGGCACAGTCCCGCCATCGGCGCGGaac                                                        >  2:1383729/1‑238 (MQ=255)
                                                                                                                                                               gATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGCAAGATAACAACTTTAGCCAGGGGCCAGTCCCGCAGTCGGCGCGGAAAg                                                       <  2:610369/239‑1 (MQ=255)
                                                                                                                                                                           aCGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACg                                                                                                                                                                                             >  1:631898/1‑93 (MQ=255)
                                                                                                                                                                           aCGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACg                                                                                                                                                                                             <  2:631898/93‑1 (MQ=255)
                                                                                                                                                                                         gCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGCAAGATAACAACTTTAGCCAGGGGCCAGTCCCGCAGTCGGCGCGGAAAGGGGTATTGGCATTGACGTTCGTCAt                              >  2:399465/1‑238 (MQ=255)
                                                                                                                                                                                              cccACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGCAAGATAACAACTTTAGCCAGGGGCCAGTCCCGCAGTCGGCGCGGAAAGGGGTATTGGCATTGAc                                       >  1:19864/1‑224 (MQ=255)
                                                                                                                                                                                              cccACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGCAAGATAACAACTTTAGCCAGGGGCCAGTCCCGCAGTCGGCGCGGAAAGGGGTATTGGCATTGAc                                       <  2:19864/224‑1 (MQ=255)
                                                                                                                                                                                               ccACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGCAAGATAACAACTTTAGCCAGGGGCCAGTCCCGCAGTCGGCGCGGAAAGGGGTATTGGCATTGACGTTCGTCATGCTGGGa                       <  1:175207/239‑1 (MQ=255)
                                                                                                                                                                                                   tGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCtt                                                                                                                                                                        >  2:1373214/1‑90 (MQ=255)
                                                                                                                                                                                                   tGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCtt                                                                                                                                                                        <  1:1373214/90‑1 (MQ=255)
                                                                                                                                                                                                   tGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGCAAGATAACAACTTTAGCCAGGGGCCAGTCCCGCAGTCGTCGCGGAAAGGGGTCTTGTCATTGACGTTCGTCATGCTGGGATTa                    >  2:1143427/1‑238 (MQ=255)
                                                                                                                                                                                                   tGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGCAAGATAACAACTTTAGCCAGGGGCCAGTCCCGCAGTCGGCGCGGAAAGGGGTATTGGCATTGACGTTCGTCATGCTggg                        >  2:860279/1‑234 (MQ=255)
                                                                                                                                                                                                   tGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGCAAGATAACAACTTTAGCCAGGGGCCAGTCCCGCAGTCGGCGCGGAAAGGGGTATTGGCATTGACGTTCGTCATGCTggg                        <  1:860279/234‑1 (MQ=255)
                                                                                                                                                                                                   tGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGCAAGATAACAACTTTAGCCAGGGGCCAGTCCCGCAGTCGGCGCGGAAAGGGGTATTGGCATTGACGTTCGTCATGCTGGGATTa                    >  1:1459077/1‑238 (MQ=255)
                                                                                                                                                                                                   tGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGCAAGATAACAACTTTAGCCAGGGGCCAGTCCCGCAGTCGGCGCGGAAAGGGGTATTGGCATTGACGTTCGTCATGCTGGGATTa                    >  1:752563/1‑238 (MQ=255)
                                                                                                                                                                                                         ttACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGCAAGATAACAACTTTAGCCAGGGGCCAGTCCCGCAGTCGGCGCGGa                                                          >  1:1367119/1‑195 (MQ=255)
                                                                                                                                                                                                         ttACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGCAAGATAACAACTTTAGCCAGGGGCCAGTCCCGCAGTCGGCGCGGa                                                          <  2:1367119/195‑1 (MQ=255)
                                                                                                                                                                                                          tACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGCAAGATAACAACTTTAGCCAggg                                                                                 >  2:553328/1‑171 (MQ=255)
                                                                                                                                                                                                          tACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGCAAGATAACAACTTTAGCCAggg                                                                                 <  1:553328/171‑1 (MQ=255)
                                                                                                                                                                                                             acCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGCAAGATAACAACTTTAGCCAGGGGCCAGTCCCGCAGTCGGCGCGGAAAGGGGTATTGGCATTGACGTTCGTCATGCTGGGATTAACCTTCtttt          <  1:1711337/239‑1 (MQ=255)
                                                                                                                                                                                                                gATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTccc                                                                                                                                                                                                 >  2:156165/1‑53 (MQ=255)
                                                                                                                                                                                                                gATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTccc                                                                                                                                                                                                 <  1:156165/53‑1 (MQ=255)
                                                                                                                                                                                                                gATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCtt                                                                                                                                                                        >  1:2106841/1‑78 (MQ=255)
                                                                                                                                                                                                                gATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCtt                                                                                                                                                                        <  2:2106841/78‑1 (MQ=255)
                                                                                                                                                                                                                gATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCaaaa                                                                                                                          >  2:1298797/1‑124 (MQ=255)
                                                                                                                                                                                                                gATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCaaaa                                                                                                                          <  1:1298797/124‑1 (MQ=255)
                                                                                                                                                                                                                gATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGCAAGATAACAACTTTAGCCAGGGGCCAGTCCCGCAGTCGGCGCGGAAAGGGGTATTGGCATTGACGTTCGTCATGCTGGGATTAACCTTCTTTTcc        >  2:649189/1‑238 (MQ=255)
                                                                                                                                                                                                                  tAATTTTTCCCCCTCCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGCAAGATAACAACTTTAGCCAGGGGCCAGTCCCGCAGTCGGCGCGGAAAGGGGTATTGGCATTGACGTTCGTCATGCTGGGATTAACCTTCTTTTCCGcc     <  1:1649879/239‑1 (MQ=255)
                                                                                                                                                                                                                  tAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCttt                                                                                                                                                                               >  2:180213/1‑69 (MQ=255)
                                                                                                                                                                                                                  tAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCttt                                                                                                                                                                               <  1:180213/69‑1 (MQ=255)
                                                                                                                                                                                                                  tAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGCAAGATAACAACTTTAGCCAGGGGCCAGACCCGCAGTCGGCGCGGAAAGGGGTATTGGCATTGAc                                       >  1:1304126/1‑205 (MQ=255)
                                                                                                                                                                                                                   aaTTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGCAAGATAACAACTTTAGCCAGGGGCCAGTCCCGCAGTCGGCGCGGAAAGGGGTATTGGCATTGACGTTCGTCATGCTGGGATTAACCTTCTTTTCCGCCa    <  1:1226696/239‑1 (MQ=255)
                                                                                                                                                                                                                   aaTTTTTCCCCCACCTTTTGGAATTCATTCATATAAAAAATATATGTCCCCACGAAAACGATTGCTGTTTGTCTTCAGATGAATAGAAGGCGGCGGATTTTTTGGGTTTCAAACGGCAAAAAGGGGGAATTTCGTGTCGAAAGATAACAACTTTAGCCAGGGGCCAGTCCCGCAGTCGGCGCGGAAAGGGGTATTGGCATTGACGTGCGTCATGCTGGGATTAACCTTCTTTTCCGcaa    <  2:1536398/239‑3 (MQ=255)
                                                                                                                                                                                                                     tttttCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGCAAGATAACAACTTTAGCCAGGGGCCAGTCCCGCAGTCGGCGCGGAAAGGGGTATTGGCATTGACGTTCGTCATGCTGGGATTAACCTTCTTTTCCGCCAGt  <  1:399465/239‑1 (MQ=255)
                                                                                                                                                                                                                               aCCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGCAAGATAACAACTTTAGCCAGGGGc                                                                               >  1:1917817/1‑152 (MQ=255)
                                                                                                                                                                                                                               aCCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGCAAGATAACAACTTTAGCCAGGGGc                                                                               <  2:1917817/152‑1 (MQ=255)
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                         
CGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCACTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGCAAGATAACAACTTTAGCCAGGGGCCAGTCCCGCAGTCGGCGCGGAAAGGGGTATTGGCATTGACGTTCGTCATGCTGGGATTAACCTTCTTTTCCGCCAGT  >  NC_000913/354579‑355029

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: