Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,175,949 +G intergenic (+82/‑76) nusG → / → rplK NusG/50S ribosomal subunit protein L11

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,175,9471.G98.7% 839.1 / 1.3 231intergenic (+80/‑78)nusG/rplKNusG/50S ribosomal subunit protein L11
Reads supporting (aligned to +/- strand):  ref base . (1/2);  new base G (102/126);  total (103/128)
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.98e-01

Only 100 of 269 total aligned reads displayed.

TGAAAAAGCCTAACCCAGCGATCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACTTGAGGAATTTATAATGGCTAAG  >  NC_000913/4175855‑4176033
                                                                                              |                                                                                      
tGAAAAAGCCTAACCCAGCGATCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTAt                                                                                         <  1:433909/93‑1 (MQ=255)
tGAAAAAGCCTAACCCAGCGATCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTAt                                                                                         <  1:5977735/93‑1 (MQ=255)
tGAAAAAGCCTAACCCAGCGATCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTAt                                                                                         <  1:1830025/93‑1 (MQ=255)
 gAAAAAGCCTAACCCAGCGATCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATg                                                                                       >  1:1500985/1‑93 (MQ=255)
 gAAAAAGCCTAACCCAGCGATCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATg                                                                                       <  1:3253081/93‑1 (MQ=255)
 gAAAAAGCCTAACCCAGCGATCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATg                                                                                       <  1:2392182/93‑1 (MQ=255)
  aaaaaGCCTAACCCAGCGATCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATgg                                                                                      <  2:4156834/93‑1 (MQ=255)
   aaaaGCCTAACCCAGCGATCAAAAGAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATggg                                                                                     >  1:6406168/1‑92 (MQ=255)
   aaaaGCCTAACCCAGCGATCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATggg                                                                                     >  1:1901459/1‑92 (MQ=255)
    aaaGCCTAACCCAGCGATCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATgggc                                                                                    <  2:4060538/93‑3 (MQ=255)
    aaaGCCTAACCCAGCGATCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATgggc                                                                                    >  1:6335774/1‑91 (MQ=255)
     aaGCCTAACCCAGCGATCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATgggcc                                                                                   <  2:6054940/93‑4 (MQ=255)
      aGCCTAACCCAGCGATCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATgggcct                                                                                  <  2:2746821/93‑5 (MQ=255)
      aGCCTAACCCAGCGATCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATgggcct                                                                                  <  1:4765128/93‑5 (MQ=255)
      aGCCTAACCCAGCGATCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATgggcct                                                                                  <  2:321273/93‑5 (MQ=255)
        ccTAACCCAGCGATCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTg                                                                                 <  2:5515811/93‑1 (MQ=255)
         cTAACCCAGCGATCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACGATTTCGCGCCTTTTGTTTTTATGGGCCTTGc                                                                                >  2:1116671/1‑93 (MQ=255)
         cTAACCCAGCGATCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGc                                                                                >  2:813996/1‑93 (MQ=255)
          tAACCCAGCGATCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGcc                                                                               >  2:2339830/1‑93 (MQ=255)
          tAACCCAGCGATCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGcc                                                                               <  1:2344613/93‑1 (MQ=255)
          tAACCCAGCGATCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGcc                                                                               <  2:3305947/93‑1 (MQ=255)
          tAACCCAGCGATCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGcc                                                                               >  1:6069815/1‑93 (MQ=255)
          tAACCCAGCGATCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGcc                                                                               >  1:4214155/1‑93 (MQ=255)
           aaCCCAGCGATCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGccc                                                                              <  2:5180667/93‑1 (MQ=255)
           aaCCCAGCGATCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGccc                                                                              >  1:3744192/1‑93 (MQ=255)
            aCCCAGCGATCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCg                                                                             >  2:1248841/1‑93 (MQ=255)
              ccAGCGATCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTa                                                                           >  1:5302975/1‑93 (MQ=255)
               cAGCGATCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTaa                                                                          >  1:478392/1‑93 (MQ=255)
                aGCGATCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTaaa                                                                         <  1:1106128/93‑1 (MQ=255)
                 gCGATCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTaaaa                                                                        >  2:235844/1‑93 (MQ=255)
                 gCGATCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTaaaa                                                                        >  2:3644576/1‑93 (MQ=255)
                 gCGATCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTaaaa                                                                        >  1:5388070/1‑93 (MQ=255)
                 gCGATCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTaaaa                                                                        >  1:2007414/1‑93 (MQ=255)
                   gATCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACg                                                                      >  2:584831/1‑93 (MQ=255)
                    aTCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGa                                                                     >  2:1983913/1‑93 (MQ=255)
                      cAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTACGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGAtt                                                                   >  2:2914046/1‑93 (MQ=255)
                       aaaaaaGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGAttt                                                                  <  2:2243035/93‑1 (MQ=255)
                       aaaaaaGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGAttt                                                                  <  2:6626134/93‑1 (MQ=255)
                       aaaaaaGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGAttt                                                                  <  1:2682219/93‑1 (MQ=255)
                          aaaGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGAtttttt                                                               >  1:46531/1‑93 (MQ=255)
                           aaGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTta                                                              >  1:1068044/1‑93 (MQ=255)
                             gcggcgATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTtata                                                            <  2:5129479/93‑1 (MQ=255)
                             gcggcgATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTtata                                                            <  2:5371561/93‑1 (MQ=255)
                             gcggcgATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTtata                                                            <  2:5739301/93‑1 (MQ=255)
                             gcggcgATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTtaca                                                            >  1:3915947/1‑91 (MQ=255)
                              cggcgATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTatat                                                           <  2:5687671/93‑1 (MQ=255)
                               ggcgATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATc                                                          >  2:4028881/1‑93 (MQ=255)
                               ggcgATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATc                                                          >  1:2293649/1‑93 (MQ=255)
                                gcgATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCa                                                         <  2:887831/93‑1 (MQ=255)
                                 cgATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCAc                                                        <  2:6586552/93‑1 (MQ=255)
                                  gATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACg                                                       >  2:1527482/1‑93 (MQ=255)
                                  gATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACg                                                       >  1:537143/1‑93 (MQ=255)
                                   aTTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACgg                                                      <  2:523488/93‑1 (MQ=255)
                                    tttAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACggg                                                     >  2:4531971/1‑93 (MQ=255)
                                     ttAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACgggg                                                    <  1:2618035/93‑1 (MQ=255)
                                      tAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGa                                                   <  2:1420950/93‑1 (MQ=255)
                                       aaTCGTTGCACAAGGCGCGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAg                                                  >  1:2203445/1‑93 (MQ=255)
                                        aTCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGc                                                 >  1:5853711/1‑93 (MQ=255)
                                         tCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGcc                                                <  2:345262/93‑1 (MQ=255)
                                         tCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGcc                                                <  2:1783894/93‑1 (MQ=255)
                                         tCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGcc                                                >  2:2529908/1‑93 (MQ=255)
                                          cGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCct                                               <  1:5028977/93‑1 (MQ=255)
                                          cGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCct                                               <  1:304716/93‑1 (MQ=255)
                                           gTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCctc                                              >  1:2974421/1‑93 (MQ=255)
                                            gtGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCctct                                             <  1:1170250/92‑1 (MQ=255)
                                             tGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCtctc                                            <  2:4814257/93‑1 (MQ=255)
                                             tGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCtctc                                            >  2:5988176/1‑93 (MQ=255)
                                             tGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCtctc                                            <  1:3563685/93‑1 (MQ=255)
                                              gCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCa                                           <  1:5678057/93‑1 (MQ=255)
                                                acaAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCaga                                         >  2:6297040/1‑93 (MQ=255)
                                                acaAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCaga                                         >  2:4579991/1‑93 (MQ=255)
                                                acaAGGCCTGAGATTGGAATTCAATTTCGCGCTTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATCTCACGGGGAGCCTCTCaga                                         <  1:6416581/93‑1 (MQ=255)
                                                 caAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAGCGATTTTTTATATCACGGGGAGCCTCTCagag                                        >  1:5077756/1‑93 (MQ=255)
                                                  aaGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTACGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAgg                                       >  2:3546531/1‑93 (MQ=255)
                                                  aaGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAgg                                       <  2:1092569/93‑1 (MQ=255)
                                                  aaGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAgg                                       <  2:1278124/93‑1 (MQ=255)
                                                   aGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGc                                      >  2:5279617/1‑93 (MQ=255)
                                                      cGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGtt                                   <  2:563203/93‑1 (MQ=255)
                                                       gTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGtta                                  >  1:5284359/1‑93 (MQ=255)
                                                       gTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGtta                                  <  1:4714820/93‑1 (MQ=255)
                                                       gTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGtta                                  >  1:1687794/1‑93 (MQ=255)
                                                              tGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCa                           <  2:3757958/93‑1 (MQ=255)
                                                               ggAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGGGGCGTTATTACCCaa                          <  1:3601318/93‑1 (MQ=255)
                                                                gAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCAAc                         >  2:1881309/1‑93 (MQ=255)
                                                                gAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCAAc                         <  2:5404703/93‑1 (MQ=255)
                                                                  aTACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACtt                       >  2:6184858/1‑93 (MQ=255)
                                                                   tACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGGGGCGTTATTACCCAACTTg                      >  2:3241776/1‑93 (MQ=255)
                                                                       aTTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACTTGAGGa                  <  1:5141175/93‑1 (MQ=255)
                                                                        tttCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACTTgaagaa                 <  1:461061/93‑1 (MQ=255)
                                                                        tttCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACTTGAGGaa                 >  1:79152/1‑93 (MQ=255)
                                                                         ttCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGGGGCGTTATTACCCAACTTGAGGAAt                <  1:3440053/93‑1 (MQ=255)
                                                                         ttCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACTTGAGGAAt                <  2:6436409/93‑1 (MQ=255)
                                                                           cgcgCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGGGGCGTTATTACCCAACTTGATGAAttt              >  2:500251/1‑93 (MQ=255)
                                                                           cgcgCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACTTGAGGAAttt              <  2:4233655/93‑1 (MQ=255)
                                                                            gcgcCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACTTGAGGAATTta             >  1:626125/1‑93 (MQ=255)
                                                                              gcCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACTTGAGGAATTtata           >  2:4958341/1‑93 (MQ=255)
                                                                               ccTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACTTGAGGAATTTATaa          <  2:729258/93‑1 (MQ=255)
                                                                                cTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACTTGAGGAATTTATAAt         >  1:2824951/1‑93 (MQ=255)
                                                                                 ttttGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACTTGAGGAATTTATAATg        >  1:154169/1‑93 (MQ=255)
                                                                                   ttGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACTTGAGGAATTTATAATGGc      <  1:1576956/93‑1 (MQ=255)
                                                                                              |                                                                                      
TGAAAAAGCCTAACCCAGCGATCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACTTGAGGAATTTATAATGGCTAAG  >  NC_000913/4175855‑4176033

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