Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,118,829 T→C 5.2% intergenic (+69/+16) 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,8290TC 5.2% 339.9 / 6.9 117intergenic (+69/+16)zapB/rraAFtsZ stabilizer, septal ring assembly factor, cell division stimulator/ribonuclease E (RNase E) inhibitor protein
Reads supporting (aligned to +/- strand):  ref base T (42/67);  new base C (3/3);  total (47/70)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.75e-01

Only 100 of 120 total aligned reads displayed.

AGAACAACCATCTGAAAGAACAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGCGGCA  >  NC_000913/4118681‑4118977
                                                                                                                                                    |                                                                                                                                                    
agaACAACCATCTGAAAGAACAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCGTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCtt                                                                                                                                                      >  2:926045/1‑149 (MQ=255)
 gaACAACCATCTGAAAGAACAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTg                                                                                                                                                     >  2:1011346/1‑149 (MQ=255)
   acaacCATCTGAAAGCACAGCAGAACGGCTGTCAGGACCGGCTTCAGTCCCTGCTGGGTCGCATGGAAGATGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGGTGAAAAAAGGCACCTTGCg                                                                                                                                                   <  2:274189/149‑1 (MQ=255)
       ccATCTGAAATAACAGCAGAACTGCTGGCAGGAACGTCTGCAGGCCCTGCTTGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGc                                                                                                                                               <  1:496134/149‑1 (MQ=255)
          tCTGAAAGAACAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCtt                                                                                                                                            >  2:558860/1‑149 (MQ=255)
           cTGAAAGAACAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGCGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCttt                                                                                                                                           <  1:796176/149‑1 (MQ=255)
            tGAAAGAACAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTc                                                                                                                                          <  2:210598/149‑1 (MQ=255)
             gAAAGTACTTCAGAACTTCTGGCAGGAACGTCTTAAGGCCATGATGGGTCTCTTGGAAGAGGTCTGATTACGCTGCATCTTCAGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCt                                                                                                                                         <  2:1420441/149‑1 (MQ=255)
               aaGAACAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTa                                                                                                                                       <  1:143431/149‑1 (MQ=255)
                aGAACAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCGCAAGGTGCCttttttct                                                                                                                                      <  1:670370/149‑3 (MQ=255)
                    cagcagAACGGCTGGCAGGAACGTCTTCAGGCCCTGCTGGGGCGCATGGAAGAGGTCTGATTACGCGTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCAtt                                                                                                                                  <  2:1227342/149‑1 (MQ=255)
                      gcagAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTGTTTAATGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTGTTATCATTCa                                                                                                                                >  2:166880/1‑149 (MQ=255)
                      gcagAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCa                                                                                                                                >  2:1130620/1‑149 (MQ=255)
                      gcagAACGGCGGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGTCACCTTGCGGTGCCTTTCTTATCATTCa                                                                                                                                <  1:486221/149‑1 (MQ=255)
                            cGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATcc                                                                                                                          >  2:53847/1‑149 (MQ=255)
                             ggCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCa                                                                                                                         >  2:834510/1‑149 (MQ=255)
                                   cAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGAt                                                                                                                   <  1:464000/149‑1 (MQ=255)
                                    aGGAACGTCTGCAGGCCCTGCTGGGTCACATGGAAGAGGTCTCAATACGCTTCAACTTCCGCCTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAAAAAAAAGGCACCTTGCGGTGCCTATCTAATCAATCAAAATCCAGcctagc                                                                                                                  >  2:1469111/1‑144 (MQ=255)
                                      gAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGAtctt                                                                                                                <  1:252625/149‑1 (MQ=255)
                                      gAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCAAGCGGAtctt                                                                                                                >  1:405424/1‑149 (MQ=255)
                                       aaCGTATGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGAAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGAtcttc                                                                                                               <  1:1078062/149‑1 (MQ=255)
                                        aCGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGAtcttct                                                                                                              <  2:888438/149‑1 (MQ=255)
                                        aCGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGAtcttct                                                                                                              <  1:1344984/149‑1 (MQ=255)
                                         cGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTg                                                                                                             <  2:250824/149‑1 (MQ=255)
                                          gTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGa                                                                                                            >  1:14949/1‑149 (MQ=255)
                                           tCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCGAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGAActgctgta                                                                                                           >  1:1444572/1‑147 (MQ=255)
                                           gcggCAGGCCCTGCGGGGTCGCATGAAGACGTCTGATTACGCGTCAGCATCCGCGTAATTTCTTTCCTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGaaa                                                                                                          <  2:1433701/146‑1 (MQ=255)
                                            cTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGaaa                                                                                                          <  1:1011346/149‑1 (MQ=255)
                                             tGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAg                                                                                                         <  1:250887/149‑1 (MQ=255)
                                             cgCAGCCCCTGCTGGGTCGCATGGAAGAGCTCTGCTTCCGCGTCATCTTCCGCGGAATTGCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAACGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAg                                                                                                         <  2:327261/148‑1 (MQ=255)
                                                 ggCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGaata                                                                                                     <  1:1130620/149‑1 (MQ=255)
                                                  gCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGaataa                                                                                                    <  1:72392/149‑1 (MQ=255)
                                                    ccTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATc                                                                                                  >  1:245507/1‑149 (MQ=255)
                                                     ctgctgGGTCGCATGGAAGATGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATcc                                                                                                 <  2:595504/149‑1 (MQ=255)
                                                     ctgctgGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATcc                                                                                                 >  1:616338/1‑149 (MQ=255)
                                                     ctgctgGGTAGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATACAGCGGATCTTCTGAAATAATAATcc                                                                                                 >  2:541817/1‑149 (MQ=255)
                                                        ctgGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCgg                                                                                              <  2:816768/149‑1 (MQ=255)
                                                         tgGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGt                                                                                             <  2:1155271/149‑1 (MQ=255)
                                                         tgGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGt                                                                                             <  1:1566657/149‑1 (MQ=255)
                                                             tCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTg                                                                                         >  2:925198/1‑149 (MQ=255)
                                                              cGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGt                                                                                        <  2:1509742/149‑1 (MQ=255)
                                                               gCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTc                                                                                       <  2:750016/149‑1 (MQ=255)
                                                                cATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCAATCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCg                                                                                      >  2:1178852/1‑149 (MQ=255)
                                                                cATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAAATTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCg                                                                                      >  2:664128/1‑149 (MQ=255)
                                                                 aTGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCgg                                                                                     <  1:568285/149‑1 (MQ=255)
                                                                   ggAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCa                                                                                   <  1:1360844/149‑1 (MQ=255)
                                                                    gAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCAt                                                                                  >  2:1342439/1‑149 (MQ=255)
                                                                      agagGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATaa                                                                                <  1:1028912/149‑1 (MQ=255)
                                                                        agGTCTGATTACGCTTCATCTGCCGCGTAATTTCATTACTGTGTAGGCCAGCAGGTAACGCAGAAAAAAGGCACCTAGAGGTGCCTTTCGTATCATTCAATATCCACCGGATCTTCTGAAAGAATAATCCCGGTATTCTCGGCATAAAg                                                                              >  1:342696/1‑149 (MQ=255)
                                                                         ggTCTGATTACGCCTCATCTTCCGCGTAATTTCTTAACTGTGTAGGCCAACAGGTAACGCAGAAAAAACGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGa                                                                             <  2:289102/149‑1 (MQ=255)
                                                                           tCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATg                                                                           <  2:1442200/149‑1 (MQ=255)
                                                                           tCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATg                                                                           <  2:445368/149‑1 (MQ=255)
                                                                            cTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCCGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATgg                                                                          <  1:1177626/149‑1 (MQ=255)
                                                                            cTGATTACGATTCATCTTCCGCGTAATTTCTTTACTGTGTATGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATgg                                                                          <  1:840315/149‑1 (MQ=255)
                                                                             tGATTACGCTTCATCTTCCGCGTAATTTCTATACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGt                                                                         <  1:103906/149‑1 (MQ=255)
                                                                              gATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCAAGCGGATCTTCTGAAAGATTAATCCCGGTATTGTCGGCATAAAGATGGTc                                                                        >  2:603819/1‑149 (MQ=255)
                                                                                  aCGCTTCATCTTCCGCGTAATTTCTTTACTGTGTCGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCg                                                                    <  2:865899/149‑1 (MQ=255)
                                                                                  aCGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCg                                                                    <  1:1349515/149‑1 (MQ=255)
                                                                                  aCGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCg                                                                    <  1:570121/149‑1 (MQ=255)
                                                                                     cTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCGTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGgag                                                                 >  2:1483574/1‑149 (MQ=255)
                                                                                     cTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGAAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGTTGGTCGCCGgag                                                                 <  1:400105/149‑1 (MQ=255)
                                                                                       tCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGagaa                                                               <  1:834510/149‑1 (MQ=255)
                                                                                         aTCTTAGGCGTAATTTATTTACTGTGTAGGCCAACAGGTAACGCAGAAAAACTGAACCTTGCGGTGCCTTTCTTATCATTCTATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGagaaga                                                             <  1:919515/149‑1 (MQ=255)
                                                                                           cTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGaagaag                                                           <  2:612935/149‑1 (MQ=255)
                                                                                           cTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGaagaag                                                           <  2:471047/149‑1 (MQ=255)
                                                                                             tCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGGAGAAAAAAGGCACCTTGCGGTTCCTTTCTTATCATTCAATATCCATCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGt                                                         <  2:1436202/149‑1 (MQ=255)
                                                                                             tCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGt                                                         <  2:598595/149‑1 (MQ=255)
                                                                                               cgcgTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGa                                                       <  2:1149101/149‑1 (MQ=255)
                                                                                                gcgTAATTTCTTTACTGTGTAGGCCAACAGGAAAAGGAGAAAAAAGGCACCTTGCGGTGCCTATCTTATCATTCAAAATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGac                                                      >  2:76284/1‑149 (MQ=255)
                                                                                                   tAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGGCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATAACGGTATTGAAGTCATAAAGATGGTCGCCGGAGAAGAAGGTGACAcc                                                   >  1:1178272/1‑149 (MQ=255)
                                                                                                   tAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACAcc                                                   >  1:1195300/1‑149 (MQ=255)
                                                                                                         cTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCaaa                                             >  2:407043/1‑149 (MQ=255)
                                                                                                         cTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCaaa                                             >  1:470400/1‑149 (MQ=255)
                                                                                                          tttACTGTTTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGCATCTTCTGAAAGAATAAACCCGTTATTGTCGGCATAAAGATGGGCGCCGGAGAAGAAGGTGACACCGCCaaaa                                            <  1:1469111/149‑1 (MQ=255)
                                                                                                          tttACTGGGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCATTATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGGGACACCGCCaaaa                                            <  1:1344658/149‑1 (MQ=255)
                                                                                                             agtgtgTAGGCCTACAGGTATCTAATAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGCATCTTCTGAAAGAATAATCCCTGAATTGGCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTg                                         <  2:1298163/147‑1 (MQ=255)
                                                                                                             aCTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTg                                         <  1:1483574/149‑1 (MQ=255)
                                                                                                              cTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGa                                        <  1:57994/149‑1 (MQ=255)
                                                                                                              cTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGa                                        <  1:926045/149‑1 (MQ=255)
                                                                                                               tgtgTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTTACTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAATAATAATCCCGGTATTGTCGGCATAATGATGGTCGACGCATAAGAATGTTACACAGCCAAAATTGAc                                       >  2:1364703/1‑149 (MQ=255)
                                                                                                                  gtAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGAcgcg                                    <  1:53847/149‑1 (MQ=255)
                                                                                                                      gCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCATCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAATATGGTCGCCGGAGAAGAACGTGACACCGCCAAAATTGACGCGGACa                                >  2:931350/1‑149 (MQ=255)
                                                                                                                        cAACAGGTAACGCAGAAACAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATCTCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATc                              <  2:406551/149‑1 (MQ=255)
                                                                                                                        cAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCAATCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATc                              <  2:412285/149‑1 (MQ=255)
                                                                                                                         aaCAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCg                             <  2:320868/149‑1 (MQ=255)
                                                                                                                         aaCAGGTAACGCAGAAAAAAGGCACCTAGCGGTGCCTTTCATATCATTCAATATCCAGCGGATCTTCTGAAAGCATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGCAGCAGGTGACCCCGCCAAAATTGACGCGGACATCg                             >  2:491591/1‑149 (MQ=255)
                                                                                                                         aaCAGGTAACGCAGAAAAAACGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTAATGTCGGAATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCAGACATCg                             >  1:504398/1‑149 (MQ=255)
                                                                                                                              gTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTc                        >  2:369909/1‑149 (MQ=255)
                                                                                                                              gTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGGGACACCGCCAAAATTGACGCGGACATCGCTTTc                        <  1:925198/149‑1 (MQ=255)
                                                                                                                               tAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCg                       <  1:1114580/149‑1 (MQ=255)
                                                                                                                                   gCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGCCACCGCCAAAATTGACGCGGACATCGCTTTCGCCaa                   >  1:778388/1‑149 (MQ=255)
                                                                                                                                   gCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCaa                   >  2:1002160/1‑149 (MQ=255)
                                                                                                                                    cAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAAt                  >  1:218897/1‑149 (MQ=255)
                                                                                                                                     aGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATg                 >  1:1238230/1‑149 (MQ=255)
                                                                                                                                       aagaaaGGCACCGCAAGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGcc               >  1:528100/4‑149 (MQ=255)
                                                                                                                                       aagaaaGGCACCGCAAGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGcc               >  2:1277383/4‑149 (MQ=255)
                                                                                                                                       aaaaaaGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGcc               >  1:1363322/1‑149 (MQ=255)
                                                                                                                                            aGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGc          <  2:1383820/149‑1 (MQ=255)
                                                                                                                                                 ccTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCTGGATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTCACACCGCCAAAATTGACGCGGACATCGCGATCGCCAAAGCCTGCGCCAGCg     >  2:880584/1‑149 (MQ=255)
                                                                                                                                                 ccTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGCg     <  2:1529845/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AGAACAACCATCTGAAAGAACAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGCGGCA  >  NC_000913/4118681‑4118977

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: