Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,118,830 G→A 8.0% intergenic (+70/+15) zapB → / ← rraA FtsZ stabilizer, septal ring assembly factor, cell division stimulator/ribonuclease E (RNase E) inhibitor protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,118,8300GA 8.0% 400.0 / 16.4 138intergenic (+70/+15)zapB/rraAFtsZ stabilizer, septal ring assembly factor, cell division stimulator/ribonuclease E (RNase E) inhibitor protein
Reads supporting (aligned to +/- strand):  ref base G (63/63);  new base A (4/7);  total (67/71)
Fisher's exact test for biased strand distribution p-value = 5.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.92e-02

Only 100 of 138 total aligned reads displayed.

GAACAACCATCTGAAAGAACAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGCGGC  >  NC_000913/4118682‑4118976
                                                                                                                                                    |                                                                                                                                                  
gaACAACCATCTGAAAGAACAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGc                                                                                                                                                   <  2:1101342/150‑1 (MQ=255)
gaACAACCATCTGAAAGAACAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCGTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAATGCACCTTGc                                                                                                                                                   >  2:1240851/1‑150 (MQ=255)
  acaacCATCTGAAAGAACAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCgg                                                                                                                                                 >  2:903324/1‑150 (MQ=255)
   caacCATCTGAAAGAACAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGt                                                                                                                                                >  2:1113159/1‑150 (MQ=255)
      ccctCTGACAGAACAGCTGAACGGCTGGCAGGATCGTCTGCAGGCCCTGCTGGGTCGCATGGAAGCGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGcc                                                                                                                                             <  1:1609699/147‑1 (MQ=255)
         tCTGAAAGACCAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCttt                                                                                                                                          <  1:745828/150‑1 (MQ=255)
         tCTGAAAGAACAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCttt                                                                                                                                          <  1:1464880/150‑1 (MQ=255)
         tCTGAAAGAACAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCGCAAGGTGCCttt                                                                                                                                          <  1:1156801/150‑1 (MQ=255)
           tGAACGAACAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCGGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCGCAAGGTGCCttttt                                                                                                                                        <  1:1119679/150‑3 (MQ=255)
           tGAAACAACAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCAGGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCt                                                                                                                                        <  2:286432/150‑1 (MQ=255)
            gAAAGAACAGCAGATCGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCATCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCtt                                                                                                                                       <  2:126527/150‑1 (MQ=255)
            gAAAGAACAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCGCAAGGTGCCtttttt                                                                                                                                       <  2:210909/150‑1 (MQ=255)
            gAAAGAACAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCGCAAGGTGCCtttttt                                                                                                                                       <  2:1311967/150‑1 (MQ=255)
            gAAAGAACAGCAGAACGGCTGGCAGGAACGGCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCtt                                                                                                                                       <  2:1142673/150‑1 (MQ=255)
              aaGAACAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTAt                                                                                                                                     <  2:1597578/150‑1 (MQ=255)
              aaGAACAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTAt                                                                                                                                     >  2:1595145/1‑150 (MQ=255)
              aaGAACAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCATCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCGCAAGGTGCCttttttct                                                                                                                                     <  1:1269319/150‑3 (MQ=255)
               aGAACAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCTTGGAAGATGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATc                                                                                                                                    <  1:1481503/150‑1 (MQ=255)
               aGAACAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATc                                                                                                                                    <  1:666875/150‑1 (MQ=255)
                 aaCAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCAt                                                                                                                                  <  2:1560898/150‑1 (MQ=255)
                  aCAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCAtt                                                                                                                                 >  1:671235/1‑150 (MQ=255)
                   cagcagACCGGCTGGCAGGAACGGCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTc                                                                                                                                <  1:1318952/150‑1 (MQ=255)
                     gcagAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAATGCACCGTGCGGTGCCGTTAGTATCATTCaa                                                                                                                              >  1:1563459/1‑150 (MQ=255)
                      cagAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCCTTCaaa                                                                                                                             >  2:1422326/1‑149 (MQ=255)
                        gAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTa                                                                                                                                      >  2:1123033/1‑139 (MQ=255)
                           cGGCTGGCAGGAACGTCTGCAGGCCCTGCGGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTGTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATAATTCAATATCCa                                                                                                                        >  1:951692/1‑150 (MQ=255)
                                ggCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGa                                                                                                                   >  2:382535/1‑150 (MQ=255)
                                    ggAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGAtctt                                                                                                               >  1:1117891/1‑150 (MQ=255)
                                       aCGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTg                                                                                                            <  1:1335020/150‑1 (MQ=255)
                                        cGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGa                                                                                                           <  1:1113159/150‑1 (MQ=255)
                                         gTCTGCAGTGCCTGCTGGGTCGCATGGAACAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCCGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGaa                                                                                                          <  2:308664/150‑1 (MQ=255)
                                          tCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGaaa                                                                                                         <  2:144069/150‑1 (MQ=255)
                                           cTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAg                                                                                                        >  1:137026/1‑150 (MQ=255)
                                           cTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAg                                                                                                        <  1:1713014/150‑1 (MQ=255)
                                           cTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTCCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAg                                                                                                        <  2:676027/150‑1 (MQ=255)
                                               aGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTACCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGaata                                                                                                    >  1:524146/1‑150 (MQ=255)
                                                ggCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACCGGGAACGCAGCAAAAAGGCCCCTGGCGGGGCCGGGCTGATCCTTCAATATCCAGCGGAGCTTCTGACAgaagaa                                                                                                   >  1:1674972/1‑150 (MQ=255)
                                                   ccTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATcc                                                                                                >  2:978890/1‑150 (MQ=255)
                                                   ccTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATcc                                                                                                >  1:1456556/1‑150 (MQ=255)
                                                   ccTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATcc                                                                                                <  2:235341/150‑1 (MQ=255)
                                                    ctgctgGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATccc                                                                                               <  2:223133/150‑1 (MQ=255)
                                                     tgctgGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCg                                                                                              <  2:1013942/150‑1 (MQ=255)
                                                      gctgGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCgg                                                                                             <  2:605882/150‑1 (MQ=255)
                                                      gctgGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGAGCTTCTGAAAGAATAATCCCgg                                                                                             <  1:1127906/150‑1 (MQ=255)
                                                           gTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTg                                                                                        >  2:1897/1‑150 (MQ=255)
                                                           gTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTg                                                                                        >  2:507022/1‑150 (MQ=255)
                                                               cATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCgg                                                                                    >  1:952170/1‑150 (MQ=255)
                                                                 tGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCa                                                                                  <  1:1338296/150‑1 (MQ=255)
                                                                        ggTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCAtccaggc                                                                           >  1:1385548/1‑144 (MQ=255)
                                                                        ggTCTGATTACGCTTCAACTTACGCGAAATTTATTTACTGTGTATGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCATATCATTCAATATCAAGCGCATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAACAt                                                                           >  2:171309/1‑150 (MQ=255)
                                                                         gTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGTTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATg                                                                          >  2:1225580/1‑150 (MQ=255)
                                                                           cTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGt                                                                        <  1:517679/150‑1 (MQ=255)
                                                                            tGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTc                                                                       <  2:1058395/150‑1 (MQ=255)
                                                                               ttACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGcc                                                                    <  1:542949/150‑1 (MQ=255)
                                                                                 aCGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCgg                                                                  <  2:1252861/150‑1 (MQ=255)
                                                                                  cGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGga                                                                 >  1:467158/1‑150 (MQ=255)
                                                                                  cGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGga                                                                 >  2:1213857/1‑150 (MQ=255)
                                                                                    cTTCATCTTCCGTGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGgaga                                                               <  1:223222/150‑1 (MQ=255)
                                                                                     ttCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGAAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATAATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGagaa                                                              <  2:100262/150‑1 (MQ=255)
                                                                                        aTCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTAACATTCAATATCCAGCGGATCTTCTGAAAGAATAATACCGGTATTTTCGGCATAAAGATGGTAGCCGGAgaagaa                                                           >  2:1318923/1‑150 (MQ=255)
                                                                                          cTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATTGTCGCCGGAGAAGAAgg                                                         <  1:1058991/150‑1 (MQ=255)
                                                                                          cTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAgg                                                         <  1:1559462/150‑1 (MQ=255)
                                                                                          cTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAgg                                                         <  1:1897/150‑1 (MQ=255)
                                                                                          cTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGACTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAgg                                                         <  2:833028/150‑1 (MQ=255)
                                                                                            tCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCATAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTg                                                       <  2:306660/150‑1 (MQ=255)
                                                                                            tCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAATGGCAACTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTg                                                       <  2:1463313/150‑1 (MQ=255)
                                                                                             ccGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGa                                                      >  2:458824/1‑150 (MQ=255)
                                                                                               gcgTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGaca                                                    >  2:1363042/1‑150 (MQ=255)
                                                                                                 gTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACAcc                                                  <  1:831737/150‑1 (MQ=255)
                                                                                                 gTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACAcc                                                  >  1:408343/1‑150 (MQ=255)
                                                                                                    aTTTCTTTAGTGTGTAGGCCAACAGGGAACGCAGAAAAACGGCACCTTGCTGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCAGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCTGAGAAGAAGGTGACACCGcc                                               <  1:128360/150‑1 (MQ=255)
                                                                                                    aTTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATAGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGcc                                               >  1:1718002/1‑150 (MQ=255)
                                                                                                       tCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAATGCAACTTGCGGTGCCTTTCTTATCATTAAATATACAGCGGATCTTCTGAAATTATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAAAAGTTGACACCTCCaaa                                            >  2:348640/1‑150 (MQ=255)
                                                                                                         tttACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTAGTATCATTAAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACCCCGCCAAAAt                                          >  1:1304720/1‑150 (MQ=255)
                                                                                                          ttACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTAATGTCGCCATAAAGATGGTCGCCGGACAAGAAGGTGACACCGCCaaaaat                                         >  2:747773/1‑148 (MQ=255)
                                                                                                          ttACTGTGTAGGCAAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCACGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAAtt                                         <  2:90439/150‑1 (MQ=255)
                                                                                                            aCTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGa                                       <  2:864625/150‑1 (MQ=255)
                                                                                                                 gtAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCgg                                  >  1:94544/1‑150 (MQ=255)
                                                                                                                 gtAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATATTCTGAAAGAATAATCCCGGTATAGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAAATGAAGCgg                                  >  1:73021/1‑150 (MQ=255)
                                                                                                                  tAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGa                                 >  1:554362/1‑150 (MQ=255)
                                                                                                                   aGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAGGGTGACACCGCCAAAATTGACGCGGAc                                <  2:551891/150‑1 (MQ=255)
                                                                                                                        aaCAGTTAACGCAGAAAAAAGGCAACTTTCTGTGCCTTTCTTATCATTCAATATCCAGCGGATATTCTGAAAGAATAATCCCGGTATTTTCGGCATAAAGATGGTCTCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGc                           <  1:1083577/150‑1 (MQ=255)
                                                                                                                        aaCAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGc                           <  1:336927/150‑1 (MQ=255)
                                                                                                                        aaCAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGc                           <  1:789058/150‑1 (MQ=255)
                                                                                                                          cAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCtt                         >  1:821345/1‑150 (MQ=255)
                                                                                                                            ggTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTc                       <  1:664024/150‑1 (MQ=255)
                                                                                                                             gTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACAACGCTTTCg                      <  2:1020992/150‑1 (MQ=255)
                                                                                                                                 cGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAAAATACAGCGGATCTTATGAAAGAATAATCCCGGTATTTTCGGCATAAAGATGGTCGCCGGATAAGAAGGTGACCCCGCCAAAATTTACGCGGACATCGCTTTCGCCaa                  >  2:1451218/1‑150 (MQ=255)
                                                                                                                                   cAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATg                <  1:619963/150‑1 (MQ=255)
                                                                                                                                    ataagaaaGGCACCGCAAGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGc               >  1:85870/6‑150 (MQ=255)
                                                                                                                                     gAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCAGACATCGCTTTCGCCAATGcc              >  1:492577/1‑150 (MQ=255)
                                                                                                                                      aaaaaaGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCt             >  2:1145615/1‑150 (MQ=255)
                                                                                                                                       agaaaGGCACCACAAGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCtt            >  2:645215/3‑150 (MQ=255)
                                                                                                                                       aaaaaGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCtt            >  1:1527022/1‑150 (MQ=255)
                                                                                                                                        aaaaGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTc           <  1:908446/150‑1 (MQ=255)
                                                                                                                                        aaaaGGCACCTTGCGGTGCCTTTCTTATCATTCAATATACAGCGGTTCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAAGGCCTTc           >  2:1059089/1‑150 (MQ=255)
                                                                                                                                            ggCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCa       >  2:1427927/1‑150 (MQ=255)
                                                                                                                                                ccTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGCgg   <  1:1020120/150‑1 (MQ=255)
                                                                                                                                                ccTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGCgg   <  2:952170/150‑1 (MQ=255)
                                                                                                                                                 cTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGCGGc  <  2:1040238/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GAACAACCATCTGAAAGAACAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGCGGC  >  NC_000913/4118682‑4118976

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