Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 801,357 G→A 6.2% intergenic (+30/‑159) rpmA → / → obg 50S ribosomal protein L27/GTPase Obg

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947801,3570GA 6.2% 239.0 / 2.7 81intergenic (+30/‑159)rpmA/obg50S ribosomal protein L27/GTPase Obg
Reads supporting (aligned to +/- strand):  ref base G (35/39);  new base A (3/2);  total (38/43)
Fisher's exact test for biased strand distribution p-value = 6.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.38e-01

CTGGCTACGGCGTTGGCATGGGCAAGGACCACACCTTGTTCGCCAAGATCGAAGGCGTGATCAAGTTCGAGAAGAAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGCTCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTGTTTGTCTGGGTTGCCGGCGTTGGTTTCTACGGTCGTACGGGGCCGCCGCGCTCATTTTTGCAAGAGCCTCATGTTTTTGAGTCATTCAGCTCGTCTGCCGAC  >  NC_002947/801209‑801501
                                                                                                                                                    |                                                                                                                                                
cTGGCTACGGCGTTGGCATGGGCAAGGACCACACCTTGTTCGCCAAGATCGAAGGCGTGATCAAGTTCGAGAAGAAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGCTCCAGAAGCCCCGTCATg                                                                                                                                                  <  8:312280/149‑1 (MQ=255)
cTGGCTACGGCGTTGGCATGGGCAAGGACCACACCTTGTTCGCCAAGATCGAAGGCGTGATCAAGTTCGAGAAGAAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGCTCCAGAAGCCCCGTCATg                                                                                                                                                  <  7:2715/149‑1 (MQ=255)
  ggCTACGGCGTTGGCATGGGCAAGGACCACACCTTGTTCGCCAAGATCGAAGGCGTGATCAAGTTCGAGAAGAAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGCTCCAGAAGCCCCGTCATGc                                                                                                                                                 >  4:40833/1‑148 (MQ=255)
   gCTACGGCGTTGGCATGGGCAAGGACCACACCTTGTTCGCCAAGATCGAAGGCGTGATCAAGTTCGAGAAGAAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGCTCCAGAAGCCCCGTCATGc                                                                                                                                                 <  8:26100/147‑1 (MQ=255)
     tACGGCGTTGGCATGGGCAAGGACCACACCTTGTTCGCCAAGATCGAAGGCGTGATCAAGTTCGAGAAGAAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGCTCCAGAAGCCCCGTCATGCGACg                                                                                                                                             >  8:240378/1‑149 (MQ=255)
     tACGGCGTTGGCATGGGCAAGGACCACACCTTGTTCGCCAAGATCGAAGGCGTGATCAAGTTCGAGAAGAAAGGCGAGTTCATGCGCCGTTACGTGAGAATCGTCGCCGCTTAATCGCGACGTCGCTCCAGAAGCCCCGTCATGCGACg                                                                                                                                             >  3:72692/1‑149 (MQ=255)
       cGGCGTTGGCATGGGCAAGGACCACACCTTGTTCGCCAAGATCGAAGGCGTGATCAAGTTCGAGAAGAAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGCTCCAGAAGCCCCGTCATGc                                                                                                                                                 <  8:204907/143‑1 (MQ=255)
       cGGCGTTGGCATGGGCAAGGACCACACCTTGTTCGCCAAGATCGAAGGCGTGATCAAGTTCGAGAAGAAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGCTCCAGAAGCCCCGTCATGc                                                                                                                                                 >  7:204907/1‑143 (MQ=255)
       cGGCGTTGGCATGGGCAAGGACCACACCTTGTTCGCCAAGATCGAAGGCGTGATCAAGTTCGAGAAGAAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGCTCCAGAAGCCCCGTCATGCGACg                                                                                                                                             <  5:174962/147‑1 (MQ=255)
           gTTGGCATGGGCAAGGACCACACCTTGTTCGCCAAGATCGAAGGCGTGATCAAGTTCGAGAAGAAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGCTCCAGAAGCCCCGTCATGCGACGGGGCtt                                                                                                                                       >  5:348963/1‑149 (MQ=255)
              ggCATGGGCAAGGACCACACCTTGTTCGCCAAGATCGAAGGCGTGATCAAGTTCGAGAAGAAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGCTCCAGAAGCCCCGTCATGCGACGGGGCtttt                                                                                                                                     <  6:80735/148‑1 (MQ=255)
               gCATGGGCAAGGACCACACCTTGTTCGCCAAGATCGAAGGCGTGATCAAGTTCGAGAAGAAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGCTCCAGAAGCCCCGTCATGCGACGGGGCTTTTTc                                                                                                                                   >  1:149095/1‑149 (MQ=255)
                cATGGGCAAGGACCACACCTTGTTCGCCAAGATCGAAGGCGTGATCAAGTTCGAGAAGAAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGCTCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCg                                                                                                                                  >  3:157408/1‑149 (MQ=255)
                cATGGGCAAGGACCACACCTTGTTCGCCAAGATCGAAGGCGTGATCAAGTTCGAGAAGAAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGCTCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCg                                                                                                                                  >  4:32606/1‑149 (MQ=255)
                cATGGGCAAGGACCACACCTTGTTCGCCAAGATCGAAGGCGTGATCAAGTTCGAGAAGAAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGCTCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCg                                                                                                                                  >  1:74209/1‑149 (MQ=255)
                cATGGGCAAGGACCACACCTTGTTCGCCAAGATCGAAGGCGTGATCAAGTTCGAGAAGAAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGCTCCAAAAGCCCCGTCGCATGACGGGGCttctgga                                                                                                                                  <  6:39355/149‑6 (MQ=255)
                cATGGGCAAGGACCACACCTTGTTCGCCAAGATCGAAGGCGTGATCAAGTTCGAGAAGAAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGCTCCAAAAGCCCCGTCGCATGACGGGGCttctg                                                                                                                                    >  5:39355/1‑144 (MQ=255)
                 aTGGGCAAGGACCACACCTTGTTCGCCAAGATCGAAGGCGTGATCAAGTTCGAGAAGAAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGCTCCAGAAGCCCCGTCGCATGACGGGGCttctggag                                                                                                                                 <  8:78613/149‑7 (MQ=255)
                  tGGGCAAGGACCACACCTTGTTCGCCAAGATCGAAGGCGTGATCAAGTTCGAGAAGAAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGCTCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGtt                                                                                                                                <  4:157408/149‑1 (MQ=255)
                    ggCAAGGACGACACCCTGTTCGGCAAGATCGAAGGCGTGATCAAGTTCGAGAAGAAAGGCGAGATCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGCTCCAGAAGCCGCGTCATGCGACGGGGCTTTTTCGttt                                                                                                                               <  6:86659/148‑1 (MQ=255)
                     gCAAGGACCACACCTTGTTCGCCAAGATCGAAGGCGTGATCAAGTTCGAGAAGAAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGCTCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTtgt                                                                                                                             >  5:59757/1‑149 (MQ=255)
                     gCAAGGACCACACCTTGTTCGCCAAGATCGAAGGAGTGATCAAGTTCGAGAACAAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATAGCGACGTCGCTCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTtgt                                                                                                                             >  5:309743/1‑149 (MQ=255)
                         ggACCACACCTTGTTCGCCAAGATCGAAGGCGTGATCAAGTTCGAGAAGAAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGCTCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTgtggtg                                                                                                                         <  3:119804/149‑1 (MQ=255)
                          gACCACACCTTGTTCGCCAAGATCGAAGGCGTGATCAAGTTCGAGAAGAAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGCTCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGa                                                                                                                        <  3:129884/149‑1 (MQ=255)
                           aCCACACCTTGTTCGCCAAGATCGAAGGCGTGATCAAGTTCGAGAAGAAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGCTCCAGAAGCCCCGTCATCCGACGGGGCTTTTTCGTTTGTGGTGAg                                                                                                                       <  1:211864/149‑1 (MQ=255)
                           aCCACACCTTGTTCGCCAAGATCGAAGGCGTGATCAAGTTCGAGAAGAAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGCTCCAGAAGCCCCGTCATCCGACGGGGCTTTTTCGTTTGTGGTGAg                                                                                                                       <  1:218104/149‑1 (MQ=255)
                             cacaCCTTGTTCGCCAAGATCGAAGGCGTGATCAAGTTCGAGAAGAAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGCTCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGcc                                                                                                                     >  2:56751/1‑149 (MQ=255)
                                 ccTTGTTCGCCAAGATCGAAGGCGTGATCAAGTTCGAGAAGAAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGCTCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCtt                                                                                                                 <  7:240378/149‑1 (MQ=255)
                                  cTTGTTCGCCAAGATCGAAGGCGTGATCAAGTTCGAGAAGAAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGCTCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCtt                                                                                                                 >  7:39909/1‑148 (MQ=255)
                                  cTTGTTCGCCAAGATCGAAGGCGTGATCAAGTTCGAGAAGAAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGCTCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCtt                                                                                                                 <  8:39909/148‑1 (MQ=255)
                                  cTTGTTCGCCAAGATCGAAGGCGTGATCAAGTTCGAGAAGAAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGCTCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTg                                                                                                                <  7:103492/149‑1 (MQ=255)
                                         gCCAAGATCGAAGGCGTGATCAAGTTCGAGAAGAAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGCTCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCt                                                                                                         <  3:32606/149‑1 (MQ=255)
                                         gCCAAGATCGAAGGCGTGATCAAGTTCGAGAAGAAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGCTCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCt                                                                                                         <  6:348963/149‑1 (MQ=255)
                                          ccAAGATCGAAGGCGTGATCAAGTTCGAGAAGAAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGCTCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTg                                                                                                        >  1:368237/1‑149 (MQ=255)
                                          ccAAGATCGAAGGCGTGATCAAGTTCGAGAAGAAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGCTCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTg                                                                                                        >  5:196912/1‑149 (MQ=255)
                                           cAAGATCGAAGGCGTGATCAAGTTCGAGAAGAAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGCTCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTGt                                                                                                       >  1:32279/1‑149 (MQ=255)
                                           cAAGATCGAAGGCGTGATCAAGTTCGAGAAGAAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGCTCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTGt                                                                                                       <  5:361327/149‑1 (MQ=255)
                                            aaGATCGAAGGCGTGATCAAGTTCGAGAAGAAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGCTCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTGtt                                                                                                      <  2:32279/149‑1 (MQ=255)
                                            aaGATCGAAGGCGTGATCAAGTTCGAGAAGAAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGCTCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTGtt                                                                                                      <  6:196912/149‑1 (MQ=255)
                                             aGATCGAAGGCGTGATCAAGTTCGAGAAGAAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGCTCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTGttt                                                                                                     <  6:207359/149‑1 (MQ=255)
                                              gATCGAAGGCGTGATCAAGTTCGAGAAGAAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGCTACAGAAGCCCCGTCATGCGACGGGGCTTTTTAGTTTGTGGTGAGCCTCTTGCAAAGCTGTTTg                                                                                                    >  7:161383/1‑149 (MQ=255)
                                                  gAAGGCGTGATCAAGTTCGAGAAGAAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGCTCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCtt                                                                                                                 <  3:90299/132‑1 (MQ=255)
                                                  gAAGGCGTGATCAAGTTCGAGAAGAAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGCTCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCtt                                                                                                                 >  4:90299/1‑132 (MQ=255)
                                                   aaGGCGTGATCAAGTTCGAGAAGAAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGCTCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTGTTTGTCTgg                                                                                               >  1:250329/1‑149 (MQ=255)
                                                   aaGGCGTGATCAAGTTCGAGAAGAAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGCTCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTGTTTGTCTgg                                                                                               <  6:300037/149‑1 (MQ=255)
                                                     ggCGTGATCAAGTTCGAGAAGAAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGCTCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTGTTTGTCTGGGt                                                                                             <  3:248109/149‑1 (MQ=255)
                                                       cGTGATCAAGTTCGAGAAGAAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGCTCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTGTTTGTCTGGGTTg                                                                                           <  2:74209/149‑1 (MQ=255)
                                                         tGATCAAGTTCGAGAAGAAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGCTCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTGTTTGTCTGGGTTGcc                                                                                         >  6:128641/1‑149 (MQ=255)
                                                              aaGTTCGAGAAGAAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGCTCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTGTTTGTCTGGGTTGCCGGCGt                                                                                    >  3:185711/1‑149 (MQ=255)
                                                                gTTCGATAAGAAAGGCCAGTTCATCCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCTCTACAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTGTTTGTCTGGGTTGCCGGCGTTg                                                                                  <  6:309743/149‑1 (MQ=255)
                                                                gTTCGAGAAGAAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGCTCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTGTTTGTCTGGGTTGCCGGCGTTg                                                                                  <  2:149095/149‑1 (MQ=255)
                                                                     agaagaAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGCTCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTGTTTGTCTGGGTTGCCGGCGTTGGTTTc                                                                             >  1:247848/1‑149 (MQ=255)
                                                                         gaaAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGCTCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTg                                                                                                        <  3:65803/118‑1 (MQ=255)
                                                                         gaaAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGCTCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTg                                                                                                        >  4:65803/1‑118 (MQ=255)
                                                                         gaaAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGCTCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTGTTTGTCTGGGTTGCCGGCGTTGGTTTCTACg                                                                         >  8:203456/1‑149 (MQ=255)
                                                                            aGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGCTCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTGTTTGTCTGGGTTGCCGGCGTTGGTTTCTACGGTc                                                                      >  2:299551/1‑149 (MQ=255)
                                                                               cGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGCTCCAGAAGCCCCGTCAAGCGACGGGGCTTTTTATTTTGTGGTGAGCCTCATGCAAAGCTGTTTGTCTGGGTTGACGGCGTTGGTTTCTACGTTCGTa                                                                   >  8:337408/1‑149 (MQ=255)
                                                                                 aGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGCTCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTGTTTGTCTGGGTTGCCGGCGTTGGTTTCTACGGTCGTACg                                                                 <  4:72692/149‑1 (MQ=255)
                                                                                       tGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGCTCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTGTTTGTCTGGGTTGCCGGCGTTGGTTTCTACGGTCGTACGGGgccg                                                           >  4:332225/1‑149 (MQ=255)
                                                                                          gcCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGCTCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTGTTTGTCTGGGTTGCCGGCGTTGGTTTCTACGGTCGTACggg                                                               >  5:77108/1‑142 (MQ=255)
                                                                                          gcCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGCTCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTGTTTGTCTGGGTTGCCGGCGTTGGTTTCTACGGTCGTACggg                                                               <  6:77108/142‑1 (MQ=255)
                                                                                            cGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGCTCCATAAGCCCCGTCATGGGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTGTTTGTCTGGGTTGCCGGCGTTGGTTTCTACGGTCGTACGGGGCCGCcgcg                                                      <  7:337408/149‑1 (MQ=255)
                                                                                             gTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGCTCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTGTTTGTCTGGGTTGCCGGCGTTGGTTTCTACGGTCGTACGGGGCCGCCgcgc                                                     >  6:239235/1‑149 (MQ=255)
                                                                                                   tGAGCATCGTCGCCGCTTAATCGCGACGTCGCTCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTGTTTGTCTGGGTTGCCGGCGTTGGTTTCTACGGTCGTACGGGGCCGCCGCGCTCAttt                                               <  2:250329/149‑1 (MQ=255)
                                                                                                        aTCGTCGCCGCTTAATCGCGACGTCGCTCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTGTTTGTCTGGGTTGCCGGCGTTGGTTTCTACGGTCGTACGGGGCCGCCGCGCTCAttttt                                             >  2:48368/1‑146 (MQ=255)
                                                                                                        aTCGTCGCCGCTTAATCGCGACGTCGCTCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTGTTTGTCTGGGTTGCCGGCGTTGGTTTCTACGGTCGTACGGGGCCGCCGCGCTCAttttt                                             <  1:48368/146‑1 (MQ=255)
                                                                                                           gtcgCCGCTTAATCGCGACGTCGCTCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTGTTTGTCTGGGTTGCCGGCGTTGGTTTCTACGGTCGTACGGGGCCGCCGCGCTCATTTTTGCaa                                         >  5:388381/1‑147 (MQ=255)
                                                                                                                       tCGCGACGTCGCTCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTGTTTGTCGGGGTTGCCGGCGTTGGTTTCTACGGTCGTACGGGGCCGCCGCGCTCATTTTTGCAAGAGCCTCATGtttt                           >  6:72307/1‑149 (MQ=255)
                                                                                                                         gcgACGTCGCTCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTGTTTGTCTGGGTTGCCGGCGTTGGTTTCTACGGTCGTACGGGGCCGCCGCGCTCATTTTTGCAAGAGCCTCATGttttt                          <  3:385220/148‑1 (MQ=255)
                                                                                                                          cgACGTCGCTCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTGTTTGTCTGGGTTGCCGGCGTTGGTTTCTACGGTCGTACGGGGCCGCCGCGCTCATTTTTGCAAGAGCCTCATGTTTTTGa                        <  3:40833/149‑1 (MQ=255)
                                                                                                                          cgACGTCGCTCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTGTTTGTCTGGGTTGCCGGCGTTGGTTTCTACGGTCGTACGGGGCCGCCGCGCTCATTTTTGCAAGAGCCTCATGTTTTTGa                        <  5:128641/149‑1 (MQ=255)
                                                                                                                          cgACGTCGCTCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTGTTTGTCTGGGTTGCCGGCGTTGGTTTCTACGGTCGTACGGGGCCGCCGCGCTCATTTTTGCAAGAGCCTCATGTTTTTGa                        <  1:231612/149‑1 (MQ=255)
                                                                                                                          cgACGGCGCTCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGATCATCTTGCAAAGCTGTTTGTCTGGGTTGCCGGCGTTGGTTGCTACGGTGGTACGGGGCCGGCGCGCTAATTTTTGCAAGAGCCTCATGTTTTTGa                        <  8:161383/149‑1 (MQ=255)
                                                                                                                                 gCTCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTGTTTGTCTGGGTTGCCGGCGTTGGTTTCTACGGTCGTATGGGGCCGCCGCGCTCATTTTTGCAAGAGCCTCATGTTTTTGAGTCAtt                  <  1:56751/148‑1 (MQ=255)
                                                                                                                                 gCTCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTGTTTGTCTGGGTTGCCGGCGTTGGTTTCTACGGTCGTACGGGGCCGCCGCGCTCATTTTTGCAAGAGCCTCATGTTTTTGAGTCAtt                  <  4:185711/148‑1 (MQ=255)
                                                                                                                                 gCTCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTGTTTGTCTGGGTTGCCGGCGTTGGTTTCTACGGTCGTACGGGGCCGCCGCGCTCATTTTTGCAAGAGCCTCATGTTTTTGAGTCATTc                 <  6:388381/149‑1 (MQ=255)
                                                                                                                                   tCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTGTTTGTCTGGGTTGCCGGCGTTGGTTTCTACGGTCGTACGGGGCCGCCGCGCTCATTTTTGCAAGAGCCTCATGTTTTTGAGTCATTCAg               >  8:360561/1‑149 (MQ=255)
                                                                                                                                   tCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTGTTTGTCTGGGTTGCCGGCGTTGGTTTCTACGGTCGTACGGGGCCGCCGCGCTCATTTTTGCAAGAGCCTCATGTTTTTGAGTAATTCAg               >  2:272599/1‑149 (MQ=255)
                                                                                                                                    gaaaaaGCCCCGTCGCATGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTGTTTGTCTGGGTTGCCGGCGTTGGTTTCTACGGTCGTACGGGGCCGCCGCGCTCATTTTTGCAAGAGCCTCATGTTTTTGAGTCATTCAGc              >  4:270516/5‑149 (MQ=255)
                                                                                                                                     aaaaaGCCCCGTCGCATGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTGTTTGTCTGGGTTGCCGGCGTTGGTTTCTACGGTCGTACGGGGCCGCCGCGCTCATTTTTGCAAGAGCCTCATGTTTTTGAGTCATTCAGCt             >  1:87972/4‑149 (MQ=255)
                                                                                                                                          gCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTGTTTGTCGGGGTTGCCGGCGTTGGTTTCTACGGTCGTACGGGGCCGCCGCGCTCATTTTTGCAAGAGCCTCATGTTTTTGAGTCATTCAGCTCGTCt        <  5:72307/149‑1 (MQ=255)
                                                                                                                                            cccGTCATGCGAAGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCTAAGCTGTTTGTCTGGTTTGCCGGCGTTGGTTTATACGGTCGTACGGGGCCGCCGCGCTCATTTTTGCAATAGCATCATGTTTTTGAGTCATTCAGATCGTCTGc      >  6:78307/1‑149 (MQ=255)
                                                                                                                                                gcctGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTGTTTGTCTGGGTGGCCGGCGTGGGTTTCTACGGTCGTACGGGGCCGCCGCGCTCATTTTTGCAAGAGCCTCATGTTTTTGAGTCATTCAGCTCGTCTGCcgac  <  6:59757/146‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
CTGGCTACGGCGTTGGCATGGGCAAGGACCACACCTTGTTCGCCAAGATCGAAGGCGTGATCAAGTTCGAGAAGAAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGCTCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTGTTTGTCTGGGTTGCCGGCGTTGGTTTCTACGGTCGTACGGGGCCGCCGCGCTCATTTTTGCAAGAGCCTCATGTTTTTGAGTCATTCAGCTCGTCTGCCGAC  >  NC_002947/801209‑801501

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