Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,600,990 C→T 6.7% intergenic (‑207/+38) 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,9900CT 6.7% 374.9 / 10.2 119intergenic (‑207/+38)rpoH/ftsXRNA polymerase, sigma 32 (sigma H) factor/putative ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base C (54/57);  new base T (6/2);  total (60/59)
Fisher's exact test for biased strand distribution p-value = 2.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.94e-01

Only 100 of 120 total aligned reads displayed.

AGCTTAAGAGCAACGAGGTTATCATTCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGC  >  NC_000913/3600841‑3601137
                                                                                                                                                       |                                                                                                                                                   
        agCAACGAGGTTATCATTCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCa                                                                                                                                             <  2:221515/150‑1 (MQ=255)
        agCAACGAGGTTATCATTCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCa                                                                                                                                             <  1:1813457/150‑1 (MQ=255)
          cAACGAGGTTATCATTCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCAtt                                                                                                                                           >  1:1069437/1‑150 (MQ=255)
           aaCGAGGTTATCATTCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACACACAAAATGTGACATAGAGATGAATTACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTg                                                                                                                                          <  1:297685/150‑1 (MQ=255)
             cGAGGTTATCATTCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCa                                                                                                                                        <  2:1787621/150‑1 (MQ=255)
                ggTTATCATTCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAggg                                                                                                                                     <  2:1157476/150‑1 (MQ=255)
                 gTTATCATTCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTAGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGa                                                                                                                                    <  2:1196783/150‑1 (MQ=255)
                   tATCATTCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGaaa                                                                                                                                  <  2:919472/150‑1 (MQ=255)
                   tATCATTCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGaaa                                                                                                                                  >  2:375832/1‑150 (MQ=255)
                   tATCATTCACTGTTTTATCAGACCGTGATTTTATCCAAAAGTTCTATGCAAGCTTGTGAATTAATTTCGCACAAAATGTGTCATAGAGATGAATTACAGGGAATAGACAACGGGGTCGCTTTCCCTGCAATGGGTTCCGTAGCAGGGaaa                                                                                                                                  <  1:346020/150‑1 (MQ=255)
                    aTCATTCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAag                                                                                                                                 <  2:1320772/150‑1 (MQ=255)
                    aTCATTCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCAATGGGTTCCGTAGCAGGGAAag                                                                                                                                 <  2:1666355/150‑1 (MQ=255)
                      cATTCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGATTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAagag                                                                                                                               <  2:895277/150‑1 (MQ=255)
                       agtCACTGTTTAATCAGACCCTGATTTTATCCACAAGTTCAATGCAAGCTTCTTAATAAATAACGCACAAAATTTGAAATAGAGATGATTATAAGGGAAGCGAAAACGGGGTCTCTTACCATTCTACGGCACCCATTGCAGGGAAagaga                                                                                                                              >  2:1223171/3‑149 (MQ=255)
                       aTTCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCGATTGCAGGGAAAGAGt                                                                                                                              >  2:1191178/1‑150 (MQ=255)
                         tCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGtat                                                                                                                            <  1:1431211/150‑1 (MQ=255)
                          cACTGTTTTATCAGACCATGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGtat                                                                                                                            >  2:1482076/1‑150 (MQ=255)
                          cACTGTTATATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGAGGAAATACCGGGAAGAGACAACGGTGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGtata                                                                                                                           >  1:286383/1‑150 (MQ=255)
                           aCTGTTTTATCAGACCGTGATTTTATCCACAAGATCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGATACGGAACACATTGCAGGGAAAGAGTATaa                                                                                                                          >  2:1191433/1‑150 (MQ=255)
                            cTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAac                                                                                                                         >  1:723189/1‑150 (MQ=255)
                             tGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAaca                                                                                                                        <  1:1514814/150‑1 (MQ=255)
                              gTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAacac                                                                                                                       >  2:450736/1‑150 (MQ=255)
                              gTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAacac                                                                                                                       >  1:1194600/1‑150 (MQ=255)
                                 ttATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCt                                                                                                                    <  1:947429/150‑1 (MQ=255)
                                  tATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCtt                                                                                                                   >  2:1273546/1‑150 (MQ=255)
                                      aGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTttat                                                                                                               >  2:701504/1‑150 (MQ=255)
                                      aGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTttat                                                                                                               >  1:176582/1‑150 (MQ=255)
                                         ccGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGACCCCATTGCAGGGAAAGAGTAAAACCCGCTTTtattat                                                                                                            >  2:1362004/1‑150 (MQ=255)
                                         ccGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTtattat                                                                                                            <  2:1351581/150‑1 (MQ=255)
                                          cGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTattatt                                                                                                           <  1:973978/150‑1 (MQ=255)
                                           gTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTattatta                                                                                                          >  1:1295521/1‑149 (MQ=255)
                                           gTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTc                                                                                                          <  1:1217484/150‑1 (MQ=255)
                                            tGATTTTATCCACAAGTTCAATGCAAGTTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGTGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCa                                                                                                         <  2:13547/150‑1 (MQ=255)
                                              aTTTTATCCACAAGTTCAATTCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAgg                                                                                                       <  1:519392/150‑1 (MQ=255)
                                              aTTTTATACACAAGTTCCATGCAAGCTTGTGAATAAATTACGCACAAAAGGTGACATAGAGATGAAATACCGTGAAGAGACAACGGGGGCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAgg                                                                                                       <  2:1268567/150‑1 (MQ=255)
                                                 ttATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGt                                                                                                    <  2:839179/150‑1 (MQ=255)
                                                  tATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTa                                                                                                   >  2:1026339/1‑150 (MQ=255)
                                                    tCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTaaa                                                                                                 >  2:1287858/1‑150 (MQ=255)
                                                     ccACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAg                                                                                                >  1:25132/1‑150 (MQ=255)
                                                      cacaAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGt                                                                                               >  2:863053/1‑150 (MQ=255)
                                                        caAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTgg                                                                                             <  1:750792/150‑1 (MQ=255)
                                                           gTTCAATGCAAGCTTGTGAATAAATTCCGCACAAAATGGGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGt                                                                                          <  2:734788/150‑1 (MQ=255)
                                                           gTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGt                                                                                          >  1:1754093/1‑150 (MQ=255)
                                                           gTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGt                                                                                          >  1:1099294/1‑150 (MQ=255)
                                                           gTTCAATGCAAGCTTGTGAATAAATTACGCACAAAAAGTGACAAAGAGATGAAATACCGGGAAGACACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTTCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGAAAAGTGGAGt                                                                                          >  1:1423955/1‑150 (MQ=255)
                                                            ttCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTa                                                                                         <  2:668913/150‑1 (MQ=255)
                                                                aTGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGTAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAAtg                                                                                     <  2:137702/150‑1 (MQ=255)
                                                                 tGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAAtgt                                                                                    >  1:1770781/1‑150 (MQ=255)
                                                                  gCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAAtgtt                                                                                   >  2:786594/1‑150 (MQ=255)
                                                                   cAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAAtgttg                                                                                  <  2:442619/150‑1 (MQ=255)
                                                                   cAAGATTGTGAATAAATAACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGTCAACGGGGTCTCTTTCCCTGCTACGGATCCCATTGCAGTGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAATGTGGCGTAAAtgttg                                                                                  <  1:1191433/150‑1 (MQ=255)
                                                                    aaGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAAtgttgt                                                                                 >  2:802411/1‑150 (MQ=255)
                                                                    aaGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAAtgttgt                                                                                 <  1:1183434/150‑1 (MQ=255)
                                                                           tgAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTgg                                                                          <  1:253390/150‑1 (MQ=255)
                                                                           tgAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTgg                                                                          >  2:1619857/1‑150 (MQ=255)
                                                                           tgAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGACACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCATTGAAATAGTAAAACACGCTTTTATTATTCAGGCGGAAAGTGGCGTAAATGTTGTACCGTgg                                                                          <  2:280389/150‑1 (MQ=255)
                                                                             aaTAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTTCGGAACCCATTGCAGGGAAAGAGTATAACAAGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTTTACCGTGGCa                                                                        >  1:1623565/1‑150 (MQ=255)
                                                                                   ttACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCAc                                                                  <  1:739841/150‑1 (MQ=255)
                                                                                     aCGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGc                                                                <  1:600995/150‑1 (MQ=255)
                                                                                     aCGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGc                                                                <  1:199339/150‑1 (MQ=255)
                                                                                      cGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCt                                                               >  1:12694/1‑150 (MQ=255)
                                                                                      cGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCt                                                               >  2:1240382/1‑150 (MQ=255)
                                                                                        cacaAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTTAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGc                                                             >  1:763208/1‑150 (MQ=255)
                                                                                        cacaAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGc                                                             >  2:1199078/1‑150 (MQ=255)
                                                                                         acaAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGcc                                                            <  2:1308100/150‑1 (MQ=255)
                                                                                             aaTGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCAccc                                                        <  1:360836/150‑1 (MQ=255)
                                                                                               tgggACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAg                                                      <  2:1031921/147‑1 (MQ=255)
                                                                                                  gACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGAGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCa                                                   >  2:1065384/1‑150 (MQ=255)
                                                                                                    cATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAat                                                 <  1:500180/150‑1 (MQ=255)
                                                                                                         agaTGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCg                                            <  2:746082/150‑1 (MQ=255)
                                                                                                           aTGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAg                                          <  1:1512804/150‑1 (MQ=255)
                                                                                                            tGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGc                                         >  1:384940/1‑150 (MQ=255)
                                                                                                            tGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGc                                         >  1:1374838/1‑150 (MQ=255)
                                                                                                            tGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGc                                         >  1:1651745/1‑150 (MQ=255)
                                                                                                            tGAAATACCGGGAAGAGACAAAGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCCAGCCACGCTGCCACCCAGCCAATCATCCAGc                                         >  2:136482/1‑150 (MQ=255)
                                                                                                             gAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCa                                        <  2:840527/150‑1 (MQ=255)
                                                                                                                    aCCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACCCGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATCTTGTACCGTGGCCAGCCACGCTGACCCCCAGCCAATCATCGAGCATACca                                   >  2:330032/1‑150 (MQ=255)
                                                                                                                    aCCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACca                                   <  1:735987/150‑1 (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)
                                                                                                                                 aaCGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAg                      <  2:204497/150‑1 (MQ=255)
                                                                                                                                     tatactctTTCCCTGCAATGGGTTCCGTAGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCAt                  >  1:1100060/5‑150 (MQ=38)
                                                                                                                                      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)
                                                                                                                                          tctTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCACGCCACGCTGCCACCCAGCCACTCATCGAGCATACCAGCAGCAATAGCAGGCATTCATc             >  1:767457/1‑150 (MQ=255)
                                                                                                                                          tctTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATc             <  2:513865/150‑1 (MQ=255)
                                                                                                                                          tctTTCCCTGCAATGGGTTCCGTAGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATc             >  1:1332484/1‑150 (MQ=39)
                                                                                                                                           ctTTCCCTGCAATGGGTTCCGTAGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCg            >  2:1508648/1‑150 (MQ=255)
                                                                                                                                            tttCCCTGCAATGGGTTCCGTAGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGa           >  1:372465/1‑150 (MQ=255)
                                                                                                                                                ccTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGa       >  2:1940/1‑150 (MQ=255)
                                                                                                                                                 cTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGAt      >  1:1703161/1‑150 (MQ=255)
                                                                                                                                                 cTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGAt      >  1:1092682/1‑150 (MQ=255)
                                                                                                                                                       |                                                                                                                                                   
AGCTTAAGAGCAACGAGGTTATCATTCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGC  >  NC_000913/3600841‑3601137

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: