Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,600,994 A→T 6.7% intergenic (‑211/+34) rpoH ← / ← ftsX RNA polymerase, sigma 32 (sigma H) factor/putative ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,600,9940AT 6.7% 339.6 / 9.6 124intergenic (‑211/+34)rpoH/ftsXRNA polymerase, sigma 32 (sigma H) factor/putative ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (56/58);  new base T (6/3);  total (63/61)
Fisher's exact test for biased strand distribution p-value = 4.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.87e-01

Only 100 of 129 total aligned reads displayed.

AGCAACGAGGTTATCATTCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGA  >  NC_000913/3600849‑3601143
                                                                                                                                                   |                                                                                                                                                     
agCAACGAGGTTATCATTCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCa                                                                                                                                                   <  1:1813457/150‑1 (MQ=255)
agCAACGAGGTTATCATTCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCa                                                                                                                                                   <  2:221515/150‑1 (MQ=255)
  cAACGAGGTTATCATTCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCAtt                                                                                                                                                 >  1:1069437/1‑150 (MQ=255)
     cGAGGTTATCATTCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCa                                                                                                                                              <  2:1787621/150‑1 (MQ=255)
        ggTTATCATTCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAggg                                                                                                                                           <  2:1157476/150‑1 (MQ=255)
         gTTATCATTCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAAACCATTGCAGGGa                                                                                                                                          >  1:1735452/1‑150 (MQ=255)
           tATCATTCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGaaa                                                                                                                                        >  2:375832/1‑150 (MQ=255)
           tATCATTCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGaaa                                                                                                                                        <  2:919472/150‑1 (MQ=255)
           tATCATTCACTGTTTTATCAGACCGTGATTTTATCCAAAAGTTCTATGCAAGCTTGTGAATTAATTTCGCACAAAATGTGTCATAGAGATG‑A‑ATTACAGGGAATAGACAACGGGGTCGCTTTCCCTGCAATGGGTTCCGTAGCAGGGaaa                                                                                                                                        <  1:346020/150‑1 (MQ=255)
            aTCATTCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAag                                                                                                                                       <  2:1320772/150‑1 (MQ=255)
            aTCATTCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCAATGGGTTCCGTAGCAGGGAAag                                                                                                                                       <  2:1666355/150‑1 (MQ=255)
              cATTCACTGTTTTATCAGACAGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGAAATAGACATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAagag                                                                                                                                     >  1:497458/1‑150 (MQ=255)
               agtCACTGTTTAATCAGACCCTGATTTTATCCACAAGTTCAATGCAAGCTTCTTAATAAATAACGCACAAAATTTGAAATAGAGATG‑ATTATA‑AGGGAAGCGAAAACGGGGTCTCTTACCATTCTACGGCACCCATTGCAGGGAAagaga                                                                                                                                    >  2:1223171/3‑149 (MQ=255)
               aTTCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCGATTGCAGGGAAAGAGt                                                                                                                                    >  2:1191178/1‑150 (MQ=255)
                ttCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGta                                                                                                                                   <  1:300849/150‑1 (MQ=255)
                 tCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGtat                                                                                                                                  <  1:1431211/150‑1 (MQ=255)
                  cACTGTTTTATCAGACCATGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAA‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGtat                                                                                                                                  >  2:1482076/1‑150 (MQ=255)
                  cACTGTTATATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGAGG‑A‑AATACCGGGAAGAGACAACGGTGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGtata                                                                                                                                 >  1:286383/1‑150 (MQ=255)
                   aCTGTTTTATCAGACCGTGATTTTATCCACAAGATCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGATACGGAACACATTGCAGGGAAAGAGTATaa                                                                                                                                >  2:1191433/1‑150 (MQ=255)
                     tGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAaca                                                                                                                              <  1:1514814/150‑1 (MQ=255)
                      gTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAacac                                                                                                                             >  2:450736/1‑150 (MQ=255)
                      gTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAacac                                                                                                                             >  1:1194600/1‑150 (MQ=255)
                          tATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCtt                                                                                                                         >  2:1273546/1‑150 (MQ=255)
                              aGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTttat                                                                                                                     >  2:701504/1‑150 (MQ=255)
                                 ccGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGACCCCATTGCAGGGAAAGAGTAAAACCCGCTTTtattat                                                                                                                  >  2:1362004/1‑150 (MQ=255)
                                 ccGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTtattat                                                                                                                  <  2:1351581/150‑1 (MQ=255)
                                  cGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTattatt                                                                                                                 <  1:973978/150‑1 (MQ=255)
                                   gTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTattatta                                                                                                                >  1:1295521/1‑149 (MQ=255)
                                   gTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTc                                                                                                                <  1:1217484/150‑1 (MQ=255)
                                     taTTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAg                                                                                                              <  1:728445/149‑1 (MQ=255)
                                      aTTTTATCCACAAGTTCAATTCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAgg                                                                                                             <  1:519392/150‑1 (MQ=255)
                                      aTTTTATACACAAGTTCCATGCAAGCTTGTGAATAAATTACGCACAAAAGGTGACATAGAGATG‑A‑AATACCGTGAAGAGACAACGGGGGCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAgg                                                                                                             <  2:1268567/150‑1 (MQ=255)
                                         ttATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGt                                                                                                          <  2:839179/150‑1 (MQ=255)
                                          tATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTa                                                                                                         >  2:1026339/1‑150 (MQ=255)
                                            tCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTaaa                                                                                                       >  2:1287858/1‑150 (MQ=255)
                                              cacaAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGt                                                                                                     >  2:863053/1‑150 (MQ=255)
                                                caAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTgg                                                                                                   <  1:750792/150‑1 (MQ=255)
                                                   gTTCAATGCAAGCTTGTGAATAAATTCCGCACAAAATGGGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGt                                                                                                <  2:734788/150‑1 (MQ=255)
                                                   gTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGt                                                                                                >  1:1754093/1‑150 (MQ=255)
                                                   gTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGt                                                                                                >  1:1099294/1‑150 (MQ=255)
                                                         tGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAAtgt                                                                                          >  1:1770781/1‑150 (MQ=255)
                                                          gCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAAtgtt                                                                                         >  2:786594/1‑150 (MQ=255)
                                                          gCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAAtgtt                                                                                         <  2:944668/150‑1 (MQ=255)
                                                           cAAGATTGTGAATAAATAACGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGTCAACGGGGTCTCTTTCCCTGCTACGGATCCCATTGCAGTGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAATGTGGCGTAAAtgttg                                                                                        <  1:1191433/150‑1 (MQ=255)
                                                            aaGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAAtgttgt                                                                                       >  2:802411/1‑150 (MQ=255)
                                                            aaGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAAtgttgt                                                                                       <  1:1183434/150‑1 (MQ=255)
                                                                ttGTGAACAAATTACGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTATTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTAAGGCGTAAAGTTTCGTAAATGTTGTACCg                                                                                   >  1:773396/1‑150 (MQ=255)
                                                                   tgAATAAATTACGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTgg                                                                                >  2:1619857/1‑150 (MQ=255)
                                                                   tgAATAAATTACGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTgg                                                                                <  1:253390/150‑1 (MQ=255)
                                                                   tgAATAAATTACGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGAAGACACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCATTGAAATAGTAAAACACGCTTTTATTATTCAGGCGGAAAGTGGCGTAAATGTTGTACCGTgg                                                                                <  2:280389/150‑1 (MQ=255)
                                                                     aaTAAATTACGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTTCGGAACCCATTGCAGGGAAAGAGTATAACAAGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTTTACCGTGGCa                                                                              >  1:1623565/1‑150 (MQ=255)
                                                                           ttACGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCAc                                                                        <  1:739841/150‑1 (MQ=255)
                                                                             aCGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGc                                                                      <  1:199339/150‑1 (MQ=255)
                                                                             aCGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGc                                                                      <  1:600995/150‑1 (MQ=255)
                                                                              cGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCt                                                                     >  1:12694/1‑150 (MQ=255)
                                                                               ccacaAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTg                                                                    <  2:463936/149‑1 (MQ=255)
                                                                                cacaAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGc                                                                   >  2:1199078/1‑150 (MQ=255)
                                                                                 acaAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGcc                                                                  <  2:1308100/150‑1 (MQ=255)
                                                                                     aaTGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCAccc                                                              <  1:360836/150‑1 (MQ=255)
                                                                                       tgtgACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAg                                                            <  2:135591/150‑1 (MQ=255)
                                                                                       tgggACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAg                                                            <  2:1031921/147‑1 (MQ=255)
                                                                                          gACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGAGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCa                                                         >  2:1065384/1‑150 (MQ=255)
                                                                                            cATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAat                                                       <  1:500180/150‑1 (MQ=255)
                                                                                                 agaTG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCg                                                  <  2:746082/150‑1 (MQ=255)
                                                                                                   aTG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAg                                                <  1:1512804/150‑1 (MQ=255)
                                                                                                    tG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGc                                               >  1:1374838/1‑150 (MQ=255)
                                                                                                    tG‑A‑AATACCGGGAAGAGACAAAGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCCAGCCACGCTGCCACCCAGCCAATCATCCAGc                                               >  2:136482/1‑150 (MQ=255)
                                                                                                            aCCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACCCGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATCTTGTACCGTGGCCAGCCACGCTGACCCCCAGCCAATCATCGAGCATACca                                         >  2:330032/1‑150 (MQ=255)
                                                                                                            aCCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACca                                         <  1:735987/150‑1 (MQ=255)
                                                                                                             ccGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACcag                                        >  2:1199300/1‑150 (MQ=255)
                                                                                                               gggAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACcagca                                      >  2:1485121/1‑150 (MQ=255)
                                                                                                                 gAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCagcagc                                    <  1:1089580/150‑1 (MQ=255)
                                                                                                                 gAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGAATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCagcagc                                    >  1:522089/1‑150 (MQ=255)
                                                                                                                   agagACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCaa                                  <  2:1304712/150‑1 (MQ=255)
                                                                                                                    gagaCAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCATGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGATGCCACCCATCCAATAAGAGAGAATACCAGCAGCAAt                                 >  2:93156/1‑150 (MQ=255)
                                                                                                                       aCAACTGGGTCTCTTTCCCTGCTACTGAACCCATTGCAGGGAAAGAGTATAACACGCTTTAATAATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGc                              <  1:294406/150‑1 (MQ=255)
                                                                                                                        cAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCa                             <  2:689066/150‑1 (MQ=255)
                                                                                                                        cAACGGGGGCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGGGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCa                             <  2:1706635/150‑1 (MQ=255)
                                                                                                                              ggTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGTCGTAAAGTTGCGTAAATGTTGTACCGTGGCAAGCCAAGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCAtt                       <  2:286383/150‑1 (MQ=255)
                                                                                                                              atactctTTCCCTGCAATGGGTTCCGTAGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCAtt                       >  2:394828/4‑150 (MQ=38)
                                                                                                                              acactctTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCAtt                       >  2:1569875/4‑150 (MQ=255)
                                                                                                                                tctcTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCa                     <  2:885453/150‑1 (MQ=255)
                                                                                                                                 ctctTTCCCTGCTACGGAACCCATTGCAGGGTAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACAGTGGCAAGCCACGCTGCCACCCAGCCACTCATCGAGCATACCAGCAGCAATAGCAGGCATTCAt                    >  1:45390/1‑150 (MQ=255)
                                                                                                                                 ctctTTCCCTGCAATGGGTTCCGTAGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCAt                    >  2:150487/1‑150 (MQ=39)
                                                                                                                                  tctTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCACGCCACGCTGCCACCCAGCCACTCATCGAGCATACCAGCAGCAATAGCAGGCATTCATc                   >  1:767457/1‑150 (MQ=255)
                                                                                                                                  tctTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATc                   <  2:513865/150‑1 (MQ=255)
                                                                                                                                  tctTTCCCTGCAATGGGTTCCGTAGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATc                   >  1:1332484/1‑150 (MQ=39)
                                                                                                                                   ctTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTCCCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCGTACCTGCAGCAATAGCAGGCATTCATCg                  >  1:1132172/1‑150 (MQ=255)
                                                                                                                                    tttCCCTGCAATGGGTTCCGTAGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGa                 >  1:372465/1‑150 (MQ=255)
                                                                                                                                        ccTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGa             >  2:862323/1‑150 (MQ=255)
                                                                                                                                        ccTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGa             >  2:1940/1‑150 (MQ=255)
                                                                                                                                         cTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGAt            >  1:1092682/1‑150 (MQ=255)
                                                                                                                                         cTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGAt            >  1:1703161/1‑150 (MQ=255)
                                                                                                                                             tACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGACACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGc        >  1:223792/1‑150 (MQ=255)
                                                                                                                                                ggAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCAt     <  2:1770781/150‑1 (MQ=255)
                                                                                                                                                ggAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCAAt     >  1:1100355/1‑150 (MQ=255)
                                                                                                                                                  aaCCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTg   <  2:367180/150‑1 (MQ=255)
                                                                                                                                                  aaCCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGACAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTg   >  2:1227422/1‑150 (MQ=255)
                                                                                                                                                   aCCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGa  <  2:666854/150‑1 (MQ=255)
                                                                                                                                                   aCCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGa  <  1:802411/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
AGCAACGAGGTTATCATTCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATG‑A‑AATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGA  >  NC_000913/3600849‑3601143

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