Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,326,925:1 (C)6→7 100% coding (817/1185 nt) atoB → acetyl‑CoA acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,326,9191.C53.3% 165.4 / 180.4 90coding (811/1185 nt)atoBacetyl‑CoA acetyltransferase
Reads supporting (aligned to +/- strand):  ref base . (24/18);  new base C (26/22);  total (50/40)
Fisher's exact test for biased strand distribution p-value = 8.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGCAGGAACAGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTA‑‑CCCC‑C‑CTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTACCTGCCACGCAAAAAGCGTTAC  >  NC_000913/2326821‑2327012
                                                                                                   |                                                                                                
aaGCAGGAACAGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTA‑‑c                                                                                                >  2:1689493/1‑100 (MQ=255)
 ngCAGGAACAGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCAGCTCTGGTGATTATGGAATAATCTGCGGCGCTGGCAGCAGGCCTTA‑‑cc                                                                                               >  2:2035549/2‑100 (MQ=255)
 aGCAGGAACAGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTA‑‑cc                                                                                               >  1:2268331/1‑100 (MQ=255)
   cAGGAACAGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTA‑‑cccc                                                                                             <  1:632289/100‑1 (MQ=255)
     ggAACAGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTA‑‑cccc‑c‑c                                                                                         >  2:884522/1‑100 (MQ=255)
       aaCAGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTA‑‑CCCC‑C‑CTg                                                                                       <  2:913500/100‑1 (MQ=255)
          aGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTCCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTA‑‑CCCC‑C‑CTGGCt                                                                                    <  2:331683/100‑1 (MQ=255)
           gTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTA‑‑CCCC‑C‑CTGGCTc                                                                                   >  2:1845267/1‑100 (MQ=255)
            tCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGTTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTA‑‑CCCC‑C‑CTGGCTCg                                                                                  <  2:462696/100‑1 (MQ=255)
             cACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTAC‑CCCC‑C‑CTGGCTCg                                                                                  <  1:440897/100‑1 (MQ=255)
              aCCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTAC‑CCCC‑C‑CTGGCTCGc                                                                                 <  2:1045731/100‑1 (MQ=255)
               ccGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTAC‑CCCC‑C‑CTGGCTCGCa                                                                                >  1:961934/1‑100 (MQ=255)
               ccGCTGGGAACGCGTCTGGTATTAACGACCGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTAC‑CCCC‑C‑CTGGCTCGCa                                                                                <  2:1412185/100‑1 (MQ=255)
                cGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTAC‑CCCC‑C‑CTGGCTCGCAt                                                                               <  1:1229087/100‑1 (MQ=255)
                 gCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTAC‑CCCC‑C‑CTGGCTCGCAtt                                                                              >  1:1482177/1‑100 (MQ=255)
                    gggAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTAC‑CCCC‑C‑CTGGCTCGCATTaaa                                                                           >  1:994988/1‑100 (MQ=255)
                      gAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTAC‑CCCC‑C‑CTGGCTCGCATTAAAAg                                                                         >  1:348982/1‑100 (MQ=255)
                      gAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTAC‑CCCC‑C‑CTGGCTCGCATTAAAAg                                                                         >  1:575519/1‑100 (MQ=255)
                       aaCGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTAC‑CCCC‑C‑CTGGCTCGCATTAAAAGt                                                                        <  1:1931096/100‑1 (MQ=255)
                       aaCGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTA‑‑CCCC‑C‑CTGGCTCGCATTAAAAGtt                                                                       >  1:171539/1‑100 (MQ=255)
                         cgcgTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTAC‑CCCC‑C‑CTGGCTCGCATTAAAAGTTa                                                                      >  1:634000/1‑100 (MQ=255)
                         cgcgTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTAC‑CCCC‑C‑CTGGCTCGCATTAAAAGTTa                                                                      <  2:1745521/100‑1 (MQ=255)
                          gcgTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTAC‑CCCC‑C‑CTGGCTCGCATTAAAAGTTAt                                                                     <  1:797752/100‑1 (MQ=255)
                            gTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAGGAATCTGCGGCGCTGGCAGCAGGCCTTA‑‑CCCCGC‑CTGGCTCGCATTAAACGTGATGc                                                                   >  2:1831604/1‑100 (MQ=255)
                             tCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTAC‑CCCC‑C‑CTGGCTCGCATTAAAAGTTATGcc                                                                  <  1:1821950/100‑1 (MQ=255)
                             tCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTAC‑CCCC‑C‑CTGGCTCGCATTAAAAGTTATGcc                                                                  >  1:2164890/1‑100 (MQ=255)
                               tGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTA‑‑CCCC‑C‑CTGGCTCGCATTAAAAGTTATGCCAg                                                                >  2:112166/1‑99 (MQ=255)
                               tGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTA‑‑CCCC‑C‑CTGGCTCGCATTAAAAGTTATGCCAGc                                                               <  2:2350224/100‑1 (MQ=255)
                               tGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTA‑‑CCCC‑C‑CTGGCTCGCATTAAAAGTTATGCCAGc                                                               >  1:2305040/1‑100 (MQ=255)
                                ggTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTAC‑CCCC‑C‑CTGGCTCGCATTAAAAGTTATGCCAGc                                                               >  1:581293/1‑100 (MQ=255)
                                ggTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTAC‑CCCC‑C‑CTGGCTCGCATTAAAAGTTATGCCAGc                                                               >  1:254265/1‑100 (MQ=255)
                                ggTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTA‑‑CCCC‑C‑CTGGCTCGCATTAAAAGTTATGCCAGCg                                                              <  2:1630680/100‑1 (MQ=255)
                                ggTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTA‑‑CCCC‑C‑CTGGCTCGCATTAAAAGTTATGCCAGCg                                                              >  2:474169/1‑100 (MQ=255)
                                ggTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTA‑‑CCCC‑C‑CTGGCTCGCATTAAAAGTTATGCCAGCg                                                              >  1:1799547/1‑100 (MQ=255)
                                  tATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTAC‑CCCC‑C‑CTGGCTCGCATTAAAAGTTATGCCAGCgg                                                             >  1:1864665/1‑100 (MQ=255)
                                   aTTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTAC‑CCCC‑C‑CTGGCTCGCATTAAAAGTTATGCCAGCGGt                                                            <  2:1580935/100‑1 (MQ=255)
                                    ttAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTAC‑CCCC‑C‑CTGGCTCGCATTAAAAGTTATGCCAGCGGTg                                                           >  1:1221457/1‑100 (MQ=255)
                                     tAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTAC‑CCCC‑C‑CTGGCTCGCATTAAAAGTTATGCCAGCGGTgg                                                          >  2:2008799/1‑100 (MQ=255)
                                      aaCGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTAC‑CCCC‑C‑CTGGCTCGCATTAAAAGTTATGCCAGCGGTGGc                                                         >  1:824158/1‑100 (MQ=255)
                                        cgacgGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTAC‑CCCC‑C‑CTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGt                                                       >  1:1290280/1‑100 (MQ=255)
                                         gacgGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTAC‑CCCC‑C‑CTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTg                                                      <  1:1341480/100‑1 (MQ=255)
                                         gacgGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTAC‑CCCC‑C‑CTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTg                                                      <  1:2217191/100‑1 (MQ=255)
                                         gacgGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTA‑‑CCCC‑C‑CTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGc                                                     >  2:118140/1‑100 (MQ=255)
                                           cgGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTAC‑CCCC‑C‑CTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGcc                                                    <  1:1378615/100‑1 (MQ=255)
                                               gctgcCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTA‑‑CCCC‑C‑CTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCg                                               <  1:2276710/100‑1 (MQ=255)
                                                 tgcCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTA‑‑CCCC‑C‑CTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCa                                             >  1:1768831/1‑100 (MQ=255)
                                                  gcCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTAC‑CCCC‑C‑CTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCa                                             <  1:2056403/100‑1 (MQ=255)
                                                  gcCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTAC‑CCCC‑C‑CTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCa                                             >  1:2185114/1‑100 (MQ=255)
                                                  gcCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTA‑‑CCCC‑C‑CTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCAt                                            >  1:2187741/1‑100 (MQ=255)
                                                  gcCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTA‑‑CCCC‑C‑CTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCAt                                            <  2:768257/100‑1 (MQ=255)
                                                   ccGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTA‑‑CCCC‑C‑CTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCAtt                                           <  1:615646/100‑1 (MQ=255)
                                                     gCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTAC‑CCCC‑C‑CTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTg                                          <  2:325207/100‑1 (MQ=255)
                                                        ctGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTAC‑CCCC‑C‑CTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATg                                       <  1:229970/100‑1 (MQ=255)
                                                        ctGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTAC‑CCCC‑C‑CTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATg                                       >  1:1789823/1‑100 (MQ=255)
                                                         tGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTA‑‑CCCC‑C‑CTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATggg                                     <  2:585022/100‑1 (MQ=255)
                                                         tGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTA‑‑CCCC‑C‑CTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATggg                                     <  2:491758/100‑1 (MQ=255)
                                                         tGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTA‑‑CCCC‑C‑CTGACTCGCATTAAAAGATATGCCAGCGGTGGCGTGCCCCCCGCATTGATggg                                     >  2:402222/1‑100 (MQ=255)
                                                          ggTGATTATGGAAGAATCTGCGGCGCTGGCAGCGGGCCTTA‑‑CCCC‑C‑CTGGCTCGCATTAAAAGTTATGCCGGCGGTGGCGTGCCCCCCGCATTGATGGGt                                    >  2:1419332/1‑100 (MQ=255)
                                                          ggTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTAC‑CCCC‑C‑CTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATggg                                     >  1:1416434/1‑100 (MQ=255)
                                                            tGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTAC‑CCCC‑C‑CTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTCATggg                                     >  2:1845973/1‑98 (MQ=255)
                                                            tGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTA‑‑CCCC‑C‑CTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTAt                                  >  1:2167513/1‑100 (MQ=255)
                                                             gATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTA‑‑CCCC‑C‑CTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATg                                 <  2:1416475/100‑1 (MQ=255)
                                                             gATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTA‑‑CCCC‑C‑CTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATg                                 >  2:1715390/1‑100 (MQ=255)
                                                              aTTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTA‑‑CCCC‑C‑CTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATgg                                <  2:289434/100‑1 (MQ=255)
                                                              aTTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTA‑‑CCCC‑C‑CTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATgg                                >  1:1190637/1‑100 (MQ=255)
                                                               ttATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTA‑‑CCCC‑C‑CTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATggg                               >  1:1111800/1‑100 (MQ=255)
                                                                tATGGAACAATCTGCGGCGCTGGCAGCAGGCCTTAC‑CCCC‑C‑CTAGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATggg                               <  2:998094/100‑1 (MQ=255)
                                                                 aTGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTAC‑CCCC‑C‑CTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATgggg                              >  2:859434/1‑100 (MQ=255)
                                                                  tGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTAC‑CCCC‑C‑CTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGc                             >  1:412526/1‑100 (MQ=255)
                                                                    gaagaaTCTGCGGCGCTGGCAGCAGGCCTTA‑‑CCCC‑C‑CTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAg                          >  2:927198/1‑100 (MQ=255)
                                                                    gaagaaTCTGCGGCGCTGGCAGCAGGCCTTA‑‑CCCC‑C‑CTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAg                          >  1:1409224/1‑100 (MQ=255)
                                                                      agaaTCTGCGGCGCTGGCAGCAGGCCTTAC‑CCCC‑C‑CTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGt                         >  1:1213826/1‑100 (MQ=255)
                                                                       gaaTCTGCGGCGCTGGCAGCAGGCCTTAC‑CCCC‑C‑CTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTa                        <  1:576319/100‑1 (MQ=255)
                                                                       gaaTCTGCGGCGCTGGCAGCAGGCCTTA‑‑CCCC‑C‑CTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTAc                       >  1:1681350/1‑100 (MQ=255)
                                                                          tCTGCGGCGCTGGCAGCAGGCCTTACCCCCC‑C‑CTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTAcc                      <  2:1303591/100‑1 (MQ=255)
                                                                          tCTGCGGCGCTGGCAGCAGGCCTTA‑‑CCCC‑C‑CTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTACCTg                    <  2:1501771/100‑1 (MQ=255)
                                                                          tCTGCGGCGCTGGCAGCAGGCCTTA‑‑CCCC‑C‑CTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTACCTg                    <  1:1572514/100‑1 (MQ=255)
                                                                           cTGCGGCGCTGGCAGCAGGCCTTAC‑CCCC‑C‑CTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTACCTg                    >  2:921483/1‑100 (MQ=255)
                                                                            tGCGGCGCTGGCAGCAGGCCTTA‑‑CCCC‑C‑CTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTCTGGGGCCAGTACCCGcc                  >  2:398390/1‑100 (MQ=255)
                                                                            tGCGGCGCTGGCAGCAGGCCTTA‑‑CCCC‑C‑CTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTACCTGcc                  >  1:319468/1‑100 (MQ=255)
                                                                               ggcgCTGGCAGCAGGCCTTA‑‑CCCC‑CGCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTACCTGc                   <  2:1351411/97‑1 (MQ=255)
                                                                               ggcgCTGGCAGCAGGCCTTA‑‑CCCC‑C‑CTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTACCTGCCACg               <  2:1799534/100‑1 (MQ=255)
                                                                                gcgcTGGCAGCAGGCCTTAC‑CCCC‑C‑CTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTACCTGCCACg               >  1:1533267/1‑100 (MQ=255)
                                                                                 cgcTGGCAGCAGGCCTTAC‑CCCC‑C‑CTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTACCTGCCACGc              >  1:1761221/1‑100 (MQ=255)
                                                                                     ggCAGCAGGCCTTAC‑CCCC‑C‑CTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTACCTGCCACGCaaaa          >  1:493885/1‑100 (MQ=255)
                                                                                       cagcagGCCTTAC‑CCCC‑C‑CTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTACCTGCCACGCAAAAAg        <  1:2263835/100‑1 (MQ=255)
                                                                                       cagcagGCCTTA‑‑CCCC‑C‑CTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTACCTGCCACGCAAAAAGc       >  2:162399/1‑100 (MQ=255)
                                                                                       cagcagGCCTTA‑‑CCCC‑C‑CTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTACCTGCCACGCAAAAAGc       <  2:1338971/100‑1 (MQ=255)
                                                                                        agcagGCCTTA‑‑CCCC‑C‑CTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTACCTGCCACGCAAAAAGCg      <  2:1768818/100‑1 (MQ=255)
                                                                                        agcagGCCTTA‑‑CCCC‑C‑CTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTACCTGCCACGCAAAAAGCg      >  1:2054263/1‑100 (MQ=255)
                                                                                         gcagGCCTTAC‑CCCC‑C‑CTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTACCTGCCACGCAAAAAGCg      <  2:309761/100‑1 (MQ=255)
                                                                                          cagGCCTTAC‑CCCC‑C‑CTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTACCTGCCACGCAAAAAGCGt     >  1:1826345/1‑100 (MQ=255)
                                                                                          cagGCCTTA‑‑CCCC‑C‑CTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTACCTGCCACGCAAAAAGCGtt    >  2:427896/1‑100 (MQ=255)
                                                                                             gCCTTAC‑CCCC‑C‑CTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTACCTGCCACGCAAAAAGCGTTAc  <  1:1485400/100‑1 (MQ=255)
                                                                                             gCCTTAC‑CCCC‑C‑CTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTACCTGCCACGCAAAAAGCGTTAc  <  2:2339800/100‑1 (MQ=255)
                                                                                                   |                                                                                                
AAGCAGGAACAGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTA‑‑CCCC‑C‑CTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTACCTGCCACGCAAAAAGCGTTAC  >  NC_000913/2326821‑2327012

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