Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,185,306 A→T 5.6% intergenic (+33/‑44) rpoB → / → rpoC RNA polymerase subunit beta/RNA polymerase subunit beta'

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,185,3060AT 5.6% 561.8 / 16.4 178intergenic (+33/‑44)rpoB/rpoCRNA polymerase subunit beta/RNA polymerase subunit beta'
Reads supporting (aligned to +/- strand):  ref base A (83/85);  new base T (2/8);  total (85/93)
Fisher's exact test for biased strand distribution p-value = 1.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.05e-01

Only 100 of 189 total aligned reads displayed.

TATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTGCTCTGGCTTCGCCAGACATGATCCGTTCATGGTCTTTC  >  NC_000913/4185157‑4185454
                                                                                                                                                     |                                                                                                                                                    
tataAAAACATCGTGGACGTCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTa                                                                                                                                                      <  1:665520/150‑1 (MQ=255)
tataAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTa                                                                                                                                                      >  1:604309/1‑150 (MQ=255)
  taAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTaaa                                                                                                                                                    >  1:345628/1‑150 (MQ=255)
  aaaaaaCATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTaaa                                                                                                                                                    >  1:393720/2‑150 (MQ=255)
   aaaaaCATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTaaaa                                                                                                                                                   >  2:1247950/1‑150 (MQ=255)
        cATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCgg                                                                                                                                              >  2:1379385/1‑150 (MQ=255)
         aTCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGc                                                                                                                                             <  1:1366519/150‑1 (MQ=255)
           cGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGCCGGGTTTTAACCCGACAg                                                                                                                                           <  2:431335/150‑1 (MQ=255)
           cGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGCCGGGTTTTAACCCGACAg                                                                                                                                           >  1:1191564/1‑150 (MQ=255)
            gTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTTTTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGCCGGGTTTTAACCCGACAGc                                                                                                                                          <  1:1324944/150‑1 (MQ=255)
            gTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGc                                                                                                                                          >  2:138469/1‑150 (MQ=255)
              ggACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGCCGGGTTTTAACCCGACAGcag                                                                                                                                        <  2:1089592/150‑3 (MQ=255)
               gACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGa                                                                                                                                       >  1:1249742/1‑150 (MQ=255)
               gACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGa                                                                                                                                       <  1:1379669/150‑1 (MQ=255)
                aCGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGCCGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGAt                                                                                                                                      >  2:1313758/1‑150 (MQ=255)
                aCGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGCCGGGTTTTAACCCGACAGcagtg                                                                                                                                      <  2:1191564/150‑5 (MQ=255)
                aCGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGCCGGGTTTTAACCCGACAGcagtg                                                                                                                                      <  2:1163745/150‑5 (MQ=255)
                   gCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATtgt                                                                                                                                   >  2:26070/1‑150 (MQ=255)
                    cAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCTGCAGCGGATtttg                                                                                                                                  <  2:345628/150‑1 (MQ=255)
                           cAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTTGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTc                                                                                                                           >  2:302910/1‑150 (MQ=255)
                            aGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTcc                                                                                                                          >  2:1208677/1‑150 (MQ=255)
                            aGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTcc                                                                                                                          >  2:742421/1‑150 (MQ=255)
                             gATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCg                                                                                                                         >  1:1380151/1‑150 (MQ=255)
                                ggAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACg                                                                                                                      <  2:351280/150‑1 (MQ=255)
                                ggAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACg                                                                                                                      <  1:1237274/150‑1 (MQ=255)
                                 gAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACgg                                                                                                                     >  1:1059153/1‑150 (MQ=255)
                                    ccGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAg                                                                                                                  <  1:904627/150‑1 (MQ=255)
                                     cGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCg                                                                                                                                               >  1:1392089/1‑120 (MQ=255)
                                     cGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGc                                                                                                                 >  2:1104724/1‑150 (MQ=255)
                                     cGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGc                                                                                                                 <  1:305866/150‑1 (MQ=255)
                                     cGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGc                                                                                                                 <  2:769725/150‑1 (MQ=255)
                                      gggCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCa                                                                                                                <  1:1161815/150‑1 (MQ=255)
                                         cATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAAt                                                                                                             <  2:1303559/150‑1 (MQ=255)
                                         cATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACACGTCACTGCTGTCGGCTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAAt                                                                                                             >  2:899784/1‑150 (MQ=255)
                                           tGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATcc                                                                                                           >  1:407470/1‑150 (MQ=255)
                                            gCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCg                                                                                                          >  2:688527/1‑150 (MQ=255)
                                            gCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCg                                                                                                          >  1:177408/1‑150 (MQ=255)
                                            gCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCg                                                                                                          >  1:578313/1‑150 (MQ=255)
                                             ccAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGt                                                                                                         <  1:276369/150‑1 (MQ=255)
                                              cAGAATCCTTCAACGTTTTGTTGAAAGAGATTCGTTCGCTGTTTATCAACATGGATCTGGAAGACGAGTAATTCTCGCTCAAACATGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCGGACGGGAGCAAATCCGTg                                                                                                        <  2:387270/150‑1 (MQ=255)
                                                gAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGaa                                                                                                      <  1:9570/150‑1 (MQ=255)
                                                 aaTCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGaaa                                                                                                     <  2:853350/150‑1 (MQ=255)
                                                  aTCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAg                                                                                                    <  2:26718/150‑1 (MQ=255)
                                                    ccTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGAt                                                                                                  >  2:296692/1‑150 (MQ=255)
                                                    ccTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACAGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGAt                                                                                                  >  1:1351231/1‑150 (MQ=255)
                                                     cTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGt                                                                                                         >  1:1093885/1‑142 (MQ=255)
                                                     cTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGt                                                                                                         <  2:1093885/142‑1 (MQ=255)
                                                       tCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATtta                                                                                               >  1:1037929/1‑150 (MQ=255)
                                                         aaCGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAATAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGt                                                                                                         <  2:1000836/138‑1 (MQ=255)
                                                         aaCGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAATAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGt                                                                                                         >  1:1000836/1‑138 (MQ=255)
                                                         aaCGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATttatt                                                                                             >  1:509724/1‑150 (MQ=255)
                                                          aCGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAg                                                                                                                                           <  1:442404/103‑1 (MQ=255)
                                                          aCGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAg                                                                                                                                           >  2:442404/1‑103 (MQ=255)
                                                          aCGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATttatta                                                                                            <  2:1362715/150‑1 (MQ=255)
                                                            gTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTaaa                                                                                          <  1:705784/150‑1 (MQ=255)
                                                             tATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAg                                                                                         >  1:216411/1‑150 (MQ=255)
                                                                tgttgAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGttt                                                                                      >  1:391826/1‑150 (MQ=255)
                                                                 gttgAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTc                                                                                     >  2:610680/1‑150 (MQ=255)
                                                                 gttgAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTc                                                                                     <  2:216411/150‑1 (MQ=255)
                                                                 gttgAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTc                                                                                     >  1:634421/1‑150 (MQ=255)
                                                                  ttgAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCt                                                                                    <  2:514788/150‑1 (MQ=255)
                                                                   tgAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTg                                                                                   >  1:128851/1‑150 (MQ=255)
                                                                   tgAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTg                                                                                   <  2:700089/150‑1 (MQ=255)
                                                                    gAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGa                                                                                  >  1:227542/1‑150 (MQ=255)
                                                                            ttcgttcgCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAg                                                                          <  2:124434/150‑1 (MQ=255)
                                                                              cgttcgCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGAc                                                                        <  2:1423327/150‑1 (MQ=255)
                                                                               gttcgCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACt                                                                       >  1:877647/1‑150 (MQ=255)
                                                                                ttcgCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTa                                                                      <  2:190644/150‑1 (MQ=255)
                                                                                ttcgCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTa                                                                      >  2:1261924/1‑150 (MQ=255)
                                                                                     tGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAAcc                                                                 <  1:1208677/150‑1 (MQ=255)
                                                                                            cAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAGACCCGGCAGTGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGt                                                          >  1:747977/1‑150 (MQ=255)
                                                                                            cAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGt                                                          >  2:1354570/1‑150 (MQ=255)
                                                                                             aaCATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGtt                                                         >  2:715524/1‑150 (MQ=255)
                                                                                             aaCATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGtt                                                         >  1:1019321/1‑150 (MQ=255)
                                                                                                 tCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGAt                                                     <  2:1249742/150‑1 (MQ=255)
                                                                                                   gAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGc                                                   <  2:1059153/150‑1 (MQ=255)
                                                                                                    aaCTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCg                                                  <  2:1142988/150‑1 (MQ=255)
                                                                                                      cTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGAt                                                <  2:1037929/150‑1 (MQ=255)
                                                                                                      cTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGAt                                                >  1:801711/1‑150 (MQ=255)
                                                                                                        ggAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCa                                              >  2:1029173/1‑150 (MQ=255)
                                                                                                        cgAATACGAGTAATTCTCGCTCAAACAGGTCACTCCTCTCGGGTTAAAAACAGGCAGCGGATTGTGCTAACTCCGACGGGAGAAAATACGTCAAAGAATTATAAAAGTTTCTGAAAGCGCAGACTAAAACAGAAGAGTTTGATGAGAACa                                              >  2:1293820/2‑150 (MQ=255)
                                                                                                         gAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCaa                                             <  1:1171540/150‑1 (MQ=255)
                                                                                                          aaGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCaaa                                            <  2:874836/150‑1 (MQ=255)
                                                                                                              cGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTg                                        >  1:775710/1‑150 (MQ=255)
                                                                                                                       ctcGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTGCTCTGGCtt                               >  1:1301235/1‑150 (MQ=255)
                                                                                                                         cGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTGCTCTGGCTTCg                             <  2:645660/150‑1 (MQ=255)
                                                                                                                         cGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTGCTCTGGCTTCg                             >  2:952667/1‑150 (MQ=255)
                                                                                                                          gCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTGCTCTGGCTTCGc                            <  2:924507/150‑1 (MQ=255)
                                                                                                                           cTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTGCTCTGGCTTCGcc                           >  2:1096783/1‑150 (MQ=255)
                                                                                                                            tCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTGCTCTGGCTTCGCCa                          <  1:899210/150‑1 (MQ=255)
                                                                                                                              aaaCAGGTAACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTGCTCTGGCTTCGCCAGa                        >  1:844616/1‑150 (MQ=255)
                                                                                                                                 cAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGGTTGATGCGATCAAAATTGCTCTGGCTTCGCCAGACAt                     <  1:969791/150‑1 (MQ=255)
                                                                                                                                        ctgctgTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTGCTCTGGCTTCGCCAGACATGATCCGt              >  1:1321192/1‑150 (MQ=255)
                                                                                                                                          gctgTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAGAACCGAAGAGTTTGATGCGATCAAAATTGCTCTGGCTTCGCCAGCCATGATCCGTTc            >  2:1217947/1‑150 (MQ=255)
                                                                                                                                          gctgTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTGCTCTGGCTTCGCCAGACATGATCCGTTc            >  1:122675/1‑150 (MQ=255)
                                                                                                                                           ctgTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTGCTCTGGCTTCGCCAGACATGATCCGTTCa           <  1:663724/150‑1 (MQ=255)
                                                                                                                                                 ggTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTGCTCTGGCTTCg                             <  1:722319/126‑1 (MQ=255)
                                                                                                                                                  gTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTGCTCTGGCTTCGCCAGACATGATCCGTTCATGGTCtt    >  1:533884/1‑150 (MQ=255)
                                                                                                                                                    tAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCTTGTAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTGCTCTGGCTTCGCCAGACATGATCCGTTCATGGTCTTTc  <  1:96059/150‑1 (MQ=255)
                                                                                                                                                    tAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTGCTCTGGCTTCGCCAGACATGATCCGTTCATGGTCTTTc  <  2:775710/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAACCCGGCAGCGGATTGTGCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTGCTCTGGCTTCGCCAGACATGATCCGTTCATGGTCTTTC  >  NC_000913/4185157‑4185454

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