Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 770,682 A→C 16.1% intergenic (+71/‑776) mngB → / → cydA alpha‑mannosidase/cytochrome d terminal oxidase, subunit I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913770,6820AC16.1% 362.2 / 62.4 141intergenic (+71/‑776)mngB/cydAalpha‑mannosidase/cytochrome d terminal oxidase, subunit I
Reads supporting (aligned to +/- strand):  ref base A (48/70);  new base C (13/10);  total (61/80)
Fisher's exact test for biased strand distribution p-value = 1.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.25e-01

Only 100 of 143 total aligned reads displayed.

ACGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGC  >  NC_000913/770584‑770781
                                                                                                  |                                                                                                   
aCGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATAc                                                                                                    <  1:3502802/100‑1 (MQ=255)
aCGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATAc                                                                                                    >  2:430769/1‑100 (MQ=255)
  ggACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACgg                                                                                                  <  2:2535590/100‑1 (MQ=255)
   gACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGc                                                                                                 <  1:566672/100‑1 (MQ=255)
   gACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGc                                                                                                 <  2:2731571/100‑1 (MQ=255)
    aCGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCa                                                                                                <  2:2634361/100‑1 (MQ=255)
     cGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCa                                                                                                <  2:2657338/99‑1 (MQ=255)
       ttCAGTTACCGGCTTGGCTGAATATCAATCAAACCGATGCCACATATGCGCGGCCAGATTGTTGACAGAGGGCGCTTTGTTCATGCCGGATACTGCATGa                                                                                             <  1:2780845/100‑1 (MQ=38)
       ttCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGa                                                                                             >  2:3333222/1‑100 (MQ=255)
       cttccgatctggCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGTATCCGGCATGa                                                                                             <  1:1431647/90‑1 (MQ=38)
         cAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGTATCCGGCATGAAc                                                                                           <  2:833601/100‑1 (MQ=255)
         cAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGTATCCGGCATGAAc                                                                                           <  2:3357859/100‑1 (MQ=255)
         cAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGTATCCGGCATGAAc                                                                                           <  2:1874300/100‑1 (MQ=255)
         cAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGTATCCGGCATGAAc                                                                                           <  2:74945/100‑1 (MQ=255)
         cAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGGTCATGCCGGATACGGCATGAAc                                                                                           <  1:2065339/100‑1 (MQ=255)
         cAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATACGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGTATCCGGCATGAAc                                                                                           <  1:770224/100‑1 (MQ=255)
          aGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACg                                                                                          >  2:421539/1‑100 (MQ=255)
           gTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGc                                                                                         >  1:3066141/1‑100 (MQ=255)
            ttACCTGCTTGCCTGAATAGCAATCAATCCTAAGCAACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGTATGAACg                                                                                          <  1:1434754/98‑1 (MQ=38)
             tACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTCGTGCATGCCGTATCCGGCATGAAcaaag                                                                                       <  1:1568113/100‑5 (MQ=38)
              aCCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATACCGGATACGGCATGAACGCttt                                                                                      <  2:3099891/100‑1 (MQ=255)
               ccGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCACGAACGCTTTa                                                                                     >  2:1645166/1‑100 (MQ=255)
                 ggCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGTATCCGGCATGaagatcggaag                                                                                   >  2:1431654/1‑91 (MQ=38)
                  gCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTAttn                                                                                  <  2:2944020/100‑2 (MQ=255)
                  gCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTc                                                                                  <  1:556181/100‑1 (MQ=255)
                  gCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTc                                                                                  <  1:168408/100‑1 (MQ=255)
                   cTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCg                                                                                 <  2:1281069/100‑1 (MQ=255)
                    ttGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCgg                                                                                <  2:2997342/100‑1 (MQ=255)
                    ttGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCgg                                                                                <  1:1831020/100‑1 (MQ=255)
                    ttGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCgg                                                                                <  1:2251297/100‑1 (MQ=255)
                     tGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCg                                                                                 >  1:2712575/1‑98 (MQ=255)
                      gCCTGAATAGCAATCAAATCGAAGCCGCATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTTTTTGGt                                                                               >  1:188648/1‑99 (MQ=38)
                       ccTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCt                                                                             <  2:1989724/100‑1 (MQ=255)
                       ccTGAATAGCAAGCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGCTACGGCATGAACGCTTTATTCGGTc                                                                              >  1:1482768/1‑99 (MQ=255)
                        cTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTa                                                                            <  1:2872294/100‑1 (MQ=255)
                          gAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACa                                                                          >  2:2360299/1‑100 (MQ=255)
                          gAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACa                                                                          <  2:233843/100‑1 (MQ=255)
                          gAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACa                                                                          <  1:109492/100‑1 (MQ=255)
                              aGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAg                                                                      >  1:700187/1‑100 (MQ=255)
                               gCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGc                                                                     >  1:2094442/1‑100 (MQ=255)
                                cAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCa                                                                    <  2:2038673/100‑1 (MQ=255)
                                cAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCa                                                                    <  2:1484623/100‑1 (MQ=255)
                                 aaTCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAg                                                                   >  2:2757032/1‑100 (MQ=255)
                                   tCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGc                                                                 <  1:2321259/100‑1 (MQ=255)
                                   tCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGc                                                                 <  1:2817144/100‑1 (MQ=255)
                                   tCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGc                                                                 <  1:1202085/100‑1 (MQ=255)
                                       aCCGAAGCCACATACGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAAt                                                             >  1:672946/1‑100 (MQ=255)
                                         cGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTc                                                           <  2:2665465/100‑1 (MQ=255)
                                             gCCACATATGCGCGGCCAGATTGTTGCCAAAGGGCGCTTTGTTCATGCCGGATACGGCATAAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAat                                                        <  1:1213318/99‑1 (MQ=255)
                                             gCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAata                                                       <  2:2684356/100‑1 (MQ=255)
                                              ccACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAatat                                                      <  2:484329/100‑1 (MQ=255)
                                                acaTATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAatatat                                                    <  2:1132036/100‑1 (MQ=255)
                                                  atatGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATAtt                                                   >  1:1451117/1‑99 (MQ=255)
                                                   tatGCGCGGCCAGATTGTTGACAAAGTGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATAtt                                                   <  1:2231433/98‑1 (MQ=255)
                                                   tatGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAATGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGc                                                 >  1:2494082/1‑100 (MQ=255)
                                                    atGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCa                                                <  2:80896/100‑1 (MQ=255)
                                                       cgcgGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCagag                                             <  2:672949/100‑1 (MQ=255)
                                                         cgGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAgat                                           >  1:2118252/1‑100 (MQ=255)
                                                           gCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAgatga                                         >  2:474518/1‑100 (MQ=255)
                                                           gCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAgatga                                         <  2:2890522/100‑1 (MQ=255)
                                                           gCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAgat                                           >  1:1121803/1‑98 (MQ=255)
                                                                 ttgttgACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAg                                    <  1:951537/99‑1 (MQ=255)
                                                                   gttgACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGaca                                 >  2:2624686/1‑100 (MQ=255)
                                                                   gttgACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGaca                                 >  2:2797675/1‑100 (MQ=255)
                                                                   gttgACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGac                                  >  1:2713471/1‑99 (MQ=255)
                                                                   gttgACAAAGGGCGCTTTGTTCATACCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGaca                                 >  1:1237102/1‑100 (MQ=255)
                                                                        cAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGc                            >  2:1309336/1‑100 (MQ=255)
                                                                        cAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGc                            >  2:3433238/1‑100 (MQ=255)
                                                                            gggCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGc                        >  1:749656/1‑100 (MQ=255)
                                                                             ggCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTCTTCGGTCTACACAAGCAGGTAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGc                        <  1:1935182/99‑1 (MQ=255)
                                                                              gcgcTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGt                      >  2:3231920/1‑100 (MQ=255)
                                                                              gcgcTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGt                      >  2:3041178/1‑100 (MQ=255)
                                                                               ccctctcttCATGCCGTATCCGGCATGAACGCTTTATTAGGTCTACCAAAGAAGCCAAATTCAATATATTGCAGGGATGATGTAGACACTGGCAAg                         <  1:294409/89‑1 (MQ=25)
                                                                                gcTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAg                    <  1:2948606/100‑1 (MQ=255)
                                                                                 cTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAgc                   >  1:1797944/1‑100 (MQ=255)
                                                                                 cTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAgc                   <  1:582380/100‑1 (MQ=255)
                                                                                 cTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAgc                   <  2:3202182/100‑1 (MQ=255)
                                                                                 cTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAgc                   >  2:1648085/1‑100 (MQ=255)
                                                                                 aagcgTTCATGCCGTATCCGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAgc                   >  2:220242/5‑100 (MQ=37)
                                                                                 aagcgTTCATGCCGTATCCGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAgc                   >  1:199340/5‑100 (MQ=37)
                                                                                  tttGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAgcg                  <  2:2854959/100‑1 (MQ=255)
                                                                                    cgTTCATGCCGTATCCGGCATGAACGCTTTATTCGGTCTACAACAGCAGGCAACTTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCa                >  1:3451008/2‑100 (MQ=25)
                                                                                    cgTTCATGCCGTATCCGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCa                >  2:197779/2‑100 (MQ=255)
                                                                                     gTTCATGCCGTATCCGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCAt               >  1:2171622/1‑100 (MQ=255)
                                                                                     gTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCAt               >  1:261479/1‑100 (MQ=255)
                                                                                     gTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCAt               >  2:463384/1‑100 (MQ=255)
                                                                                        cATGCCGTATCCGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATGTATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAg            >  2:3034009/1‑100 (MQ=37)
                                                                                        cATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCa             <  1:321997/99‑1 (MQ=255)
                                                                                        cATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAg            >  2:146295/1‑100 (MQ=255)
                                                                                         aTGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCa             <  1:1131292/98‑1 (MQ=255)
                                                                                         aTGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAgg           >  2:2683893/1‑100 (MQ=255)
                                                                                          tGCCGTATCCGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGt          >  2:3282414/1‑100 (MQ=255)
                                                                                           gCCGTATCCGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTa         >  1:967572/1‑100 (MQ=255)
                                                                                           gCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTa         <  1:3333198/100‑1 (MQ=255)
                                                                                           gCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTa         >  1:3388116/1‑100 (MQ=255)
                                                                                            ccGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTaa        >  2:1895412/1‑100 (MQ=255)
                                                                                             cGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAAt       <  1:3544576/100‑1 (MQ=255)
                                                                                                aTACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAAtttt    <  2:1385438/100‑1 (MQ=255)
                                                                                                  aCGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTGATTTTgn  <  2:2115286/100‑2 (MQ=255)
                                                                                                  aCGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGc  <  2:1570951/100‑1 (MQ=255)
                                                                                                  |                                                                                                   
ACGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGC  >  NC_000913/770584‑770781

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