Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,083,596 G→A 6.9% intergenic (‑61/+207) purE ← / ← PP_5337 5‑(carboxyamino)imidazole ribonucleotide mutase/LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,083,5960GA 6.9% 214.5 / 3.9 72intergenic (‑61/+207)purE/PP_53375‑(carboxyamino)imidazole ribonucleotide mutase/LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (36/31);  new base A (2/3);  total (38/34)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.26e-01

CCGAGTTTTTCCAGCATATCGGCGGTGTGGCTAAGGGTGGACCAATCGGACTTGGAGCCCATGATCACGCCAACCAGTGCACTCATCGTCGAGCCTCTCCTGCAAGCGCCTTCTGGCGCGCTAAAAGCAACAAGCCACGCAGGTGGGCCGGCGTGGCTTGTCATACGGATTCTGAATCGACCAGGTCGGCCGAAGGCGCGGGATTGTAACGTATTTAAAGGCTGGCCGCCCACCCCTACAAGCACATGCCGGATCAAAGGGCAAGCTGTCTCGTGGGACAACTTTTTCCTGCTG  >  NC_002947/6083450‑6083743
                                                                                                                                                  |                                                                                                                                                   
ccGAGTTTTTCCAGCATATCGGCGGTGTGGCTAAGGGTGGACCAATCGGACTTGGAGCCCATGATCACGCCAACCAGTGCACTCATCGTCGAGCCTCTCCTGCAAGCGCCTTCTGGCGCGCTAAAAGCAACAAGCCACGCAGGTGGGcc                                                                                                                                                   >  5:156118/1‑149 (MQ=255)
 cGAGTTTTTCCAGCATATCGGCGGTGTGGCTAAGGGTGGACCAATCGGACTTGGAGCCCATGATCACGCCAACCAGTGCACTCATCGTCGAGCCTCTCCTGCAAGCGCCTTCTGGCGCGCTAAAAGCAACAAGCCACGCAGGTGGGCCg                                                                                                                                                  >  1:282454/1‑149 (MQ=255)
 cGAGTTTTTCCAGCATATCGGCGGTGTGGCTAAGGGTGGACCAATCGGACTTGGAGCCCATGATCACGCCAACCAGTGCACTCATCGTCGAGCCTCTCCTGCAAGCGCCTTCTGGCGCGCTAAAAGCAACAAGCCACGCAGGTGGGCCg                                                                                                                                                  >  7:340194/1‑149 (MQ=255)
        ttCCAGCATATCGGCGGTGTGGCTAAGGGTGGACCAATCGGACTTGGAGCCCATGATCACGCCACCCAGTCCACTCATCGTCGAGCCTCTCCTGCAAGCGCCTTCTGGCGCGCTAAAAGCAACAAGCCACGCAGGTGGGCCGGCGTgg                                                                                                                                            <  6:381221/148‑1 (MQ=255)
         tCCAGCATATCGGCGGTGTGGCTAAGGGTGGACCAATCGGACTTGGAGCCCATGATCACGCCAACCAGTGCACTCATCGTCGAGCCTCTCCTGCAAGCGCCTTCTGGCGCGCTAAAAGCAACAAGCCACGCAGGTGGGCCGGCGTGGCt                                                                                                                                          >  5:302130/1‑149 (MQ=255)
           cAGCATATCGGCGGTGTGGCTAAGGGTGGACCAATCGGACTTGGAGCCCATGATCACGCCAACCAGTGCACTCATCGTCGAGCCTCTCCTGCAAGCGCCTTCTGGCGCGCTAAAAGCAACAAGCCACGCAGGTGGGCCGGCGTGGCTTg                                                                                                                                        >  8:360515/1‑149 (MQ=255)
                  tCGGCGGTGTGGCTAAGGGTGGACCAATCGGACTTGGAGCCCATGATCACGCCAACCAGTGCACTCATCGTCGAGCCTCTCCTGCAAGCGCCTTCTGGCGCGCTAAAAGCAACAACCCACGCCGGCCCACCTGCGTGGCTTGttgcttt                                                                                                                                 <  5:41654/149‑7 (MQ=38)
                   cggcggTGTGGCTAAGGGTGGACCAATCGGACTTGGAGCCCATGATCACGCCAACCAGTGCACTCATCGTCGAGCCTCTCCTGCAAGCGCCTTCTGGCGCGCTAAAAGCAACAAGCCACGCAGGTGGGCCGGCGTGGCTTGTCATACg                                                                                                                                 >  5:220381/1‑148 (MQ=255)
                          gtgGCTAAGGGTGGACCAATCGGACTTGGAGCCCATGATCACGCCAACCAGTGCACTCATCGTCGAGCCTCTCCTGCAAGCGCCTTCTGGCGCGCTAAAAGCAACAAGCCACGCAGGTGGACCGGCGTGGc                                                                                                                                           >  8:355403/1‑131 (MQ=255)
                          gtgGCTAAGGGTGGACCAATCGGACTTGGAGCCCATGATCACGCCAACCAGTGCACTCATCGTCGAGCCTCTCCTGCAAGCGCCTTCTGGCGCGCTAAAAGCAACAAGCCACGCAGGTGGACCGGCGTGGc                                                                                                                                           <  7:355403/131‑1 (MQ=255)
                           tgGCTAAGGGTGGACCAATCGGACTTGGAGCCCATGATCACGCCAACCAGTGCACTCATCGTCGAGCCTCTCCTGCAAGCGCCTTCTGGCGCGCTAAAAGCAACAAGCCACGCAGGTGGGCCGGCGTGGCTTGTCATACGGATTCTGa                                                                                                                         <  6:156118/148‑1 (MQ=255)
                               tAAGGGTGGACCAATCGGACTTGGAGCCCATGATCACGCCAACCAGTGCACTCATCGTCGAGCCTCTCCTGCAAGCGCCTTCTGGCGCGCTAAAAGCAACAAGCCACGCAGGTGGGCCGGCGTGGCTTGTCATACGGATTCTGAAt                                                                                                                       <  4:211380/146‑1 (MQ=255)
                               tAAGGGTGGACCAATCGGACTTGGAGCCCATGATCACGCCAACCAGTGCACTCATCGTCGAGCCTCTCCTGCAAGCGCCTTCTGGCGCGCTAAAAGCAACAAGCCACGCAGGTGGGCCGGCGTGGCTTGTCATACGGATTCTGAAt                                                                                                                       >  3:211380/1‑146 (MQ=255)
                                 aGGGTGGACCAATCGGACTTGGAGCCCATGATCACGCCAACCAGTGCACTCATCGTCGAGCCTCTCCTGCAAGCGCCTTCTGGCGCGCTAAAAGCAACAAGCCACGCAGGTGGGCCGGCGTGGCTTGTCATACGGATTCTGAATCGAcc                                                                                                                  <  2:17119/149‑1 (MQ=255)
                                   ggTGGACCAATTGGACTTGGAGCCCATGATCACGCCAACCAGTGCACTCATCGTCGAGCCTCTCCTGCAAGCGCCTTCTGGCGCGCTAAAAGCAACAAGCCACGCAGGTGGGCCGGCGTGGCTTGTCATACGGATTCTGAATCGACCAg                                                                                                                >  2:124444/1‑149 (MQ=255)
                                       gACCAATCGGACTTGGAGCCCATGATCACGCCAACCAGTGCACTCATCGTCGAGCCTCTCCTGCAAGCGCCTTCTGGCGCGCTAAAAGCAACAAGCCACGCAGGTGGGCCGGCGTGGCTTGTCATACGGATTCTGAATCGACCAGGTCg                                                                                                            >  2:118062/1‑149 (MQ=255)
                                        aCCAATCGGACTTGGAGCCCATGATCACGCCAACCAGTGCACTCATCGTCGAGCCTCTCCTGCAAGCGCCTTCTGGCGCGCTAAAAGCAACAAGCCACGCAGGTGGGCCGGCGTGGCTTGTCATACGGATTCTGAATCGACCAGGTCgg                                                                                                           <  4:52073/149‑1 (MQ=255)
                                         ccAATCGGACTTGGAGCCCATGATCACGCCAACCAGTGCACTCATCGTCGAGCCTCTCCTGCAAGCGCCTTCTGGCGCGCTAAAAGCAACAAGCCACGCAGGTGGGCCGGCGTGGCTTGTCATACGGATTCTGAATCGACCAGGTCGGc                                                                                                          <  2:282454/149‑1 (MQ=255)
                                          cAATCGGACTTGGAGCCCATGATCACGCCAACCAGTGCACTCATCGTCGAGCCTCTCCTGCAAGCGCCTTCTGGCGCGCTAAAAGCAACAAGCCACGCAGGTGGGCCGGCGTGGCTTGTCATACGGATTCTGAATCGACCAGGTCGGcc                                                                                                         <  4:272672/149‑1 (MQ=255)
                                            aTCGGACTTGGAGCCCATGATCACGCCAACCAGTGCACTCATCGTCGAGCCTCTCCTGCAAGCGCCTTCTGGCGCGCTAAAAGCAACAAGCCACGCAGGTGGGCCGGCGt                                                                                                                                              <  1:378597/110‑1 (MQ=255)
                                            aTCGGACTTGGAGCCCATGATCACGCCAACCAGTGCACTCATCGTCGAGCCTCTCCTGCAAGCGCCTTCTGGCGCGCTAAAAGCAACAAGCCACGCAGGTGGGCCGGCGt                                                                                                                                              >  2:378597/1‑110 (MQ=255)
                                             tCGGACTTGGAGCCCATGATCACGCCAACCAGTGCACTCATCGTCGAGCCTCTCCTGCAAGCGCCTTCTGGCGCGCTAAAAGCAACAAGCCACGCAGGTGGGCCGGCGTGGCTTGTCATACGGATTCTGAATCGACCAGGTCGGCCg                                                                                                        >  4:136280/1‑147 (MQ=255)
                                             tCGGACTTGGAGCCCATGATCACGCCAACCAGTGCACTCATCGTCGAGCCTCTCCTGCAAGCGCCTTCTGGCGCGCTAAAAGCAACAAGCCACGCAGGTGGGCCGGCGTGGCTTGTCATACGGATTCTGAATCGACCAGGTCGGCCGaa                                                                                                      <  8:332979/149‑1 (MQ=255)
                                               ggACTTGGAGCCCATGATCACGCCAACCAGTGCACTCATCGTCGAGCCTCTCCTGCAAGCGCCTTCTGGCGCGCTAAAAGCAACAAGCCACGCAGGTGGGCCGGCGTGGCTTGTCATACGGATTCTGAATCGACCAGGTCGGCCGAAg                                                                                                     <  4:135227/148‑1 (MQ=255)
                                                   ttGGAGCCCATGATCACGCCAACCAGTGCACTCATCGTCGAGCCTCTCCTGCAAGCGCCTTCTGGCGCGCTAAAAGCAACAAGCCACGCAGGTGGGCCGGCGTGGCTTGTCATACGGATTCTGAATCGACCAGGTCGGCCGAAGGcgcg                                                                                                >  7:166386/1‑149 (MQ=255)
                                                     ggAGCCCATGATCACGCCAACCAGTGCACTCATCGTCGAGCCTCTCCTGCAAGCGCCTTCTGGCGCGCTAAAAGCAACAAGCCACGCAGGTGGGCCGGCGTGGCTTGTCATACGGATTCTGAATCGACCAGGTCGGCCGAAGGcgcg                                                                                                <  4:135185/147‑1 (MQ=255)
                                                     ggAGCCCATGATCACGCCAACCAGTGCACTCATCGTCGAGCCTCTCCTGCAAGCGCCTTCTGGCGCGCTAAAAGCAACAAGCCACGCAGGTGGGCCGGCGTGGCTTGTCATACGGATTCTGAATCGACCAGGTCGGCCGAAGGCGCgg                                                                                               >  1:325651/1‑148 (MQ=255)
                                                       aGCCCATGATCACGCCAACCAGTGCACTCATCGTCGAGCCTCTCCTTCAAGCGCCTTCTTGCGCGCTAAAAGCAACAAGCCACGCAGGTGGGCCGGCGTGGCTTGTCATACGGATTCTGAATCGACCAGGTCGGCCGAAGGCGCGGGAt                                                                                            <  6:402111/149‑1 (MQ=255)
                                                         cccATGATCACGCCAACCAGTGCACTCATCGTCGAGCCTCTCCTGCAAGCGCCTTCTGGCGCGCTAAAAGCAACAAGCCACGCAGGTGGGCCGGCGTGGCTTGTCATACGGATTCTGAATCGACCAGGTCGGCCGAAGGCGCGGGATTg                                                                                          >  3:341050/1‑149 (MQ=255)
                                                             tGATCACGCCAACCAGTGCACTCATCGTCGAGCCTCTCCTGCAAGCGCCTTCTGGCGCGCTAAAAGCAACAAGCCACGCAGGTGGGCCGGCGTTGCTTGTCATACGGATTCTGAATCGACCAGGTCGGCCGAAGGCGCGGGATTGTAAc                                                                                      <  7:360515/149‑1 (MQ=255)
                                                             tGATCACGCCAACCAGTGCACTCATCGTCGAGCCTCTCCTGCAAGCGCCTTCTGGCGCGCTAAAAGCAACAAGCCACGCAGGTGGGCCGGCGTGGCTTGTCATACGGATTCTGAATCGACCAGGTCGGCCGAAGGCGCGGGATTGTAAc                                                                                      <  8:166386/149‑1 (MQ=255)
                                                              gATCACGCCAACCAGTGCACTCATCGTCGAGCCTCTCCTGCAAGCGCCTTCTGGCGAGCTAAAAGCAACAAGCCACGCAGGTGGGCCGGCGTGGCTTGTCATACGGATTCTGAATCGACCAGGt                                                                                                              >  1:111481/1‑124 (MQ=255)
                                                                  aCGCCAACCAGTGCAGTCATCGTCGTGCCTCTCCTGCAAGCGCCTTCTTGCGCGCTAAAAGCAAGAAGCCACGCAGGTTGGCCGGCGTGGCTTGTCATACGGATTCTGAATCGACCAGGGCGGCCGAAGGCGCGGGATTGTAACGTAt                                                                                  <  7:347451/148‑1 (MQ=255)
                                                                          cattGCACTCATCGTCGAGCCTCTCCTGCAAGCGCCTTCTGGCGCGCTAAAAGCAACAAGCCACGCAGGTGGGCCGGCGTGGCTTGTCATACGGATTCTGAATCGACCAGGTCGGCCGAAGGCGCGGGATTGTAACGTATTTAAAGGCt                                                                         >  3:105501/4‑149 (MQ=255)
                                                                           aGTGCACTCATCGTCGAGCCTCTCCTGCAAGCGCCTTCTGGCGCGCTAAAAGCAACAAGCCACGCAGGTGGGCCGGCGTGGCTTGTCATAc                                                                                                                                  <  7:272686/91‑1 (MQ=255)
                                                                           aGTGCACTCATCGTCGAGCCTCTCCTGCAAGCGCCTTCTGGCGCGCTAAAAGCAACAAGCCACGCAGGTGGGCCGGCGTGGCTTGTCATAc                                                                                                                                  >  8:272686/1‑91 (MQ=255)
                                                                                 cTCATCGTCGAGCCTCTCCTGCAAGCGCCTTCTGGCGCGCTAAAAGCAACAAGCCACGCAGGTGGGCCGGCGTGGCTTGTCATACGGATTCTGAATCGACCAGGTCGGCCGAAGGCGCGGGATTGTAACGTATTTAAAGGCTGgccgcc                                                                  <  1:118062/149‑1 (MQ=255)
                                                                                    aTCGTCGAGCCTCTCCTGCAAGCGCCTTCTGGCGCGCTAAAAGCAACAAGCCACGCAGGTGGGCCGGCGTGGCTTGTCATACGGATTCTGAATCGACCAGGTCGGCCGAAGGCGCGGGATTGTAACGTATTTAAAGGCTGGCCGCCCa                                                                <  6:302130/148‑1 (MQ=255)
                                                                                    aTCGTCGAGCCTCTCCTGCAAGCGCCTTCTGGCGCGCTAAAAGCAACAAGCCACGCAGGTGGGCCGGCGTGGCTTGTCATACGGATTCTGAATCGACCAGGTCGGCCGAAGGCGCGGGATTGGAACGTATTTAAAGGCTGGCCGCCCAc                                                               >  3:396172/1‑149 (MQ=255)
                                                                                          gAGCCTCTCCTGCAAGCGCCTTCTGGCGCGCTAAAAGCAACAAGCCACGCAGGTGGGCCGGCGTGGCTTGTCATACGGATTCTGAATCGACCAGGTCGGCCGAAGGCGCGGGATTGTAACGTATTTAAAGGCTGGCCGCCCACCCCTAc                                                         <  6:220381/149‑1 (MQ=255)
                                                                                             ccTCTCCTGCAAGCGGCTTCTGGCGCGCTATAAGCAACAAGCGACGCAGGTGGGCAGGCGTGCATTGTCATACGGATTCTGAATCGACCAGGTCGGCCGAAGGCGCGGGATTGTAACGTATTTAAAGGCTGGCCGACCTCCCCTCCAAg                                                      >  5:293487/1‑149 (MQ=255)
                                                                                               tctcCTGCATGCTCCTTCTGGGGCGCTAAAAGGAACAAGCCACGCAGGTGGGCCGGCGTGGCTTGTCATACGGATTCTGAATCGACCAGGTCGGCCGAAGGCGCGGGATTGTAACGTATTTaaa                                                                             <  8:347136/124‑1 (MQ=255)
                                                                                               tctcCTGCAAGCGCCTTCTGGCGCGCTAAAAGCAACAAGCCACGCAGGTGGGCCGGCGTGGCTTGTCATACGGATTCTGAATCGACCAGGTCGGCCGAAGGCGCGGGATTGTAACGTATTTaaa                                                                             >  7:347136/1‑124 (MQ=255)
                                                                                                   cTGCAAGCGCCTTCTGGCGCGCTAAAAGCAACAAGCCACGCAGGTGGGCCGGCGTGGCTTGTCATACGGATTCTGAATCGACCAGGTCGGCCGAAGGCGCGGGATTGTAACGTATTTaaa                                                                             <  6:210762/120‑1 (MQ=255)
                                                                                                   cTGCAAGCGCCTTCTGGCGCGCTAAAAGCAACAAGCCACGCAGGTGGGCCGGCGTGGCTTGTCATACGGATTCTGAATCGACCAGGTCGGCCGAAGGCGCGGGATTGTAACGTATTTaaa                                                                             >  5:210762/1‑120 (MQ=255)
                                                                                                     gCAAGCGCCTTCTGGCGCGCTAAAAGCAACAAGCCACGCAGGTGGGCCGGCGTGGCTTGTCATACGGATTCTGAATCGACCAGGTCGGCCGAAGGCGCGGGATTGTAACGTATTTAAAGGCTGGCCGCCCACCCCTACAAGCACATGcc                                              >  4:125669/1‑149 (MQ=255)
                                                                                                         gcgcCTTCTGGCGCGCTAAAAGCAACAAGCCACGCAGGTGGGCCGGCGTGGCTTGTCATACGGATTCTGAATCGACCAGGTCGGCCGAAGGCGCGGGATTGTAACGTATTTAAAGGCTGGCCGCCCACCCCTACAAGCACATGCCGGa                                           >  7:143515/1‑148 (MQ=255)
                                                                                                         gcgcCTTCTGGCGCGCTAAAAGCAACAAGCCACGCAGGTGGGCCGGCGTGGCTTGTCATACGGATTCTGAATCGACCAGGTCGGCCGAAGGCGCGGGATTGTAACGTATTTAAAGGCTGGCCGCCCACCCCTACAAGCACATGCCGGAt                                          <  8:340194/149‑1 (MQ=255)
                                                                                                              ttCTGGCGCGCTAAAAGCAACAAGCCACGCAGGTGGGCCGGCGTGGCTTGTCATACGGATTCTGAATCGACCAGGTCGGCCGAAGGCGCGGGATTGTAACGTATTTAAAGGCTGGCCGCCCACCCCTACAAGCACATGCCGGATCAAAg                                     >  8:44853/1‑149 (MQ=255)
                                                                                                                   gcgcgcTAAAAGCAACAAGCCACGCAGGTGGGCCGGCGTGGCTTGTCATACGGATTCTGAATCGACCAGGTCGGCCGAAGGCGCGGGATTGTAACGTATTTAAAGGCTGGCCGCCCACCCCTACAAGCACATGCCg                                             <  6:104678/136‑1 (MQ=255)
                                                                                                                   gcgcgcTAAAAGCAACAAGCCACGCAGGTGGGCCGGCGTGGCTTGTCATACGGATTCTGAATCGACCAGGTCGGCCGAAGGCGCGGGATTGTAACGTATTTAAAGGCTGGCCGCCCACCCCTACAAGCACATGCCg                                             >  5:104678/1‑136 (MQ=255)
                                                                                                                   gcgcgcTAAAAGCAACAAGCCACGCAGGTGGGCCGGCGTGGCTTGTCATACGGATTCTGAATCGACCAGGTCGGCCGAAGGCGCGGGATTGTAACGTATTTAAAGGCTGGCCGCCCACCCCTACAAGCACATGCCGGATCa                                        <  1:265035/141‑1 (MQ=255)
                                                                                                                   gcgcgcTAAAAGCAACAAGCCACGCAGGTGGGCCGGCGTGGCTTGTCATACGGATTCTGAATCGACCAGGTCGGCCGAAGGCGCGGGATTGTAACGTATTTAAAGGCTGGCCGCCCACCCCTACAAGCACATGCCGGATCa                                        >  2:265035/1‑141 (MQ=255)
                                                                                                                   gcgcgcTAAAAGCAACAAGCCACGCAGGTGGGCCGGCGTGGCTTGTCATACGGATTCTGAATCGACCAGGTCGGCCGAAGGCGCGGGATTGTAACGTATTTAAAGGCTGGCCGCCCACCCCTACAAGCACATGCCGGATCAAAg                                     <  7:251240/144‑1 (MQ=255)
                                                                                                                   gcgcgcTAAAAGCAACAAGCCACGCAGGTGGGCCGGCGTGGCTTGTCATACGGATTCTGAATCGACCAGGTCGGCCGAAGGCGCGGGATTGTAACGTATTTAAAGGCTGGCCGCCCACCCCTACAAGCACATGCCGGATCAAAg                                     >  8:251240/1‑144 (MQ=255)
                                                                                                                         tAAAAGCAACAAGCCACGCAGGTGGGCCGGCGTGGCTTGTCATACGGATTCTGAATCGACCAGGTCGGCCGAAGGCGCGGGATTGTAACGTATTTAAAGGCTGGCCGCCCACCCCTACAAGCACATGCCGGATCAAAGGGCAAGCt                             <  5:198807/146‑1 (MQ=255)
                                                                                                                         tAAAAGCAACAAGCCACGCAGGTGGGCCGGCGTGGCTTGTCATACGGATTCTGAATCGACCAGGTCGGCCGAAGGCGCGGGATTGTAACGTATTTAAAGGCTGGCCGCCCACCCCTACAAGCACATGCCGGATCAAAGGGCAAGCt                             >  6:198807/1‑146 (MQ=255)
                                                                                                                         tAAAAGCAACAAGCCACGCAGGTGGGCCGGCGTGGCTTGTCATACGGATTCTGAATCGACCAGGTCGGCCGAAGGCGCGGGATTGTAACGTATTTAAAGGCTGGCCGCCCACCCCTACAAGCACATGCCGGATCAAAGGGCAAGCTGtc                          >  1:315884/1‑149 (MQ=255)
                                                                                                                          aaaaGCAACAAGCCACGCAGGTGGGCCGGCGTGGCTTGTCATACGGATTCTGAATCGACCAGGTCGGCCGAAGGCGCGGGATTGTAACGTATTTAAAGGCTGGCCGCCCACCCCTACAAGCACATGCCGGATCAAAGGGCAAGCTGtct                         <  8:143515/149‑1 (MQ=255)
                                                                                                                             aGCAACAAGCCACGCAGGTGGGCCGGCGTGGCTTGTCATACGGATTCTGAATCGACCAGGTCGGCCGAAGGCGCGGGATTGTAACGTATTTAAAGGCTGGCCGCCCACCCCTACAAGCACATGCCGGATCAAAGGGCAAGCTGTCTCGt                      >  8:6367/1‑149 (MQ=255)
                                                                                                                             aGCAACAAGCCACGCAGGTGGGCCGGCGTGGCTTGTCATACGGATTCTGAATCGACCAGGTCGGCCGAAGGCGCGGGATTGTAACGTATTTAAAGGCTGGCCGCCCACCCCTACAAGCACATGCCGGATCAAAGGGCAAGCTGTCTCGt                      <  4:341050/149‑1 (MQ=255)
                                                                                                                             aGCAACAAGCCACGCAGGTGGGCCGGCGTGGCTTGTCATACGGATTCTGAATCGACCAGGTCGGCCGAAGGCGCGGGATTGTAACGTATTTAAAGGCTGGCCGCCCACCCCTACAAGCACATGCCGGATCAAAGGGCAAGCTGTCTCGt                      >  1:216506/1‑149 (MQ=255)
                                                                                                                                  caaGCCACGCCGGCCCACCTGCGTGGCTTGTCATACGGATTCTGAATCGACCAGGTCGGCCGAAGGCGCGGGATTGTAACGTATTTAAAGGCTGGCCGCCCACCCCTACAAGCACATGCCGGa                                           >  3:181272/1‑123 (MQ=38)
                                                                                                                                  caaGCCACGCCGGCCCACCTGCGTGGCTTGTCATACGGATTCTGAATCGACCAGGTCGGCCGAAGGCGCGGGATTGTAACGTATTTAAAGGCTGGCCGCCCACCCCTACAAGCACATGCCGGa                                           <  4:181272/123‑1 (MQ=38)
                                                                                                                                       cACGCAGGTGGGCCGGCGTGGCTTGTCATACGGATTCTGAATCGACCAGGTCGGCCGAAGGCGCGGGATTGTAACGTATTTAAAGGCTGGCCGCCCACCCCTACAAGCACATGCCGGATCAAAGGGCAAGCTGTCTCGTg                     >  7:54737/1‑140 (MQ=255)
                                                                                                                                       cACGCAGGTGGGCCGGCGTGGCTTGTCATACGGATTCTGAATCGACCAGGTCGGCCGAAGGCGCGGGATTGTAACGTATTTAAAGGCTGGCCGCCCACCCCTACAAGCACATGCCGGATCAAAGGGCAAGCTGTCTCGTg                     <  8:54737/140‑1 (MQ=255)
                                                                                                                                        aCGCAGGTGGGCCGGCGTGGCTTGTCATACGGATTCTGAATCGACCAGGTCGGCCGAAGGCGCGGGATTGTAACGTATTTAAAGGCTGGCCGCCCACCCCTACAAGCACATGCCGGATCAAAGGGCAAGCTGTCTCGTGGGACAACttt           >  3:19046/1‑149 (MQ=255)
                                                                                                                                           cAGGTGGGCCGGCGTGGCTTGTCATACGGATTCTGAATCGACCAGGTCGGCCGAAGGCGCGGGATTGTAACGTATTTAAAGGCTGGCCGCCCACCCCTACAAGCACATGCCGGATCAAAGGGCAAGCTGTCTCGTGGGAAAACTTTTTc        >  6:295464/1‑149 (MQ=255)
                                                                                                                                            aGGTGGGCCGGCGTGGCTTGTCATACGGATTCTGAATCGACCAGGTCGGCCGAAGGCGCGGGATTGTAACGTATTTAAAGGCTGGCCGCCCACCCCTACAAGCACATGCCGGATCAAAg                                     >  5:418214/1‑119 (MQ=255)
                                                                                                                                            aGGTGGGCCGGCGTGGCTTGTCATACGGATTCTGAATCGACCAGGTCGGCCGAAGGCGCGGGATTGTAACGTATTTAAAGGCTGGCCGCCCACCCCTACAAGCACATGCCGGATCAAAg                                     <  6:418214/119‑1 (MQ=255)
                                                                                                                                             ggTGGGCCGGCGTGGCTTGTCATACGGATTCTGAATCGACCAGGTCGGCCGAAGGCGCGGGATTGTAACGTATTTAAAGGCTGGCCGCCCACCCCTACAAGCACATGCCGGATCAAAGGGCAAGCTGTCTCGTGGGACAACTTTTTCct      <  1:124444/149‑1 (MQ=255)
                                                                                                                                              gTGGGCCGGCGTGGCTTGTCATACGGATTCTGAATCGACCAGGTCGGCCGAAGGCGCGGGATTGTAACGTATTTAAAGGCTGGCCGCCCACCCCTACAAGCACATGCCGGATCAAAGGGCAAGCTGTCTCGTGGGACAACTTTTTCctg     >  3:253495/1‑149 (MQ=255)
                                                                                                                                                gggCCGGCGTGGCTTGTCATACGGATTCTGAATCGACCAGGTCGGCCGAAGGCGCGGGATTGTAACGTATTTAAAGGCTGGCCGCCCACCCCTACAAGCACATGCCGGATCAAAGGGCAAGCTGTCTCGTGGGACAACTTTTTCctgat   >  5:315426/1‑147 (MQ=255)
                                                                                                                                                 ggCCGGCGTGGCTTGTCATACGGATTCTGAATCGACCAGGTCGGCCGAAGGCGCGGGATTGTAACGTATTTAAAGGCTGGCCGCCCACCCCTACAAGCACATGCCGGATCAAAGGGCAAGCTGTCTCGTGGGACAACTTTTTCctgctg  <  3:136280/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CCGAGTTTTTCCAGCATATCGGCGGTGTGGCTAAGGGTGGACCAATCGGACTTGGAGCCCATGATCACGCCAACCAGTGCACTCATCGTCGAGCCTCTCCTGCAAGCGCCTTCTGGCGCGCTAAAAGCAACAAGCCACGCAGGTGGGCCGGCGTGGCTTGTCATACGGATTCTGAATCGACCAGGTCGGCCGAAGGCGCGGGATTGTAACGTATTTAAAGGCTGGCCGCCCACCCCTACAAGCACATGCCGGATCAAAGGGCAAGCTGTCTCGTGGGACAACTTTTTCCTGCTG  >  NC_002947/6083450‑6083743

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: ‑