Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,452,057 A→T 100% I77N (ATC→AAC)  rplD ← 50S ribosomal subunit protein L4

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,452,0570AT93.3% 286.4 / 8.2 92I77N (ATC→AAC) rplD50S ribosomal subunit protein L4
Reads supporting (aligned to +/- strand):  ref base A (4/2);  new base T (46/39);  total (51/41)
Fisher's exact test for biased strand distribution p-value = 6.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.29e-01

CAGGATGCTTTTCAGCGCGCCGCGGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGATCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGCACGAGTCTTCTGAG  >  NC_000913/3451960‑3452153
                                                                                                 |                                                                                                
cAGGATGCTTTTCAGCGCGCCGCTGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTc                                                                                                <  2:620105/100‑1 (MQ=255)
  ggATGCTTTTCAGCGCGCCGCGGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCg                                                                                               >  1:225587/1‑99 (MQ=255)
  atctGCTTTTCAGCGCGCCGCGGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGATGAGCAGCAAAGGTCACGCCACCAGAACGCCAGATCgg                                                                                              <  2:1366161/97‑1 (MQ=255)
    aTGCTTTTCAGCGCGCCGCGGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGc                                                                                            <  1:2310796/100‑1 (MQ=255)
     tGCTTTTCAGCGCGCCGCGGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGct                                                                                           >  1:78197/1‑100 (MQ=255)
     tGCTTTTCAGCGCGCCGCGGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGct                                                                                           >  1:1671711/1‑100 (MQ=255)
      gCTTTTCAGCGCGCCGCGGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGATCggagat                                                                                          >  1:1366170/1‑96 (MQ=255)
       cTTTTCAGCGCGCCGCGGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGctct                                                                                         <  2:73922/100‑1 (MQ=255)
       cTTTTCAGCGCGCCGCGGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGATCGGGctct                                                                                         >  2:2110178/1‑100 (MQ=255)
        ttttCAGCGCGCCGCGGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGGCACGCCACCAGAACGCCAGTTCGGGCTCtt                                                                                        >  1:2324762/1‑100 (MQ=255)
           tCAGCGCGCCGCGGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGAt                                                                                     <  1:2003246/100‑1 (MQ=255)
           tCAGCGCGCCGCGGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGATCGGGCTCTTGAt                                                                                     <  1:2110192/100‑1 (MQ=255)
             aGCGCGCCGCGGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAg                                                                                   >  1:2105605/1‑100 (MQ=255)
              gcgcgcCGCGGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCACAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGa                                                                                  >  2:333551/1‑100 (MQ=255)
                   ccgcGGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAg                                                                             <  2:1482354/100‑1 (MQ=255)
                    cgcgGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGa                                                                            >  1:414503/1‑100 (MQ=255)
                    cgcgGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGa                                                                            >  1:2117498/1‑100 (MQ=255)
                      cgGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAAc                                                                          <  2:2105591/100‑1 (MQ=255)
                       ggTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAAcg                                                                         >  1:466302/1‑100 (MQ=255)
                         tACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAAcgcg                                                                       >  2:735045/1‑100 (MQ=255)
                          aCATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACgcgc                                                                      <  2:2229238/100‑1 (MQ=255)
                          aCATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACgcgc                                                                      <  2:465144/100‑1 (MQ=255)
                           cATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCa                                                                     <  2:355150/100‑1 (MQ=255)
                           cATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCa                                                                     >  2:2364897/1‑100 (MQ=255)
                             tcttctTGTTAACTTTTTTTTTTTGTTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACg                                                                   >  2:1712831/1‑100 (MQ=255)
                             tcttctTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACg                                                                   <  2:2292438/100‑1 (MQ=255)
                              cttcttGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACg                                                                   <  1:2448466/99‑1 (MQ=255)
                                  ttGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCgg                                                              >  1:2348965/1‑100 (MQ=255)
                                   tGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGATCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGt                                                             >  1:1005519/1‑100 (MQ=255)
                                    gTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTg                                                            <  2:1287874/100‑1 (MQ=255)
                                    gTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTg                                                            >  2:386412/1‑100 (MQ=255)
                                      tAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGcc                                                          <  1:1033461/100‑1 (MQ=255)
                                      tAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGcc                                                          >  2:1604222/1‑100 (MQ=255)
                                      tAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGATCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGcc                                                          >  1:1448133/1‑100 (MQ=255)
                                       aaCTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCt                                                         >  2:2146558/1‑100 (MQ=255)
                                         cTTTTTGACTTTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCtt                                                        <  2:2303951/99‑1 (MQ=255)
                                         cTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCttt                                                       >  1:1206539/1‑100 (MQ=255)
                                         cTTTTTGACTCTGGTCCTGCCGACGCTCAGCCAATGTCACGCCACCACACCGCCAGCTCGGGCCGTTGATAGAACCAGAACGAACacggacgggcact                                                         >  2:1359793/1‑93 (MQ=255)
                                          tttttGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTc                                                      <  2:308845/100‑1 (MQ=255)
                                               gACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCAACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGgcg                                                 <  2:448611/100‑1 (MQ=255)
                                                 cTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGcc                                               >  2:1597412/1‑100 (MQ=255)
                                                   gtgGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCAc                                             >  1:1432108/1‑100 (MQ=255)
                                                   gtgGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCAc                                             <  1:1307327/100‑1 (MQ=255)
                                                      gTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGt                                          >  2:1024336/1‑100 (MQ=255)
                                                       tCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCCCGGCCGGTGCCTTTCTGGCGCCACGGtt                                         >  2:530110/1‑100 (MQ=255)
                                                        ccTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGttt                                        <  1:2468388/100‑1 (MQ=255)
                                                         cTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGtttt                                       <  2:1465496/100‑1 (MQ=255)
                                                           gCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGtttttt                                     >  2:506414/1‑100 (MQ=255)
                                                            cGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTa                                    <  1:1308098/100‑1 (MQ=255)
                                                            cGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTa                                    >  1:2484097/1‑100 (MQ=255)
                                                              gACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTTATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTAcc                                  <  1:2444492/100‑1 (MQ=255)
                                                               aCGAGCCTCAAAGGTCACGGCACCAGATCGCGAGTTCGGGCTCTGGATAGAGCCGGAACGCGCACGGCCGGTGCGTTTCTGGCGCCACGGTTTTTTACCg                                 <  2:1649831/100‑1 (MQ=255)
                                                               aCGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCg                                 <  1:1994922/100‑1 (MQ=255)
                                                               aCGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCg                                 <  2:1651706/100‑1 (MQ=255)
                                                               aCGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCg                                 <  1:456548/100‑1 (MQ=255)
                                                               aCGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTGCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCg                                 <  2:901392/100‑1 (MQ=255)
                                                                  agcagcAAAGGTCACGCCACCAGAACACCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGaa                              >  1:2482114/1‑100 (MQ=255)
                                                                   gcagcaAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAAc                             >  1:415591/1‑100 (MQ=255)
                                                                   gcagcaAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAAc                             <  1:1502513/100‑1 (MQ=255)
                                                                    cagcaAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAAcc                            >  2:2433692/1‑100 (MQ=255)
                                                                    cagcaAAGGTCACGCCACCAGAACACCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAAcc                            <  2:2482100/100‑1 (MQ=255)
                                                                      gcaAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAg                          >  1:1173237/1‑100 (MQ=255)
                                                                      gcaAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGCTTTTTACCGGAACCAg                          >  2:1777074/1‑100 (MQ=255)
                                                                       caAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGTGCACGGCCGGTGCCTTTCTGGCGTCACGGTTTTTTACCGGAACCAGt                         <  2:6385/100‑1 (MQ=255)
                                                                       caAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGt                         >  1:435980/1‑100 (MQ=255)
                                                                       caAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTACTGGCGCCACGGTTTTTTACCGGAACCAGt                         <  2:1992550/100‑1 (MQ=255)
                                                                        aaaGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGtt                        <  1:334337/100‑1 (MQ=255)
                                                                        aaaGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACACGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGtt                        >  1:1760350/1‑100 (MQ=255)
                                                                         aaGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTa                       >  1:1330747/1‑100 (MQ=255)
                                                                          aGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTAc                      <  1:1712843/100‑1 (MQ=255)
                                                                           ggTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACt                     >  1:1329663/1‑100 (MQ=255)
                                                                           ggTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACt                     >  1:2506506/1‑100 (MQ=255)
                                                                              cACGCCACCAGAACGCCAGTACGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCa                  >  2:214977/1‑100 (MQ=255)
                                                                                 gCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGCa               <  2:1175410/100‑1 (MQ=255)
                                                                                 gCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGCa               <  2:691023/100‑1 (MQ=255)
                                                                                  ccaccaGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGCAc              >  1:200123/1‑100 (MQ=255)
                                                                                  ccaccaGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGCAc              <  2:1952283/100‑1 (MQ=255)
                                                                                   caccaGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGCACg             >  2:1646015/1‑100 (MQ=255)
                                                                                   caccaGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGCACg             >  1:1128203/1‑100 (MQ=255)
                                                                                   caccaGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGCACg             >  2:846954/1‑100 (MQ=255)
                                                                                    accaGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGCACGa            >  2:2454090/1‑100 (MQ=255)
                                                                                    accaGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGCACGa            >  2:2126829/1‑100 (MQ=255)
                                                                                       aGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTGAGCACGAGtc         >  2:1680511/1‑100 (MQ=255)
                                                                                          aCGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGGGCCTTTCTGGCGCCACGGATTTTTACCGGAACCAGTTACTTCAGCACGAgattt       >  2:1647568/1‑95 (MQ=255)
                                                                                            gCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACTTTATCCTGTTACTTCAGCACGAGtcttct     >  2:593904/1‑99 (MQ=255)
                                                                                            gCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGCACGAGtcttc      <  2:1549121/98‑1 (MQ=255)
                                                                                            gCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGCACGAGTCTTCTg    >  2:2176617/1‑100 (MQ=255)
                                                                                            gCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGCACGAGTCTTCTg    <  2:132360/100‑1 (MQ=255)
                                                                                             ccAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGACCCAGTTACTTCAGCACGAGtcttc      >  2:2262873/1‑97 (MQ=255)
                                                                                             ccAGTTCGGGCTCTGGATAGAACCAGAACGCGCACGGCCGGTGCCGTGCTTGCGCCACGGTCGTTTACCGGCACCAGTTACTTCAGCACGAGTCTTCTg    >  2:169921/1‑99 (MQ=255)
                                                                                              cAGTTCGGGCTCTTGATGGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGCACGAGTCTTCTGAg  <  2:837777/100‑1 (MQ=255)
                                                                                              cAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGCACGAGTCTTCTGAg  <  1:827085/100‑1 (MQ=255)
                                                                                                 |                                                                                                
CAGGATGCTTTTCAGCGCGCCGCGGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGATCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGCACGAGTCTTCTGAG  >  NC_000913/3451960‑3452153

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: