Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,376,321 A→T 8.3% intergenic (+30/+19) kup → / ← PP_1201 potassium transport system protein Kup/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,376,3210AT 8.3% 208.4 / 7.2 72intergenic (+30/+19)kup/PP_1201potassium transport system protein Kup/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (37/29);  new base T (2/4);  total (39/33)
Fisher's exact test for biased strand distribution p-value = 4.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

TGTCGCGCTGGCGGGGCAACCTGTTCGCGTTCTTGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGCAGTTTCAGGCGCTCACCTACTGCC  >  NC_002947/1376173‑1376469
                                                                                                                                                    |                                                                                                                                                    
tGTCGCGCTGGCGGGGCAACCTGTTCGCGTTCTTGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCa                                                                                                                                                      <  1:35866/149‑1 (MQ=255)
      gcTGGCGGGGCAACCTGTTCGCGTTCTTGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCgg                                                                                                                                                >  2:134024/1‑149 (MQ=255)
          gCGGGGCAACCTGTTCGCGTTCTTGCTTAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCTCTGCGGCCCt                                                                                                                                            <  3:25636/149‑1 (MQ=255)
            ggggCAACCTGTTCGCGTTCTTGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCTCTGCGGCCCTGa                                                                                                                                          <  3:89236/149‑1 (MQ=255)
            ggggCAACCTGTTCGCGTTCTTGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGa                                                                                                                                          <  6:126314/149‑1 (MQ=255)
            ggggCAACCTGTTCGCGTTCTTGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGa                                                                                                                                          <  7:278269/149‑1 (MQ=255)
            ggggCAACCTGTTCGCGTTCTTGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGa                                                                                                                                          <  7:272548/149‑1 (MQ=255)
                   ccTGTTCGCGTTCTTGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTg                                                                                                                                    >  5:158941/1‑148 (MQ=255)
                   ccTGTTCGCGTTCTTGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTg                                                                                                                                    <  6:158941/148‑1 (MQ=255)
                   ccTGTTCGCGTTCTTGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGt                                                                                                                                   <  6:116687/149‑1 (MQ=255)
                             ttCTTGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCt                                                                                                                         <  5:24/149‑1 (MQ=255)
                                  gCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGcc                                                                                                                                              <  1:160566/123‑1 (MQ=255)
                                  gCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGcc                                                                                                                                              >  2:160566/1‑123 (MQ=255)
                                   cTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGTGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCt                                                                                                                    <  5:186683/148‑1 (MQ=255)
                                    tGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCTCTGCGGCCCTGa                                                                                                                                          <  4:81909/125‑1 (MQ=255)
                                    tGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCTCTGCGGCCCTGa                                                                                                                                          >  3:81909/1‑125 (MQ=255)
                                    tGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTAc                                                                                                                  >  4:130198/1‑149 (MQ=255)
                                     gaagaaTGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGa                                                                                                                                          <  3:189456/124‑1 (MQ=255)
                                     gaagaaTGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGa                                                                                                                                          >  4:189456/1‑124 (MQ=255)
                                             ccAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACtt                                                                                                         >  2:175856/1‑149 (MQ=255)
                                              cAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGa                                                                                                                                          <  5:108702/115‑1 (MQ=255)
                                              cAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGa                                                                                                                                          >  6:108702/1‑115 (MQ=255)
                                                    cAACTTGCGCTTCTTCAACCTGCCGTTGATCCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTgg                                                                                                  <  7:260116/149‑1 (MQ=255)
                                                    cAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTc                                                                                                                        <  4:25942/127‑1 (MQ=255)
                                                    cAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTc                                                                                                                        >  3:25942/1‑127 (MQ=255)
                                                      aCTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGa                                                                                                                                          <  4:230518/107‑1 (MQ=255)
                                                      aCTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGa                                                                                                                                          >  3:230518/1‑107 (MQ=255)
                                                                cttcAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGGCTTTTGAGTTCTGCTAAGCTACGCTGGACGTGCCTTGGCGCGTTTCGAta                                                                                      >  3:6529/1‑148 (MQ=255)
                                                                    aaCCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGc                                                                                  >  2:33867/1‑149 (MQ=255)
                                                                       cTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCaa                                                                               >  8:104131/1‑149 (MQ=255)
                                                                                  ccGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTg                                                                    <  7:104131/149‑1 (MQ=255)
                                                                                   cGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTgg                                                                   <  3:130198/149‑1 (MQ=255)
                                                                                    gggTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGc                                                                  >  7:291203/1‑149 (MQ=255)
                                                                                     ggTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGc                                                                  >  6:167547/1‑148 (MQ=255)
                                                                                         aTCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTg                                                             >  6:103831/1‑149 (MQ=255)
                                                                                         aTCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTAGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTg                                                             >  5:236297/1‑149 (MQ=255)
                                                                                          tCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGc                                                            >  1:112210/1‑149 (MQ=255)
                                                                                             aGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCgg                                                         >  7:146223/1‑149 (MQ=255)
                                                                                             aGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCgg                                                         >  3:136977/1‑149 (MQ=255)
                                                                                               cTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGt                                                       <  4:136977/149‑1 (MQ=255)
                                                                                                  gggACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTa                                                      <  1:112081/147‑1 (MQ=255)
                                                                                                      cccAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTtctc                                                <  1:175856/149‑1 (MQ=255)
                                                                                                       ccAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCtcg                                               >  8:36015/1‑149 (MQ=255)
                                                                                                       ccAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCtcg                                               >  8:37227/1‑149 (MQ=255)
                                                                                                         aGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTa                                                      <  3:254186/140‑1 (MQ=255)
                                                                                                         aGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTa                                                      >  4:254186/1‑140 (MQ=255)
                                                                                                         aGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCtcgtc                                             <  1:130025/149‑1 (MQ=255)
                                                                                                         aGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCtcgtc                                             <  1:134024/149‑1 (MQ=255)
                                                                                                          ggTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCtcgtcg                                            <  1:33867/149‑1 (MQ=255)
                                                                                                                gaTCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCAAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAgtgg                                      >  2:41214/1‑149 (MQ=255)
                                                                                                                    tGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTtctc                                                >  8:286663/1‑135 (MQ=255)
                                                                                                                    tGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTtctc                                                <  7:286663/135‑1 (MQ=255)
                                                                                                                    tGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGGAGTTCTCGTCGAAGTGGTGcc                                  >  3:264370/1‑149 (MQ=255)
                                                                                                                          gggTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAg                            >  3:263362/1‑149 (MQ=255)
                                                                                                                           ggTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAgg                           >  7:303415/1‑149 (MQ=255)
                                                                                                                                taCAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGCAGtt                      >  7:77313/1‑149 (MQ=255)
                                                                                                                                taCAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGCAGt                       <  5:103831/148‑1 (MQ=255)
                                                                                                                                  cagcagGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGCAGttt                     <  2:112210/148‑1 (MQ=255)
                                                                                                                                   cgatctcTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGCAGttt                     <  7:273270/141‑1 (MQ=255)
                                                                                                                                   agcagGCTCAGGGCCGCAGAGCCGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGCAGttt                     >  8:207343/1‑147 (MQ=255)
                                                                                                                                      agGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGCAGTTTCAGgc                <  8:303415/149‑1 (MQ=255)
                                                                                                                                        gCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGCAGTTTCAGgcgc              >  7:202203/1‑149 (MQ=255)
                                                                                                                                        gCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGCAGTTTCAGgcgc              >  5:264333/1‑149 (MQ=255)
                                                                                                                                         cTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGCAGTTTCAGGCGCt             >  4:92801/1‑149 (MQ=255)
                                                                                                                                          tCAGGGCCGCAGAGCGGCCCTGGGCCTTGTCAGTTATGCTCAGCTACGCTGGACTTGCCGTGGCGCGTTTCGATCACGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGCAGTTTAAGGCGCTc            >  7:205965/1‑149 (MQ=255)
                                                                                                                                          tCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGCAGTTTCAGGCGCt             >  5:233384/1‑148 (MQ=255)
                                                                                                                                            aGGGCCGCTCTGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGaa                                          <  6:219941/117‑1 (MQ=255)
                                                                                                                                            aGGGCCGCTCTGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGaa                                          >  5:219941/1‑117 (MQ=255)
                                                                                                                                                 cGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGCAGTTTCAGGCGCTCACCTACt     >  4:44222/1‑149 (MQ=255)
                                                                                                                                                  gCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGCAGTTTc                    >  7:215061/1‑133 (MQ=255)
                                                                                                                                                  gCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGCAGTTTc                    <  8:215061/133‑1 (MQ=255)
                                                                                                                                                   cAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGCAGTTTc                    <  3:316558/132‑1 (MQ=255)
                                                                                                                                                   cAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGCAGTTTc                    >  4:316558/1‑132 (MQ=255)
                                                                                                                                                    agagCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGCAGTTTCAGGCGCTCACCTACTGcc  >  1:326539/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TGTCGCGCTGGCGGGGCAACCTGTTCGCGTTCTTGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGCAGTTTCAGGCGCTCACCTACTGCC  >  NC_002947/1376173‑1376469

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