Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,600,993 A→G 6.5% intergenic (‑210/+35) 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,9930AG 6.5% 363.5 / 11.7 124intergenic (‑210/+35)rpoH/ftsXRNA polymerase, sigma 32 (sigma H) factor/putative ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (56/59);  new base G (6/2);  total (63/61)
Fisher's exact test for biased strand distribution p-value = 2.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.38e-01

Only 100 of 127 total aligned reads displayed.

AGCAACGAGGTTATCATTCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTG  >  NC_000913/3600849‑3601142
                                                                                                                                                 |                                                                                                                                                     
agCAACGAGGTTATCATTCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCa                                                                                                                                                  <  2:221515/150‑1 (MQ=255)
  cAACGAGGTTATCATTCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCAtt                                                                                                                                                >  1:1069437/1‑150 (MQ=255)
   aaCGAGGTTATCATTCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACACACAAAATGTGACATAGAGATGAATTACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTg                                                                                                                                               <  1:297685/150‑1 (MQ=255)
        ggTTATCATTCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAggg                                                                                                                                          <  2:1157476/150‑1 (MQ=255)
         gTTATCATTCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAAACCATTGCAGGGa                                                                                                                                         >  1:1735452/1‑150 (MQ=255)
           tATCATTCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGaaa                                                                                                                                       <  2:919472/150‑1 (MQ=255)
            aTCATTCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAag                                                                                                                                      <  2:1320772/150‑1 (MQ=255)
            aTCATTCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCAATGGGTTCCGTAGCAGGGAAag                                                                                                                                      <  2:1666355/150‑1 (MQ=255)
               agtCACTGTTTAATCAGACCCTGATTTTATCCACAAGTTCAATGCAAGCTTCTTAATAAATAACGCACAAAATTTGAAATAGAGATGATTATAAGGGAAGCGAAAACGGGGTCTCTTACCATTCTACGGCACCCATTGCAGGGAAagaga                                                                                                                                   >  2:1223171/3‑149 (MQ=255)
               aTTCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCGATTGCAGGGAAAGAGt                                                                                                                                   >  2:1191178/1‑150 (MQ=255)
                ttCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGta                                                                                                                                  <  1:300849/150‑1 (MQ=255)
                  cACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGtata                                                                                                                                >  2:735987/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)
                          tATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCtt                                                                                                                        >  2:1273546/1‑150 (MQ=255)
                              aGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTttat                                                                                                                    >  1:176582/1‑150 (MQ=255)
                              aGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTttat                                                                                                                    >  2:701504/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)
                                   gTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTc                                                                                                               <  1:1217484/150‑1 (MQ=255)
                                     taTTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAg                                                                                                             <  1:728445/149‑1 (MQ=255)
                                      aTTTTATCCACAAGTTCAATTCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAgg                                                                                                            <  1:519392/150‑1 (MQ=255)
                                      aTTTTATACACAAGTTCCATGCAAGCTTGTGAATAAATTACGCACAAAAGGTGACATAGAGATGAAATACCGTGAAGAGACAACGGGGGCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAgg                                                                                                            <  2:1268567/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)
                                                   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)
                                                          gCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAAtgtt                                                                                        <  2:944668/150‑1 (MQ=255)
                                                           cAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAAtgttg                                                                                       <  2:442619/150‑1 (MQ=255)
                                                            aaGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAAtgttgt                                                                                      >  2:802411/1‑150 (MQ=255)
                                                            aaGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAAtgttgt                                                                                      <  1:1183434/150‑1 (MQ=255)
                                                                ttGTGAACAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTATTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTAAGGCGTAAAGTTTCGTAAATGTTGTACCg                                                                                  >  1:773396/1‑150 (MQ=255)
                                                                     aaTAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTTCGGAACCCATTGCAGGGAAAGAGTATAACAAGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTTTACCGTGGCa                                                                             >  1:1623565/1‑150 (MQ=255)
                                                                             aCGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGc                                                                     <  1:199339/150‑1 (MQ=255)
                                                                              cGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCt                                                                    >  1:12694/1‑150 (MQ=255)
                                                                               gCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTg                                                                   <  2:1303780/150‑1 (MQ=255)
                                                                               ccacaAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTg                                                                   <  2:463936/149‑1 (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)
                                                                                       tgtgACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAg                                                           <  2:135591/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:1374838/1‑150 (MQ=255)
                                                                                                    tGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGc                                              >  1:1651745/1‑150 (MQ=255)
                                                                                                    tGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGc                                              >  1:384940/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)
                                                                                                            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)
                                                                                                                   gagaCAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCATGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGATGCCACCCATCCAATAAGAGAGAATACCAGCAGCAAt                                >  2:93156/1‑150 (MQ=255)
                                                                                                                    agaCAACGGTGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATa                               <  2:1736762/150‑1 (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)
                                                                                                                         aCGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAgg                          <  1:1273546/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)
                                                                                                                                 tctTTCCCTGCAATGGGTTCCGTAGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATc                  >  1:1332484/1‑150 (MQ=39)
                                                                                                                                  ctTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTCCCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCGTACCTGCAGCAATAGCAGGCATTCATCg                 >  1:1132172/1‑150 (MQ=255)
                                                                                                                                  ctTTCCCTGCAATGGGTTCCGTAGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCg                 >  2:1508648/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)
                                                                                                                                            tACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGACACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGc       >  1:223792/1‑150 (MQ=255)
                                                                                                                                               ggAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCAt    >  2:360129/1‑150 (MQ=255)
                                                                                                                                               ggAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCAt    <  2:1770781/150‑1 (MQ=255)
                                                                                                                                               ggAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCAAt    >  1:1100355/1‑150 (MQ=255)
                                                                                                                                                gAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCAtt   <  1:1191178/150‑1 (MQ=255)
                                                                                                                                                gAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGAGTAAATGTTGTACCGTGGCAAGCCACGCAGCCACACAGCCAATCATCGAGCATACCAGCAGCAAAATCATGCATTCATCCAATGATAAGCCaat   >  1:516182/1‑148 (MQ=255)
                                                                                                                                                 aaCCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTg  <  2:367180/150‑1 (MQ=255)
                                                                                                                                                 aaCCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGACAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTg  >  2:1227422/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
AGCAACGAGGTTATCATTCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTG  >  NC_000913/3600849‑3601142

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: