Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 933,387 C→T 68.0% 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,3870CT68.0% 151.0 / 112.4 131A55V (GCC→GTC) ftsKDNA translocase at septal ring sorting daughter chromsomes
Reads supporting (aligned to +/- strand):  ref base C (24/18);  new base T (53/36);  total (77/54)
Fisher's exact test for biased strand distribution p-value = 8.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

Only 100 of 133 total aligned reads displayed.

TCTGGAGGCGTTGCTGATCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGCCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGATGGCTTACACCATTCCCGT  >  NC_000913/933289‑933486
                                                                                                  |                                                                                                   
 cTGGAGGCGTTGCTGATCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCt                                                                                                   >  2:488532/1‑100 (MQ=255)
  tGGAGGCGTTGCTGATCCTTATTGTACTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTg                                                                                                  >  1:3025946/1‑100 (MQ=255)
   ggAGGCGTTGCTGATCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTgg                                                                                                 <  2:105358/100‑1 (MQ=255)
   ggAGGCGTTGCTGATCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGCCTgg                                                                                                 <  1:2499320/100‑1 (MQ=255)
   ggAGGCGTTGCTGATCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGCCTgg                                                                                                 >  2:430577/1‑100 (MQ=255)
   ggAGGCGTTGCTGATCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGCCTgg                                                                                                 <  1:151331/100‑1 (MQ=255)
    gAGGCGTTGCTGATCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGCCTGGc                                                                                                <  2:87274/100‑1 (MQ=255)
     aGGCGTTGCTGATCCTTGTTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGCCTGGCa                                                                                               <  1:2572023/100‑1 (MQ=255)
     aGGCGTTGCTGATCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGCCTGGCa                                                                                               <  2:3046803/100‑1 (MQ=255)
      ggCGTTGCTGATCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGCCTGGCAt                                                                                              >  1:1074226/1‑100 (MQ=255)
       gCGTTGCTGATCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATg                                                                                             <  2:2885996/100‑1 (MQ=255)
        cGTTGCTGATCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGa                                                                                            >  2:2977156/1‑100 (MQ=255)
           tGCTGATCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAAcc                                                                                         >  1:2142705/1‑100 (MQ=255)
           tGCTGATCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGCCTGGCATGAAcc                                                                                         <  2:1811610/100‑1 (MQ=255)
            gCTGATCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCt                                                                                        >  2:1369431/1‑100 (MQ=255)
            gCTGATCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGCCTGGCATGAACCt                                                                                        >  2:1789406/1‑100 (MQ=255)
             cTGATCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTa                                                                                       >  2:2361647/1‑100 (MQ=255)
              tGATCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTAt                                                                                      >  1:1332634/1‑100 (MQ=255)
               gATCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATc                                                                                     >  2:410119/1‑100 (MQ=255)
               gATCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATc                                                                                     >  1:1881174/1‑100 (MQ=255)
               gATCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGCCTGGCATGAACCTATc                                                                                     >  1:2975790/1‑100 (MQ=255)
                 tCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCa                                                                                   >  2:2213390/1‑100 (MQ=255)
                 tCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGCCTGGCATGAACCTATCCa                                                                                   <  2:1836629/100‑1 (MQ=255)
                     tATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAAt                                                                               <  1:413360/100‑1 (MQ=255)
                     tATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAAt                                                                               <  2:1988394/100‑1 (MQ=255)
                      aTTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAAtt                                                                              <  1:2300902/100‑1 (MQ=255)
                      aTTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAAtt                                                                              >  2:2926892/1‑100 (MQ=255)
                       ttGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAAttt                                                                             <  1:2301807/100‑1 (MQ=255)
                        tGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTa                                                                            <  2:1420055/100‑1 (MQ=255)
                         gTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGCCTGGCATGAACCTATCCATAATTTAg                                                                           >  2:1932412/1‑100 (MQ=255)
                          tCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAgg                                                                          >  2:629998/1‑100 (MQ=255)
                           ccTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGt                                                                         <  2:1182396/100‑1 (MQ=255)
                            cTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCCGCTGGTCGCAAACGGCCTGGCATGAACCTATCAATAATTTAGGTg                                                                        >  2:1974379/1‑100 (MQ=255)
                            cTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTg                                                                        <  1:2470741/100‑1 (MQ=255)
                            cTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTg                                                                        >  2:177458/1‑100 (MQ=255)
                            cTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTg                                                                        <  2:3128539/100‑1 (MQ=255)
                            cTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGCCTGGCATGAACCTATCCATAATTTAGGTg                                                                        >  1:165538/1‑100 (MQ=255)
                              gTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTggg                                                                      <  1:1082252/100‑1 (MQ=255)
                                ttGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGCCTGGCATGAACCTATCCATAATTTAGGTGGGAt                                                                    <  1:924817/100‑1 (MQ=255)
                                  gCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGCCTGGCATGAACCTATCCATAATTTAGGTGGGATGc                                                                  >  2:2079626/1‑100 (MQ=255)
                                   ccGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGcc                                                                 <  1:346467/100‑1 (MQ=255)
                                   ccGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGcc                                                                 >  2:2114970/1‑100 (MQ=255)
                                   ccGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGcc                                                                 >  2:2174547/1‑100 (MQ=255)
                                   ccGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGCCTGGCATGAACCTATCCATAATTTAGGTGGGATGcc                                                                 >  1:2788961/1‑100 (MQ=255)
                                   ccGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGCCTGGCATGAACCTATCCATAATTTAGGTGGGATGcc                                                                 >  1:2788542/1‑100 (MQ=255)
                                     gTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCg                                                               <  2:1538451/100‑1 (MQ=255)
                                     gTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGCCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCg                                                               >  2:697405/1‑100 (MQ=255)
                                       cTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGCCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGt                                                             >  1:2392280/1‑100 (MQ=255)
                                        tGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTg                                                            <  2:1006822/100‑1 (MQ=255)
                                         ggTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGtcc                                                           >  1:1528685/1‑98 (MQ=255)
                                         ggTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGc                                                           >  2:1556081/1‑100 (MQ=255)
                                         ggTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGCCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGc                                                           >  1:1730549/1‑100 (MQ=255)
                                            tGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTg                                                        >  1:1460338/1‑100 (MQ=255)
                                            tGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTCGTCGCAAACGGCCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCg                                                          >  1:567574/1‑98 (MQ=255)
                                             gATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTgg                                                       >  1:129237/1‑100 (MQ=255)
                                               tGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGCCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGtt                                                     <  2:905055/100‑1 (MQ=255)
                                               tGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGCCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGtt                                                     <  2:1244051/100‑1 (MQ=255)
                                                ggCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTg                                                    <  1:3171817/100‑1 (MQ=255)
                                                ggCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGCCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTg                                                    <  2:2392297/100‑1 (MQ=255)
                                                          cTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGc                                          <  2:50895/100‑1 (MQ=255)
                                                          cTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGc                                          >  2:507033/1‑100 (MQ=255)
                                                          cTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGc                                          >  2:3245011/1‑100 (MQ=255)
                                                          cTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGCCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGc                                          >  2:256501/1‑100 (MQ=255)
                                                           tAAGCTTTAACCCTTCGGACCCGAGCTGGTCGTAAACGGTCTCGCATGAACCTATCCATCATTTAGGTGGGATGCCCGGGGCGTGGTTGGCATAGACg                                           >  1:3170874/1‑98 (MQ=255)
                                                            aaGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTg                                        >  1:1751164/1‑100 (MQ=255)
                                                              gCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGCCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGtt                                      >  2:1653427/1‑100 (MQ=255)
                                                               cTTTAACCCTTCGGACCCCAGCTTGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGGGCGTGGTTGGCAGATACGCTGTTc                                     <  1:323961/100‑1 (MQ=255)
                                                               cTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTc                                     >  2:407234/1‑100 (MQ=255)
                                                               cTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTc                                     >  2:1303100/1‑100 (MQ=255)
                                                               cTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGCCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTc                                     >  2:2546121/1‑100 (MQ=255)
                                                                tttAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCt                                    <  1:2180498/100‑1 (MQ=255)
                                                                 ttAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAGCCTATCCATAATTTAGGTGGGATGCGCGGTGCGTGGTTGGCAGATACGCTGTTCtt                                   >  1:670568/1‑100 (MQ=255)
                                                                  tAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCttt                                  >  1:2583039/1‑100 (MQ=255)
                                                                    aCCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTAt                                >  1:2551050/1‑100 (MQ=255)
                                                                     cccTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTAtt                               >  1:1726076/1‑100 (MQ=255)
                                                                      ccTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAGTTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTAttt                              >  1:2660154/1‑100 (MQ=255)
                                                                      ccTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTAttt                              >  2:2999233/1‑100 (MQ=255)
                                                                      ccTTCGGACCCCAGCTGGTCGCAAACGGCCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTAttt                              >  2:2822064/1‑100 (MQ=255)
                                                                      ccTTCGGACCCCAGCTGGTCGCAAACGGCCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTAttt                              >  1:1010400/1‑100 (MQ=255)
                                                                        ttCGGACCCCAGCTGGTCGCAAACGGTCTTGCATGAACCTATCCATAATTTAGGTGGGATGTCCGGTGCGTGGTTGGCAGATACGCTCTTCTTTAttttt                            <  1:3176291/100‑1 (MQ=255)
                                                                         tCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTg                           <  2:150219/100‑1 (MQ=255)
                                                                         tCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTg                           <  2:150559/100‑1 (MQ=255)
                                                                          cGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTgg                          <  2:2540570/100‑1 (MQ=255)
                                                                          cGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTgg                          <  2:1366111/100‑1 (MQ=255)
                                                                            gACCCCAGCTGGTCGCAAACGGCCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCg                        <  2:2411348/100‑1 (MQ=255)
                                                                              ccccAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGt                       <  1:900447/99‑1 (MQ=255)
                                                                              ccccAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTg                      >  2:1911172/1‑100 (MQ=255)
                                                                              ccccAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTg                      <  1:541089/100‑1 (MQ=255)
                                                                                ccAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTg                      >  1:1635271/1‑98 (MQ=255)
                                                                                    cTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGATGGCt                >  2:853780/1‑100 (MQ=255)
                                                                                      ggTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGATGGCTTa              >  2:1078005/1‑100 (MQ=255)
                                                                                      ggTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGATGGCTTa              <  1:1471796/100‑1 (MQ=255)
                                                                                       gTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGATGGCTTac             <  2:1356896/100‑1 (MQ=255)
                                                                                       gTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGATGGCTTac             >  2:72220/1‑100 (MQ=255)
                                                                                       gTCGCAAACGGCCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGATGGCTTac             >  1:302219/1‑100 (MQ=255)
                                                                                             aaCGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGATGGCTTACACCAtt       >  1:1184622/1‑100 (MQ=255)
                                                                                             aaCGGCCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGATGGCTTACACCAtt       >  1:1410969/1‑100 (MQ=255)
                                                                                              aCGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGATGGCTTACACCATTc      <  2:348234/100‑1 (MQ=255)
                                                                                               cGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGATGGCTTACACCATTcc     <  1:885133/100‑1 (MQ=255)
                                                                                                  ccTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGATGGCTTACACCATTCCCGt  >  1:1774561/1‑100 (MQ=255)
                                                                                                  |                                                                                                   
TCTGGAGGCGTTGCTGATCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGCCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGATGGCTTACACCATTCCCGT  >  NC_000913/933289‑933486

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