Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,214,770 A→C E575A (GAA→GCA)  rpoD → RNA polymerase, sigma 70 (sigma D) factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,214,7700AC98.6% 268.2 / ‑4.2 69E575A (GAA→GCA) rpoDRNA polymerase, sigma 70 (sigma D) factor
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base C (38/30);  minor base T (1/0);  total (39/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.97e-01

CCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCG  >  NC_000913/3214623‑3214916
                                                                                                                                                   |                                                                                                                                                  
ccaccCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGaagcag                                                                                                                                                  >  2:917660/1‑150 (MQ=255)
   cccTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGCAGTggg                                                                                                                                              >  2:1288818/1‑151 (MQ=255)
   cccTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGCAGTgg                                                                                                                                               <  1:145752/150‑1 (MQ=255)
     cTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGCAGTGGGt                                                                                                                                             >  2:1165052/1‑150 (MQ=255)
       cGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGCAGTGGGTaa                                                                                                                                           <  1:96290/150‑1 (MQ=255)
        gAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGCAGTGGGTAAAc                                                                                                                                         >  1:888609/1‑151 (MQ=255)
          gCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGCAGTGGGTAAACAg                                                                                                                                       <  2:336274/151‑1 (MQ=255)
           cTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGCAGTGGGTAAACAGt                                                                                                                                      >  1:446470/1‑151 (MQ=255)
            tGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGCAGTGGGTAAACAGtt                                                                                                                                     >  1:634948/1‑151 (MQ=255)
             gCCGCTGGATTCTGCGACCACCGAAGGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGCAGTGGGTAAACAGTTc                                                                                                                                    <  2:1240286/151‑1 (MQ=255)
               cGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGCAGTGGGTAAACAGTTCg                                                                                                                                   >  1:720223/1‑150 (MQ=255)
                   ggATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGCAGTGGGTAAACAGTTCGACGtt                                                                                                                              <  1:1197980/151‑1 (MQ=255)
                       tCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGCAGTGGGTAAACAGTTCGACGTTAccc                                                                                                                          >  1:560783/1‑151 (MQ=255)
                           cGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGCAGTGGGTAAACAGTTCGACGTTACCcgcg                                                                                                                       >  2:450944/1‑150 (MQ=255)
                             accaccGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGCAGTGGGTAAACAGTTCGACGTTACCCGCGAAc                                                                                                                    <  2:504926/151‑1 (MQ=255)
                              ccaccGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGCAGTGGGTAAACAGTTCGACGTTACCCGCGAACg                                                                                                                   <  2:1420387/151‑1 (MQ=255)
                                accGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGCAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTa                                                                                                                 <  1:712373/151‑1 (MQ=255)
                                  cGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGCAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATc                                                                                                               >  1:27358/1‑151 (MQ=255)
                                      agactgcctgcGCCAACGCACCACGTGATGCCTGGACTGACTCCGCGTGAACCAAAAGCTCTGGGCAGGCGTTTCGGGATAGATATGAACACCGACTACACGCTGGAAGCAGTGGGGAAACAGTTCGACGCTACCCGCGCAAGTATGCGTc                                                                                                           >  2:113835/4‑151 (MQ=255)
                                       gCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGCAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCa                                                                                                          <  2:937590/151‑1 (MQ=255)
                                           gcgtgcgGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGCAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGa                                                                                                        <  1:277063/149‑1 (MQ=255)
                                           gcgtgcgGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGCAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATc                                                                                                      <  2:390090/151‑1 (MQ=255)
                                            cgtgcgGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGCAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATcg                                                                                                     <  1:797707/151‑1 (MQ=255)
                                                cgGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGCAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATcgaagc                                                                                                 <  1:831523/151‑1 (MQ=255)
                                                    aaCGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGCAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATcgaagcga                                                                                               >  2:754946/1‑149 (MQ=255)
                                                      cGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGCAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGgc                                                                                           <  1:1378556/151‑1 (MQ=255)
                                                            acgTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGCAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTgcg                                                                                     <  2:265986/151‑1 (MQ=255)
                                                                gctggctgGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGCAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCaa                                                                                  <  2:720210/150‑1 (MQ=255)
                                                                 ctggctggCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGCAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAAc                                                                                >  1:932614/1‑151 (MQ=255)
                                                                  tggctggcCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGCAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACt                                                                               >  2:877774/1‑151 (MQ=255)
                                                                   ggctggcCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGCAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTg                                                                              <  1:691705/151‑1 (MQ=255)
                                                                       ggcCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGCAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGt                                                                           <  2:1295505/150‑1 (MQ=255)
                                                                          cTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGCAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCAcc                                                                       >  2:1316652/1‑151 (MQ=255)
                                                                           tGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGCAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCAccc                                                                      >  1:420700/1‑151 (MQ=255)
                                                                             aCCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGCAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCg                                                                     >  1:205044/1‑150 (MQ=255)
                                                                                 cgcgcgAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGCAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCg                                                                >  2:810520/1‑151 (MQ=255)
                                                                                 cgcgTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGCAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCg                                                                >  1:77423/1‑151 (MQ=255)
                                                                                    gTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGCAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTc                                                             >  2:433442/1‑151 (MQ=255)
                                                                                      gAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGCAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTg                                                           >  2:275526/1‑151 (MQ=255)
                                                                                          cAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGCAGTGGGTAAACACTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGt                                                       <  1:529065/151‑1 (MQ=255)
                                                                                            aaaGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGCAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGtgc                                                     >  1:183185/1‑151 (MQ=255)
                                                                                               gTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGCAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGtgctgc                                                  >  1:320640/1‑151 (MQ=255)
                                                                                               gTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGCAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGtgctgc                                                  <  2:938489/151‑1 (MQ=255)
                                                                                                ttCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGCAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCg                                                 >  1:882557/1‑151 (MQ=255)
                                                                                                          gCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGCAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTg                                       <  2:41535/151‑1 (MQ=255)
                                                                                                           cGTTTCGGTATCGATATGAACACCAACTACACGCTGGAAGCAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTgg                                      >  2:782641/1‑151 (MQ=255)
                                                                                                            gTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGCAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTgg                                      >  1:95766/1‑150 (MQ=255)
                                                                                                              ttCGGTTTCGATATAATCTCCGCCTACACGCTGAAAGCAGTGGGTAAACAGTTCGACGTTCCCTGCGAACGTACCCGACAGAGCGAGGCGAAGGCGCTGCGCAAACAGCCTCACCCGCGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGAc                                    <  1:1250489/150‑1 (MQ=255)
                                                                                                                cGGTATCGATATGAACACCGACTACACGCTGGAAGCAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGAt                                 >  2:889076/1‑151 (MQ=255)
                                                                                                                 ggTATCGATATGAACACCGACTACACGCTGGAAGCAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGAtt                                >  1:729388/1‑151 (MQ=255)
                                                                                                                  gTATCGATATGAACACCGACTACACGCTGGAAGCAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTa                               <  1:362620/151‑1 (MQ=255)
                                                                                                                    aTCGATATGAACACCGACTACACGCTGGAAGCAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAAt                             >  2:783207/1‑151 (MQ=255)
                                                                                                                      cGATATGAACACCGACTACACGCTGGAAGCAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCg                           >  2:96305/1‑151 (MQ=255)
                                                                                                                       gATATGAACGCCGGCTACGCGCGGGAAGCGGTGGGTAAACAGTTCGGCGTTACCGGCGAGCGTGTCCGTCAGATCGAGGGGAGGGCGCGGCGCAAGCTGGGTCGCCCGAGCCGTTCTGAAGTGCTGCGGGGCGTCCTGGACGATTaagc                            <  2:438268/149‑3 (MQ=255)
                                                                                                                       gATATGAACACCGACTACACGCTGGAAGCAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCgg                          >  1:1211531/1‑151 (MQ=255)
                                                                                                                          atGAACACCGACTACACGCTGGAAGCAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTa                        <  1:130176/150‑1 (MQ=255)
                                                                                                                            gAACACCGACTACACGCTGGAAGCAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAgg                      <  1:450951/150‑1 (MQ=255)
                                                                                                                               cacCGACTACACGCTGGAAGCAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCgg                  <  2:634940/151‑1 (MQ=255)
                                                                                                                                acCGACTACACGCTGGAAGCAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCgg                  >  1:949633/1‑150 (MQ=255)
                                                                                                                                 ccGACTACACGCTGGAAGCAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGAt                <  2:321129/151‑1 (MQ=255)
                                                                                                                                   gACTACACGCTGGAAGCAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATc               >  1:616292/1‑150 (MQ=255)
                                                                                                                                   gACTACACGCTGGAAGCAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCa              <  1:864599/151‑1 (MQ=255)
                                                                                                                                   gACTACACGCTGGAAGCAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCa              <  2:484610/151‑1 (MQ=255)
                                                                                                                                      tACACGCTGGAAGCAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGc           >  2:52056/1‑151 (MQ=255)
                                                                                                                                       acacGCTGGAAGCAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCg          >  2:101859/1‑151 (MQ=255)
                                                                                                                                        cacGCTGGAAGCAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGt         >  1:1276201/1‑151 (MQ=255)
                                                                                                                                         acGCTGGAAGTAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGtt        >  1:1359798/1‑151 (MQ=255)
                                                                                                                                            cTGGAAGCAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTAcg     >  2:192755/1‑151 (MQ=255)
                                                                                                                                              ggAAGCAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTAcgcc   >  2:515950/1‑151 (MQ=255)
                                                                                                                                              ggAAGCAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTAcgcc   >  1:33484/1‑151 (MQ=255)
                                                                                                                                               gaagcaGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTAcgccg  >  2:310137/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCG  >  NC_000913/3214623‑3214916

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑