Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,861,172 A→G 7.2% intergenic (‑78/+159) ligA ← / ← PP_4275 DNA ligase/cell division protein ZipA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,861,1720AG 7.2% 406.6 / 13.6 139intergenic (‑78/+159)ligA/PP_4275DNA ligase/cell division protein ZipA
Reads supporting (aligned to +/- strand):  ref base A (59/70);  new base G (5/5);  total (64/75)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.31e-01

Only 100 of 142 total aligned reads displayed.

GAGCACGTAGTAGCGATAATTGTGCTGATCGAGCTCGGCACGCAGTTCGAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGT  >  NC_002947/4861023‑4861321
                                                                                                                                                     |                                                                                                                                                     
gAGCACGTAGTAGCGATAATTGTGCTGATCGAGCTCGGCACGCAGTTCGAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCa                                                                                                                                                       <  2:2259440/150‑1 (MQ=255)
  gCACGTAGTAGCGATAATTGTGCTGATCGAGCTGGGCACGCAGTTCGAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGTGTTCCAGGTGTCTGTTGACTGTGCCAGC                                                                                                                                                     <  2:213416‑M1/150‑1 (MQ=255)
  gCACGTAGTAGCGATAATTGTGCTGATCGAGCTCGGCACGCAGTTCGAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGC                                                                                                                                                     >  1:865194‑M1/1‑150 (MQ=255)
          tagCGATAATTGTGCTGATCGAGCTCGGCACGCAGTTCGAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGC                                                                                                                                             <  1:2078771‑M1/150‑1 (MQ=255)
          tagCGATAATTGTGCTGATCGAGCTCGGCACGCAGTTCGAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCAAGCCTCTTCGC                                                                                                                                             >  1:960630‑M1/1‑150 (MQ=255)
            gCGATAATTGTGCTGATCGAGCTCGGCACGCAGTTCGAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGG                                                                                                                                           <  1:1036990‑M1/150‑1 (MQ=255)
            gCGATAATTGTGCTGATCGAGCTCGGCACGCAGTTCGAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGG                                                                                                                                           <  2:960630‑M1/150‑1 (MQ=255)
            gCGATAATTGTGCTGATCGAGCTCGGCACGCAGTTCGAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGG                                                                                                                                           <  1:2220741‑M1/150‑1 (MQ=255)
            gCGATAATTGTGCTGATCGAGCTCGGCACGCAGTTCGAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGG                                                                                                                                           <  2:156304‑M1/150‑1 (MQ=255)
              gATAATTGTGCTGATCGAGCTCGGCACGCAGTTCGAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGT                                                                                                                                         >  1:2312456‑M1/1‑150 (MQ=255)
                  aTTGTGCTGATCGAGCTCGGCACGCAGTTCGAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAAC                                                                                                                                     <  2:1705710‑M1/150‑1 (MQ=255)
                   ttGTGCTGATCGAGCTCGGCACGCAGTTCGAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACC                                                                                                                                    >  2:316770‑M1/1‑150 (MQ=255)
                    tgtgCTGATCGAGCTCGGCACGCAGTTCGAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCC                                                                                                                                   <  2:2797637‑M1/150‑1 (MQ=255)
                          gATCGAGCTCGGTCCGACGTTCGCGGTTACTAGTTGCGGCTGTCATGGTATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTccc                                                                                                                             <  2:2959614/150‑1 (MQ=255)
                          gATCGAGCTCGGCACGCAGTTCGAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTccc                                                                                                                             >  1:2615589/1‑150 (MQ=255)
                           aTCGAGCTCGGCACGCAGTTCGAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCca                                                                                                                            >  1:1043134/1‑150 (MQ=255)
                               aGCTCGGCACGCAGTTCGAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCtacagg                                                                                                                        >  1:2120157‑M1/1‑144 (MQ=255)
                                 cTCGGCACGCAGTTCGAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGTGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAggg                                                                                                                      <  2:2487919/150‑1 (MQ=255)
                                  tCGGCACGCAGTTCGAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGc                                                                                                                     <  2:2558007/150‑1 (MQ=255)
                                  tCGGCACGCAGTTCGAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGc                                                                                                                     <  1:2815898/150‑1 (MQ=255)
                                   cGGCACGCAGTTCGAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCg                                                                                                                    <  2:1043134/150‑1 (MQ=255)
                                    ggCACGCAGTTCGAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGt                                                                                                                   <  1:316770/150‑1 (MQ=255)
                                       aCGCAGTTCGAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCtacaggtacatcac                                                                                                                >  1:175967‑M1/1‑136 (MQ=255)
                                       aCGCAGTTCGAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTcac                                                                                                                >  2:1767653/1‑150 (MQ=255)
                                       aCGCAGTTCGAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTcac                                                                                                                <  1:1812215/150‑1 (MQ=255)
                                         gCAGTTCGAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTTCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCtacaggtacatcaccg                                                                                                              >  1:1979218‑M1/1‑134 (MQ=255)
                                         gCAGTTCGAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGCCTGTGCCGGCCTCTTCGCGGG                                                                                                                                          <  1:1508564‑M1/122‑1 (MQ=255)
                                                gAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAg                                                                                                       >  1:1440118/1‑150 (MQ=255)
                                                  ggATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAggg                                                                                                     >  2:568032/1‑150 (MQ=255)
                                                  ggATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAggg                                                                                                     <  2:963893/150‑1 (MQ=255)
                                                     tACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACAt                                                                                                              <  2:240405/138‑1 (MQ=255)
                                                     tACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACAt                                                                                                              >  1:240405/1‑138 (MQ=255)
                                                     tACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACC                                                                                                                                    >  2:532104‑M1/1‑116 (MQ=255)
                                                     tACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACC                                                                                                                                    <  1:532104‑M1/116‑1 (MQ=255)
                                                       cGGGTTTCGGCGGTCATGGGATTTGTCTTATCTTGTAGCGGAGAAGTTGTGCTTGGTCTGGTGCGGGGTGTTCCATGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCAGCTCCCACACGGGCGTCACATGCCCTAGGGCGGTg                                                                                                <  2:1440118/150‑1 (MQ=255)
                                                        gggTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGa                                                                                               <  2:2615589/150‑1 (MQ=255)
                                                        gggTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGa                                                                                               >  1:1904470/1‑150 (MQ=255)
                                                        gggTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGa                                                                                               >  2:1522830/1‑150 (MQ=255)
                                                            ttCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATgag                                                                                           >  1:658661/1‑148 (MQ=255)
                                                               ggcggTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAg                                                                                                                        <  1:1685242/118‑1 (MQ=255)
                                                               ggcggTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCt                                                                                        <  2:2683254/150‑1 (MQ=255)
                                                                  ggTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTAc                                                                                          <  2:899649/145‑1 (MQ=255)
                                                                   gTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAg                                                                                    >  2:1094976‑M2/1‑149 (MQ=255)
                                                                    tCATGGGTTTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTa                                                                                     <  1:1484935/148‑1 (MQ=255)
                                                                    tCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTa                                                                                     >  2:1484935/1‑148 (MQ=255)
                                                                      aTGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGa                                                                                               >  2:2221101/1‑136 (MQ=255)
                                                                      aTGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGa                                                                                               <  1:2221101/136‑1 (MQ=255)
                                                                       tGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACaa                                                                                                                         >  1:2642466/1‑109 (MQ=255)
                                                                       tGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACaa                                                                                                                         <  2:2642466/109‑1 (MQ=255)
                                                                       tGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGc                                                                                                                     >  1:949615/1‑113 (MQ=255)
                                                                         tgaTTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCtacaggtacatcaccgccctagggcatgtgacgcccttgtgggagcgg                                                                              <  2:192258‑M1/149‑49 (MQ=255)
                                                                         ggATTTGTCTTCTCTTTCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCtacaggtacatcaccgccctagggcatgtgacgcccttgtgggagcgg                                                                              <  2:1979218‑M1/150‑49 (MQ=255)
                                                                         ggATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCGGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAggagcgg                                                                              >  2:1342499‑M2/1‑143 (MQ=255)
                                                                         ggATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGAGCGCCGCTCGCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCtacaggtacatcaccgccctagggcatgtgacgcccttgtgggagcgg                                                                              <  2:81898‑M1/150‑49 (MQ=255)
                                                                          gATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCtacaggtacatcaccgccctagggcatgtgacgcccttgtgggagcggg                                                                             <  1:1685629‑M1/150‑50 (MQ=255)
                                                                          gATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAggagcggg                                                                             >  1:1767653‑M2/1‑142 (MQ=255)
                                                                           aTTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCtacaggtacatcaccgccctagggcatgtgacgcccttgtgggagc                                                                                <  2:2402428‑M1/146‑47 (MQ=255)
                                                                           aTTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCtacaggtacatcaccgccctagggcatgtgacgcccttgtgggagc                                                                                >  1:2402428‑M1/1‑100 (MQ=255)
                                                                            tttGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAggagcgggtt                                                                           >  1:1522830‑M2/1‑140 (MQ=255)
                                                                                tcttctCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTACGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCTCAGTACGGCGTCACATGCCCTAGGGCGGCGATGTACCTGTAggagcgggtttacc                                                                       <  1:2490792‑M2/150‑15 (MQ=255)
                                                                                 cttctCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAggagcgggtttaccc                                                                      >  2:2294965‑M2/1‑135 (MQ=255)
                                                                                   tctcTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCtacaggtacatcaccgccctagggcatgtgacgcccttgtgggagcggg                                                                             >  2:450337‑M1/1‑92 (MQ=255)
                                                                                   tctcTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCtacaggtacatcaccgccctagggcatgtgacgcccttgtgggagcggg                                                                             <  1:450337‑M1/141‑50 (MQ=255)
                                                                                   tctcTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACCTGCCTTAGGGCGGGGATGTACCTGTAggagcgggtttacccgc                                                                    <  2:337608‑M2/150‑18 (MQ=255)
                                                                                      ctTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAg                                                                                                                        >  2:1719531/1‑95 (MQ=255)
                                                                                      ctTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAggagcgggtttacccgcgaa                                                                 >  2:1904470‑M2/1‑130 (MQ=255)
                                                                                         gcaaagcaaaAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGt                                                                                      <  2:1796403/126‑1 (MQ=255)
                                                                                         gcaaagcaaaAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGt                                                                                      >  1:1796403/1‑126 (MQ=255)
                                                                                         gcaaagcaaaAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGC                                                                                                                                             <  2:94164‑M1/71‑1 (MQ=255)
                                                                                         gcaaagcaaaAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGC                                                                                                                                             >  1:94164‑M1/1‑71 (MQ=255)
                                                                                           aaagcaaaAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCtacaggtacatcaccgccctagggcatgtgacgcccttgtgggagcgggtttacccgcgaagaggc                                                            <  2:2312456‑M1/150‑67 (MQ=255)
                                                                                                aaaaGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCtacaggtacatcaccgccctagggcatgtgacgcccttgtgggagcgggt                                                                            >  1:1121973‑M1/1‑79 (MQ=255)
                                                                                                aaaaGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCtacaggtacatcaccgccctagggcatgtgacgcccttgtgggagcgggt                                                                            <  2:1121973‑M1/129‑51 (MQ=255)
                                                                                                aaaaGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAggagcggg                                                                             <  1:1157918‑M2/128‑9 (MQ=255)
                                                                                                aaaaGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAggagcgg                                                                              <  1:1431140‑M2/127‑8 (MQ=255)
                                                                                                 aaaGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCt                                                                                        <  2:1970505/116‑1 (MQ=255)
                                                                                                 aaaGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCt                                                                                        >  1:1970505/1‑116 (MQ=255)
                                                                                                 aaaGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGcac                                                      >  1:872332/1‑150 (MQ=255)
                                                                                                   agagCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTGGACCGTGCCAGCCTCGGCGCGTGTAAACCCGCCGCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTACACGCGGGTTTACCCGCGAAGAGGCCGGCACAg                                                    <  2:2490792/150‑1 (MQ=255)
                                                                                                   agagCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGGCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAg                                                    <  1:413198/150‑1 (MQ=255)
                                                                                                     agcagcTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCtacaggtacatcaccgccctagggcatgtgacgcccttgtgggagcgggtttacccgcgaagaggctggcacagtc                                                  <  2:2120157‑M1/150‑77 (MQ=255)
                                                                                                          cTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCg                                                                                                   <  2:122494/96‑1 (MQ=255)
                                                                                                               ggtggtGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCtacaggtacatcaccgccctagggcatgtgacgcccttgtgggagcggg                                                                             >  2:774068‑M1/1‑64 (MQ=255)
                                                                                                                 tggtgCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATAcc                                      >  2:2796063/1‑150 (MQ=255)
                                                                                                                         ggggTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCTGTGATGTa                                                                                           <  2:1528083/89‑1 (MQ=255)
                                                                                                                          gggTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTTATCCCGCCCCCACAAGGGCGGCACATGCCCTAGGGCGGTAATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCTGGCACAGGCAACGACGATACCCGAACCGGc                             >  1:2651765/1‑150 (MQ=38)
                                                                                                                          gggTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGc                             <  2:2434291/150‑1 (MQ=255)
                                                                                                                            gTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCAGCTCCca                                                                                                                            <  1:2587360/53‑1 (MQ=255)
                                                                                                                            gTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCAGCTCCca                                                                                                                            >  2:2587360/1‑53 (MQ=255)
                                                                                                                            gTTCCAGGTGTCTGTGGACTGTGCCAGCCTCTTCCCGGGTAGACCCGCTCCTACAAGGGCTTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTCCCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCa                           <  2:1935951/150‑1 (MQ=38)
                                                                                                                             ttCCAGGTGTCTGTTGACTGTGCCAGCCTATTCGACGGGAGAGCGACTCGCTCGAGGGCTACACATGCCCTAGCGAGCTAAAGGCGCTGCATGGTCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACAAAGATACCCGAACCGGCCAg                          >  1:2994793/1‑150 (MQ=12)
                                                                                                                              tCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAGCCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGc                         <  2:1313507/150‑1 (MQ=255)
                                                                                                                                        tGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAg               >  1:2453148/1‑150 (MQ=255)
                                                                                                                                         gTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCtacaggtacatcaccgccctagggcatgtgacgcccttgtgggagcgggtttacccgcgaa                                                                 >  1:1960740‑M1/1‑38 (MQ=255)
                                                                                                                                          ttGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTAc                                                                                          >  1:65654/1‑73 (MQ=255)
                                                                                                                                          ttGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTAc                                                                                          <  2:65654/73‑1 (MQ=255)
                                                                                                                                           tGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCAc                       <  1:1058149/139‑1 (MQ=255)
                                                                                                                                           tGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCAc                       >  2:1058149/1‑139 (MQ=255)
                                                                                                                                            gACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAggagcgggtttacccgcgaagaggccggcacaggcaacagacacctggaaccccacgcaccaccccaagctgct           <  2:2120157‑M2/150‑75 (MQ=255)
                                                                                                                                                     aGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGt  >  2:424897/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GAGCACGTAGTAGCGATAATTGTGCTGATCGAGCTCGGCACGCAGTTCGAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGT  >  NC_002947/4861023‑4861321

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: