Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 933,387 C→T 100% A55V (GCC→GTC)  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_000913933,3870CT98.9% 303.3 / ‑4.7 95A55V (GCC→GTC) ftsKDNA translocase at septal ring sorting daughter chromsomes
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (56/37);  total (57/38)
Fisher's exact test for biased strand distribution p-value = 4.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.13e-01

GGAGGCGTTGCTGATCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGCCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGATGGCTTACACCATTCC  >  NC_000913/933292‑933483
                                                                                               |                                                                                                
ggAGGCGCTGCTGATCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGGCCCCAGCTGGTCGCAAAAGGTCTgg                                                                                              >  2:717185/1‑100 (MQ=255)
  aGGCGTTGCTGATCCTTATTGTCCTGTTTGCCGTCTGGCTGATGGCTGCCTTACTAAGTTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCa                                                                                            <  2:1394373/100‑1 (MQ=255)
     cGTTGCTGATCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGa                                                                                         >  1:2448165/1‑100 (MQ=255)
     cGTTGCTGATCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGa                                                                                         >  1:1878488/1‑100 (MQ=255)
        tGCTGATCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAAcc                                                                                      >  1:1242981/1‑100 (MQ=255)
        tGCTGATCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAAcc                                                                                      <  2:27340/100‑1 (MQ=255)
        tGCTGATCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAAcc                                                                                      <  2:35012/100‑1 (MQ=255)
         gCTGATCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCt                                                                                     >  1:2349221/1‑100 (MQ=255)
           tGATCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTAt                                                                                   >  2:356484/1‑100 (MQ=255)
            gATCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATc                                                                                  <  2:1930371/100‑1 (MQ=255)
             aTCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCATACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATcc                                                                                 <  2:294100/100‑1 (MQ=255)
               ccTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCAt                                                                               <  1:2429715/100‑1 (MQ=255)
                cTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCAt                                                                               >  2:1038540/1‑99 (MQ=255)
                    ttGTCCTGTTTGCCGTCTGGTTGATGGCTGTCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAAttt                                                                          <  1:1278153/100‑1 (MQ=255)
                     tGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTa                                                                         <  1:594501/100‑1 (MQ=255)
                      gTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAg                                                                        >  2:2351661/1‑100 (MQ=255)
                        ccTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGt                                                                      >  2:1510720/1‑100 (MQ=255)
                         cTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTg                                                                     <  2:1651130/100‑1 (MQ=255)
                         cTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTg                                                                     <  1:2351676/100‑1 (MQ=255)
                          tGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTgg                                                                    >  2:1437770/1‑100 (MQ=255)
                           tttcgCCGTCTGGTTGATAGATGCCTTAACAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTgg                                                                    >  2:442339/5‑99 (MQ=255)
                           gTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTggg                                                                   <  2:1384345/100‑1 (MQ=255)
                            tttGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGa                                                                  >  1:1549279/1‑100 (MQ=255)
                             ttGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCACACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGAt                                                                 >  1:437226/1‑100 (MQ=255)
                              tGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATg                                                                >  1:2500305/1‑100 (MQ=255)
                               gCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGc                                                               >  1:1969060/1‑100 (MQ=255)
                                 cGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGccc                                                             <  1:2211388/100‑1 (MQ=255)
                                 cGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGccc                                                             <  1:2162652/100‑1 (MQ=255)
                                     tGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCCATCCATAATTTAGGTGGGCTGCCCGGt                                                          >  2:1461286/1‑99 (MQ=255)
                                     tGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGATCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTg                                                         >  2:427433/1‑100 (MQ=255)
                                      ggTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGc                                                        >  2:1266731/1‑100 (MQ=255)
                                      ggTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGc                                                        >  1:36840/1‑100 (MQ=255)
                                       gTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGc                                                        <  2:877006/99‑1 (MQ=255)
                                       gTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGCACCCCAGCTGGTCGCCCACGGTCTGGCATGCACCTATCCCTAATTTAGGTGGGATGCCCgg                                                           >  2:1604767/1‑96 (MQ=255)
                                         tGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTg                                                     <  1:1203699/100‑1 (MQ=255)
                                         tGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTg                                                     <  2:616837/100‑1 (MQ=255)
                                         tGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCAGAATTTAGGTGGGATGCCCGGAGCGt                                                      >  1:1648617/1‑99 (MQ=255)
                                          gATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTgg                                                    >  2:1311183/1‑100 (MQ=255)
                                            tGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGtt                                                  >  2:1755685/1‑100 (MQ=255)
                                             ggCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTg                                                 <  2:898057/100‑1 (MQ=255)
                                             ggCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTg                                                 >  2:311656/1‑100 (MQ=255)
                                             ggCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTg                                                 <  1:1549628/100‑1 (MQ=255)
                                             ggCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTg                                                 <  1:42971/100‑1 (MQ=255)
                                              gCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTgg                                                >  2:244731/1‑100 (MQ=255)
                                              gCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTgg                                                >  2:2063424/1‑100 (MQ=255)
                                              gCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTg                                                 <  2:2295704/99‑1 (MQ=255)
                                                tGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCa                                              <  1:317400/100‑1 (MQ=255)
                                                tGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCa                                              <  1:1583377/100‑1 (MQ=255)
                                                tGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCa                                              >  2:1400214/1‑100 (MQ=255)
                                                tGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCa                                              <  1:2319768/100‑1 (MQ=255)
                                                 gCCTTACTAAGCTTTAACCATTCGGACCGCAGCTGGTCACAAACGGTCTGCCATGATCCACTCCCTCATTTAGCTGGGAGGCCCGGTGTGTGGTTCGCa                                              >  2:216101/1‑99 (MQ=255)
                                                  ccTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGa                                            >  1:808524/1‑100 (MQ=255)
                                                  ccTTACTAAACTTTAACCCTTCGGACCCCAGCTGGTCGCAAACCGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGa                                            <  1:2063029/100‑1 (MQ=255)
                                                     tACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATAc                                         >  2:1714437/1‑100 (MQ=255)
                                                     tACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATAc                                         >  1:774192/1‑100 (MQ=255)
                                                     tACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATAc                                         >  1:2244788/1‑100 (MQ=255)
                                                       cTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCCTAATCCAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGc                                       >  2:2157773/1‑100 (MQ=255)
                                                       cTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGc                                       <  2:1889478/100‑1 (MQ=255)
                                                       cTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGc                                       >  1:232542/1‑100 (MQ=255)
                                                          aGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGt                                    <  2:853676/100‑1 (MQ=255)
                                                           gCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGtt                                   >  1:9216/1‑100 (MQ=255)
                                                           gCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGtt                                   >  1:638779/1‑100 (MQ=255)
                                                            cTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATCCGCTGTTc                                  >  2:570392/1‑100 (MQ=255)
                                                            cTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTc                                  <  1:164742/100‑1 (MQ=255)
                                                            cTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTc                                  >  1:2303106/1‑100 (MQ=255)
                                                             tttAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCt                                 >  2:2463233/1‑100 (MQ=255)
                                                             tttAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCt                                 <  1:2044889/100‑1 (MQ=255)
                                                               tAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCttt                               <  2:2115012/100‑1 (MQ=255)
                                                                aaCCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTa                              <  2:586658/100‑1 (MQ=255)
                                                                 aCCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCACGGTGCGTGGTTGGCAGATACGCTGTTCTTTAt                             >  1:1644976/1‑100 (MQ=255)
                                                                  cgcTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTAtt                            >  1:1725592/3‑100 (MQ=255)
                                                                    cTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTAtttt                          <  2:1781517/100‑1 (MQ=255)
                                                                    cTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTAtttt                          >  2:2535544/1‑100 (MQ=255)
                                                                    cTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTAtttt                          >  1:1595842/1‑100 (MQ=255)
                                                                     ttCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTAtttt                          <  2:2509589/99‑1 (MQ=255)
                                                                       cGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTgg                       >  1:785224/1‑100 (MQ=255)
                                                                       cGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTgg                       <  2:91883/100‑1 (MQ=255)
                                                                       cGGACCCCAGCTGGTCACAAACGGACTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGAGGTTGGCAGATACGCTGTTCTTTATTTTTg                        >  2:271790/1‑99 (MQ=255)
                                                                        ggACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGc                      >  1:2195827/1‑100 (MQ=255)
                                                                           ccccAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTg                   >  1:1177123/1‑100 (MQ=255)
                                                                           ccccAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTg                   >  1:2422887/1‑100 (MQ=255)
                                                                            cccAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGATGTTCTTTATTTTTGGCGTg                   >  2:92442/1‑99 (MQ=255)
                                                                             ccAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGAt                 >  2:255825/1‑100 (MQ=255)
                                                                             ccAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGAt                 >  2:1063301/1‑100 (MQ=255)
                                                                             ccAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGAt                 >  1:2552529/1‑100 (MQ=255)
                                                                             ccAGCTGGTCGCAAACGGCCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGAt                 <  1:138143/100‑1 (MQ=255)
                                                                             ccAGCTGGTCGAAAACGGTCTGGCATGAGCCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATGCGCTGTGCTTTATTTTTGGCGTGAt                 <  2:642021/100‑1 (MQ=255)
                                                                              cAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGATg                <  2:1241122/100‑1 (MQ=255)
                                                                                   ggTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGATGGCTTa           >  2:171870/1‑100 (MQ=255)
                                                                                      cGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGATGGCTTacac        <  2:2552513/100‑1 (MQ=255)
                                                                                      cGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGATGGCTTacac        >  1:639622/1‑100 (MQ=255)
                                                                                      cGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGATGGCTTacac        >  1:2228499/1‑100 (MQ=255)
                                                                                        cAAACGGTCTAGCATGAACCTCTCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGATGGCTTacac        >  2:749991/1‑98 (MQ=255)
                                                                                           aCGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGATGGCTTACACCATTc   <  1:1063309/100‑1 (MQ=255)
                                                                                            cGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGATGGCTTACACCATTcc  >  2:1529045/1‑100 (MQ=255)
                                                                                               |                                                                                                
GGAGGCGTTGCTGATCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGCCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGATGGCTTACACCATTCC  >  NC_000913/933292‑933483

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