Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,962,589 T→A 6.6% K356* (AAG→TAG)  fleS ← two‑component system sensor histidine kinase FleS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,962,5890TA 6.6% 178.5 / 3.0 61K356* (AAG→TAG) fleStwo‑component system sensor histidine kinase FleS
Reads supporting (aligned to +/- strand):  ref base T (23/34);  new base A (2/2);  total (25/36)
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.63e-01

ACCCCCTCCGCCAATTGCCCGTCGATCAACGGCAGCTCTACCCGCACGCAGGTACCACGGCCGACTTTCGAGCGCAGGTTGATCGTGCCCTGATGCGCACGCACCACCGCCTTGACCACGGCCAGGCCCAGGCCGGTACCGGTGGCCTTGGTGGTAAGGAACGGTTCGCCCAAACGACCGAGCAAGCCTGCATCGATGCCACTGCCGGCATCGCTGACGCACAGGTGCAAGCTACGTTCGCGGCGGTACAGGTGCACCTTCAGCCGCGCCGGGCCATCGCTGGCCTGCAGGGC  >  NC_002947/4962441‑4962733
                                                                                                                                                    |                                                                                                                                                
acCCCCTCCGCCAATTGCCCGTCGATCAACGGCAGCTCTACCCGCACGCAGGTACCACGGCCGACTTTCGAGCGCAGGTTGATCGTGCCCT‑ATGCGCACGCACCACCGCCTTGACCACGGCCAGGCCCAGGCCGGTACCGGTGGCCTtg                                                                                                                                                 >  2:108871/1‑149 (MQ=255)
  ccccTCCGCCAATTGCCCGTCGATCAACGGCAGCTCTACCCGCACGCAGGTACCACGGCCGACTTTCGAGCGCAGGTTGATCGTGCCCTGATGCGCACGCACCACCGCCTTGACCACGGCCAGGCCCAGGCCGGTACCGGTGGCCTtgg                                                                                                                                                <  1:217498/149‑1 (MQ=255)
   cccTCCGCCAATTGCCCGTCGATCAACGGCAGCTCTACCCGCACGCAGGTACCACGGCCGACTTTCGAGCGCAGGTTGATCGTGCCCTGATGCGCACGCACCACCGCCTTGACCACGGCCAGGCCCAGGCCGGTACCGGTGGCCTtggt                                                                                                                                               <  5:217873/149‑1 (MQ=255)
    ccTCCGCCAATTGCCCGTCGATCAACGGCAGCTCTACCCGCACGCAGGTACCACGGCCGACTTTCGAGCGCAGGTTGATCGTGCCCTGATGCGCACGCACCACCGCCTTGACCACGGCCAGGCCCAGGCCGGTACCGGTGGCCTtggtg                                                                                                                                              >  3:51801/1‑149 (MQ=255)
    ccTCCGCCAATTGCCCGTCGATCAACGGCAGCTCTACCCGCACGCAGGTACCACGGCCGACTTTCGAGCGCAGGTTGATCGTGCCCTGATGCGCACGCACCACCGCCTTGACCACGGCCAGGCCCAGGCCGGTACCGGTGGCCTtggtg                                                                                                                                              <  2:189744/149‑1 (MQ=255)
        cGCCAATTGCCCGTCGATCAACGGCAGCTCTACCCGCACGCAGGTACCACGGCCGACTTTCGAGCGCAGGTTGATCGTGCCCTGATGCGCACGCACCACCGCCTTGACCACGGCCAGGCCCAGGCCGGTACCGGTGGCCTTggtggt                                                                                                                                            >  3:311887/1‑147 (MQ=255)
        cGCCAATTGCCCGTCGATCAACGGCAGCTCTACCCGCACGCAGGTACCACGGCCGACTTTCGAGCGCAGGTTGATCGTGCCCTGATGCGCACGCACCACCGCCTTGACCACGGCCAGGCCCAGGCCGGTACCGGTGGCCTTGGTGGTaa                                                                                                                                          <  1:17488/149‑1 (MQ=255)
         gCCAATTGCCCGTCGATCAACGGCAGCTCTACCCGCACGCAGGTACCACGGCCGACTTTCGAGCGCAGGTTGATCGTGCCCTGATGCGCACGCACCACCGCCTTGACCACGGCCAGGCCCAGGCCGGTACCGGTGGCCTTGGTGGTAAg                                                                                                                                         <  6:325069/149‑1 (MQ=255)
            aaTTGCCCGTCGATCAACGGCAGCTCTACCCGCACGCAGGTACCACGGCCGACTTTCGAGCGCAGGTTGATCGTGCCCTGATGCGCACGCACCACCGCCTTGACCACGGCCAGGCCCAGGCCGGTACCGGTGGCCTTGGTGGTAAGGaa                                                                                                                                      >  5:292898/1‑149 (MQ=255)
             aTTGCCCGTCGATCAACGGCAGCTCTACCCGCACGCAGGTACCACGGCCGACTTTCGAGCGCATGTTGATCGTGCCCTGATGCGCACGCACCACCGCCTTGACCACGGCCAGGCCCAGGCCGGTACCGGTGGCCTTGGTGGTAAGGAAc                                                                                                                                     <  8:299581/149‑1 (MQ=255)
              ttGCCCGTCGATCAACGGCAGCTCTACCCGCACGCAGGTACCACGGCCGACTTTCGAGCGCAGGTTGATCGTGCCCTGATGCGCACGCACCACCGCCTTGACCACGGCCAGGCCCAGGCCGGTACCGGTGGCCTTGGTGGTAAGGAACg                                                                                                                                    <  6:292898/149‑1 (MQ=255)
                 cccGTCGATCAACGGCAGCTCTACCCGCACGCAGGTACCACGGCCGACTTTCGAGCGCAGGTTGATCGTGCCCTGATGCGCACGCACCACCGCCTTGACCACGGCCAGGCCCAGGCCGGTACCGGTGGCCTTGGTGGTAAGGAACGGtt                                                                                                                                 <  2:19077/149‑1 (MQ=255)
                  ccGTCGATCAACGGCAGCTCTACCCGCACGCAGGTACCACGGCCGACTTTCGAGCGCAGGTTGATCGTGCCCTGATGCGCACGCACCACCGCCTTGACCACGGCCAGGCCCAGGCCGGTACCGGTGGCCTTGGTGGTAAGGAACGGTTc                                                                                                                                <  2:64389/149‑1 (MQ=255)
                    gTCGATCAACGGCAGCTCTACCCGCACGCAGGTACCACGGCCGACTTTCGAGCGCAGGTTGATCGTGCCCTGATGCGCACGCACCACCGCCTTTACCACGGCCAGGCCCAGGCCGGTACCGGTGGCCTTGGTTGTAAGGAACGGTTCGc                                                                                                                              >  6:304670/1‑149 (MQ=255)
                     tCGATCAACGGCAGCTCTACCCGCACGCAGGTACCACGGCCGACTTTCGAGCGCAGGTTGATCGTGCCCTGATGCGCACGCACCACCGCCTTGACCACGGCCAGGCCCAGGCCGGTACCGGTGGCCTTGGTGGTAAGGAACGGTTCGcc                                                                                                                             >  7:260990/1‑149 (MQ=255)
                        aTCAACGGCAGCTCTACCCGCACGCAGGTACCACGGCCGACTTTCGAGCGCAGGTTGATCGTGCCCTGATGCGCACGCACCACCGCCTTGACCACGGCCAGGCCCAGGCCGGTACCGGTGGCCTTGGTGGTAAGGAACGGTTCGCCCa                                                                                                                           >  8:254631/1‑148 (MQ=255)
                        aTCAACGGCAGCTCTACCCGCACGCAGGTACCACGGCCGACTTTCGAGCGCAGGTTGATCGTGCCCT‑ATGCCCACGCACCACCGCCTTGACCACGGCCAGGCCCAGGCCGGTACCGGTGGCCTTGGTGGTAAGGAACGGTTCGCCCaaa                                                                                                                         <  1:108871/149‑1 (MQ=255)
                         tCAACGGCAGCTCTACCCGCACGCAGGTACCACGGCCGACTTTCGAGCGCAGGTTGATCGTGCCCTGATGCGCACGCACCACCGCCTTGACCACGGCCAGGCCCAGGCCGGTACCGGTGGCCTTGGTGGTAAGGAACGGTTCGCCCaaa                                                                                                                         >  1:33989/1‑149 (MQ=255)
                          catcGGCAGCTCTACCCGCACGGAGGTACCAGGGCCGACTTTCGAGCGCAGGTTGATCGTGCCCTTATGGGCACGCACCACCGCCTTGACCACGGCCAGGCCCAGGCCGGTACCGGTGGCCTTGGTGGTAAGGAACGGTTCGCCCaaa                                                                                                                         <  3:31954/145‑1 (MQ=255)
                                cAGCTCCACCCGCACGCAGGTACCACGGCCGACTTTCGAGCGCAGGTTGATCGTGCCCTGATGCGCACGCACCACCGCCTTGACCACGGCCAGGCCCAGGCCGGTACCGGTGGCCTTGGTGGTAAGGAACGGTTCGCCCAAACGACCGa                                                                                                                  <  1:159941/149‑1 (MQ=255)
                                  gCTCTACCCGCACGCAGGTACCACGGCCGACTTTCGAGCGCAGGTTGATCGTGCCCTGATGCGCACGCACCACCGCCTTGACCACGGCCAGGCCCAGGCCGGTACCGGTGGCCTTGGTGGTAAGGAACGGTTCGCCCAAACGACCGAg                                                                                                                 >  3:63768/1‑148 (MQ=255)
                                  gCTCTACCCGCACGCAGGTACCACGGCCGACTTTCGAGCGCAGGTTGATCGTGCCCTGATGCGCACGCACCACCGCCTTGACCACGGCCAGGCCCAGGCCGGTACCGGTGGCCTTGGTGGTAAGGAACGGTTCGCCCAAACGACCGAGc                                                                                                                <  4:311887/149‑1 (MQ=255)
                                   ctctACCCGCACGCAGGTACCACGGCCGACTTTCGAGCGCAGGTTGATCGTGCCCTGATGCGCACGCACCACCGCCTTGACCACGGCCAGGCCCAGGCCGGTACCGGTGGCCTTGGTGGTAAGGAACGGTTCGCCCAAACGACCGAGc                                                                                                                >  5:329649/1‑148 (MQ=255)
                                    tctACCCGCACGCAGGTACCACGGCCGACTTTCGAGCGCAGGTTGATCGTGCCCTGATGCGCACGCACCACCGCCTTGACCACGGCCAGGCCCAGGCCGGTACCGGTGGCCTTGGTGGTAAGGAACGGTTCGCCCAAACGACCGAGc                                                                                                                >  3:112185/1‑147 (MQ=255)
                                      tACCCGCACGCAGGTACCACGGCCGACTTTCGAGCGCAGGTTGATCGTGCCCTGATGCGCACGCACCACCGCCTTGACCACGGCCAGGCCCAGGCCGGTACCGGTGGCCTTGGTGGTAAGGAACGGTTCGCCCAAACGACCGAGCaa                                                                                                              <  1:163438/147‑1 (MQ=255)
                                        cccGCACGCAGGTACCACGGCCGACTTTCGAGCGCAGGTTGATCGTGCCCTGATGCGCACGCACCACCGCCTTGACCACGGCCAGGCCCAGGCCGGTACCGGTGGCCTTGGTGGTAAGGAACGGTTCGCCCAAACGACCGAGCAAGCCt                                                                                                          <  8:198914/149‑1 (MQ=255)
                                         ccGCACGCAGGTACCACGGCCGACTTTCGAGCGCAGGTTGATCGTGCCCTGATGCGCACGCACCACCGCCTTGACCACGGCCAGGCCCAGGCCGGTACCGGTGGCCTAGGTGGTAAGGAACGGTTCGCCCAAACGACCGAGCAAGCCTg                                                                                                         >  7:161147/1‑149 (MQ=255)
                                         ccGCACGCAGGTACCACGGCCGACTTTCGAGCGCAGGTTGATCGTGCCCTGATGCGCACGCACCACCGCCTTGACCACGGCCAGGCCCAGGCCGGTACCGGTGGCCTAGGTGGTAAGGAACGGTTCGCCCAAACGACCGAGCAAGCCTg                                                                                                         >  3:323391/1‑149 (MQ=255)
                                           gcacgcaGGTACCACGGCCGACTTTCGAGCGCAGGTTGATCGTGCCCTGATGCGCACGCACCACCGCCTTGACCACGGCCAGGCCCAGGCCGGTACCGGTGGCCTTGGTGGTAAGGAACGGTTCGCCCAAACGACCGAGCAAGCCTGCa                                                                                                       <  8:260990/149‑1 (MQ=255)
                                                caGGTACCACGGCCGACTTTCGAGCGCAGGTTGATCGTGCCCTGATGCGCACGCACCACCGCCTTGACCACGGCCAGGCCCAGGCCGGTACCGGTGGCCTTGGTGGTAAGGAACGGTTCGCCCAAACGACCGAGCAAGCCTGCATCGAt                                                                                                  <  4:112185/149‑1 (MQ=255)
                                                   gTACCACGGCCGACTTTCGAGCGCAGGTTGATCGTTCCCTGATGCGCACGCACCACCGCCTTGACCACGGCCAGGCCCAGGCCGGTACCGGGGGCCTTGGTGGTAAGGAACGGTTCGCCCAAACGACCGAGCAAGCCTGCATCGATGcc                                                                                               <  2:267210/149‑1 (MQ=255)
                                                     aCCACGGCCGACTTTCGAGCGCAGGTTGATCGTGCCCTGATGCGCACGCACCACCGCCTTGACCACGGCCAGGCCCAGGCCGGTACCGGTGGCCTTggtggt                                                                                                                                            <  8:154335/102‑1 (MQ=255)
                                                         cGGCCGACTTTCGAGCGCAGGTTGATCGTGCCCTGATGCGCACGCACCACCGCCTTGACCACGGCCAGGCCCAGGCCGGTACCGGTGGCCTTGGTGGTAAGGAACGGTTCGCCCAAACGACCGAGCAAGCCt                                                                                                          >  1:44458/1‑132 (MQ=255)
                                                         cGGCCGACTTTCGAGCGCAGGTTGATCGTGCCCTGATGCGCACGCACCACCGCCTTGACCACGGCCAGGCCCAGGCCGGTACCGGTGGCCTTGGTGGTAAGGAACGGTTCGCCCAAACGACCGAGCAAGCCt                                                                                                          <  2:44458/132‑1 (MQ=255)
                                                          ggCCGACTTTCGAGCGCAGGTTGATCGTGCCCTGATGCGCACGCACCACCGCCTTGACCACGGCCAGGCCCAGGCCGGTACCGGTGGCCTTGGTGGTAAGGAACGGTTCGCCCAAACGACCGAGCAAGCCTGCATCGATGCCACTGCCg                                                                                        <  4:51801/149‑1 (MQ=255)
                                                            ccGACTTTCGAGCGCAGGTTGATCGTGCCCTGATGCGCACGCACCACCGCCTTGACCACGGCCAGGCCCAGGCCGGTACCGGTGGCCTTGGTGGTAAGGAACGGTTCGCCCAAACGACCGAGCAAGCCTGCATCGATGCCACTGCCgg                                                                                       >  8:161752/1‑148 (MQ=255)
                                                            ccGACTTTCGAGCGCAGGTTGATCGTGCCCTGATGCGCACGCACCACCGCCTTGACCACGGCCAGGCCCAGGCCGGTACCGGTGGCCTTGGTGGTAAGGAACGGTTCGCCCAAACGACCGAGCAAGCCTGCATCGATGCCACTGCCGGc                                                                                      <  7:306238/149‑1 (MQ=255)
                                                              gACTTTCGAGCGCAGGTTGATCGTGCCCTGATGTGCACGCACCACCGCCTTGACCACGGCCAGGCCCAGGCCGGTACCGGTGGCCTTGGTGGTAAGGAACGGTTCGCCCAAACGACCGAGCAAGCCTGCATCGATGCCACTGCCGGCAt                                                                                    <  7:100238/149‑1 (MQ=255)
                                                                cTTTCGAGCGCAGGTTGATCGTGCCCTGATGCGCACGCACCACCGCCTTGACCACGGCCAGGCCCAGGCCGGTACCGGTGGCCTTGGTGGTAAGGAACGGTTCGCCCAAACGACCGAGCAAGCCTGCATCGATGCCACTGCCGGCATCg                                                                                  >  8:244316/1‑149 (MQ=255)
                                                                         gcAGGTTGATCGTGCCCTGATGCGCACGCACCCCCGCCTTGACCACGGCCAGGCCCAGGCCGGTACCGGTGGCCTTGGTGGTAAGGAACGGTTCGCCCAAACGACCGAGCAAGCCTGCATCGATGCCACTGCCGGCATCGCTGACGca                                                                          >  4:250215/1‑148 (MQ=255)
                                                                           aGGTTGATCGTGCCCTGATGCGCACGCACCACCGCCTTGACCACGGCCAGGCCCAGGCCGGTACCGGTGGCCTTGGTGGTAAGGAACGGTTCGCCCAAACGACCGAGCAAGCCTGCATCGATGCCACTGCCGGCATCGCTGACGCACAg                                                                       <  5:59317/149‑1 (MQ=255)
                                                                                  tCGTGCCCTGATTCGCACTCACCACCCCCTTGACCACGGCCATGCCCACGCCGGTACCGGTGGCCTTGGTGGTAAGGAACGGATCCCCCAAACGACCGAGCAAGCCTGCATCGATGCCACTGCCGGCATCGCTGACGCACAGGAGCAAg                                                                <  5:304670/149‑1 (MQ=255)
                                                                                   cGTGCCCTGATGCGCACGCACCACCGCCTTGACCACGGCCAGGCCCAGGCCGGTACCGGTGGCCTTGGTGGTAAGGAACGGTTCGCCCAAACGACCGAGCAAGCCTGCATCGAt                                                                                                  <  4:118119/114‑1 (MQ=255)
                                                                                   cGTGCCCTGATGCGCACGCACCACCGCCTTGACCACGGCCAGGCCCAGGCCGGTACCGGTGGCCTTGGTGGTAAGGAACGGTTCGCCCAAACGACCGAGCAAGCCTGCATCGAt                                                                                                  >  3:118119/1‑114 (MQ=255)
                                                                                       cccTGATGCGCACGCACCACCGCCTTGACCACGGCCAGGCCCAGGCCGGTACCGGTGGCCTTGGTGGTAAGGAACGGTTCGCCCAAACGACCGAGCAAGCCTGCATCGATGGCACTGCCGGCATCGCTGACGCACAGGTGCAAGCTACg                                                           >  4:198884/1‑149 (MQ=255)
                                                                                          tGATGCGCACGCACCACCGCCTTGACCACGGCCAGGCCCAGGCCGGTACCGGTGGCCTTGGTGGT‑AGGAACGGTTCGCCCAAACGACCGAGCAAGCCTGCATCGATGCCACTGCCGCCATCGCTGACGCACAGGTGCAAGCTACGTTcg                                                       >  8:305420/1‑149 (MQ=255)
                                                                                               cgcacgcaCCACCGCCTTGACCACGGCCAGGCCCAGGCCGGTACCGGTGGCCTTGGTGGTAAGGAACGGTTCGCCCAAACGACCGAGCAAGCCTGCATCGATGCCACTGCCGGCATCGCTGACGCACAGGTGCAAGCTACGTTCgcgg                                                    <  7:161752/148‑1 (MQ=255)
                                                                                                gcacgcacCACCGCCTTGACCACGGCCAGGCCCAGGCCGGTACCGGTGGCCTTGGTGGTAAGGAACGGTTCGCCCAAACGACCGAGCAAGCCTGCATCGATGCCACTGCCGGCATCGCTGACGCACAGGTGCAAGCTACGTTCgcggcg                                                  >  1:91307/1‑149 (MQ=255)
                                                                                                gcacgcacCACCGCCTTGACCACGGCCAGGCCCAGGCCGGTACCGGTGGCCTTGGTGGTAAGGAACGGTTCGCCCAAACGACCGAGCAAGCCTGCATCGATGCCACTGCCGGCATCGCTGACGCACAGGTGCAAGCTACGTTCgcggcg                                                  >  7:34085/1‑149 (MQ=255)
                                                                                                    gcacCACCGCCTTGACCACGGCCAGGCCCAGGCCGGTACCGGTGGCCTTGGTGGTAAGGAACGGTTCGCCCAAACGACCGAGCAAGCCTGCATCGATGCCACTGCCGGCATCGCTGACGCACAGGTGCAAGCTACGTTCGCGGCGGTa                                               >  4:94420/1‑148 (MQ=255)
                                                                                                      accaccGCCTTGACCACGGCCAGGCCCAGGCCGGTACCGGTGGCCTTGGTGGTAAGGAACGGTTCGCCCAAACGACCGAGCAAGCCTGCATCGATGGCACTGCCGGCATCGCTGACGCACAGGTGCAAGCTACGTTCGCGGCGGTACAg                                            <  3:198884/149‑1 (MQ=255)
                                                                                                      accaccGCCTTGACCACGGCCAGGCCCAGGCCGGTACCGGTGGCCTTGGTGGTAAGGAACGGTTCGCCCAAACGACCGAGCAAGCCTGCATCGATGCCACTGCCGGCATCGCTGACGCACAGGTGCAAGCTACGTTCgcggcg                                                  >  6:39682/1‑143 (MQ=255)
                                                                                                      accaccGCCTTGACCACGGCCAGGCCCAGGCCGGTACCGGTGGCCTTGGTGGTAAGGAACGGTTCGCCCAAACGACCGAGCAAGCCTGCATCGATGCCACTGCCGGCATCGCTGACGCACAGGTGCAAGCTACGTTCgcggcg                                                  <  5:39682/143‑1 (MQ=255)
                                                                                                      accaccGCCTTGACCACGGCCAGGCCCAGGCCGGTACCGGTGGCCTTGGTGGTAAGGAACGGTTCGCCCAAACGACCGAGCAAGCCTGCATCGATGCCACTGCCGGCATCGCTGACGCACAGGTGCAAGCTACGTTCGCGGCGGTACAg                                            >  2:9806/1‑149 (MQ=255)
                                                                                                       ccaccGCCTTGACCACGGCCAGGCCCAGGCCGGTACCGGTGGCCTTGGTGGTAAGGAACGGTTCGCCCAAACGACCGAGCAAGCCTGCATCGATGCCACTGCCGGCATCGCTGACGCACAGGTGCAAGCTACGTTCGCGGCGGTACAgg                                           <  2:334446/149‑1 (MQ=255)
                                                                                                         accGCCTTGACCACGGCCAGGCCCAGGCCGGTACCGGTGGCCTTGGTGGTAAGGAACGGTTCGCCCAAACGACCGAGCAAGCCTGCATCGATGCCACTGCCGGCATCGCTGACGCACAGGTGCAAGCTACGTTCGCGGCGGTACAGGTg                                         <  6:329649/149‑1 (MQ=255)
                                                                                                            gCCTTGACCACGGCCAGGCCCAGGCCGGTACCGGTGGCCTTGGTGGTAAGGAACGGTTCGCCCAAACGACCGAGCAAGCCTGCATCGATGCCACTGCCGGCATCGCTGACGCACAGGTGCAAGCTACGTTCGCGGCGGTACAGGTGCAc                                      <  8:34085/149‑1 (MQ=255)
                                                                                                                cggtaccGGCCTGGGCCTGGCCGGTACCGGTGGCCTTGGTGGTAAGGAACGGTTCGCCCAAACGACCGAGCAAGCCTGCATCGATGCCACTGCCGGCATCGCTGACGCACAGGTGCAAGCTACGTTCGCGGCGGTACAGGTGCACCtt                                   >  4:116518/7‑148 (MQ=255)
                                                                                                                    cACGGCCAGGCCCAGGCCGGTACCGGTGGCCTTGGTGGTAAGGAACGGTTCGCCCAAACGACCGAGCAAGCCTGCATCGATGCCACTGCCGGCATCGCTGACGCACAGGTGCAAGCTACGTTCGCGGCGGTACAGGTGCACCTTCAGcc                              >  7:304795/1‑149 (MQ=255)
                                                                                                                                  ggccGGTACCGGTGGCCTTGGTGGTAAGGAACGGTTCGCCCAAACGACCGAGCAAGCCTGCATCGATGCCACTGCCGGCATCGCTGACGCACAGGTGCAAGCTACGTTCGCGGCGGTACAGGTGCACCTTCAGCCGCGCCGGGCCATCg                <  7:254631/149‑1 (MQ=255)
                                                                                                                                   gccGGTACCGGTGGCCTTGGTGGTAAGGAACGGTTCGCCCAAACGACCGAGCAAGCCTGCATCGATGCCACTGCCGGCATCGCTGACGCACAGGTGCAAGCTACGTTCGCGGCGGTACAGGTGCACCTTCAGCCGCGCCGGGCCATc                 <  3:116518/147‑1 (MQ=255)
                                                                                                                                    ccGGTACCGGTGGCCTAGGTGGTAAGGAACGGTTCGCCCAAACGACCGAGCAAGCCTGCATCGATGCCACTGCCGGCATCGCTGACGCACAGGTGCAAGCTACGTTCGCGGCGGTACAGGTGCACCTTCAGCCGCGCCGGGCCATCGCt              <  8:161147/149‑1 (MQ=255)
                                                                                                                                    ccGGTACCGGTGGCCTAGGTGGTAAGGAACGGTTCGCCCAAACGACCGAGCAAGCCTGCATCGATGCCACTGCCGGCATCGCTGACGCACAGGTGCAAGCTACGTTCGCGGCGGTACAGGTGCACCTTCAGCCGCGCCGGGCCATCGCt              <  4:323391/149‑1 (MQ=255)
                                                                                                                                     cGGTACCGGTGGCCTTGGTGGTAAGGAACGGTTCGCCCAAACGACCGAGCAAGCCTGCATCGATGCCACTGCCGGCATCGCTGACGCACAGGTGCAAGCTACGTTCGCGGCGGTACAGGTGCACCTTCAGCCGCGCCGGGCCATCGCTg             <  2:234859/149‑1 (MQ=255)
                                                                                                                                          ccGGTGGCCTTGGTGGTAAGGAACGGTTCGCCCAAACGACCGAGCAAGCCTGCATCGATGCCACTGCCGGCATCGCTGACGCACAGGTGCAAGCTACGTTCGCGGCGGTACAGGTGCACCTTCAGCCGCGCCGGGCCATCGCTGGCCTg        <  3:250215/149‑1 (MQ=255)
                                                                                                                                                gCCTTGGTGGTAAGGAACGGTTCGCCCAAACGACCGAGCAAGCCTGCATCGATGCCACTGCCGGCATCGCTGACGCACAGGTGCAAGCTACGTTCGCGGCGGTACAGGTGCACCTTCAGCCGCGCCGGGCCATCGCTGGCCTGCAGGGc  <  4:63768/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
ACCCCCTCCGCCAATTGCCCGTCGATCAACGGCAGCTCTACCCGCACGCAGGTACCACGGCCGACTTTCGAGCGCAGGTTGATCGTGCCCTGATGCGCACGCACCACCGCCTTGACCACGGCCAGGCCCAGGCCGGTACCGGTGGCCTTGGTGGTAAGGAACGGTTCGCCCAAACGACCGAGCAAGCCTGCATCGATGCCACTGCCGGCATCGCTGACGCACAGGTGCAAGCTACGTTCGCGGCGGTACAGGTGCACCTTCAGCCGCGCCGGGCCATCGCTGGCCTGCAGGGC  >  NC_002947/4962441‑4962733

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