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,0570AT99.2% 428.6 / ‑4.2 128I77N (ATC→AAC) rplD50S ribosomal subunit protein L4
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base T (61/66);  total (61/67)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.08e-01

Only 100 of 128 total aligned reads displayed.

GGATGCTTTTCAGCGCGCCGCGGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGATCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGCACGAGTCTTCTGA  >  NC_000913/3451962‑3452152
                                                                                               |                                                                                               
  aTGCTTTTCAGCGCGCCGCGGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGc                                                                                           <  1:1747138/100‑1 (MQ=255)
  aTGCTTTTCAGCGCGCCGCGGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGc                                                                                           <  1:1388985/100‑1 (MQ=255)
   tGCTTTTCAGCGCGCCGCGGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGct                                                                                          <  1:602274/100‑1 (MQ=255)
   tGCTTTTCAGCGCGCCGCGGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGATCGGGct                                                                                          <  1:2293540/100‑1 (MQ=255)
     cTTTTCAGCGCGCCGCGGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGctct                                                                                        >  2:403727/1‑100 (MQ=255)
     cTTTTCAGCGCGCCGCGGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGctct                                                                                        <  2:1812087/100‑1 (MQ=255)
      ttttCAGCGCGCCGCGGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCtt                                                                                       <  1:1475064/100‑1 (MQ=255)
      ttttCAGCGCGCCGCGGTACATCTTCTTGTTAACTTTTTGACTGTGATCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCtt                                                                                       >  2:1327862/1‑100 (MQ=255)
          cAGCGCGCCGCGGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATa                                                                                   <  2:845686/100‑1 (MQ=255)
           aGCGCGCCGCGGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAg                                                                                  >  1:2008263/1‑100 (MQ=255)
             cgcgcCGCGGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGaa                                                                                >  1:69274/1‑100 (MQ=255)
              gcgcCGCGGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAAC                                                                               <  2:1264796‑M1/100‑1 (MQ=255)
                 ccgcGGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAG                                                                            >  1:47791‑M1/1‑100 (MQ=255)
                 ccgcGGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCA                                                                             <  2:256290‑M1/99‑1 (MQ=255)
                   gcgGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAA                                                                          >  2:1121602‑M1/1‑100 (MQ=255)
                     ggTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCATCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAAcg                                                                        <  2:1278371/100‑1 (MQ=255)
                     ggTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAAcg                                                                        >  2:2341580/1‑100 (MQ=255)
                       tACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAAcgcg                                                                      <  1:930785/100‑1 (MQ=255)
                        aCATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACgcgc                                                                     <  2:2382334/100‑1 (MQ=255)
                          aTCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCa                                                                    >  2:522292/1‑99 (MQ=255)
                            cttcttGTTAACTTTTTGACTATGATACTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCa                                                                    >  2:1028198/1‑97 (MQ=255)
                             ttcttGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGc                                                                >  1:1664651/1‑100 (MQ=255)
                              tcttGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGcc                                                               >  2:865407/1‑100 (MQ=255)
                               cttGTTAACTTTTTGACTGTAATCCTGCGGACGAGCAGCAAAGGTAACGCCACCAGAAAGCCAGTTCGGGCTCTTGATAGAACCAGAAAGCGCACGGcc                                                               >  2:2117417/1‑99 (MQ=255)
                                 tGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCATGGCCGGt                                                            <  1:1258608/100‑1 (MQ=255)
                                 tGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGt                                                            >  1:1139155/1‑100 (MQ=255)
                                 tGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGt                                                            >  1:1283369/1‑100 (MQ=255)
                                  gTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTg                                                           <  1:55173/100‑1 (MQ=255)
                                  gTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTg                                                           >  1:2234791/1‑100 (MQ=255)
                                    tctcTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGcc                                                         <  2:771835/97‑1 (MQ=255)
                                     aaCTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGc                                                          >  1:111956/1‑98 (MQ=255)
                                      aCTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTt                                                       <  1:815804‑M2/100‑2 (MQ=255)
                                      aCTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTt                                                       <  1:231019‑M2/100‑2 (MQ=255)
                                       cTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGccaga                                                      >  1:771844/1‑97 (MQ=255)
                                       cTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTtt                                                      <  1:1453511‑M2/100‑3 (MQ=255)
                                       cTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTtt                                                      <  2:118545‑M2/100‑3 (MQ=255)
                                           ttGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTttctgg                                                  >  2:1704702‑M2/1‑94 (MQ=255)
                                            tGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGgc                                                 >  1:1271733/1‑100 (MQ=255)
                                             gACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGCAaggcaccggccgtgcgcgttctgg                                                  <  2:350560‑M1/98‑25 (MQ=255)
                                             gACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGgc                                                 >  1:147577/1‑99 (MQ=255)
                                              aCTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCTGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGgcg                                                <  2:332341/99‑1 (MQ=255)
                                               cTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTTTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGcc                                              <  2:378326/100‑1 (MQ=255)
                                               cTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAAaggcaccggccgtgcgcgttctggttct                                              <  2:2191742‑M1/100‑29 (MQ=255)
                                               cTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGcc                                              <  2:1776934/100‑1 (MQ=255)
                                               cTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGcc                                              <  2:1867941/100‑1 (MQ=255)
                                               cTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGcc                                              >  2:2386312/1‑100 (MQ=255)
                                                 gtgGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAAaggcaccggccgtgcgcgttctggttctat                                            <  1:1796842‑M1/100‑31 (MQ=255)
                                                 gtgGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCAc                                            >  1:1112230/1‑100 (MQ=255)
                                                    gTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAAaggcaccggccgtgcgcgttctggttctatcaa                                         <  1:1543640‑M1/100‑34 (MQ=255)
                                                     tCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAAaggcaccggccgtgcgcgttctggttctatcaag                                        <  1:1013106‑M1/100‑35 (MQ=255)
                                                     aagaagcgccGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGtt                                        <  2:1992249/91‑1 (MQ=255)
                                                       cTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGtttt                                      >  2:1963254/1‑100 (MQ=255)
                                                        tGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGttttt                                     >  2:2303010/1‑100 (MQ=255)
                                                        tGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGtttt                                      <  1:656478/99‑1 (MQ=255)
                                                        tGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGtttt                                      >  1:2117658/1‑99 (MQ=255)
                                                         gCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGtttttt                                    <  1:801231/100‑1 (MQ=255)
                                                          cGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTa                                   <  2:254348/100‑1 (MQ=255)
                                                          cGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTa                                   >  1:1362240/1‑100 (MQ=255)
                                                           ggtcGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTAc                                  <  2:1453495/97‑1 (MQ=255)
                                                            gACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTAcc                                 >  1:204414/1‑100 (MQ=255)
                                                             aCGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCg                                <  1:2222720/100‑1 (MQ=255)
                                                              ctagcagcAAAGGTCACGCCACCAGAACGCGAGTTGCGGCTCTTTATAGAACCAGAACTCGCACGGCCGTTGCCTTTCTGGCGCCACGGTTTTTTACCgg                               <  2:1612734/98‑1 (MQ=255)
                                                              cGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCgg                               >  1:1816197/1‑100 (MQ=255)
                                                              cGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCgg                               >  1:2120552/1‑100 (MQ=255)
                                                              cGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCgg                               <  1:788000/100‑1 (MQ=255)
                                                               gAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGa                              <  2:813186/100‑1 (MQ=255)
                                                                agcagcAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGACACGGTATTTTAcc                                 >  2:1468495/1‑96 (MQ=255)
                                                                 gcagcaAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAAc                            >  2:718880/1‑100 (MQ=255)
                                                                 gcagcaAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAAc                            >  1:400885/1‑100 (MQ=255)
                                                                   agcaAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAAcc                           >  2:1583518/1‑99 (MQ=255)
                                                                   agcaAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCa                          >  2:1089538/1‑100 (MQ=255)
                                                                   agcaAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCa                          >  1:1191926/1‑100 (MQ=255)
                                                                   agcaAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCAGGTGCATTTCTGGCTCCACGGTTTTTTACCGGaa                             >  2:1613154/1‑97 (MQ=255)
                                                                   agcaAAGCTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCa                          >  1:1257686/1‑100 (MQ=255)
                                                                   aacaAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGGGCACGGCCGGTGCCTTTCTGGCTCCACGGTTTTTTACCGGAACCa                          <  2:1923822/98‑1 (MQ=255)
                                                                    gcaAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAg                         >  1:580187/1‑100 (MQ=255)
                                                                    gcaAAGGTCACGCCACCAGAACGCCAGTTCGAGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTTGCGCCATGGTGTTTTACCGGAACCAg                         >  2:228946/1‑100 (MQ=255)
                                                                      aaaGGTCACGCCACCAGGACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGtt                       <  2:400882/100‑1 (MQ=255)
                                                                      aaaGGTCACACCACCACAACGCCAGTTCGCGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCTCCACGGTTTTTTACCGGAACCAGtt                       <  2:1442490/100‑1 (MQ=255)
                                                                       aaGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTa                      <  2:192165/100‑1 (MQ=255)
                                                                        aGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAATCAGTTAc                     >  2:205057/1‑100 (MQ=255)
                                                                         ggTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACt                    <  1:1781041/100‑1 (MQ=255)
                                                                         ggTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACt                    >  1:552550/1‑100 (MQ=255)
                                                                           tCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTc                  <  2:2039779/100‑1 (MQ=255)
                                                                             aCGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAg                <  1:1193838/100‑1 (MQ=255)
                                                                              cGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGc               >  1:30250/1‑100 (MQ=255)
                                                                              cGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTCTTTACCGGAACCAGTTACTTCAGc               >  1:1791789/1‑100 (MQ=255)
                                                                               gCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGCa              <  1:454769/100‑1 (MQ=255)
                                                                                 caccaGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGCACg            >  2:1079134/1‑100 (MQ=255)
                                                                                 caccaGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGCACg            <  1:648278/100‑1 (MQ=255)
                                                                                  accaGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGCACGa           >  2:1221811/1‑100 (MQ=255)
                                                                                   ccaGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGCACGAg          <  2:204411/100‑1 (MQ=255)
                                                                                   ccaGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGCACGAg          <  1:273221/100‑1 (MQ=255)
                                                                                       aaCGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAACACTAGtcta      >  1:2169530/1‑99 (MQ=255)
                                                                                        aCGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGCACGAGtcttc     <  1:75564/100‑1 (MQ=255)
                                                                                         cGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGCACGAGtcttct    <  2:1386896/100‑1 (MQ=255)
                                                                                         cGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGCACGAGtcttct    >  1:164273/1‑100 (MQ=255)
                                                                                           ccAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGCACGAGTCTTCTGa  >  1:1888587/1‑100 (MQ=255)
                                                                                            cAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGCACGAGTCTTCTg   <  2:1455062/98‑1 (MQ=255)
                                                                                            cAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGCACGAGTCTTCTGa  >  2:490879/1‑99 (MQ=255)
                                                                                               |                                                                                               
GGATGCTTTTCAGCGCGCCGCGGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGATCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGCACGAGTCTTCTGA  >  NC_000913/3451962‑3452152

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: