Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,452,057 A→T 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,0570AT95.6% 313.1 / 3.4 91I77N (ATC→AAC) rplD50S ribosomal subunit protein L4
Reads supporting (aligned to +/- strand):  ref base A (2/2);  new base T (38/49);  total (40/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

GGATGCTTTTCAGCGCGCCGCGGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGATCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGCACGAGTCTTCTGAG  >  NC_000913/3451962‑3452153
                                                                                               |                                                                                                
ggATGCTTTTCAGCGCGCCGCGGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCgg                                                                                              >  1:342845/1‑100 (MQ=255)
ggATGCTTTTCAGCGCGCCGCGGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCgg                                                                                              >  1:2081574/1‑100 (MQ=255)
  aTGCTTTTCAGCGCGCCGCGGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGc                                                                                            <  2:893899/100‑1 (MQ=255)
      ttttCAGCGCGCCGCGGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCtt                                                                                        <  1:1161148/100‑1 (MQ=255)
       tttCAGCGCGCCGCGGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTg                                                                                       <  1:108711/100‑1 (MQ=255)
       tttCAGCGCGCCGCGGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTg                                                                                       <  1:753687/100‑1 (MQ=255)
       tttCAGCGCGCCGCGGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTg                                                                                       <  2:118058/100‑1 (MQ=255)
        ttCAGCGCGCCGCGGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGa                                                                                      <  1:2016453/100‑1 (MQ=255)
         tCAGCGCGCCGCGGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGAt                                                                                     <  2:342843/100‑1 (MQ=255)
          cAGCGCGCCGCGGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATa                                                                                    >  2:1735864/1‑100 (MQ=255)
           aGCGCGCCGCGGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAg                                                                                   <  2:38904/100‑1 (MQ=255)
           aGCGCGCCGCGGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAg                                                                                   >  2:822098/1‑100 (MQ=255)
            gcgcgcCGCGGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGa                                                                                  >  2:2145780/1‑100 (MQ=255)
              gcgcCGCGGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAAc                                                                                <  1:2189461/100‑1 (MQ=255)
               cgccgcGGTACATCTTCTTGTTAACTTTTTTACTGTTGTCCTGCGGACTAGCAGCAAAGGTCACGCCACCAGAACTCCAGTTCGGGCTCTTGATAGAAcc                                                                               <  2:804773/100‑1 (MQ=255)
               cgccgcGGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAAcc                                                                               <  1:947746/100‑1 (MQ=255)
                 ccgcGGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGATCGGGCTCTTGATAGAACCAg                                                                             >  1:334993/1‑100 (MQ=255)
                   gcgGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGaa                                                                           <  2:108710/100‑1 (MQ=255)
                     ggTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAAcg                                                                         >  1:1176673/1‑100 (MQ=255)
                         cATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCa                                                                     <  2:401346/100‑1 (MQ=255)
                         cATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCa                                                                     >  1:1660865/1‑100 (MQ=255)
                           tcttctTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACg                                                                   >  2:1514855/1‑100 (MQ=255)
                           tcttctTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACg                                                                   >  2:348584/1‑100 (MQ=255)
                            cttcttGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACgg                                                                  >  1:1836805/1‑100 (MQ=255)
                            cttcttGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACgg                                                                  >  2:1270457/1‑100 (MQ=255)
                            cttcttGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACgg                                                                  >  2:2063244/1‑100 (MQ=255)
                            cttcttGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACgg                                                                  >  1:1791393/1‑100 (MQ=255)
                               cttGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCg                                                               >  2:1814367/1‑100 (MQ=255)
                                  gTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTg                                                            >  1:2173200/1‑100 (MQ=255)
                                  gTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTg                                                            >  1:352541/1‑100 (MQ=255)
                                  gTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTg                                                            <  2:422334/100‑1 (MQ=255)
                                    tAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGcc                                                          <  2:1614547/100‑1 (MQ=255)
                                    tAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGcc                                                          >  2:1606277/1‑100 (MQ=255)
                                     aaCTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCt                                                         >  1:1113246/1‑100 (MQ=255)
                                       cTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCttt                                                       >  1:644100/1‑100 (MQ=255)
                                       cTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCACTTGATAGAACCAGAACGCGCACGGCCGGTGCCttt                                                       >  1:1684048/1‑100 (MQ=255)
                                         ttttGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCt                                                     <  1:717183/100‑1 (MQ=255)
                                          tttGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCtn                                                    <  2:860279/100‑2 (MQ=255)
                                            tGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGgc                                                  <  1:1814372/100‑1 (MQ=255)
                                             gACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGgcg                                                 <  1:1923664/100‑1 (MQ=255)
                                              accgtggtgCTCCGTCCGAGCAGCAAAGGTCACGCCACTAGAACGCCGGTTCGGTCTCTTGATAGAACCAGAGCGCGCACGGCCGGTGcc                                                          <  2:726035/87‑1 (MQ=255)
                                               cTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGcc                                               >  1:2080352/1‑100 (MQ=255)
                                               cTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGcc                                               <  1:1514859/100‑1 (MQ=255)
                                               cTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGcc                                               >  2:1100915/1‑100 (MQ=255)
                                               cTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGAGCCTTTCTGGCGcc                                               <  2:1287851/100‑1 (MQ=255)
                                                tgtgGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCa                                              >  1:1463715/1‑100 (MQ=255)
                                                   ggTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACgg                                           >  2:1473607/1‑100 (MQ=255)
                                                     tCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGtt                                         <  1:2145787/100‑1 (MQ=255)
                                                      ccTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGttt                                        <  2:352539/100‑1 (MQ=255)
                                                        tGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGttttt                                      <  1:181728/100‑1 (MQ=255)
                                                          cGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTa                                    >  1:1530527/1‑100 (MQ=255)
                                                            gACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTAcc                                  <  1:823925/100‑1 (MQ=255)
                                                            gACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTAcc                                  <  1:348586/100‑1 (MQ=255)
                                                             aCGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCg                                 <  2:820289/100‑1 (MQ=255)
                                                              cGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCgg                                <  1:2185429/100‑1 (MQ=255)
                                                               gAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGa                               <  2:2173193/100‑1 (MQ=255)
                                                                agcagcAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGaa                              <  2:1546010/100‑1 (MQ=255)
                                                                agcagcAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGaa                              <  2:1545850/100‑1 (MQ=255)
                                                                agcagcAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGaa                              >  2:1345427/1‑100 (MQ=255)
                                                                  cagcaAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAAcc                            <  1:1270460/100‑1 (MQ=255)
                                                                  cagcaAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAAcc                            >  2:389244/1‑100 (MQ=255)
                                                                  cagcaAAGGTCACGCCACCAGAACGCCAGATCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAAcc                            >  1:1793267/1‑100 (MQ=255)
                                                                   agcaAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCa                           >  2:1187703/1‑100 (MQ=255)
                                                                     caAAGGTCACGCCACCAGGACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGt                         <  2:235345/100‑1 (MQ=255)
                                                                        aGGTCACGCCACCAGAACGCCAGATCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTAc                      <  2:334991/100‑1 (MQ=255)
                                                                         ggTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACt                     <  1:500970/100‑1 (MQ=255)
                                                                         ggTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACt                     >  1:820112/1‑100 (MQ=255)
                                                                         ggTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACt                     >  2:769875/1‑100 (MQ=255)
                                                                           tctcGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTc                   <  1:12621/97‑1 (MQ=255)
                                                                           tCACGCCCCCAGCACGCCAGTTCGGGCTCTTGATAGAACCAGCACGCGCACGGCCGGTGCCTTTCTGGCGCCGCGGGTTTTTACCGGAACCAGTTACTTc                   <  2:1756278/100‑1 (MQ=255)
                                                                           tCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTc                   <  2:1756277/100‑1 (MQ=255)
                                                                           tCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTc                   <  1:994439/100‑1 (MQ=255)
                                                                            cACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCa                  >  2:967091/1‑100 (MQ=255)
                                                                            cACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCa                  >  2:965788/1‑100 (MQ=255)
                                                                             aCGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAg                 <  2:109391/100‑1 (MQ=255)
                                                                              cGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAga                >  2:12621/1‑99 (MQ=255)
                                                                                ccaccaGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGCAc              <  1:161970/100‑1 (MQ=255)
                                                                                 caccaGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGCACg             <  2:1463711/100‑1 (MQ=255)
                                                                                 caccaGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGCACg             >  2:541463/1‑100 (MQ=255)
                                                                                 caccaGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGCACg             >  1:217739/1‑100 (MQ=255)
                                                                                    caGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGCACGAGt          >  2:1550677/1‑100 (MQ=255)
                                                                                     aGAACGCCAGATCGGTCTCTTGATTGAACCAGAACGCGCTCGGCCGGTGCCGTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGCACGAg           <  2:1613716/98‑1 (MQ=255)
                                                                                        aCGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGCACGAGtcttc      <  1:1735869/100‑1 (MQ=255)
                                                                                        aCGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGCACGAGtcttc      <  1:965790/100‑1 (MQ=255)
                                                                                        aCGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGCACGAGtcttc      <  1:967093/100‑1 (MQ=255)
                                                                                        aCGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGCACGAGtcttc      <  2:1791388/100‑1 (MQ=255)
                                                                                          gCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGCACGAGTCTTCTg    <  2:569465/100‑1 (MQ=255)
                                                                                           ccAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGCACGAGTCTTCTGa   <  2:1684043/100‑1 (MQ=255)
                                                                                           ccAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGCACGAGTCTTCTGa   >  2:2287425/1‑100 (MQ=255)
                                                                                            cAGTTCGGGCTCTTGATAGAACCAGAACTCGCACGGCCGGTGCCTTTCTGGCGCCACGGTGTTTTACCGGAACCAGTTACTTCAGCAGGAGTCGTCGGAg  >  1:1277421/1‑100 (MQ=255)
                                                                                            cAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGCACGAGTCTTCTGAg  <  2:644099/100‑1 (MQ=255)
                                                                                               |                                                                                                
GGATGCTTTTCAGCGCGCCGCGGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGATCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGCACGAGTCTTCTGAG  >  NC_000913/3451962‑3452153

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