Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 937,130 G→T 6.3% E1303* (GAA→TAA)  ftsK → DNA translocase at septal ring sorting daughter chromsomes

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913937,1300GT 6.3% 213.4 / 3.7 64E1303* (GAA→TAA) ftsKDNA translocase at septal ring sorting daughter chromsomes
Reads supporting (aligned to +/- strand):  ref base G (33/27);  new base T (2/2);  total (35/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.87e-01

TGGTGCGGGTGGTTTCGATGGCGCTGAAGAACTGGATCCGTTGTTCGATCAGGCGGTGCAGTTTGTCACTGAAAAACGCAAAGCGTCAATTTCTGGCGTACAGCGTCAGTTCCGCATTGGTTATAACCGTGCAGCGCGTATTATCGAACAGATGGAAGCGCAGGGGATTGTCAGCGAACAGGGGCACAACGGTAATCGTGAAGTGCTGGCCCCACCGCCGTTTGACTAACTAATGCATCGTATGCCGGATAAGGCGCGGTAGCGTCGCATCCGGCACTCTATCAACTGAAAA  >  NC_000913/936985‑937276
                                                                                                                                                 |                                                                                                                                                  
tggtgCGGGTGGTTTCGATGGCGCTGAAGAACTGGATCCGTTGTTCGATCAGGCGGTGCAGTTTGTCACTGAAAAACGCAAAGCGTCAATTTCTGGCGTACAGCGTCAGTTCCGCATTGGTTATAACCGTGCAGCGCGTATTATCGAAc                                                                                                                                                 <  8:95171/149‑1 (MQ=255)
  gtgCGGGTGGTTTCGATGGCGCTGAAGAACTGGATCCGTTGTTCGATCAGGCGGTGCAGTTTGTCACTGAAAAACGCAAAGCGTCAATTTCTGGCGTACAGCGTCAGTTCCACATTGGTTATAACCGTGCAGCGCGTATTATCGaa                                                                                                                                                  >  2:326803/1‑146 (MQ=255)
   tgCGGGTGGTTTCGATGGCGCTGAAGAACTGGATCCGTTGTTCGATCAGGCGGTGCAGTTTGTCACTGAAAAACGCAAAGCGTCAATTTCTGGCGTACAGCGTCAGTTCCGCATTGGTTATAACCGTGCAGCGCGTATTATCGAACAGa                                                                                                                                              <  8:176709/149‑1 (MQ=255)
                gATGGCGCTGAAGAACTGGATCCGTTGTTCGATCAGGCGGTGCAGTTTGTCACTGAAAAACGCAAAGCGTCAATTTCTGGCGTACAGCGTCAGTTCCGCATTGGTTATAACCGTGCAGCGCGTATTATCGAACAGATGGAAGCGCAggg                                                                                                                                 >  8:69864/1‑149 (MQ=255)
                gATGGCGCTGAAGAACTGGATCCGTTGTTCGATCAGGCGGTGCAGTTTGTCACTGAAAAACGCAAAGCGTCAATTTCTGGCGTACAGCGTCAGTTCCGCATTGGTTATAACCGTGCAGCGCGTATTATCGAACAGATGGAAGCGCAgg                                                                                                                                  <  1:281766/148‑1 (MQ=255)
                  tGGCGCTGAAGAACTGGATCCGTTGTTCGATCAGGCGGTGCAGGTTGTCACTGAAAAACGCAAAGCGTCAATTTCTGGCGTACAGCGTCAGTTCCGCATTGGTTATAACCGTGCAGCGCGTATTATCGAACAGATGGAAGCGCAGGGGa                                                                                                                               <  3:336446/149‑1 (MQ=255)
                      gcTGAAGAACTGGATCCGTTGTTCGATCAGGCGGTGCAGTTTGTCACTGAAAAACGCAAAGCGTCAATTTCTGGCGTACAGCGTCAGTTCCGCATTGGTTATAACCGTGCAGCGCGTATTATCGAACAGATGGAAGCGCAGGGGATTGt                                                                                                                           >  1:63172/1‑149 (MQ=255)
                       cTGAAGAACTGGATCCGTTGTTCGATCAGGCGGTGCAGTTTGTCACTGAAAAACGCAAAGCGTCAATTTCTGGCGTACAGCGTCAGTTCCGCATTGGTTATAACCGTGCAGCGCGTATTATCGAACAGATGGAAGCGCAGGGGATTGTc                                                                                                                          <  2:63172/149‑1 (MQ=255)
                         gaagaaCTGGATCCGTTGTTCGATCAGGCGGTGCAGTTTGTCACTGAAAAACGCAAAGCGTCAATTTCTGGCGTACAGCGTCAGTTCCGCATTGGTTATAACCGTGCAGCGCGTATTATCGAACAGATGGAAGCGCAGGGGATTGTCAg                                                                                                                        >  6:355257/1‑149 (MQ=255)
                         gaagaaCTGGATCCGTTGTTCGATCAGGCGGTGCAGTTTGTCACTGAAAAACGCAAAGCGTCAATTTCTGGCGTACAGCGTCAGTTCCGCATTGGTTATAACCGTGCAGCGCGTATTATCGAACAGATGGAAGCGCAGGGGATTGTCAg                                                                                                                        >  8:150824/1‑149 (MQ=255)
                               cTGGATCCGTTGTTCGATCAGGCGGTGCAGTTTGTCACTGAAAAACGCAAAGCGTCAATTTCTGGCGTACAGCGTCAGTTCCGCATTGGTTATAACCGTGCAGCGCGTATTATCGAACAGATGGAAGCGCAGGGGATTGTCAGCGaa                                                                                                                    >  4:299737/1‑147 (MQ=255)
                               cTGGATCCGTTGTTCGATCAGGCGGTGCAGTTTGTCACTGAAAAACGCAAAGCGTCAATTTCTGGCGTACAGCGTCAGTTCCGCATTGGTTATAACCGTGCAGCGCGTATTATCGAACAGATGGAAGCGCAGGGGATTGTCAGCGAACa                                                                                                                  <  2:240026/149‑1 (MQ=255)
                                tGGATCCGTTGTTCGATCAGGCGGTGCAGTTTGTCACTGAAAAACGCAAAGCGTCAATTTCTGGCGTACAGCGTCAGTTCCGCATTGGTTATAACCGTGCAGCGCGTATTATCGAACAGATGGAAGCGCAGGGGATTGTCAGCGAACAg                                                                                                                 >  6:328956/1‑149 (MQ=255)
                                  gATCCGTTGTTCGATCAGGCGGTGCAGTTTGTCACTGAAAAACGCAAAGCGTCAATTTCTGGCGTACAGCGTCAGTTCCGCATTGGTTATAACCGTGCAGCGCGTATTATCGAACAGATGGAAGCGCAGGGGATTGTCAGCGAACAggg                                                                                                               <  8:347113/149‑1 (MQ=255)
                                  gATCCGTTGTTCGATCAGGCGGTGCAGTTTGTCACTGAAAAACGCAAAGCGTCAATTTCTGGCGTACAGCGTCAGTTCCGCATTGGTTATAACCGTGCAGCGCGTATTATCGAACAGATGGAAGCGCAGGGGATTGTCAGCGAACAgg                                                                                                                <  3:299737/148‑1 (MQ=255)
                                    tCCGTTGTTCGATCAGGCGGTGCAGTTTGTCACTGAAAAACGCAAAGCGTCAATTTCTGGCGTACAGCGTCAGTTCCGCATTGGTTATAACCGTGCAGCGCGTATTATCGAACAGATGGAAGCGCAGGGGATTGTCAGCGAACAgggg                                                                                                              >  1:203693/1‑148 (MQ=255)
                                               aTCAGGCGGTGCAGTTTGTCACTGAAAAACGCAAAGCGTCAATTTCTGGCGTACAGCGTCAGTTCCGCATTGGTTATAACCGTGCAGCGCGTATTATCGAACAGATGGAAGCTCAGGGGATTGTCAGCGAACAGGGGCACAACGGTAAt                                                                                                  >  5:300970/1‑149 (MQ=255)
                                                tCAGGCGGTGCAGTTTGTCACTGAAAAACGCAAAGCGTCAATTTCTGGCGTACAGCGTCAGTTCCGCATTGGTTATAACCGTGCAGCGCGTATTATCGAACAGATGGAAGCGCAGGGGATTGTCAGCGAACAGGGGCACAACGGTAATc                                                                                                 <  7:150824/149‑1 (MQ=255)
                                                     cGGTGCAGTTTGTCACTGAAAAACGCAAAGCGTCAATTTCTGGCGTACAGCGTCAGTTCCGCATTGGTTATAACCGTGCAGCGCGTATTATCGAACAGATGGAAGCGCAGGGGATTGTCAGCGAACAGGGGCACAACGGTAATCGTGaa                                                                                            <  7:69864/149‑1 (MQ=255)
                                                        tGCAGTTTGTCACTGAAAAACGCAAAGCGTCAATTTCTGGCGTACAGCGTCAGTTCCGCATTGGTTATAACCGTGCAGCGCGTATTATCGAACAGATGGAAGCGCAGGGGATTGTCAGCGAACAGGGGCACAACGGTAATCGTGAAGTg                                                                                         >  6:22054/1‑149 (MQ=255)
                                                                  cACTGAAAAACGCAAAGCGTCAATTTCTGGCGTACAGCGTCAGTTCCGCATTGGTTATAACCGTGCAGCGCGTATTATCGAACAGATGGAAGCGCAGGGGATTGTCAGCGAACAGGGGCACAACGGTAATCGTGAAGTGCTGGCCCCAc                                                                               >  1:220093/1‑149 (MQ=255)
                                                                  cACTGAAAAACGCAAAGCGTCAATTTCTGGCGTACAGCGTCAGTTCCGCATTGGTTATAACCGTGCAGCGCGTATTATCGAACAGATGGAAGCGCAGGGGATTGTCAGCGAACAGGGGCACAACGGTAATCGTGAAGTGCTGGCCCCAc                                                                               >  6:46724/1‑149 (MQ=255)
                                                                   aCTGAAAAACGCAAAGCGTCAATTTCTGGCGTACAGCGTCAGTTCCGCATTGGTTATAACCGTGCAGCGCGTATTATCGAACAGATGGAAGCGCAGGGGATTGTCAGCGAACAGGGGCACAACGGTAATCGTGaa                                                                                            <  5:75574/135‑1 (MQ=255)
                                                                   aCTGAAAAACGCAAAGCGTCAATTTCTGGCGTACAGCGTCAGTTCCGCATTGGTTATAACCGTGCAGCGCGTATTATCGAACAGATGGAAGCGCAGGGGATTGTCAGCGAACAGGGGCACAACGGTAATCGTGaa                                                                                            >  6:75574/1‑135 (MQ=255)
                                                                   aCTGAAAAACGCAAAGCGTCAATTTCTGGCGTACAGCGTCAGTTCCGCATTGGTTATAACCGTGCAGCGCGTATTATCGAACAGATGGAAGCGCAGGGGATTGTCAGCGAACAGGGGCACAACGGTAATCGTGAAGTGCTGGCCCCAcc                                                                              >  3:72354/1‑149 (MQ=255)
                                                                    cTGAAAAACGCAAAGCGTCAATTTCTGGCGTACAGCGTCAGTTCCGCATTGGTTATAACCGTGCAGCGCGTATTATCGAACAGATGGAAGCGCAGGGGATTGTCAGCGAACAGGGGCCCAACGGTAATCGTGAAGTGCTGGCCCCAccg                                                                             >  4:62932/1‑149 (MQ=255)
                                                                    cTGAAAAACGCAAAGCGTCAATTTCTGGCGTACAGCGTCAGTTCCGCATTGGTTATAACCGTGCAGCGCGTATTATCGAACAGATGGAAGCGCAGGGGATTGTCAGCGAACAGGGGCACAACGGTAATCGTGAAGTGCTGGCCCCAccg                                                                             >  8:99085/1‑149 (MQ=255)
                                                                     tGAAAAACGCAAAGCGTCAATTTCTGGCGTACAGCGTCAGTTCCGCATTGGTTATAACCGTGCAGCGCGTATTATCGAACAGATGGAAGCGCAGGGGATTGTCAGCGAACAGGGGCACAACGGTAATCGTGAAGTGCTGGCCCCAccgc                                                                            >  1:224209/1‑149 (MQ=255)
                                                                     tGAAAAACGCAAAGCGTCAATTTCTGGCGTACAGCGTCAGTTCCGCATTGGTTATAACCGTGCAGCGCGTATTATCGAACAGATGGAAGCGCAGGGGATTGTCAGCGAACAGGGGCACAACGGTAATCGTGAAGTGCTGGCCCCAccgc                                                                            >  7:240804/1‑149 (MQ=255)
                                                                       aaaaaCGCAAAGCGTCAATTTCTGGCGTACAGCGTCAGTTCCGCATTGGTTATAACCGTGCAGCGCGTATTATCGAACAGATGGAAGCGCAGGGGATTGTCAGCGAACAGGGGCACAACGGTAATCGTGAAGTGCTGGCCCCAccgccg                                                                          >  6:262774/1‑149 (MQ=255)
                                                                        aaaaCGCAAAGCGTCAATTTCTGGCGTACAGCGTCAGTTCCGCATTGGTTATAACCGTGCAGCGCGTATTATCGAACAGATGGAAGCGCAGGGGATTGTCAGCGAACAGGGGCACAACGGTAATCGTGAAGTGCTGGCCCCACCGCCGt                                                                         <  2:224209/149‑1 (MQ=255)
                                                                         aaaCGCAAAGCGTCAATTTCTGGCGTACAGCGTCAGTTCCGCATTGGTTATAACCGTGCAGCGCGTATTATCGAACAGATGGAAGCGCAGGGGATTGTCAGCGAACAGGGGCACAACGGTAATCGTGAAGTGCTGGCCCCACCGCCGtt                                                                        >  5:368389/1‑149 (MQ=255)
                                                                          aaCGCAAAGCGTCAATTTCTGGCGTACAGCGTCAGTTCCGCATTGGTTATAACCGTGCAGCGCGTATTATCGAACAGATGGAAGCGCAGGGGATTGTCAGCGAACAGGGGCACAACGGTAATCGTGAAGTGCTGGCCCCACCGCGGttt                                                                       <  5:355257/149‑1 (MQ=255)
                                                                          aaCGCAAAGCGTCAATTTCTGGCGTACAGCGTCAGTTCCGCATTGGTTATAACCGTGCAGCGCGTATTATCGAACAGATGGAAGCGCAGGGGATTGTCAGCGAACAGGGGCACAACGGTAATCGTGAAGTGCTGGCCCCACCGCCGttt                                                                       <  7:99085/149‑1 (MQ=255)
                                                                           aCGCAAAGCGTCAATTTCTGGCGTACAGCGTCAGTTCCGCATTGGTTATAACCGTGCAGCGCGTATTATCGAACAGATGGAAGCGCAGGGGATTGTCAGCGAACAGGGGCACAACGGTAATCGTGAAGTGCTGGCCCCACCGCCGttt                                                                       >  3:269848/1‑148 (MQ=255)
                                                                               aaaGCGTCAATTTCTGGCGTACAGCGTCAGTTCCGCATTGGTTATAACCGTGCAGCGCGTATTATCGAACAGATGGAAGCGCAGGGGATTGTCAGCGAACAGGGGCACAACGGTAATCGTGAAGTGCTGGCCCCACCGCCGTTTGact                                                                   >  8:1241/1‑148 (MQ=255)
                                                                               aaaGCGTCAATTTCTGGCGTACAGCGTCAGTTCCGCATTGGTTATAACCGTGCAGCGCGTATTATCGAACAGATGGAAGCGCAGGGGATTGTCAGCGAACAGGGGCACAACGGTAATCGTGAAGTGCTGGCCCCACCGCCGTTTGact                                                                   >  2:304126/1‑148 (MQ=255)
                                                                               aaaGCGTCAATTTCTGGCGTACAGCGTCAGTTCCGCATTGGTTATAACCGTGCAGCGCGTATTATCGAACAGATGGAAGCGCAGGGGATTGTCAGCGAACAGGGGCACAACGGTAATCGTGAAGTGCTGGCCCCACCGCCGTGTGact                                                                   >  3:238194/1‑148 (MQ=255)
                                                                                aaGCGTCAATTTCTGGCGTACAGCGTCAGTTCCGCATTGGTTATAACCGTGCAGCGCGTATTATCGAACAGATGGAAGCGCAGGGGATTGTCAGCGAACAGGGGCACAACGGTAATCGTGAAGTGCTGGCCCCACCGCCGTTTGactaa                                                                 <  5:328956/149‑1 (MQ=255)
                                                                                aaGCGTCAATTTCTGGCGTACAGCGTCAGTTCCACATTGGTTATAACCGTGCAGCGCGTATTATCGAACAGATGGAAGCGCAGGGGATTGTCAGCGAACAGGGGCACAACGGTAATCGTGAAGTGCTGGCCCCACCGCCGTTTGactaa                                                                 <  1:326803/149‑1 (MQ=255)
                                                                                  gCGTCAATTTCTGGCGTACAGCGTCAGTTCCGCATTGGTTATAACCGTGCAGCGCGTATTATCGAACAGATGGAAGCGCAgg                                                                                                                                  <  1:91385/82‑1 (MQ=255)
                                                                                  gCGTCAATTTCTGGCGTACAGCGTCAGTTCCGCATTGGTTATAACCGTGCAGCGCGTATTATCGAACAGATGGAAGCGCAgg                                                                                                                                  >  2:91385/1‑82 (MQ=255)
                                                                                         tttCTGGCGTACAGCGTCAGTTCCGCATTGGTTATAACCGTGCAGCGCGTATTATCGAACAGATGGAAGCGCAGGGGATTGTCAGCGAACAGGGGCACAACGGTAATCGTGAAGTGCTGGCCCCACCGCCGTTTGACTAACTAAt                                                            <  3:185829/145‑1 (MQ=255)
                                                                                         tttCTGGCGTACAGCGTCAGTTCCGCATTGGTTATAACCGTGCAGCGCGTATTATCGAACAGATGGAAGCGCAGGGGATTGTCAGCGAACAGGGGCACAACGGTAATCGTGAAGTGCTGGCCCCACCGCCGTTTGACTAACTAAt                                                            >  4:185829/1‑145 (MQ=255)
                                                                                                cGTACAGCGTCAGTTCCGCATTGGTTATAACCGTGCAGCGCGTATTATCGAACAGATGGAAGCGCAGGGGATTGTCAGCGAACAGGGGCACAACGGTAATCGTGAAGTGCTGGCCCCACCGCCGTTTGACTAACTAATGCATCGTATGc                                                 <  6:368389/149‑1 (MQ=255)
                                                                                                    cAGCGTCAGTTCCGCATTGGTTATAACCGTGCAGCGCGTATTATCGAACAGATGGAAGCGCAGGGGATTGTCAGCGAACAGGGGCACAACGGTAATCGTGAAGTGCTGGCCCCACCGCCGTTTGACTAACTAATGCATCGTATGCCgg                                              <  7:1241/148‑1 (MQ=255)
                                                                                                     aGCGTCAGTTCCGCATTGGTTATAACCGTGCAGCGCGTATTATCGAACAGATGGAAGCGCAGGGGATTGTCAGCGAACAGGGGCACAACGGTAATCGTGAAGTGCTGGCCCCACCGCCGTTTGACTAACTAATGCATCGTATGCCGGAt                                            <  8:240804/149‑1 (MQ=255)
                                                                                                       cGTCAGTTCCGCATTGGTTATAACCGTGCAGCGCGTATTATCTAACAGATGGAAGCGCAGGGGATTGTCAGCGAACAGGGGCACAACGGTAATCGTGAAGTGCTGGCCCCACCGCCGTTTGac                                                                    >  4:302511/1‑123 (MQ=255)
                                                                                                       cGTCAGTTCCGCATTGGTTATAACCGTGCAGCGCGTATTATCTAACAGATGGAAGCGCAGGGGATTGTCAGCGAACAGGGGCACAACGGTAATCGTGAAGTGCTGGCCCCACCGCCGTTTGac                                                                    <  7:137970/123‑1 (MQ=255)
                                                                                                       cGTCAGTTCCGCATTGGTTATAACCGTGCAGCGCGTATTATCTAACAGATGGAAGCGCAGGGGATTGTCAGCGAACAGGGGCACAACGGTAATCGTGAAGTGCTGGCCCCACCGCCGTTTGac                                                                    >  8:137970/1‑123 (MQ=255)
                                                                                                       cGTCAGTTCCGCATTGGTTATAACCGTGCAGCGCGTATTATCTAACAGATGGAAGCGCAGGGGATTGTCAGCGAACAGGGGCACAACGGTAATCGTGAAGTGCTGGCCCCACCGCCGTTTGac                                                                    <  3:302511/123‑1 (MQ=255)
                                                                                                             ttCCGCATTGGTTATAACCGTGCAGCGCGTATTATCGAACAGATGGAAGCGCAGGGGATTGTCAGCGAACAGGGGCACAACGGTAATCGTGAAGTGCTGGCCCCACCGCCGTTTGacta                                                                  <  8:281578/119‑1 (MQ=255)
                                                                                                             ttCCGCATTGGTTATAACCGTGCAGCGCGTATTATCGAACAGATGGAAGCGCAGGGGATTGTCAGCGAACAGGGGCACAACGGTAATCGTGAAGTGCTGGCCCCACCGCCGTTTGacta                                                                  >  7:281578/1‑119 (MQ=255)
                                                                                                                  cATTGGTTATAACCGTGCAGCGCGTATTATCGAACAGATGGAAGCGCAGGGGATTGTCAGCGAACAGGGGCACAACGGTAATCGTGAAGTGCTGGCCCCACCGCCGTTTGACTAACTAATGCATCGTATGCCGGATAAGGCGCGGTAGc                               <  5:262774/149‑1 (MQ=255)
                                                                                                                       gTTATAACCGTGCAGCGCGTATTATCGAACAGATGGAAGCGCAGGGGATTGTCAGCGAACAGGGGCACAACGGTAATCGTGAAGTGCTGGCCCCACCGCCGTTTGACTAACTAATGCATCGTATGCCGGATAAGGCGCGGTAGCGTCg                           >  3:20648/1‑148 (MQ=255)
                                                                                                                       gTTATAACCGTGCAGCGCGTATTATCGAACAGATGGAAGCGCAGGGGATTGTCAGCGAACAGGGGCACAACGGTAATCGTGAAGTGCTGGCCCCACCGCCGTTTGACTAACTAATGCATCGTATGCCGGATAAGGCGCGGTAGCGTCGc                          <  5:46724/149‑1 (MQ=255)
                                                                                                                        ttATAACCGTGCAGCGCGTATTATCGAACAGATGGAAGCGCAGGGGATTGTCAGCGAACAGGGGCACAACGGTAATCGTGAAGTGCTGGCCCAACCGCCGTTTGACTAACTAATGCATCGTATGCCGGATAAGGCGCGGTAGCGTCGCa                         >  1:74298/1‑149 (MQ=255)
                                                                                                                             aCCGTGCAGCGCGTATTATCGAACAGATGGAAGCGCAGGGGATTGTCAGCGAACAGGGGCACAACGGTAATCGTGAAGTGCTGGCCCCCCCGCCGTTTGACTAACTAATGCATCGTATGCCGGATAAGGCGCGGTAGCGTCGCATCCgg                    >  3:285927/1‑149 (MQ=255)
                                                                                                                               cGTGCAGCGCGTATTATCGAACAGATGGAAGCGCAGGGGATTGTCAGCGAACAGGGGCACAACGGTAATCGTGAAGTGCTGGCCCCACCGCCGTTTGACTAACTAATGCATCGTATGCCGGATAAGGCGCGGTAGCGTCGCATCCGGCa                  <  2:220093/149‑1 (MQ=255)
                                                                                                                                     gcgcGTATTATCGAACAGATGGAAGCGCAGGGGATTGTCTGCGAACAGGGGCACAACGGTAATAGTGAAGTGCTGGCGCCACCGCCGTTTGACTGAGGAATGCATCGTATGCCGGATAAGGCGGGGTTGGGTCGCTTGCGGCACTCTAt            >  6:352511/1‑149 (MQ=255)
                                                                                                                                           attatCGAACAGATGGAAGCGCAGGGGATTGTCAGCGAACAGGGGCACAACGGTAATCGTGAAGTGCTGGCCCCACCGCCGTTTGACTAACTAATGCATCGTATGCCGGATAAGGCGCGGTAGCGTCGCATCCGGCACTCTATCAACt       >  3:333920/1‑148 (MQ=255)
                                                                                                                                             tatCGAACAGATGGAAGCGCAGGGGATTGTCAGCGAACAGGGGCACAACGGTAATCGTGAAGTGCTGGCCCCACCGCCGTTTGACTAACTAATGCATCGTATGCCGGATAAGGCGCGGTAGCGTCGCATCCGGCACTCTATCAACTGa     >  2:246116/1‑148 (MQ=255)
                                                                                                                                               tCGAACAGATGGAAGCGCAGGGGATTGTCAGCGAACAGGGGCACAACGGTAATCGTGAAGTGCTGGCCCCACCGCCGTTTGACTAACTAATTCATCGTATGCCGGATAAGGCGCGGTAGCGTCGCATCCGGCACTCTATCAACTGaaaa  <  3:62932/149‑1 (MQ=255)
                                                                                                                                               tCGAACAGATGGAAGCGCAGGGGATTGTCAGCGAACAGGGGCACAACGGTAATCGTGAAGTGCTGGCCCCACCGCCGTTTGACTAACTAATGCATCGTATGCCGGATAAGGCGCGGTAGCGTCGCATCCGGCACTCTATCAACTGaaaa  <  5:22054/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
TGGTGCGGGTGGTTTCGATGGCGCTGAAGAACTGGATCCGTTGTTCGATCAGGCGGTGCAGTTTGTCACTGAAAAACGCAAAGCGTCAATTTCTGGCGTACAGCGTCAGTTCCGCATTGGTTATAACCGTGCAGCGCGTATTATCGAACAGATGGAAGCGCAGGGGATTGTCAGCGAACAGGGGCACAACGGTAATCGTGAAGTGCTGGCCCCACCGCCGTTTGACTAACTAATGCATCGTATGCCGGATAAGGCGCGGTAGCGTCGCATCCGGCACTCTATCAACTGAAAA  >  NC_000913/936985‑937276

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