Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 933,387 C→T 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,3870CT97.0% 422.1 / 1.4 135A55V (GCC→GTC) ftsKDNA translocase at septal ring sorting daughter chromsomes
Reads supporting (aligned to +/- strand):  ref base C (2/2);  new base T (76/54);  total (79/56)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.76e-01

Only 100 of 135 total aligned reads displayed.

TGGAGGCGTTGCTGATCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGCCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGATGGCTTACACCATTCC  >  NC_000913/933291‑933483
                                                                                                |                                                                                                
tGGAGGCGTTGCTGATCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCt                                                                                                >  1:195272/1‑99 (MQ=255)
tGGAGGCGTTGCTGATCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTg                                                                                               <  2:2029078/100‑1 (MQ=255)
 ggAGGCGTTGCTGATCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGTTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTgg                                                                                              <  1:93983/100‑1 (MQ=255)
 ggAGGCGTTGCTGATCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTgg                                                                                              <  2:292534/100‑1 (MQ=255)
 ggAGGCGTTGCTGATCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTgg                                                                                              <  2:1469088/100‑1 (MQ=255)
 ggAGGCGTTGCTGATCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTgg                                                                                              <  1:3020123/100‑1 (MQ=255)
   aGGCGTTGCTGATCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCa                                                                                            <  1:595794/100‑1 (MQ=255)
      cGTTGCTGATCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGa                                                                                         <  1:2420286/100‑1 (MQ=255)
      cGTTGCTGATCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGa                                                                                         >  2:261012/1‑100 (MQ=255)
       gTTGCTGATCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGaa                                                                                        >  1:660428/1‑100 (MQ=255)
       gTTGCTGATCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGaa                                                                                        >  2:3440097/1‑100 (MQ=255)
         tGCTGATCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAAcc                                                                                      <  1:2900738/100‑1 (MQ=255)
         tGCTGATCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAAcc                                                                                      <  1:3441734/100‑1 (MQ=255)
         tGCTGATCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAAcc                                                                                      <  1:1814588/100‑1 (MQ=255)
           cTGATCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTa                                                                                    >  2:1466290/1‑100 (MQ=255)
           cTGATCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTa                                                                                    >  2:865838/1‑100 (MQ=255)
            tGATCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTAt                                                                                   <  1:1666741/100‑1 (MQ=255)
             gATCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATc                                                                                  <  2:394730/100‑1 (MQ=255)
             gATCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATc                                                                                  >  1:1532754/1‑100 (MQ=255)
               tCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATcc                                                                                 >  1:3452363/1‑99 (MQ=255)
                ccgatctGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCAt                                                                               <  2:701055/94‑1 (MQ=255)
                ccTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCAt                                                                               >  1:309435/1‑100 (MQ=255)
                 cTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATa                                                                              <  2:1366140/100‑1 (MQ=255)
                 cTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATa                                                                              >  1:2783870/1‑100 (MQ=255)
                  ttATTGTCCGGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGATCCCAGCGGGTCGCAAACGGTCTGGCATGAACCTATCCATaa                                                                             >  1:1425241/1‑100 (MQ=255)
                   tATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAAt                                                                            <  1:1879885/100‑1 (MQ=255)
                       gTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATagatcgg                                                                        >  1:701046/1‑94 (MQ=255)
                       gTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAg                                                                        <  2:464045/100‑1 (MQ=255)
                           tGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAgg                                                                       >  1:662051/1‑97 (MQ=255)
                           tGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTgg                                                                    >  1:1346576/1‑100 (MQ=255)
                           tGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGCCATGAACCTATCCATAATTTAGGTg                                                                     >  1:2698073/1‑99 (MQ=255)
                            gTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTggg                                                                   >  2:3357470/1‑100 (MQ=255)
                            gTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTggg                                                                   <  1:2000425/100‑1 (MQ=255)
                               tGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCGGGTCGCAAAGGGTCTCGCATGAACGTATCCATAATTTAGGTGGGATg                                                                >  1:2341980/1‑100 (MQ=255)
                                 ccGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGcc                                                              >  2:2095652/1‑100 (MQ=255)
                                  cGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGccc                                                             >  2:2887125/1‑100 (MQ=255)
                                     cTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGt                                                          <  2:1584926/100‑1 (MQ=255)
                                      tGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTg                                                         <  1:997758/100‑1 (MQ=255)
                                       ggTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGc                                                        <  1:933528/100‑1 (MQ=255)
                                       ggTTGATGGCTGACTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTAACAATAATTTAGGTGGGATACCCgct                                                          >  1:556612/1‑96 (MQ=255)
                                          tGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTg                                                     <  1:3347849/100‑1 (MQ=255)
                                          tGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTg                                                     >  1:787511/1‑100 (MQ=255)
                                           gATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTgg                                                    <  1:412576/100‑1 (MQ=255)
                                           gATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTgg                                                    >  2:2668036/1‑100 (MQ=255)
                                            aTGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGt                                                   <  2:2938431/100‑1 (MQ=255)
                                             tGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGtt                                                  >  1:484390/1‑100 (MQ=255)
                                              ggCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTg                                                 <  2:448135/100‑1 (MQ=255)
                                              ggCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTg                                                 <  2:2095219/100‑1 (MQ=255)
                                                cTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGc                                               >  2:1654893/1‑100 (MQ=255)
                                                cTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGc                                               >  2:1834531/1‑100 (MQ=255)
                                                  gCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGCCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAg                                             <  2:1674881/100‑1 (MQ=255)
                                                   ccTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGa                                            >  1:1061915/1‑100 (MQ=255)
                                                   ccTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGa                                            >  1:1009037/1‑100 (MQ=255)
                                                    tttACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGAt                                           <  2:832349/99‑1 (MQ=255)
                                                     ttACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATa                                          >  2:932414/1‑100 (MQ=255)
                                                     ttACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATa                                          >  1:918783/1‑100 (MQ=255)
                                                     ttACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATa                                          <  1:3323497/100‑1 (MQ=255)
                                                      tACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATAc                                         >  1:1756317/1‑100 (MQ=255)
                                                        ttccgatctAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGc                                       <  2:3478861/92‑1 (MQ=255)
                                                        cTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTCGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGc                                       >  2:472983/1‑100 (MQ=255)
                                                        cTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGc                                       >  1:104541/1‑100 (MQ=255)
                                                        cTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGc                                       <  1:1746716/100‑1 (MQ=255)
                                                         tAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCt                                      <  1:2564842/100‑1 (MQ=255)
                                                          aaGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTg                                     >  2:2706787/1‑100 (MQ=255)
                                                           aGCTTTAACCCTTCGGACCCCAGCTGGTCTCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTCGGATGCCCGGTGCGTGGTTGGCAGATACGCTGt                                    <  1:2095621/100‑1 (MQ=255)
                                                           aGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGt                                    <  1:3308472/100‑1 (MQ=255)
                                                             cTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTc                                  >  1:667802/1‑100 (MQ=255)
                                                                 aaCCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGcagatcggaa                              >  1:3478805/1‑91 (MQ=255)
                                                                 aaCCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTa                              >  1:2310713/1‑100 (MQ=255)
                                                                  aCCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTAt                             <  2:1343101/100‑1 (MQ=255)
                                                                  aCCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTAt                             >  2:812903/1‑100 (MQ=255)
                                                                      ttCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTAttttt                         >  2:742765/1‑100 (MQ=255)
                                                                        cGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGCTCTTTATTTTTgg                       <  2:2359701/100‑1 (MQ=255)
                                                                          gACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCg                     >  2:325488/1‑100 (MQ=255)
                                                                          gACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGAGGGATGCCCGGTGCGCGGTTGTCAGATACGCTGTTCTTTATTTTTGGCg                     >  1:1507789/1‑100 (MQ=255)
                                                                          ctccccAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCg                     <  1:3263724/98‑1 (MQ=255)
                                                                            ccccAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTg                   >  2:1799292/1‑100 (MQ=255)
                                                                            ccccAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTg                   >  2:1925754/1‑100 (MQ=255)
                                                                            ccccAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTg                   >  1:502843/1‑100 (MQ=255)
                                                                            ccccAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTg                   >  2:3187428/1‑100 (MQ=255)
                                                                              ccAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATGATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGAt                 >  1:743414/1‑100 (MQ=255)
                                                                               cttcTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCTCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTGGGCGTGATg                <  1:3163980/97‑1 (MQ=255)
                                                                               cAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGATg                <  2:727863/100‑1 (MQ=255)
                                                                               cAGCTGGTCGCAAACGGCCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGATg                >  2:3219857/1‑100 (MQ=255)
                                                                                aGCTGGTCGCAAACGGCCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGATgg               <  2:2402626/100‑1 (MQ=255)
                                                                                 gCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGATGGc              <  2:3391601/100‑1 (MQ=255)
                                                                                  cTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGATGGCt             >  2:53069/1‑100 (MQ=255)
                                                                                   tGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGATGGCtt            >  2:387156/1‑100 (MQ=255)
                                                                                    ggTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGATGGCTTa           >  1:1548536/1‑100 (MQ=255)
                                                                                    ggTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGATGGCTTa           >  2:2424125/1‑100 (MQ=255)
                                                                                      tCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGATGGCTTaca         <  2:3406406/100‑1 (MQ=255)
                                                                                      tCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGATGGCTTaca         <  2:779227/100‑1 (MQ=255)
                                                                                       cGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGATGGCTTacac        >  2:3438124/1‑100 (MQ=255)
                                                                                       cGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCCGTTCTTTATTTTTGGCGTGATGGCTTacac        >  1:2592087/1‑100 (MQ=255)
                                                                                        gCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGGGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGATGGCTTACAcc       >  1:1312908/1‑100 (MQ=255)
                                                                                         cAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGATGGCTTACACCa      >  1:2048440/1‑100 (MQ=255)
                                                                                          aaaCGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGATGGCTTACACCAt     >  1:2653216/1‑100 (MQ=255)
                                                                                           aaCGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGATGGCTTACACCAtt    <  1:3187376/100‑1 (MQ=255)
                                                                                            aCGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGATGGCTTACACCATTc   <  2:3320093/100‑1 (MQ=255)
                                                                                             cGGCCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGATGGCTTACACCATTcc  >  2:340194/1‑100 (MQ=255)
                                                                                                |                                                                                                
TGGAGGCGTTGCTGATCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGCCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGATGGCTTACACCATTCC  >  NC_000913/933291‑933483

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