Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,454,977 T→C 7.3% S48G (AGC→GGC)  mfd ← transcription‑repair coupling factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,454,9770TC 7.3% 274.9 / 8.6 96S48G (AGC→GGC) mfdtranscription‑repair coupling factor
Reads supporting (aligned to +/- strand):  ref base T (41/48);  new base C (4/3);  total (45/51)
Fisher's exact test for biased strand distribution p-value = 7.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

CGGTACAGGCTGGCGATGCGCTGGGAGATGATGTCCTGGTGCGGCGAGAACAGGTCGTAGGGCAGGGTTTCCCAATCGGGGAACGGCAGCACGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTGCCGGCCGTGGCCGACATTTGCGGTAGGCGCAGAACTGACACGGGAGATT  >  NC_002947/2454829‑2455126
                                                                                                                                                    |                                                                                                                                                     
cGGTACAGGCTGGCGATGCGCTGGGAGATGATGTCCTGGTGCGGCGAGAACAGGTCGTAGGGCAGGGTTTCCCAATCGGGGAACGGCAGCACGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCt                                                                                                                                                       <  2:191557/149‑1 (MQ=255)
   tACAGGCTGGCGATGCGCTGGGAGATGATGTCCTGGTGCGGCGAGAACAGGTCGTAGGGCAGGGTTTCCCAATCGGGGAACGGCAGCACGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGCCTGTc                                                                                                                                                    >  3:21257/1‑149 (MQ=255)
            gCGATGCGCTGGGAGATGATGTCCTGGTGCGGCGAGAACAGGTCGTAGGGCAGGGTTTCCCAATCGGGGAACGGCAGCACGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTcag                                                                                                                                           <  3:335313/149‑1 (MQ=255)
            gCGATGCGCTGGGAGATGATGTCCTGGTGCGGCGAGAACAGGTCGTAGGGCAGGGTTTCCCAATCGGGGAACGGCAGCACGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTca                                                                                                                                            <  4:48500/148‑1 (MQ=255)
            gCGATGCGCTGGGAGATGATGTCCTGGTGCGGCGAGAACAGGTCGTAGGGCAGGGTGTCCCAATCGGGGAACGGCAGCACGGGCAGGTCTGCGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTc                                                                                                                                             <  2:263271/147‑1 (MQ=255)
              gATGCGCTGGGAGATGATGTCCTGGTGCGGCGAGAACAGGTCGTAGGGCAGGGTTTCCCAATCGGGGAACGGCAGCACGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTcagca                                                                                                                                         <  7:30625/149‑1 (MQ=255)
                   ggtGGGTGATGATGTCCTGTTGCGGCGAGAACAGGTCGTAGGGGAGGTTTTCGCAATCGGGGTACGGCAGCACGGGCAGTTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAgg                                                                                                                                    <  4:42436/147‑1 (MQ=255)
                   gcTGGGAGATGATGTCCTGGTGCGGCGAGAACAGGTCGTAGGGCAGGGTTTCCCAATCGGGGAACGGCAGCACGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAgg                                                                                                                                    <  2:76687/149‑1 (MQ=255)
                    cTGGGAGATGATGTCCTGGTGCGGCGAGAACAGGTCGTAGGGCAGGGTTTCCCAATCGGGGAACGGCAGCACGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTc                                                                                                                                             <  5:72886/139‑1 (MQ=255)
                    cTGGGAGATGATGTCCTGGTGCGGCGAGAACAGGTCGTAGGGCAGGGTTTCCCAATCGGGGAACGGCAGCACGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTc                                                                                                                                             >  6:72886/1‑139 (MQ=255)
                     tGGGAGATGATGTCCTGGTGCGGCGAGAACAGGTCGTAGGGCAGGGTTTCCCAATCGGGGAACGGCAGCACGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAgg                                                                                                                                    >  4:338841/1‑147 (MQ=255)
                     tGGGAGATGATGTCCTGGTGCGGCGAGAACAGGTCGTAGGGCAGGGTTTCCCAATCGGGGAACGGCAGCACGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAgg                                                                                                                                    <  3:338841/147‑1 (MQ=255)
                     tGGGAGATGATGTCCTGGTGCGGCGAGAACAGGTCGTAGGGCAGGGTTTCCCAATCGGGGAACGGCAGCACGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGt                                                                                                                                  >  1:138524/1‑149 (MQ=255)
                      gggAGATGATGTCCTGGTGCGGCGAGAACAGGTCGTAGGGCAGGGTTTCCCAATCGGGGAACGGCAGCACGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGt                                                                                                                                  >  1:312463/1‑148 (MQ=255)
                       ggAGATGATGTCCTGGTGCGGCGAGAACAGGTCGTAGGGCAGGGTTTCCCAATCGGGGAACGGCAGCACGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTg                                                                                                                                 <  6:108751/148‑1 (MQ=255)
                         agaTGATGTCCTGGTGCGGCGAGAACAGGTCGTAGGGCAGGGTTTCCCAATCGGGGAACGGCAGCACGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAg                                                                                                                              <  3:256388/149‑1 (MQ=255)
                         agaTGATGTCCTGGTGCGGCGAGAACAGGTCGTAGGGCAGGGTTTCCCAATCGGGGAACGGCAGCACGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAg                                                                                                                              <  7:233508/149‑1 (MQ=255)
                          gatgatGTCCTGGTGCGGCGAGAACAGGTCGTAGGGCAGGGTTTCCCAATCGGGGAACGGCAGCACGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGaa                                                                                                                               <  5:33127/147‑1 (MQ=255)
                          gatgatGTCCTGGTGCGGCGAGAACAGGTCGTAGGGCAGGGTTTCCCAATCGGGGAACGGCAGCACGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGaa                                                                                                                               >  6:33127/1‑147 (MQ=255)
                             gatgTCCTGGTGCGGCGAGAACAGGTCGTAGGGCAGGGTTTCCCAATCGGGGAACGGCAGCACGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTcagca                                                                                                                                         >  8:307201/1‑134 (MQ=255)
                             gatgTCCTGGTGCGGCGAGAACAGGTCGTAGGGCAGGGTTTCCCAATCGGGGAACGGCAGCACGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTcagca                                                                                                                                         <  7:307201/134‑1 (MQ=255)
                             gatgTCCTGGTGCGGCGAGAACAGGTCGTAGGGCAGGGTTTCCCAATCGGGGAACGGCAGCACGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGc                                                                                                                          >  7:299207/1‑149 (MQ=255)
                                gTCCTGGTGCGGCGAGAACAGGTCGTAGGGCAGGGTTTCCCAATCGGGGAACGGCAGCACGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAg                                                                                                                       >  5:37083/1‑149 (MQ=255)
                                 tCCTGGTGCGGCGAGAACAGGTCGTAGGGCAGGGTTTCCCAATCGGGGAACGGCAGCACGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAgc                                                                                                                      <  4:49279/149‑1 (MQ=255)
                                      gTGCGGCGAGAACAGGTCGTAGGGCAGGGTTTCCCAATCGGTGAACGGCAGCACGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCgctgc                                                                                                                 >  7:6754/1‑149 (MQ=255)
                                            cgAGAACAGGTCGTAGGGCAGGGTTTCCCAATCGGGGAACGGCAGCACGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTcag                                                                                                                                           >  2:326413/1‑117 (MQ=255)
                                            cgAGAACAGGTCGTAGGGCAGGGTTTCCCAATCGGGGAACGGCAGCACGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTcag                                                                                                                                           <  1:326413/117‑1 (MQ=255)
                                             gagaACAGGTCGTAGGGCAGGGTTTCCCAATCGGGGAACGGCAGCACGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTggcggc                                                                                                          >  3:65364/1‑149 (MQ=255)
                                             gagaACAGGTCGTAGGGCAGGGTTTCCCAATCGGGGAACGGCAGCACGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCCGCCAGCAGCGCTggccgc                                                                                                          <  8:100046/149‑1 (MQ=255)
                                                 aCAGGTCGTAGGGCAGGGTTTCCCAATCGGGGAACGGCAGCACGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTcag                                                                                                                                           >  1:34626/1‑112 (MQ=255)
                                                 aCAGGTCGTAGGGCAGGGTTTCCCAATCGGGGAACGGCAGCACGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTcag                                                                                                                                           <  2:34626/112‑1 (MQ=255)
                                                 aCAGGTCGTAGGGCAGGGTTTCCCAATCGGGGAACGGCAGCACGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCg                                                                                                      >  1:125554/1‑149 (MQ=255)
                                                 aCAGGTCGTAGGGCAGGGTTTCCCAATCGGGGAACGGCAGCACGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGACTTCg                                                                                                      <  2:312463/149‑1 (MQ=255)
                                                  cAGGTCGTAGGGCAGGGTTTCCCAATCGGGGAACGGCAGCACGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCgg                                                                                                     <  4:229350/149‑1 (MQ=255)
                                                     gTCGTAGGGCAGGGTTTCCCAATCGGGGAACGGCAGCACGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCg                                                                                                                                                   >  8:209283/1‑100 (MQ=255)
                                                     gTCGTAGGGCAGGGTTTCCCAATCGGGGAACGGCAGCACGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCg                                                                                                                                                   <  7:209283/100‑1 (MQ=255)
                                                      tCGTAGGGCAGGGTTTCCCAATCGGTGAACGGCAGCACGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGAt                                                                                                 <  8:6754/149‑1 (MQ=255)
                                                       cGTAGGGCAGGGTTTCCCAATCGGGGAACGGCAGCACGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATg                                                                                                <  1:329631/149‑1 (MQ=255)
                                                           gggCAGGGTTTCCCAATCGGGGAACGGCAGCACGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCCGCCAGCAGCGCTggccgc                                                                                                          <  1:115742/135‑1 (MQ=255)
                                                           gggCAGGGTTTCCCAATCGGGGAACGGCAGCACGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCCGCCAGCAGCGCTggccgc                                                                                                          >  2:115742/1‑135 (MQ=255)
                                                            ggCAGGGTTTCCCAATCGGGGAACGGCAGCACGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCgctgc                                                                                                                 >  2:35115/1‑127 (MQ=255)
                                                            ggCAGGGTTTCCCAATCGGGGAACGGCAGCACGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCgctgc                                                                                                                 <  1:35115/127‑1 (MQ=255)
                                                             gCAGGGTTTCCCAATCGGGGAACGGCAGCACGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCg                                                                                                                            <  4:246194/115‑1 (MQ=255)
                                                                 ggTTTCCCAATCGGGGAACGGCAGCACGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGa                                                                                      <  2:38325/149‑1 (MQ=255)
                                                                 ggTGTCCCAATCGGGGAACGGCAGCACGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGa                                                                                      <  6:232790/149‑1 (MQ=255)
                                                                  gTTTCCCAATCGGGGAACGGCAGCACGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTggc                                                                                                             >  4:196056/1‑125 (MQ=255)
                                                                  gTTTCCCAATCGGGGAACGGCAGCACGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTggc                                                                                                             <  3:196056/125‑1 (MQ=255)
                                                                      cccAATCGGGGAACGGCAGCACGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATg                                                                                                <  4:202998/134‑1 (MQ=255)
                                                                      cccAATCGGGGAACGGCAGCACGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATg                                                                                                >  3:202998/1‑134 (MQ=255)
                                                                      cccAATCGGGGAACGGCAGCACGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGggc                                                                                  <  6:37083/148‑1 (MQ=255)
                                                                       ccAATCGGGGAACGGCAGCACGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGggcgg                                                                                <  4:285103/149‑1 (MQ=255)
                                                                       ccAATCGGGGAACGGCAGCACGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGCTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGggcgg                                                                                >  8:3936/1‑149 (MQ=255)
                                                                         aaTCGGGGAACGGCAGCACGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATgg                                                                                               <  7:306191/132‑1 (MQ=255)
                                                                         aaTCGGGGAACGGCAGCACGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATgg                                                                                               >  8:306191/1‑132 (MQ=255)
                                                                             ggggAACGGCAGCACGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCg                                                                           <  8:299207/148‑1 (MQ=255)
                                                                              gggAACGGCAGCACGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCggg                                                                         >  3:320122/1‑149 (MQ=255)
                                                                               ggAACGGCAGCACGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCggg                                                                         >  5:230046/1‑148 (MQ=255)
                                                                                     gcagcaCGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCgg                                                                          >  6:72738/1‑141 (MQ=255)
                                                                                     gcagcaCGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCgg                                                                          <  5:72738/141‑1 (MQ=255)
                                                                                      cagcaCGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGc                                                                 >  5:245196/1‑149 (MQ=255)
                                                                                       agcaCGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGcc                                                                <  2:125554/149‑1 (MQ=255)
                                                                                         caCGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGccc                                                               >  1:151573/1‑148 (MQ=255)
                                                                                               cAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGgttt                                                        <  1:260542/149‑1 (MQ=255)
                                                                                               cAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGgttt                                                        <  6:245196/149‑1 (MQ=255)
                                                                                                aGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGgtttg                                                       <  6:230046/149‑1 (MQ=255)
                                                                                                  gTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGgtttgtt                                                     >  7:89176/1‑149 (MQ=255)
                                                                                                  gTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGgtttgtt                                                     <  4:320122/149‑1 (MQ=255)
                                                                                                     tGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTgc                                                  >  6:115727/1‑149 (MQ=255)
                                                                                                       gggCAAAGAAGCCCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCTGCACCGGGCAGGTTGCCCCAGGtttgtttg                                                   <  4:65364/146‑1 (MQ=255)
                                                                                                         gCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTgccggc                                              >  2:88002/1‑149 (MQ=255)
                                                                                                         gCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTgccggc                                              >  3:156354/1‑149 (MQ=255)
                                                                                                            aagaagCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGCGGGCGGCACCGGGCAGGTTGCCCCAGGTGTGTTTGACGTCCGt                                           >  4:267875/1‑149 (MQ=255)
                                                                                                             agaagCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTGCCGGCAGTg                                          >  8:314153/1‑149 (MQ=255)
                                                                                                                cgcgcAGCTCCTGCTGCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTGCCGGCCGTGGcc                                       <  4:21257/148‑1 (MQ=255)
                                                                                                                agCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTGCCGGCCGTGGcc                                       >  8:58261/1‑149 (MQ=255)
                                                                                                                    cAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTGCCGGCCGTGGCCGACa                                   <  8:89176/149‑1 (MQ=255)
                                                                                                                       cTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTgccg                                                <  1:292359/133‑1 (MQ=255)
                                                                                                                       cTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTgccg                                                >  2:292359/1‑133 (MQ=255)
                                                                                                                        tCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTg                                                                                       <  8:99931/93‑1 (MQ=255)
                                                                                                                        tCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTg                                                                                       >  7:99931/1‑93 (MQ=255)
                                                                                                                        tCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTGCCGGCCGTGGCCGACATTTg                               >  8:55954/1‑149 (MQ=255)
                                                                                                                        cccTGCTGCTGACGGCCGACAGCCAGGCCGCCGACCGTCTGGAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTGCCGGCCGTGGCCGACATTTg                               >  6:71013/2‑149 (MQ=25)
                                                                                                                        cccTGCTGCTGACGGCCGACAGCCAGGCCGCCGACCGTCTGGAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTGCCGGCCGTGGCCGACATTTg                               >  6:335939/2‑149 (MQ=25)
                                                                                                                         ccTGCTGCTGACGGCCGACAGCCAGGCCGCCGACCGTCTGGAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTGCCGGCCGTGGCCGACATTTGc                              >  3:54745/1‑149 (MQ=37)
                                                                                                                         ccTGCTGCTGACGGCCGACAGCCAGGCCGCCGACCGTCTGGAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTGCCGGCCGTGGCCGACATTTGc                              >  6:71628/1‑149 (MQ=37)
                                                                                                                          ctgctgctgACGGCCGACAGCCAGGCCGCCGACCGTCTGGAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTGCCGGCCGTGGCCGACATTTGCg                             <  5:71628/149‑1 (MQ=37)
                                                                                                                             cTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAACCGCCAGCAGCGCTggc                                                                                                             >  6:75877/1‑66 (MQ=255)
                                                                                                                             cTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAACCGCCAGCAGCGCTggc                                                                                                             <  5:75877/66‑1 (MQ=255)
                                                                                                                             cTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGACAGCAGCAGGGTGAAGCCGCCAGCGCTCCTGGCGGCTTCGGCCATGGCAAGCCTGAGGGCCCCACCTGGCAGGTTGCCCCAGCTTAGTTTGCCGGCCCTGTCCGACATTAGCCGTa                          >  7:302447/1‑149 (MQ=255)
                                                                                                                              tCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTGCCGGCCGTGGCCGACATTCGCGGTAg                         >  6:319079/1‑149 (MQ=255)
                                                                                                                                ctgACGGCCGACAGCCAGGCCGCCGACCGTCTGGAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTGCCGGCCGTGGCCGACATTTGCGGTAGgc                       <  4:54745/147‑1 (MQ=38)
                                                                                                                                cAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTGCCGGCCGTGGCCGACATTTGCGGTAGgc                       <  7:58261/149‑1 (MQ=255)
                                                                                                                                      gTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTGCCGGCCGTGGCCGACATTCGCGGTAGGCGCAGaa                 <  5:319079/149‑1 (MQ=255)
                                                                                                                                       ccgacagcCAGGCCGCCGACCGTCTGGAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTGCCGGCCGTGGCCGACATTTGCGGTAGGCGCAGAAc                <  5:335939/143‑1 (MQ=255)
                                                                                                                                        cggcggCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTGCCGGCCGTGGCCGACATTTGCGGTAGGCGCAGAACt               <  5:115727/149‑1 (MQ=255)
                                                                                                                                                tGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTGCCGGCCGTGGCCGACATTTGCGGTAGGCGCAGAACTGACACgg        >  2:72280/1‑148 (MQ=255)
                                                                                                                                                  gCTGTCGGCCGTCAGCACCACGCTGAAGCGGCCAGCGCTGCTGGCGGCTACGGCGATGGCAAGGCTGAGGGCGGCACCCGGCAGGTTGCCCCACGTTTGTTTGCCGCCCGTGGCCGACATTTGCGGTAGCCGCAGAACTGACACgg        >  4:321893/1‑146 (MQ=255)
                                                                                                                                                    tGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTGCCGGCCGTGGCCGACATTTGCGGTAGGCGCAGAACTGACACGGGAGAtt  >  4:303931/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CGGTACAGGCTGGCGATGCGCTGGGAGATGATGTCCTGGTGCGGCGAGAACAGGTCGTAGGGCAGGGTTTCCCAATCGGGGAACGGCAGCACGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTGCCGGCCGTGGCCGACATTTGCGGTAGGCGCAGAACTGACACGGGAGATT  >  NC_002947/2454829‑2455126

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: