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,0570AT97.5% 531.7 / 0.7 157I77N (ATC→AAC) rplD50S ribosomal subunit protein L4
Reads supporting (aligned to +/- strand):  ref base A (1/3);  new base T (74/79);  total (75/82)
Fisher's exact test for biased strand distribution p-value = 6.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

Only 100 of 157 total aligned reads displayed.

AGGATGCTTTTCAGCGCGCCGCGGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGATCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGCACGAGTCTTCTG  >  NC_000913/3451961‑3452151
                                                                                                |                                                                                              
aGGATGCTTTTCAGCGCGCCGCGGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCg                                                                                             >  1:1740178/1‑100 (MQ=255)
aGGATGCTTTTCAACGCGCCGCGGTACATCTTCTTGTTAACTGTTTGACTGTGGGCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCg                                                                                             >  2:2296134/1‑100 (MQ=255)
 ggATGCTTTTCAGCGCGCCGCGGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCgg                                                                                            <  1:84736/100‑1 (MQ=255)
 ggATGCTTTTCAGCGCGCCGCGGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCgg                                                                                            >  2:3053396/1‑100 (MQ=255)
   aTGCTTTTCAGCGCGCCGCGGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAATGTCACGCCACCAGAACGCCAGTTCGGGc                                                                                          >  1:828229/1‑100 (MQ=255)
    tGCTTTTCAGCGCGCCGCGGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGct                                                                                         <  2:544735/100‑1 (MQ=255)
      cTTTTCAGCGCGCCGCGGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGctct                                                                                       >  2:1348731/1‑100 (MQ=255)
      cTTTTCAGCGCGCCGCGGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGctct                                                                                       <  1:2894067/100‑1 (MQ=255)
      cTTTTCAGCGCGCCGCGGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGATCGGGctct                                                                                       >  1:2806292/1‑100 (MQ=255)
        tttCAGCGCGCCGCGGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTg                                                                                     <  1:450211/100‑1 (MQ=255)
        tttCAGCGCGCCGCGGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTg                                                                                     >  1:2660982/1‑100 (MQ=255)
        tttCAGCGCGCCGCGGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTg                                                                                     <  2:2175735/100‑1 (MQ=255)
          tCAGCGCGCCGCGGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGAt                                                                                   <  2:1941057/100‑1 (MQ=255)
          tCAGCGCGCCGCGGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGAt                                                                                   <  2:215865/100‑1 (MQ=255)
          tCAGCGCGCCGCGGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGAt                                                                                   <  1:1761805/100‑1 (MQ=255)
          tCAGCGCGCCGCGGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGAt                                                                                   >  2:664466/1‑100 (MQ=255)
          tCAGCGCGCCGCGGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGATCGGGCTCTTGAt                                                                                   <  1:3096012/100‑1 (MQ=255)
           cAGCGCGCCGCGGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATa                                                                                  <  2:3380836/100‑1 (MQ=255)
               gcgcCGCGGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCGCGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAAc                                                                              <  2:1490162/100‑1 (MQ=255)
                   cgcgGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGa                                                                          <  2:3263470/100‑1 (MQ=255)
                     cgGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAAc                                                                        <  2:727905/100‑1 (MQ=255)
                      ggTACATCTTCTTGTTAACTTTTTGAGTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAAcg                                                                       <  1:246772/100‑1 (MQ=255)
                      ggTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAAcg                                                                       <  2:415247/100‑1 (MQ=255)
                         aCATCTTCTTGTTAACTTTTTTTTTTTTTTTTTTCGGACGAGCAGCAAAGGTCACGTTACTAGAACGCCAGTTCGGGTTCTTGATTGAAc                                                                              >  1:2470538/1‑90 (MQ=255)
                          cATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCa                                                                   <  1:2265681/100‑1 (MQ=255)
                          cATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCa                                                                   <  1:204865/100‑1 (MQ=255)
                          cATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCa                                                                   >  2:2303097/1‑100 (MQ=255)
                           aTCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCAc                                                                  >  1:229008/1‑100 (MQ=255)
                             cttcttGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACgg                                                                >  2:1253314/1‑100 (MQ=255)
                               tcttGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGcc                                                              >  1:2389017/1‑100 (MQ=255)
                                  tGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGt                                                           >  2:2586043/1‑100 (MQ=255)
                                     tAACTTTTTGACTGTGGTCCTGCGGACGTCCCGCACATCTCCCGACCCTAGAACACAAATTCTGGCTCTTGATAGAACCAGAACGCGCACAGCCGGTAcc                                                        >  1:940544/1‑100 (MQ=255)
                                     tAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGcc                                                        >  2:1226709/1‑100 (MQ=255)
                                     tAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGcc                                                        <  1:1238904/100‑1 (MQ=255)
                                        cTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCttt                                                     >  2:1870759/1‑100 (MQ=255)
                                        cTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCtt                                                      >  1:1029313/1‑99 (MQ=255)
                                        cTTTTTGACTGTGGTCCTGCGGACGAGCAACAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCttt                                                     <  1:1939097/100‑1 (MQ=255)
                                         tttttGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTc                                                    >  2:1207129/1‑100 (MQ=255)
                                          ttttGACTGTGGTGCGGCGGACGAGCAGCAAAGGTCACGCGACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCtt                                                      <  1:2717016/97‑1 (MQ=255)
                                          ttttGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCt                                                   <  2:252682/100‑1 (MQ=255)
                                          ttttGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCt                                                   >  2:847582/1‑100 (MQ=255)
                                            ttGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTgg                                                 >  2:1744821/1‑100 (MQ=255)
                                              gACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGgcg                                               >  1:782851/1‑100 (MQ=255)
                                              gACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGATCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGgcg                                               <  2:3314919/100‑1 (MQ=255)
                                               aCTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGgcgc                                              <  2:2576218/100‑1 (MQ=255)
                                                cTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGcc                                             >  2:1887175/1‑100 (MQ=255)
                                                cTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGcc                                             >  1:2451039/1‑100 (MQ=255)
                                                    ggTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACgg                                         >  1:1145021/1‑100 (MQ=255)
                                                     gTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGt                                        >  2:3278407/1‑100 (MQ=255)
                                                     gTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGt                                        >  1:1586043/1‑100 (MQ=255)
                                                      tCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGtt                                       >  2:1544009/1‑100 (MQ=255)
                                                       ccTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCAGGTGCCTTTCTGGCGCCACGGttt                                      >  2:808326/1‑100 (MQ=255)
                                                        cTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGtttt                                     <  1:3069195/100‑1 (MQ=255)
                                                        cTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGtttt                                     <  1:1331236/100‑1 (MQ=255)
                                                        cTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGttt                                      >  1:406471/1‑99 (MQ=255)
                                                         tGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGttttt                                    <  2:2096643/100‑1 (MQ=255)
                                                          gCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGtttttt                                   >  2:2869597/1‑100 (MQ=255)
                                                          gCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGtttttt                                   <  1:2563608/100‑1 (MQ=255)
                                                            ggAAGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTTGCGCCACGGTTTTTTa                                  >  1:2014737/1‑99 (MQ=255)
                                                             gACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTAcc                                >  2:2720370/1‑100 (MQ=255)
                                                             gACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTAcc                                <  2:1586050/100‑1 (MQ=255)
                                                              aCGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCg                               <  2:782854/100‑1 (MQ=255)
                                                              aCGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCg                               <  2:598078/100‑1 (MQ=255)
                                                               cGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCgg                              >  2:2062937/1‑100 (MQ=255)
                                                                gAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGa                             >  1:2342164/1‑100 (MQ=255)
                                                                  gcagcaAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAAc                           <  1:432985/100‑1 (MQ=255)
                                                                  gcagcaAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAAc                           <  1:1362545/100‑1 (MQ=255)
                                                                   cagcaAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAAcc                          <  2:1147615/100‑1 (MQ=255)
                                                                    agcaAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCa                         >  2:1600005/1‑100 (MQ=255)
                                                                     gcaAAGGTCACGCCACCAGCACGCCAGTTCGGGCGCTTGATAGAACCAGAACGCGCACGGCCGGTGCCGTTCTGGCGCCACGGTTTTTTACCGGAACCAg                        <  1:2805540/100‑1 (MQ=255)
                                                                     gcaAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAg                        >  2:1314001/1‑100 (MQ=255)
                                                                     gcaAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAg                        <  1:3342813/100‑1 (MQ=255)
                                                                      caAAGGTCACGCCACCGGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGt                       <  1:1455232/100‑1 (MQ=255)
                                                                       aaaTGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGtt                      <  1:2881874/100‑1 (MQ=255)
                                                                       aaaGGTCACGCCACCAGAACGCCAGTTCGTGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGtt                      <  2:1259025/100‑1 (MQ=255)
                                                                       aaaGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGtt                      >  2:135669/1‑100 (MQ=255)
                                                                       aaaGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGtt                      >  2:1416698/1‑100 (MQ=255)
                                                                         aGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTAc                    <  2:1969385/100‑1 (MQ=255)
                                                                          ggTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACt                   <  1:785758/100‑1 (MQ=255)
                                                                          ggTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACt                   >  2:3108175/1‑100 (MQ=255)
                                                                            tCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTc                 <  2:1636736/100‑1 (MQ=255)
                                                                             cACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCa                <  1:2931527/100‑1 (MQ=255)
                                                                             cACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCa                >  1:927308/1‑100 (MQ=255)
                                                                              aCGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAg               <  2:2849255/100‑1 (MQ=255)
                                                                               cGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGc              >  2:1030783/1‑100 (MQ=255)
                                                                               cGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGc              >  2:141587/1‑100 (MQ=255)
                                                                               cGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGc              >  2:757986/1‑100 (MQ=255)
                                                                                 ccaccaGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGCAc            >  2:387662/1‑100 (MQ=255)
                                                                                  caccaGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGCACg           >  1:1871410/1‑100 (MQ=255)
                                                                                  caccaGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGCACg           <  2:1096440/100‑1 (MQ=255)
                                                                                  caccaGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGCACg           <  2:2327099/100‑1 (MQ=255)
                                                                                     caGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGCACGAGt        >  2:1851631/1‑100 (MQ=255)
                                                                                     caGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGCACGAGt        >  1:619057/1‑100 (MQ=255)
                                                                                      aGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGCACGAGtc       >  2:341236/1‑100 (MQ=255)
                                                                                       gtactccAGTTCGGGCTCTTGATAGAAACAGACCGCGGACGGCCGGTGCCTTTCTGGCGCCACTGTTTTTTACCGGAACCAGTTACTTCAGCACGAGtct      <  2:2389034/95‑1 (MQ=255)
                                                                                       gAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGCACGAGtct      >  2:1495600/1‑100 (MQ=255)
                                                                                          cGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGCACGAGtcttct   >  2:951435/1‑100 (MQ=255)
                                                                                          cGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGCACGAGtcttct   <  2:530562/100‑1 (MQ=255)
                                                                                          cGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGATTTTTTACCGGACCCAGTTACTTCAGCACGaaactt     >  1:1079100/1‑93 (MQ=255)
                                                                                           gCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGCACGAGTCTTCTg  <  1:1544002/100‑1 (MQ=255)
                                                                                                |                                                                                              
AGGATGCTTTTCAGCGCGCCGCGGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGATCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGCACGAGTCTTCTG  >  NC_000913/3451961‑3452151

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