Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,118,831 C→A 7.1% intergenic (+71/+14) 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,8310CA 7.1% 461.0 / 17.9 159intergenic (+71/+14)zapB/rraAFtsZ stabilizer, septal ring assembly factor, cell division stimulator/ribonuclease E (RNase E) inhibitor protein
Reads supporting (aligned to +/- strand):  ref base C (64/83);  new base A (4/8);  total (68/91)
Fisher's exact test for biased strand distribution p-value = 5.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.16e-01

Only 100 of 162 total aligned reads displayed.

GAACAACCATCTGAAAGAACAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGCGGCACCA  >  NC_000913/4118682‑4118980
                                                                                                                                                     |                                                                                                                                                     
gaACAACCATCTGAAAGAACATCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGTAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCATAAAAAAGGCACCTTGc                                                                                                                                                       <  1:394743/150‑1 (MQ=255)
  acaacCATCTGAAAGAACAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCgg                                                                                                                                                     >  1:693226/1‑150 (MQ=255)
   caacCATCTGAAAGAACAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGt                                                                                                                                                    >  1:69042/1‑150 (MQ=255)
     acCATCTGAAAGAACAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCTTGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGc                                                                                                                                                  <  1:350269/150‑1 (MQ=255)
      ccATCTGAAAGAACAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTTACGCAGAAAAAAGGCACCTTGCGGTGcc                                                                                                                                                 <  2:159070/150‑1 (MQ=255)
      ccATCTGAAAGAACAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGcc                                                                                                                                                 <  2:1037175/150‑1 (MQ=255)
      ccATCTGAAAGAACAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGcc                                                                                                                                                 <  1:1084464/150‑1 (MQ=255)
         tCTGAAAGAACAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCGCAAGGTGCCttt                                                                                                                                              <  2:1126549/150‑1 (MQ=255)
           tGAAAGAACAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCt                                                                                                                                            >  2:432598/1‑150 (MQ=255)
           tGAAAGAACAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCCTAATTTCATTACTCTGAAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGccatact                                                                                                                                            >  2:292860/1‑145 (MQ=255)
           tGAAAGAACAGCAGAACGGCTGGCACGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCttttt                                                                                                                                            <  1:1100570/150‑3 (MQ=255)
            gAAAGAACAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCtt                                                                                                                                           >  1:1698466/1‑150 (MQ=255)
            gAAAGAACAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCtt                                                                                                                                           <  2:1192330/150‑1 (MQ=255)
            gAAAGAACAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCtt                                                                                                                                           <  2:888811/150‑1 (MQ=255)
            gAAAGAACAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCtt                                                                                                                                           <  1:530497/150‑1 (MQ=255)
            gAAAGAACAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCtt                                                                                                                                           >  1:47123/1‑150 (MQ=255)
             aaaGAACAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTa                                                                                                                                          <  1:1125301/150‑1 (MQ=255)
               aGAACAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCGCAAGGTGCCttttttctg                                                                                                                                        <  2:20795/150‑4 (MQ=255)
               aGAACAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACACTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATc                                                                                                                                        <  2:467220/150‑1 (MQ=255)
                gAACAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCa                                                                                                                                       >  2:839157/1‑150 (MQ=255)
                gAACAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCAGGGAAGAGTTCTGATTACGCGTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCACCCGGTAACGCAGAAAAAAGGCACCGCAAGGTGCCttttttctgc                                                                                                                                       <  2:108054/150‑5 (MQ=255)
                   cagcagAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTc                                                                                                                                    >  2:551537/1‑150 (MQ=255)
                     gcagAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCaa                                                                                                                                  >  2:51678/1‑150 (MQ=255)
                      cagAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAat                                                                                                                                 <  2:1253193/150‑1 (MQ=255)
                       agAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAata                                                                                                                                <  2:311091/150‑1 (MQ=255)
                       agAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAata                                                                                                                                >  2:727577/1‑150 (MQ=255)
                        gAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAatat                                                                                                                               >  2:1668753/1‑150 (MQ=255)
                          aCGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATcc                                                                                                                             <  2:923801/150‑1 (MQ=255)
                          aCGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATcc                                                                                                                             <  2:1108039/150‑1 (MQ=255)
                            ggatGGGAGGAACGTCTGCAGGCCCTGCTTGGTCGCATGGTAGAGGTATGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAg                                                                                                                           <  1:158902/147‑1 (MQ=255)
                              cTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCg                                                                                                                         <  1:161933/150‑1 (MQ=255)
                              cTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCg                                                                                                                         <  2:450804/150‑1 (MQ=255)
                               tGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGGTTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCgg                                                                                                                        <  1:1600984/150‑1 (MQ=255)
                               tGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCgg                                                                                                                        <  1:1040135/150‑1 (MQ=255)
                                ggCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGGGCCTTTCGTATCATTAAATATCCAGCGGa                                                                                                                       >  2:781879/1‑150 (MQ=255)
                                  cAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGAtc                                                                                                                     <  2:876482/150‑1 (MQ=255)
                                    ggAACGTCTGCAGGCCCTGCTGGTTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCGAAAGGTGCCttttttc                                                                                                                                          <  1:1511063/127‑2 (MQ=255)
                                    ggAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGTTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCGCAAGGTGCCttttttc                                                                                                                                          >  2:1511066/1‑126 (MQ=255)
                                    ggAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTATGCCAACAGGTAACGCAGAAAAAAGGCACCGCAAGGTGCCttttttc                                                                                                                                          >  2:1511063/1‑126 (MQ=255)
                                    ggAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGAtctt                                                                                                                   >  2:369808/1‑150 (MQ=255)
                                              cAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGaat                                                                                                         >  1:118691/1‑150 (MQ=255)
                                               aGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGGGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGaata                                                                                                        <  1:587823/150‑1 (MQ=255)
                                                gccccTGCTGGCTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTCCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGaataa                                                                                                       >  1:275110/3‑150 (MQ=255)
                                                 gCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGaataat                                                                                                      >  1:1239155/1‑150 (MQ=255)
                                                 gCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGaataat                                                                                                      >  1:787146/1‑150 (MQ=255)
                                                      gctgGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCgg                                                                                                 >  1:413166/1‑150 (MQ=255)
                                                      gcggggTCGCATGGAAGATGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCgg                                                                                                 <  1:871472/147‑1 (MQ=255)
                                                          ggTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTAtt                                                                                             >  1:1293879/1‑150 (MQ=255)
                                                             cGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTc                                                                                          <  2:1698466/150‑1 (MQ=255)
                                                             cGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTc                                                                                          <  2:1072070/150‑1 (MQ=255)
                                                             cGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTc                                                                                          >  2:1187900/1‑150 (MQ=255)
                                                               cATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCgg                                                                                        <  1:1027899/150‑1 (MQ=255)
                                                               cATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCgg                                                                                        <  1:468788/150‑1 (MQ=255)
                                                                 tGGTAGAGGTATGATTACGCTTCATCTTCCGCTTTATTTCTTTACTTTGTAGGCCAACAGTTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTAGCATTCAATATCCAGCGTATCTTCTGAAAGAATAATCCCGGTATTGTCGGCa                                                                                      <  1:292860/150‑1 (MQ=255)
                                                                 tGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTACCGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCa                                                                                      >  1:1359752/1‑150 (MQ=255)
                                                                 tGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCa                                                                                      <  1:432598/150‑1 (MQ=255)
                                                                   gAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATa                                                                                    >  1:724738/1‑150 (MQ=255)
                                                                   gAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATa                                                                                    >  2:1665976/1‑150 (MQ=255)
                                                                     agagGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATaaa                                                                                  <  2:705543/150‑1 (MQ=255)
                                                                          tCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATgg                                                                             <  2:1638046/150‑1 (MQ=255)
                                                                             gaatACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCg                                                                          >  2:229047/4‑150 (MQ=255)
                                                                                 aCGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCgg                                                                      <  1:1419003/150‑1 (MQ=255)
                                                                                 aCGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCgg                                                                      <  1:974722/150‑1 (MQ=255)
                                                                                   gCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCATAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATTGTCGCCGgag                                                                    <  2:1142004/150‑1 (MQ=255)
                                                                                          cTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCAAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAACATGGTCGCCGGAGAAGAAgg                                                             >  2:1402299/1‑150 (MQ=255)
                                                                                           ttCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGt                                                            <  2:1159952/150‑1 (MQ=255)
                                                                                            tCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTg                                                           >  2:1558227/1‑150 (MQ=255)
                                                                                             ccGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTAAATATCAAGAGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGa                                                          >  1:972245/1‑150 (MQ=255)
                                                                                                cgTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGacac                                                       >  1:1132990/1‑150 (MQ=255)
                                                                                                cgTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGacac                                                       >  1:994237/1‑150 (MQ=255)
                                                                                                   aaTTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGc                                                    >  1:619877/1‑150 (MQ=255)
                                                                                                        cTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCaaaa                                               <  1:347308/150‑1 (MQ=255)
                                                                                                        cTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGACATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCaaaa                                               >  1:275483/1‑150 (MQ=255)
                                                                                                         tttACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAAt                                              >  1:1625461/1‑150 (MQ=255)
                                                                                                           tACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTg                                            >  2:1353952/1‑150 (MQ=255)
                                                                                                           tACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTg                                            >  2:127502/1‑150 (MQ=255)
                                                                                                                      ccAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATc                                 <  2:1625461/150‑1 (MQ=255)
                                                                                                                       cAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCg                                <  2:1239155/150‑1 (MQ=255)
                                                                                                                       cAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCg                                <  1:78262/150‑1 (MQ=255)
                                                                                                                        aaCAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGc                               >  1:1580553/1‑150 (MQ=255)
                                                                                                                        aaCAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGc                               >  1:1586967/1‑150 (MQ=255)
                                                                                                                          cAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCtt                             <  1:1076141/150‑1 (MQ=255)
                                                                                                                            ggTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTc                           <  1:229047/150‑1 (MQ=255)
                                                                                                                            ggTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTc                           >  2:578762/1‑150 (MQ=255)
                                                                                                                            ggTAACGCAGAAAAAAGACACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTc                           <  1:1441801/150‑1 (MQ=255)
                                                                                                                             gTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCg                          <  1:154547/150‑1 (MQ=255)
                                                                                                                             gTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCg                          >  2:1290807/1‑150 (MQ=255)
                                                                                                                              tAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGc                         <  1:1558227/150‑1 (MQ=255)
                                                                                                                               aaCGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGcc                        <  1:166767/150‑1 (MQ=255)
                                                                                                                               aaCGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGcc                        <  1:60559/150‑1 (MQ=255)
                                                                                                                                 cGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCaa                      >  2:1589217/1‑150 (MQ=255)
                                                                                                                                   cAGAAAAAGGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCTCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATg                    <  1:209314/150‑1 (MQ=255)
                                                                                                                                     gAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGcc                  >  1:705711/1‑150 (MQ=255)
                                                                                                                                         aaaGGCACCGCAAGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCg              >  1:1652476/1‑150 (MQ=255)
                                                                                                                                           aGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGcc            >  1:765539/1‑150 (MQ=255)
                                                                                                                                             gCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAg          >  2:253470/1‑150 (MQ=255)
                                                                                                                                                ccTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGCgg       >  1:1272344/1‑150 (MQ=255)
                                                                                                                                                    gCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGACGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGCGGCAcc   <  2:47123/150‑1 (MQ=255)
                                                                                                                                                     cGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGCGGCACCa  <  2:724738/150‑1 (MQ=255)
                                                                                                                                                     cGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGCGGCACCa  <  2:787146/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GAACAACCATCTGAAAGAACAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGCGGCACCA  >  NC_000913/4118682‑4118980

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: