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‑‑CCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTACCTGCCACGCAAAAAGCGTTAC  >  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‑‑cccccc                                                                                         >  2:884522/1‑100 (MQ=255)
       aaCAGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTA‑‑CCCCCCTg                                                                                       <  2:913500/100‑1 (MQ=255)
          aGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTCCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTA‑‑CCCCCCTGGCt                                                                                    <  2:331683/100‑1 (MQ=255)
           gTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTA‑‑CCCCCCTGGCTc                                                                                   >  2:1845267/1‑100 (MQ=255)
            tCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGTTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTA‑‑CCCCCCTGGCTCg                                                                                  <  2:462696/100‑1 (MQ=255)
             cACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTACCCCCCCTGGCTCg                                                                                  <  1:440897/100‑1 (MQ=255)
              aCCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTACCCCCCCTGGCTCGc                                                                                 <  2:1045731/100‑1 (MQ=255)
               ccGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTACCCCCCCTGGCTCGCa                                                                                >  1:961934/1‑100 (MQ=255)
               ccGCTGGGAACGCGTCTGGTATTAACGACCGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTACCCCCCCTGGCTCGCa                                                                                <  2:1412185/100‑1 (MQ=255)
                cGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTACCCCCCCTGGCTCGCAt                                                                               <  1:1229087/100‑1 (MQ=255)
                 gCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTACCCCCCCTGGCTCGCAtt                                                                              >  1:1482177/1‑100 (MQ=255)
                    gggAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTACCCCCCCTGGCTCGCATTaaa                                                                           >  1:994988/1‑100 (MQ=255)
                      gAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTACCCCCCCTGGCTCGCATTAAAAg                                                                         >  1:348982/1‑100 (MQ=255)
                      gAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTACCCCCCCTGGCTCGCATTAAAAg                                                                         >  1:575519/1‑100 (MQ=255)
                       aaCGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTACCCCCCCTGGCTCGCATTAAAAGt                                                                        <  1:1931096/100‑1 (MQ=255)
                       aaCGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTA‑‑CCCCCCTGGCTCGCATTAAAAGtt                                                                       >  1:171539/1‑100 (MQ=255)
                         cgcgTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTACCCCCCCTGGCTCGCATTAAAAGTTa                                                                      >  1:634000/1‑100 (MQ=255)
                         cgcgTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTACCCCCCCTGGCTCGCATTAAAAGTTa                                                                      <  2:1745521/100‑1 (MQ=255)
                          gcgTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTACCCCCCCTGGCTCGCATTAAAAGTTAt                                                                     <  1:797752/100‑1 (MQ=255)
                            gTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAGGAATCTGCGGCGCTGGCAGCAGGCCTTA‑‑CCCCGCCTGGCTCGCATTAAACGTGATGc                                                                   >  2:1831604/1‑100 (MQ=255)
                             tCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTACCCCCCCTGGCTCGCATTAAAAGTTATGcc                                                                  <  1:1821950/100‑1 (MQ=255)
                             tCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTACCCCCCCTGGCTCGCATTAAAAGTTATGcc                                                                  >  1:2164890/1‑100 (MQ=255)
                               tGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTA‑‑CCCCCCTGGCTCGCATTAAAAGTTATGCCAg                                                                >  2:112166/1‑99 (MQ=255)
                               tGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTA‑‑CCCCCCTGGCTCGCATTAAAAGTTATGCCAGc                                                               <  2:2350224/100‑1 (MQ=255)
                               tGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTA‑‑CCCCCCTGGCTCGCATTAAAAGTTATGCCAGc                                                               >  1:2305040/1‑100 (MQ=255)
                                ggTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTACCCCCCCTGGCTCGCATTAAAAGTTATGCCAGc                                                               >  1:581293/1‑100 (MQ=255)
                                ggTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTACCCCCCCTGGCTCGCATTAAAAGTTATGCCAGc                                                               >  1:254265/1‑100 (MQ=255)
                                ggTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTA‑‑CCCCCCTGGCTCGCATTAAAAGTTATGCCAGCg                                                              <  2:1630680/100‑1 (MQ=255)
                                ggTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTA‑‑CCCCCCTGGCTCGCATTAAAAGTTATGCCAGCg                                                              >  2:474169/1‑100 (MQ=255)
                                ggTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTA‑‑CCCCCCTGGCTCGCATTAAAAGTTATGCCAGCg                                                              >  1:1799547/1‑100 (MQ=255)
                                  tATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTACCCCCCCTGGCTCGCATTAAAAGTTATGCCAGCgg                                                             >  1:1864665/1‑100 (MQ=255)
                                   aTTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTACCCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGt                                                            <  2:1580935/100‑1 (MQ=255)
                                    ttAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTACCCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTg                                                           >  1:1221457/1‑100 (MQ=255)
                                     tAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTACCCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTgg                                                          >  2:2008799/1‑100 (MQ=255)
                                      aaCGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTACCCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGc                                                         >  1:824158/1‑100 (MQ=255)
                                        cgacgGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTACCCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGt                                                       >  1:1290280/1‑100 (MQ=255)
                                         gacgGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTACCCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTg                                                      <  1:1341480/100‑1 (MQ=255)
                                         gacgGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTACCCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTg                                                      <  1:2217191/100‑1 (MQ=255)
                                         gacgGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTA‑‑CCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGc                                                     >  2:118140/1‑100 (MQ=255)
                                           cgGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTACCCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGcc                                                    <  1:1378615/100‑1 (MQ=255)
                                               gctgcCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTA‑‑CCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCg                                               <  1:2276710/100‑1 (MQ=255)
                                                 tgcCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTA‑‑CCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCa                                             >  1:1768831/1‑100 (MQ=255)
                                                  gcCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTACCCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCa                                             <  1:2056403/100‑1 (MQ=255)
                                                  gcCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTACCCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCa                                             >  1:2185114/1‑100 (MQ=255)
                                                  gcCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTA‑‑CCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCAt                                            >  1:2187741/1‑100 (MQ=255)
                                                  gcCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTA‑‑CCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCAt                                            <  2:768257/100‑1 (MQ=255)
                                                   ccGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTA‑‑CCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCAtt                                           <  1:615646/100‑1 (MQ=255)
                                                     gCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTACCCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTg                                          <  2:325207/100‑1 (MQ=255)
                                                        ctGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTACCCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATg                                       <  1:229970/100‑1 (MQ=255)
                                                        ctGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTACCCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATg                                       >  1:1789823/1‑100 (MQ=255)
                                                         tGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTA‑‑CCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATggg                                     <  2:585022/100‑1 (MQ=255)
                                                         tGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTA‑‑CCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATggg                                     <  2:491758/100‑1 (MQ=255)
                                                         tGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTA‑‑CCCCCCTGACTCGCATTAAAAGATATGCCAGCGGTGGCGTGCCCCCCGCATTGATggg                                     >  2:402222/1‑100 (MQ=255)
                                                          ggTGATTATGGAAGAATCTGCGGCGCTGGCAGCGGGCCTTA‑‑CCCCCCTGGCTCGCATTAAAAGTTATGCCGGCGGTGGCGTGCCCCCCGCATTGATGGGt                                    >  2:1419332/1‑100 (MQ=255)
                                                          ggTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTACCCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATggg                                     >  1:1416434/1‑100 (MQ=255)
                                                            tGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTACCCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTCATggg                                     >  2:1845973/1‑98 (MQ=255)
                                                            tGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTA‑‑CCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTAt                                  >  1:2167513/1‑100 (MQ=255)
                                                             gATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTA‑‑CCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATg                                 <  2:1416475/100‑1 (MQ=255)
                                                             gATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTA‑‑CCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATg                                 >  2:1715390/1‑100 (MQ=255)
                                                              aTTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTA‑‑CCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATgg                                <  2:289434/100‑1 (MQ=255)
                                                              aTTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTA‑‑CCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATgg                                >  1:1190637/1‑100 (MQ=255)
                                                               ttATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTA‑‑CCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATggg                               >  1:1111800/1‑100 (MQ=255)
                                                                tATGGAACAATCTGCGGCGCTGGCAGCAGGCCTTACCCCCCCTAGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATggg                               <  2:998094/100‑1 (MQ=255)
                                                                 aTGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTACCCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATgggg                              >  2:859434/1‑100 (MQ=255)
                                                                  tGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTACCCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGc                             >  1:412526/1‑100 (MQ=255)
                                                                    gaagaaTCTGCGGCGCTGGCAGCAGGCCTTA‑‑CCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAg                          >  2:927198/1‑100 (MQ=255)
                                                                    gaagaaTCTGCGGCGCTGGCAGCAGGCCTTA‑‑CCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAg                          >  1:1409224/1‑100 (MQ=255)
                                                                      agaaTCTGCGGCGCTGGCAGCAGGCCTTACCCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGt                         >  1:1213826/1‑100 (MQ=255)
                                                                       gaaTCTGCGGCGCTGGCAGCAGGCCTTACCCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTa                        <  1:576319/100‑1 (MQ=255)
                                                                       gaaTCTGCGGCGCTGGCAGCAGGCCTTA‑‑CCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTAc                       >  1:1681350/1‑100 (MQ=255)
                                                                          tCTGCGGCGCTGGCAGCAGGCCTTACCCCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTAcc                      <  2:1303591/100‑1 (MQ=255)
                                                                          tCTGCGGCGCTGGCAGCAGGCCTTA‑‑CCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTACCTg                    <  2:1501771/100‑1 (MQ=255)
                                                                          tCTGCGGCGCTGGCAGCAGGCCTTA‑‑CCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTACCTg                    <  1:1572514/100‑1 (MQ=255)
                                                                           cTGCGGCGCTGGCAGCAGGCCTTACCCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTACCTg                    >  2:921483/1‑100 (MQ=255)
                                                                            tGCGGCGCTGGCAGCAGGCCTTA‑‑CCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTCTGGGGCCAGTACCCGcc                  >  2:398390/1‑100 (MQ=255)
                                                                            tGCGGCGCTGGCAGCAGGCCTTA‑‑CCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTACCTGcc                  >  1:319468/1‑100 (MQ=255)
                                                                               ggcgCTGGCAGCAGGCCTTA‑‑CCCCCGCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTACCTGc                   <  2:1351411/97‑1 (MQ=255)
                                                                               ggcgCTGGCAGCAGGCCTTA‑‑CCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTACCTGCCACg               <  2:1799534/100‑1 (MQ=255)
                                                                                gcgcTGGCAGCAGGCCTTACCCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTACCTGCCACg               >  1:1533267/1‑100 (MQ=255)
                                                                                 cgcTGGCAGCAGGCCTTACCCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTACCTGCCACGc              >  1:1761221/1‑100 (MQ=255)
                                                                                     ggCAGCAGGCCTTACCCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTACCTGCCACGCaaaa          >  1:493885/1‑100 (MQ=255)
                                                                                       cagcagGCCTTACCCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTACCTGCCACGCAAAAAg        <  1:2263835/100‑1 (MQ=255)
                                                                                       cagcagGCCTTA‑‑CCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTACCTGCCACGCAAAAAGc       >  2:162399/1‑100 (MQ=255)
                                                                                       cagcagGCCTTA‑‑CCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTACCTGCCACGCAAAAAGc       <  2:1338971/100‑1 (MQ=255)
                                                                                        agcagGCCTTA‑‑CCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTACCTGCCACGCAAAAAGCg      <  2:1768818/100‑1 (MQ=255)
                                                                                        agcagGCCTTA‑‑CCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTACCTGCCACGCAAAAAGCg      >  1:2054263/1‑100 (MQ=255)
                                                                                         gcagGCCTTACCCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTACCTGCCACGCAAAAAGCg      <  2:309761/100‑1 (MQ=255)
                                                                                          cagGCCTTACCCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTACCTGCCACGCAAAAAGCGt     >  1:1826345/1‑100 (MQ=255)
                                                                                          cagGCCTTA‑‑CCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTACCTGCCACGCAAAAAGCGtt    >  2:427896/1‑100 (MQ=255)
                                                                                             gCCTTACCCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTACCTGCCACGCAAAAAGCGTTAc  <  1:1485400/100‑1 (MQ=255)
                                                                                             gCCTTACCCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTACCTGCCACGCAAAAAGCGTTAc  <  2:2339800/100‑1 (MQ=255)
                                                                                                   |                                                                                                
AAGCAGGAACAGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTA‑‑CCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTACCTGCCACGCAAAAAGCGTTAC  >  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: