Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,922,483 G→A intergenic (‑43/+574) atpI ← / ← rsmG ATP synthase, membrane‑bound accessory factor/16S rRNA m(7)G527 methyltransferase, SAM‑dependent; glucose‑inhibited cell‑division protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,922,4830GA97.0% 245.3 / ‑3.4 67intergenic (‑43/+574)atpI/rsmGATP synthase, membrane‑bound accessory factor/16S rRNA m(7)G527 methyltransferase, SAM‑dependent; glucose‑inhibited cell‑division protein
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (31/34);  total (32/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTCTTTGAGGCTGAACAGCAATCCACTTGCTATCACCACCAGTAACTGAACGAGCAGAAGCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCT  >  NC_000913/3922341‑3922627
                                                                                                                                              |                                                                                                                                                
ggTCTTTGAGGCTGAACAGCAATCCACTTGCTATCACCACCAGTAACTGAACGAGCAGAAGCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAAc                                                                                                                                          >  1:503172/1‑151 (MQ=255)
 gTCTTTGAGGCTGAACAGCAATCCACTTGCTATCACCACCAGTAACTGAACGAGCAGAAGCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAAct                                                                                                                                         <  1:1170190/151‑1 (MQ=255)
    tttGAGGCTGAACAGCAATCCACTTGCTATCACCACCAGTAACTGAACGAGCAGAAGCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAActttc                                                                                                                                      >  2:788598/1‑151 (MQ=255)
      tGAGGCTGAACAGCAATCCACTTGCTATCACCACCAGTAACTGAACGAGCAGAAGCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAActttctt                                                                                                                                    >  2:397574/1‑151 (MQ=255)
       gAGGCTGAACAGCAATCCACTTGCTATCACCACCAGTAACTGAACGAGCAGAAGCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAActttctt                                                                                                                                    <  2:1159247/150‑1 (MQ=255)
        aGGCTGAACAGCAATCCACTTGCTATCACCACCAGTAACTGAACGAGCAGAAGCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTa                                                                                                                                  <  2:484598/151‑1 (MQ=255)
          gCTGAACAGCAATCCACTTGCTATCACCACCAGTAACTGAACGAGCAGAAGCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAg                                                                                                                                <  1:442543/151‑1 (MQ=255)
           cTGAACAGCAATCCACTTGCTATCACCACCAGTAACTGAACGAGCAGAAGCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAgg                                                                                                                               >  1:735861/1‑151 (MQ=255)
             gAACAGCAATCCACTTGCTATCACCACCAGTAACTGAACGAGCAGAAGCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGc                                                                                                                              <  1:762607/150‑1 (MQ=255)
                    aaTCCACTTGCTATCACCACCAGTAACTGAACGAGCAGAAGCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTc                                                                                                                      <  2:336588/151‑1 (MQ=255)
                    aaTCCAATTGCTATCACCACCAGTAACTGAACGAGCAGAAGCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTc                                                                                                                      >  2:338837/1‑151 (MQ=255)
                     aTCCACTTGCTATCACCACCAGTAACTGAACGAGCAGAAGCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTc                                                                                                                      >  2:391182/1‑150 (MQ=255)
                                aTCACCACCAGTAACTGAACGAGCAGAAGCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCaaaa                                                                                                           >  2:120918/1‑150 (MQ=255)
                                aTCACCACCAGTAACTGAACGAGCAGAAGCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCaaaa                                                                                                           >  1:1116814/1‑150 (MQ=255)
                                  caccacCAGTAACTGAACGAGCAGAAGCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGc                                                                                                        <  1:721438/151‑1 (MQ=255)
                                  caccacCAGTAACTGAACGAGCAGAAGCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGc                                                                                                        <  2:304632/151‑1 (MQ=255)
                                   accaccAGTAACTGAACGAGCAGAAGCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCg                                                                                                       <  1:178780/151‑1 (MQ=255)
                                      accaGTAACTGAACGAGCAGAAGCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTc                                                                                                    <  1:884390/151‑1 (MQ=255)
                                      accaGTAACTGAACGAGCAGAAGCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTc                                                                                                    >  1:670508/1‑151 (MQ=255)
                                       ccaGTAACTGAACGAGCAGAAGCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTc                                                                                                    >  2:176951/1‑150 (MQ=255)
                                        caGTAACTGAACGAGCAGAAGCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCaa                                                                                                  <  1:480483/151‑1 (MQ=255)
                                            aaCTGAACGAGCAGAAGCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAAtt                                                                                               >  2:331134/1‑150 (MQ=255)
                                            aaCTGAACGAGCAGAAGCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAAtt                                                                                               <  2:633095/150‑1 (MQ=255)
                                            aaCTGAACGAGCAGAAGCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATta                                                                                              <  1:913413/151‑1 (MQ=255)
                                             aCTGAACGAGCAGAAGCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATtat                                                                                             <  1:669529/151‑1 (MQ=255)
                                                gAACGAGCAGAAGCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACg                                                                                          <  2:314002/151‑1 (MQ=255)
                                                       cagattcTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATCTAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGcc                                                                                   >  2:1009260/7‑151 (MQ=255)
                                                         gAAGCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGcccc                                                                                 >  1:336088/1‑151 (MQ=255)
                                                          aaGCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGccccc                                                                                >  1:662598/1‑151 (MQ=255)
                                                          aaGCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGccccc                                                                                <  1:1064280/151‑1 (MQ=255)
                                                               tCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGa                                                                            <  1:397681/150‑1 (MQ=255)
                                                                ccGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAAGTATACGGTGCGCCCCCGTGAt                                                                           <  2:1082555/150‑1 (MQ=255)
                                                                  gAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTc                                                                        <  1:3624/151‑1 (MQ=255)
                                                                  gAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTc                                                                        >  1:826418/1‑151 (MQ=255)
                                                                    gCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCa                                                                       <  2:178005/150‑1 (MQ=255)
                                                                         gTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAAt                                                                 >  1:146991/1‑151 (MQ=255)
                                                                            tCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATa                                                                >  1:1011790/1‑149 (MQ=255)
                                                                             cGACTCACGAGCGACACAGACATCACGTTTTTCAATCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATACAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGCTTCGCCACCGTGATTTCAAACAACAAGt                                                             >  2:882838/1‑151 (MQ=255)
                                                                                 tCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGcc                                                         >  1:693677/1‑151 (MQ=255)
                                                                                  cACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCa                                                        <  1:429319/151‑1 (MQ=255)
                                                                                      aGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCaaa                                                      >  2:799261/1‑149 (MQ=255)
                                                                                              aGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATa                                             >  1:983818/1‑150 (MQ=255)
                                                                                                aCATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAAt                                          <  1:1127184/151‑1 (MQ=255)
                                                                                                aCATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGAGAATAAAt                                          >  2:874989/1‑151 (MQ=255)
                                                                                                  aTCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGt                                        <  2:1037918/151‑1 (MQ=255)
                                                                                                    cACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGttt                                      <  1:485883/151‑1 (MQ=255)
                                                                                                    cACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGtt                                       <  2:724450/150‑1 (MQ=255)
                                                                                                      cGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTaa                                    <  1:1003538/151‑1 (MQ=255)
                                                                                                      cGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTaa                                    <  1:24107/151‑1 (MQ=255)
                                                                                                       gTTATTCAGTACTGCTCCCTGCGAGGGATGCCGCGTGTCATATAAAACTCTCTTTAAGGCTTAGAGTCAAGCATCCAAAAGAGGTCAAATTATACGGTTCGTACCCGTGATTTCAAACAATAAGCAGCCAAAAGGCGAATAAATGTGTaa                                    >  2:711565/1‑150 (MQ=255)
                                                                                                         ttttCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAat                                  <  1:710003/150‑1 (MQ=255)
                                                                                                          tttCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAatat                                >  1:628181/1‑151 (MQ=255)
                                                                                                             cACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATAtttt                             <  1:788814/151‑1 (MQ=255)
                                                                                                              aCTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATAttttt                            <  1:944254/151‑1 (MQ=255)
                                                                                                               cTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTc                           >  1:1018640/1‑151 (MQ=255)
                                                                                                                   tGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGt                       >  1:26406/1‑151 (MQ=255)
                                                                                                                    gCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGt                       <  2:117251/150‑1 (MQ=255)
                                                                                                                    gCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTg                      <  1:438692/151‑1 (MQ=255)
                                                                                                                    gCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTg                      >  1:538685/1‑151 (MQ=255)
                                                                                                                       ccattCGAGGTATGACGCGTGTCGTATAAAACTTTCTTTTAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGAGCCCCCGTGATTTCAAACAATAAGACGCCATAAGGTGATTAAATGTTTAAATATTTTTCCAGTGCAt                   >  2:110958/4‑151 (MQ=255)
                                                                                                                                gTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCg          <  1:405043/151‑1 (MQ=255)
                                                                                                                                    gCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACtt      >  1:134249/1‑151 (MQ=255)
                                                                                                                                      cgcgTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACtttt    <  2:133793/151‑1 (MQ=255)
                                                                                                                                       gcgTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACtttt    >  2:139363/1‑150 (MQ=255)
                                                                                                                                        cgTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCt  >  2:392938/1‑151 (MQ=255)
                                                                                                                                        cgTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCt  >  2:481346/1‑151 (MQ=255)
                                                                                                                                        cgTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCt  <  1:1055193/151‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                
GGTCTTTGAGGCTGAACAGCAATCCACTTGCTATCACCACCAGTAACTGAACGAGCAGAAGCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCT  >  NC_000913/3922341‑3922627

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