Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,971,955 C→G 7.3% intergenic (‑390/+27) PP_4379 ← / ← flgL 3‑oxoacyl‑ACP synthase/flagellar hook‑associated protein FlgL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,971,9550CG 7.3% 217.3 / 5.9 68intergenic (‑390/+27)PP_4379/flgL3‑oxoacyl‑ACP synthase/flagellar hook‑associated protein FlgL
Reads supporting (aligned to +/- strand):  ref base C (30/33);  new base G (3/2);  total (33/35)
Fisher's exact test for biased strand distribution p-value = 6.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.09e-01

GCACTGGACCCCGAATGATCAGCTTCGCAGCAGATATCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAAGCATGGCCTGCTGCAGGGTCAGCTGAATGGTGGCCGAGGCCATGTCAGTATCGCCAATGG  >  NC_002947/4971809‑4972103
                                                                                                                                                    |                                                                                                                                                     
gcACTGGACCCCGAATGATCAGCTTCGCAGCAGATATCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCgcag                                                                                                                                                     <  3:250134/149‑1 (MQ=255)
gcACTGGACCCCGAATGATCAGCTTCGCAGCAGATATCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCgcag                                                                                                                                                     <  6:44202/149‑1 (MQ=255)
  aCTGGACCCCGAATGATCAGCTTCGCAGCAGATATCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCgcagca                                                                                                                                                   <  4:68055/149‑1 (MQ=255)
   cTGGACCCCGAATGATCAGCTTCGCAGCAGATATCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGcagcag                                                                                                                                                  >  6:342469/1‑149 (MQ=255)
      gACCCCGAATGATCAGCTTCGCAGCAGATATCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGccc                                                                                                                                               >  5:263290/1‑149 (MQ=255)
        ccccGAATGATCAGCTTCGCAGCAGATATCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGCTGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGcccc                                                                                                                                              <  1:96771/149‑1 (MQ=255)
        ccccGAATGATCAGCTTCGCAGCAGATATCGATGGGGCCGCATCACGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCa                                                                                                                                             <  6:273511/149‑1 (MQ=255)
            gAATGATCAGCTTCGCAGCAGATATCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGa                                                                                                                                          >  6:292111/1‑148 (MQ=255)
            gAATGATCAGCTTCGCAGCAGATATCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGAt                                                                                                                                         >  4:353302/1‑149 (MQ=255)
            gAATGATCAGCTTCGCAGCAGATATCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCGTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGAt                                                                                                                                         >  2:325139/1‑149 (MQ=255)
            gAATGATCAGCTTCGCAGCAGATATCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACAAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCTGCGCAGccctatcgc                                                                                                                                         <  2:131393/149‑7 (MQ=38)
              aTGATCAGCTTCGCAGCAGATATCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGATca                                                                                                                                       <  5:190569/149‑1 (MQ=255)
              aTGATCAGCTTCGCAGCAGATATCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGATca                                                                                                                                       <  3:353302/149‑1 (MQ=255)
              aTGATCACCTTCGCAGCAGATATCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGATca                                                                                                                                       <  6:167451/149‑1 (MQ=255)
                  tCAGCTTCGCAGCAGATATCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGATca                                                                                                                                       <  5:211722/145‑1 (MQ=255)
                  tCAGCTTCGCAGCAGATATCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGATca                                                                                                                                       >  6:211722/1‑145 (MQ=255)
                    aGCTTCGCAGCAGATATCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGATcaggcag                                                                                                                                  <  3:78390/148‑1 (MQ=255)
                     gCTTCGCAGCAGATATCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGATCaggcaggc                                                                                                                                <  1:226145/149‑1 (MQ=255)
                     gCTTCGCAGCAGATATCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCAACCACCAGAAGACAGGACGACATGTGGGGGCGTTGTGTCCCGATAGGCCTGCGCAGCAGCCCAAGGATaaggcagc                                                                                                                                 >  8:36619/1‑147 (MQ=37)
                          gcagcaGATATCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGAt                                                                                                                            >  6:8752/1‑148 (MQ=255)
                          gcagcaGATATCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGAt                                                                                                                            <  4:92744/148‑1 (MQ=255)
                           cagcagATATCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGATTa                                                                                                                          <  3:57736/149‑1 (MQ=255)
                              cagATATCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGATTAGAt                                                                                                                       <  4:119325/149‑1 (MQ=255)
                               agATTTCGATGGGGCCGCATCGGGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCGTTGTGTCGCTATAGGGCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGATTAGATc                                                                                                                      <  5:33978/149‑1 (MQ=255)
                                      gATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCGACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTc                                                                                                                 >  1:262970/1‑148 (MQ=255)
                                          gggCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGGAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAg                                                                                                          >  5:322918/1‑149 (MQ=255)
                                           ggCCGCATGGCGGGCCGTTCGCGACGCAATGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAgg                                                                                                         <  7:36619/149‑1 (MQ=255)
                                               gCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCt                                                                                                       >  7:1649/1‑147 (MQ=255)
                                                 gatgcgGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGc                                                                                                   >  3:28201/4‑149 (MQ=255)
                                                 aTCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGc                                                                                                   >  5:209360/1‑149 (MQ=255)
                                                    gcgGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTGCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTga                                                                                                <  5:292111/149‑1 (MQ=255)
                                                    gcgGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTga                                                                                                <  7:52062/149‑1 (MQ=255)
                                                     cgGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGATcaggc                                                                                                                                    <  8:321552/113‑1 (MQ=255)
                                                     cgGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGATcaggc                                                                                                                                    >  7:321552/1‑113 (MQ=255)
                                                      ggCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTgaga                                                                                              <  6:20506/149‑1 (MQ=255)
                                                           tttCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGt                                                                                                                  >  8:322146/1‑125 (MQ=255)
                                                           tttCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGt                                                                                                                  <  7:322146/125‑1 (MQ=255)
                                                             tCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATgc                                                                                         <  6:322918/147‑1 (MQ=255)
                                                               gcgACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATgcgcg                                                                                      >  5:150177/1‑148 (MQ=255)
                                                               gcgACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATgcgcg                                                                                      >  7:186127/1‑148 (MQ=255)
                                                                cgACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCa                                                                                    <  2:262970/149‑1 (MQ=255)
                                                                     cAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGc                                                                               <  6:263290/149‑1 (MQ=255)
                                                                             cTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGACCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGc                                                                        >  1:241881/1‑148 (MQ=255)
                                                                                ccACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGACCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGaa                                                                     <  2:241881/148‑1 (MQ=255)
                                                                                 cACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGc                                                                   <  5:342469/149‑1 (MQ=255)
                                                                                 cACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGc                                                                   >  6:370984/1‑149 (MQ=255)
                                                                                     gAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTc                                                                <  3:16267/149‑1 (MQ=255)
                                                                                       aTACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCaa                                                              <  5:8752/149‑1 (MQ=255)
                                                                                              aCCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAAGCATGGc                                                       <  6:209360/149‑1 (MQ=255)
                                                                                              aCCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAAGCATGGc                                                       >  8:219857/1‑149 (MQ=255)
                                                                                                 tCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAAGCATGGCctg                                                    >  6:194819/1‑149 (MQ=255)
                                                                                                   cAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAAGCATGGCctgct                                                  >  2:272449/1‑149 (MQ=255)
                                                                                                   cAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAAGCATGGCctg                                                    >  4:368434/1‑147 (MQ=255)
                                                                                                   cAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAAGCATGGCctg                                                    <  3:368434/147‑1 (MQ=255)
                                                                                                      gagagaGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAAGCATGGCctgctg                                                 >  1:95257/1‑147 (MQ=255)
                                                                                                         gagagGACGACGTGTGGGGGCGTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAAGCATGGCCTGCTGCAgg                                             <  1:325139/149‑1 (MQ=255)
                                                                                                               acgacCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAAGCATGGCCTGCTGCAGGGTCAGc                                       >  7:37507/1‑149 (MQ=255)
                                                                                                                   ccTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAAGCATGGCCTGCTGCAGGGTCAGCTGaa                                   <  6:150177/149‑1 (MQ=255)
                                                                                                                       tgGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAAGCATGGCCTGCTGCAGGGTCAGCTGAAtg                                 <  6:17481/147‑1 (MQ=255)
                                                                                                                         ggggCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAAGCATGGCCTGCTGCAGGGTCAGCTGAAtggtg                              <  1:272449/148‑1 (MQ=255)
                                                                                                                            gCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAAGCATGGCCTGCTGCAGGGTCAGCTGAATGGTGGcc                           >  7:119053/1‑148 (MQ=255)
                                                                                                                            gCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAAGCATGGCCTGCTGCAGGGTAAGCTGAATGGTGGcc                           >  3:15161/1‑148 (MQ=255)
                                                                                                                               ttGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAAGCATGGCCTGCTGCAGGGTCAGCTGAATGGTGGCCGAgg                       >  6:218379/1‑149 (MQ=255)
                                                                                                                                       cgATAGGGCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAAGCATGGCCTGCTGCAGGGTCAGCTGAATGGTGGCCGAGGCCATGTCa               <  8:37507/149‑1 (MQ=255)
                                                                                                                                         ctggggCTGCTGCGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAAGCATGGCCTGCTGCAGGGTCAGCTGAATGGTGGCCGAGGCCATGTCAg              >  7:37948/4‑148 (MQ=255)
                                                                                                                                         ctggggCTGCTGCGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAAGCATGGCCTGCTGCAGGGTCAGCTGAATGGTGGCCGAGGCCATGTCAg              >  5:95202/4‑148 (MQ=255)
                                                                                                                                          tggggCTGCTGCGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAAGCATGGCCTGCTGCAGGGTCAGCTGAAtggtg                              >  5:245890/3‑131 (MQ=255)
                                                                                                                                          tggggCTGCTGCGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAAGCATGGCCTGCTGCAGGGTCAGCTGAAtggtg                              <  6:245890/129‑1 (MQ=255)
                                                                                                                                             ggCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAAGCAt                                                          <  2:45114/100‑1 (MQ=255)
                                                                                                                                             ggCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAAGCAt                                                          >  1:45114/1‑100 (MQ=255)
                                                                                                                                              gCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAAGCATGGCCTGCTGCAGGGTCAGCTGAATGGTGGCCGAGGCCATGTCAGTATCGc        >  2:80319/1‑149 (MQ=255)
                                                                                                                                                tGCGCAGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAAGCAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAAGCATGGCCTGCTGCAGGGTCAGCTGAATGGTGGCCGAGGCCATGTCAGTATCGc        >  1:268391/1‑148 (MQ=255)
                                                                                                                                                  cgcAGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAAGCATGGCCTGCTGCAGGGTCAGCTGAATGGTGGGCGAGGCCATGTCAGTATCGCCAAt    <  8:1649/149‑1 (MQ=255)
                                                                                                                                                  cgcAGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAAGCATGGCCTGCTGCAGGGTCAGCTGAATGGTGGCCGAGGCCATGTCAGTATCGCCAAt    <  3:154630/149‑1 (MQ=255)
                                                                                                                                                   gcagcaGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAAGCATGGCCTGCTGCAGGGTCAGCTGAATGGTGGCCGTGGCCATGTCAGTATCGCCAATg   >  8:244739/1‑149 (MQ=255)
                                                                                                                                                   gcagcaGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAAGCATGGCCTGCTGCAGGGTCAGCTGAATGGTGGCCGAGGCCATGTCAGTATCGCCAATg   <  8:186127/149‑1 (MQ=255)
                                                                                                                                                   gcagcaGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAAGCATGGCCTGCTGCAGGGTCAGCTGAATGGTGGCCGAGGCCATGTCAGTATCGCCAATg   <  4:28201/149‑1 (MQ=255)
                                                                                                                                                    cagcagCCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAAGCATGGCCTGCTGCAGGGTCAGCTGAATGGTGGCCGTGGCCATGTCAGTATCGCCAATgg  <  7:244739/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GCACTGGACCCCGAATGATCAGCTTCGCAGCAGATATCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAAGCATGGCCTGCTGCAGGGTCAGCTGAATGGTGGCCGAGGCCATGTCAGTATCGCCAATGG  >  NC_002947/4971809‑4972103

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: