Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,054,478 G→T 7.2% intergenic (‑15/+22) PP_0912 ← / ← PP_0913 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,054,4780GT 7.2% 314.6 / 10.0 97intergenic (‑15/+22)PP_0912/PP_0913hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (46/44);  new base T (3/4);  total (49/48)
Fisher's exact test for biased strand distribution p-value = 7.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.71e-01

CCCCAGACATCAAGAATCTGCAGCACTGTGCCTCCTGTTAGGTTTGATCGAAGGCCTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGCAATGCCCAGCTTGCCCGCAGCCTGTTCGATCGACCCGGTCTGCTTGACCAGGGCGGCGATGGCGTCCCGCACAACCTGGTT  >  NC_002947/1054333‑1054619
                                                                                                                                                  |                                                                                                                                             
ccccAGACATCAAGAATCTGCAGCACTGTGCCTCCTGTTAGGTTTGATCGAAGGCCTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACTAAAACGCCCGGCACAAGGCCGGGCGTTTGCgg                                                                                                                                            <  8:232496/149‑1 (MQ=255)
 cccAGACATCAAGAAGCTGCTGCACTGTGCCTCATGTTAGGTTTGATCGAAGGCATACAGTATGCGCATCTGCGTGCACGCCCGGCTGCGGGAAGCGTTTTACCCGGCTGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGt                                                                                                                                           <  6:276629/149‑1 (MQ=255)
  ccAGACATCAAGAATCTGCAGCACTGTGCCTCCTGTTAGGTTTGATCGAAGGCCTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTTCgttt                                                                                                                                          <  7:320135/149‑4 (MQ=255)
  ccAGACATCAAGAATCTGCAGCACTGTGCCTCCTGTTAGGTTTGATCGAAGGCCTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTg                                                                                                                                          <  3:337458/149‑1 (MQ=255)
   cAGACATCAAGAATCTGCAGCACTGTGCCTCCTGTTAGGTTTGATCGAAGGCCTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTg                                                                                                                                          >  6:91761/1‑148 (MQ=255)
    aGACATCAAGAATCTGCAGCACTGTGCCTCCTGTTAGGTTTGATCGAAGGCCTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGaa                                                                                                                                        <  5:409664/149‑1 (MQ=255)
      aCATCAAGAATCTGCAGCACTGTGCCTCCTGTTAGGTTTGATCGAAGGCCTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAgc                                                                                                                                      <  4:237525/149‑1 (MQ=255)
         tCAAGAATCTGCAGCACTGTGCCTCCTGTTAGGTTTGATCGAAGGCCTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTg                                                                                                                                               >  8:111878/1‑137 (MQ=255)
         tCAAGAATCTGCAGCACTGTGCCTCCTGTTAGGTTTGATCGAAGGCCTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTg                                                                                                                                               <  7:111878/137‑1 (MQ=255)
           aaGAATCTGCAGCACTGTGCCTCCTGTTAGGTTTGATCGAAGGCCTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTa                                                                                                                                  >  8:346648/1‑148 (MQ=255)
           aaGAATCTGCAGCACTGTGCCTCCTGTTAGGTTTGATCGAAGGCCTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTa                                                                                                                                  >  6:287552/1‑148 (MQ=255)
                   gcagcaCTGTGCCTCCTGTTAGGTTTGATCGAAGGCCTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCCTTGTGCCGGGCGTTTTCgtt                                                                                                                                           >  3:74662/1‑129 (MQ=38)
                   gcagcaCTGTGCCTCCTGTTAGGTTTGATCGAAGGCCTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCCTTGTGCCGGGCGTTTTCgtt                                                                                                                                           <  4:74662/131‑3 (MQ=38)
                       caCTGTGCCTCCTGTTAGGTTTGATCGAAGGCCTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGTGCTACAGCAGGATCAg                                                                                                                      >  4:347712/1‑148 (MQ=255)
                       caCTGTGCCTCCTGTTAGGTTTGATCGAAGGCCTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAg                                                                                                                      >  1:88456/1‑148 (MQ=255)
                              cctcctGTTAGGTTTGATCGAAGGCCTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCg                                                                                                                                             >  4:163610/1‑118 (MQ=255)
                              cctcctGTTAGGTTTGATCGAAGGCCTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCg                                                                                                                                             <  3:163610/118‑1 (MQ=255)
                               ctcctGTTAGGTTTGATCGAAGGCCTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCa                                                                                                                       <  2:58701/139‑1 (MQ=255)
                               ctcctGTTAGGTTTGATCGAAGGCCTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCa                                                                                                                       >  1:58701/1‑139 (MQ=255)
                                tcctGTTAGGTTTGATCGAAGGCCTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGt                                                                                                                                           >  6:49845/1‑118 (MQ=255)
                                tcctGTTAGGTTTGATCGAAGGCCTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGt                                                                                                                                           <  5:49845/118‑1 (MQ=255)
                                 cctGTTAGGTTTGATCGAAGGCCTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATc                                                                                                           >  2:291238/1‑149 (MQ=255)
                                  ctGTTAGGTTTGATCGAAGGCCTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGAt                                                                                                            <  8:338248/147‑1 (MQ=255)
                                              aTCGAAGGCCTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGTGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCgg                                                                                               <  3:347712/148‑1 (MQ=255)
                                              aTCGAAGGCCTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCgg                                                                                               <  8:328299/148‑1 (MQ=255)
                                               tCGAAGGCCTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGt                                                                                             >  4:52004/1‑149 (MQ=255)
                                                 gAAGGCCTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGtt                                                                                            >  6:145234/1‑148 (MQ=255)
                                                 gAAGGCCTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACCCGTGATCCGCAGTGTCCGGGTTc                                                                                           <  2:88456/149‑1 (MQ=255)
                                                      ccTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGtt                                                                                      <  3:52004/149‑1 (MQ=255)
                                                       cTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTg                                                                                     >  8:147019/1‑149 (MQ=255)
                                                       cTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGACACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTg                                                                                     >  8:237027/1‑149 (MQ=255)
                                                           aGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGCaa                                                                                 >  8:5668/1‑149 (MQ=255)
                                                             gATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGTCACATGGCCGGGCGTTTGCGGTGAAGCGCTACAGCATGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAg                                                                                         <  7:226542/139‑1 (MQ=255)
                                                             gATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGCAATg                                                                               >  6:242184/1‑149 (MQ=255)
                                                             gATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGCAATg                                                                               <  7:121099/149‑1 (MQ=255)
                                                             gATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGCAATg                                                                               >  1:132896/1‑149 (MQ=255)
                                                             gATGAGCATCTGCGTGCAAGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGAACAAGGCCGGGAGTTTGCGGTGAAGCGCTACAGCATGATCAGAACTCGTGATACGCAGTTTCAGGGTTCAg                                                                                         >  8:226542/1‑139 (MQ=255)
                                                               tGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTa                                                                                                                                  >  7:104636/1‑96 (MQ=255)
                                                               tGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTa                                                                                                                                  <  8:104636/96‑1 (MQ=255)
                                                                gAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGCAATGcc                                                                             <  6:317430/148‑1 (MQ=255)
                                                                    aTCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGCAATGCCCAGCt                                                                        >  8:242064/1‑149 (MQ=255)
                                                                     tCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGCAATGCCCAGCtt                                                                       >  8:380926/1‑149 (MQ=255)
                                                                      cTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGCAATGCCCAGCTTg                                                                      >  3:47660/1‑149 (MQ=255)
                                                                          gTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAgg                                                                                        >  6:91464/1‑127 (MQ=255)
                                                                          gTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAgg                                                                                        <  5:91464/127‑1 (MQ=255)
                                                                            gCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCg                                                                                                                <  7:325790/101‑1 (MQ=255)
                                                                            gCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCg                                                                                                                >  8:325790/1‑101 (MQ=255)
                                                                            gCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGCAATGCCCAGCTTGCCCg                                                                  >  3:292355/1‑147 (MQ=255)
                                                                            gCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGCAATGCCCAGCTTGCCCGCa                                                                >  2:189759/1‑149 (MQ=255)
                                                                             cACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGc                                                                                   <  3:231102/129‑1 (MQ=255)
                                                                             cACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGc                                                                                   >  4:231102/1‑129 (MQ=255)
                                                                             cACGCCAGGCAGCGGGAAGCGTTTTACCCGGCACCAAACTTAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGCAATGCCCAGCTTGCCCGCa                                                                <  1:135827/149‑1 (MQ=255)
                                                                              aCGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGCAATGCCCAGCTTGCCCGCAGc                                                              <  2:11337/149‑1 (MQ=255)
                                                                                  cAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGCAATGcc                                                                             >  4:68180/1‑130 (MQ=255)
                                                                                  cAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGCAATGcc                                                                             <  3:68180/130‑1 (MQ=255)
                                                                                   aGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGCAATGCCCAGCTTGCCCGCAGCCTGtt                                                         <  7:242064/149‑1 (MQ=255)
                                                                                    ggCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCa                                                                                                       <  7:276827/102‑1 (MQ=255)
                                                                                    ggCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCa                                                                                                       >  8:276827/1‑102 (MQ=255)
                                                                                      cAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGCAATGCCCAGCTTGCCCGCAGCCTGTtcga                                                      <  5:364423/149‑1 (MQ=255)
                                                                                        gCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGCAATGCCCAGCTTGCCCGCAGCCTGTtcgatc                                                    <  5:287552/149‑1 (MQ=255)
                                                                                           ggAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGCAATGCCCAGCTTGCCCGCAGCCt                                                            >  1:316293/1‑138 (MQ=255)
                                                                                           ggAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGCAATGCCCAGCTTGCCCGCAGCCt                                                            <  2:316293/138‑1 (MQ=255)
                                                                                             aaGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGCAATGCCCAGCTTGCCCGCAGCCTGTTCGATCGAccc                                               >  2:395547/1‑149 (MQ=255)
                                                                                                 gTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGCAATGCCCAGCTTGCCCGCAGCCTGTTCGATCGACCCGGTc                                           <  5:145234/149‑1 (MQ=255)
                                                                                                     tACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGCAATGCCCAGCTTGCCCGCAGCCTGTTCGATCGACCCGGTCTGc                                        <  4:47660/148‑1 (MQ=255)
                                                                                                        ccGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGCAATGCCCAGCTTGCCCGCAGCCTGTTCGATCGACCCGGTCTGCTTGa                                    >  6:335885/1‑149 (MQ=255)
                                                                                                        ccGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGCAATGCCCAGCTTGCCCGCAGCCTGTTCGATCGACCCGGTCTGCTTGa                                    >  6:143609/1‑149 (MQ=255)
                                                                                                             agcaAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGCAATGCCCAGCTTGCCCGCAGCCTGTTCGATCGACCCGGTCTGCTTGACCAgg                               <  7:147019/149‑1 (MQ=255)
                                                                                                                gcccGAAAACGCCCGGCACCAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGCAATGCCCAGCTTGc                                                                     <  6:105764/105‑1 (MQ=255)
                                                                                                                aaaCGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGCAATGCCCAGCTTGCCCGCAGCCTGTTCGATCGACCCGGTCTGCTTGACCAGggcg                            >  1:319522/1‑149 (MQ=255)
                                                                                                                  aCGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGCAATGCCCAGCTTGCCCGCAGCCTGTTCGATCGACCCGGTCTGCTTGACCAGggcggc                          <  1:189759/149‑1 (MQ=255)
                                                                                                                  aCGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCACAGTGTCCGGGTTCAGGTTGGCAATGCCCAGCTTGCCCGCAGCCTGTTCGATCGACCCGGTCTGCTTGACCAGggcggc                          >  3:133022/1‑149 (MQ=255)
                                                                                                                   ccgcaaaCGCCCGGCCTTGTGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATACGCAGTGTCCGGGTTCAGGTTGGCAATGCCCAGCTTGCCCGCAGCCTGTTCGATCGACCCGGTCTGCTTGACCAGggcggc                          >  7:62045/5‑149 (MQ=255)
                                                                                                                    cgcaaaCGCCCGGCCTTGTGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTc                                                                                           >  1:354369/4‑83 (MQ=25)
                                                                                                                    cgcaaaCGCCCGGCCTTGTGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTc                                                                                           <  2:354369/80‑1 (MQ=25)
                                                                                                                     gcaaaCGCCCGGCCTTGTGCCGGGCGTTTTCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGCAATGCCCAGCTTGCCCGCAGCCTGTTCGATCGACCCGGTCTGCTTGACCAGGGCGGCGa                        >  2:90276/3‑149 (MQ=255)
                                                                                                                     gcaaaCGCCCGGCCTTGTGCCGGGCGTTTTCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGCAATGCCCAGCTTGCCCGCAGCCTGTTCGATCGACCCGGTCTGCTTGACCAGGGCGGCGa                        >  2:84610/3‑149 (MQ=255)
                                                                                                                       aaaCGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGCAATGCCCAGCTTGCCCGCAGCCTGTTCGATCGACCCGGTCTGCTTGACCAGGGCGGCGATg                      >  6:199269/1‑149 (MQ=255)
                                                                                                                         aCGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGCAATGCCCAGCTTGCCCGCAGCCTGTTCGATCGAcc                                                >  4:67183/1‑121 (MQ=255)
                                                                                                                         aCGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGCAATGCCCAGCTTGCCCGCAGCCTGTTCGATCGAcc                                                <  3:67183/121‑1 (MQ=255)
                                                                                                                           gCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGCAATGCCCAGCTTGCCCGCAGCCTGTTCGATCGACCCGGTCTGCTTGACCAGGGCGGCGATGGCGt                  >  7:4049/1‑149 (MQ=255)
                                                                                                                           gCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGCAATGCCCAGCTTGCCCGCAGCCTGTTCGATCGACCCGGTCTGCTTGACCAGGGCGGCGATGGCGt                  >  2:356661/1‑149 (MQ=255)
                                                                                                                            cccGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGCAATGCCCAGCTTGCCCGCAGCCTGTTCGATCGACCCGGTCTGCTTGACCAGGGCGGCGATGGCGTc                 <  8:4049/149‑1 (MQ=255)
                                                                                                                            cccGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGCAATGCCCAGCTTGCCCGCAGCCTGTTCGATCGACCCGGTCTGCTTGACCAGGGCGGCGATGGCGTc                 <  7:380926/149‑1 (MQ=255)
                                                                                                                             ccGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGCAATGCCCAGCTTGCCCGCAGCCTGTTCGATCGACCCGGTCTGCTTGACCAGGGCGGCGATGGCGTcc                >  3:348829/1‑149 (MQ=255)
                                                                                                                               ggCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGCAATGCCCAGCTTGCCCGCAGCCTGTTCGATCGACCCGGTCTGCTTGACCAGGGCGGCGATGGCGTCCCg              <  7:346648/149‑1 (MQ=255)
                                                                                                                                  acaAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGCAATGCCCAGCTTGCCCGCAGCCTGTTCGATCGACCCGGTCTGCTTGACCAGGGCGGCGa                        <  5:391532/136‑1 (MQ=255)
                                                                                                                                  acaAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGCAATGCCCAGCTTGCCCGCAGCCTGTTCGATCGACCCGGTCTGCTTGACCAGGGCGGCGa                        >  6:391532/1‑136 (MQ=255)
                                                                                                                                  acaAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGCAATGCCCAGCTTGCCCGCAGCCTGTTCGATCGACCCGGTCTGCTTGACCAGGGCGGCGATGGCGTCCCGca            >  1:115564/1‑148 (MQ=255)
                                                                                                                                      ggCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGCAATGCCCAGCTTGCCCGCAGCCTGTTCGATCGACCCGGTCTGCTTGACCAGGGCGGCGATGGCGTCCCGCACAAcc       <  5:143609/149‑1 (MQ=255)
                                                                                                                                        ccGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGCAATGCCCAGCTTGCCCGCAGCCTTTTCGATCGACCCGGTCTGCTTGACCAGGGCGGCGATGGCGTCCCGCACAACCTg     <  8:218121/149‑1 (MQ=255)
                                                                                                                                        ccGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGCAATGCCCAGCTTGCCCGCAGCCTGTTCGATCGACCCGGTCTGCTTGACCAGGGCGGCGATGGCGTCCCGCACAACCTg     <  7:237027/149‑1 (MQ=255)
                                                                                                                                         cGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGCAATGCCCAGCTTGCCCGCAGCCTGTTCGATCGACCCGGTCTGCTTGACCAgg                               >  3:223592/1‑122 (MQ=255)
                                                                                                                                         cGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGCAATGCCCAGCTTGCCCGCAGCCTGTTCGATCGACCCGGTCTGCTTGACCAgg                               <  4:223592/122‑1 (MQ=255)
                                                                                                                                         cGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGCAATGCCCAGCTTGCCCGCAGCCTGTTCGATCGACCCGGTCTGCTTGACCAGGGCGGCGATGGCGTCCCGCACAACCTgg    >  2:183678/1‑149 (MQ=255)
                                                                                                                                         cGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGCAATGCCCAGCTTGCCCGCAGCCTGTTCGATCGACCCGGTCTGCTTGACCAGGGCGGCGATGGCGTCCCGCACAACCTgg    >  6:259842/1‑149 (MQ=255)
                                                                                                                                         cGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGCAATGCCCAGCTTGCCCGCAGCCTGTTCGATCGACCCGGTCTGCTTGACCAGGGCGGCGATGGCGTCCCGCACAACCTgg    <  5:335885/149‑1 (MQ=255)
                                                                                                                                           ggCGTTTTCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGCAATGCCCAGCTTGCCCGCAGCCTGTTCGATCGACCTGGTCTGCTTGACCAGGGCGGCGATGGCGTCCCGCACAACCTGGtt  <  1:90276/149‑1 (MQ=255)
                                                                                                                                           ggCGTTTTCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGCAATGCCCAGCTTGCCCGCAGCCTGTTCGATCGACCCGGTCTGCTTGACCAGGGCGGCGATGGCGTCCCGCACAACCTGGtt  <  1:84610/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                             
CCCCAGACATCAAGAATCTGCAGCACTGTGCCTCCTGTTAGGTTTGATCGAAGGCCTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGCAAACGAAAACGCCCGGCACAAGGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGCAATGCCCAGCTTGCCCGCAGCCTGTTCGATCGACCCGGTCTGCTTGACCAGGGCGGCGATGGCGTCCCGCACAACCTGGTT  >  NC_002947/1054333‑1054619

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: