Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,454,990 C→G 7.5% L43L (CTG→CTC)  mfd ← transcription‑repair coupling factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,454,9900CG 7.5% 301.5 / 11.4 93L43L (CTG→CTC) mfdtranscription‑repair coupling factor
Reads supporting (aligned to +/- strand):  ref base C (42/44);  new base G (4/3);  total (46/47)
Fisher's exact test for biased strand distribution p-value = 7.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

GATGCGCTGGGAGATGATGTCCTGGTGCGGCGAGAACAGGTCGTAGGGCAGGGTTTCCCAATCGGGGAACGGCAGCACGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTGCCGGCCG‑TGGCCGACATTTGCGGTAGGCGCAGAACTGACACGGGAGATTGCACTC  >  NC_002947/2454843‑2455132
                                                                                                                                                   |                                                                                                                                               
gATGCGCTGGGAGATGATGTCCTGGTGCGGCGAGAACAGGTCGTAGGGCAGGGTTTCCCAATCGGGGAACGGCAGCACGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTcagca                                                                                                                                                <  7:30625/149‑1 (MQ=255)
     ggtGGGTGATGATGTCCTGTTGCGGCGAGAACAGGTCGTAGGGGAGGTTTTCGCAATCGGGGTACGGCAGCACGGGCAGTTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAgg                                                                                                                                           <  4:42436/147‑1 (MQ=255)
     gcTGGGAGATGATGTCCTGGTGCGGCGAGAACAGGTCGTAGGGCAGGGTTTCCCAATCGGGGAACGGCAGCACGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAgg                                                                                                                                           <  2:76687/149‑1 (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                                                                                                                                     <  7:233508/149‑1 (MQ=255)
           agaTGATGTCCTGGTGCGGCGAGAACAGGTCGTAGGGCAGGGTTTCCCAATCGGGGAACGGCAGCACGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAg                                                                                                                                     <  3:256388/149‑1 (MQ=255)
            gatgatGTCCTGGTGCGGCGAGAACAGGTCGTAGGGCAGGGTTTCCCAATCGGGGAACGGCAGCACGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGaa                                                                                                                                      <  5:33127/147‑1 (MQ=255)
            gatgatGTCCTGGTGCGGCGAGAACAGGTCGTAGGGCAGGGTTTCCCAATCGGGGAACGGCAGCACGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGaa                                                                                                                                      >  6:33127/1‑147 (MQ=255)
               gatgTCCTGGTGCGGCGAGAACAGGTCGTAGGGCAGGGTTTCCCAATCGGGGAACGGCAGCACGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTcagca                                                                                                                                                <  7:307201/134‑1 (MQ=255)
               gatgTCCTGGTGCGGCGAGAACAGGTCGTAGGGCAGGGTTTCCCAATCGGGGAACGGCAGCACGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTcagca                                                                                                                                                >  8:307201/1‑134 (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)
                               gagaACAGGTCGTAGGGCAGGGTTTCCCAATCGGGGAACGGCAGCACGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTggcggc                                                                                                                 >  3:65364/1‑149 (MQ=255)
                               gagaACAGGTCGTAGGGCAGGGTTTCCCAATCGGGGAACGGCAGCACGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCCGCCAGCAGCGCTggccgc                                                                                                                 <  8:100046/149‑1 (MQ=255)
                                   aCAGGTCGTAGGGCAGGGTTTCCCAATCGGGGAACGGCAGCACGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCg                                                                                                             >  1:125554/1‑149 (MQ=255)
                                   aCAGGTCGTAGGGCAGGGTTTCCCAATCGGGGAACGGCAGCACGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGACTTCg                                                                                                             <  2:312463/149‑1 (MQ=255)
                                    cAGGTCGTAGGGCAGGGTTTCCCAATCGGGGAACGGCAGCACGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCgg                                                                                                            <  4:229350/149‑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                                                                                                                    <  3:196056/125‑1 (MQ=255)
                                                    gTTTCCCAATCGGGGAACGGCAGCACGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTggc                                                                                                                    >  4:196056/1‑125 (MQ=255)
                                                        cccAATCGGGGAACGGCAGCACGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATg                                                                                                       >  3:202998/1‑134 (MQ=255)
                                                        cccAATCGGGGAACGGCAGCACGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATg                                                                                                       <  4:202998/134‑1 (MQ=255)
                                                        cccAATCGGGGAACGGCAGCACGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGggc                                                                                         <  6:37083/148‑1 (MQ=255)
                                                         ccAATCGGGGAACGGCAGCACGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGggcgg                                                                                       <  4:285103/149‑1 (MQ=255)
                                                         ccAATCGGGGAACGGCAGCACGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGCTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGggcgg                                                                                       >  8:3936/1‑149 (MQ=255)
                                                           aaTCGGGGAACGGCAGCACGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATgg                                                                                                      >  8:306191/1‑132 (MQ=255)
                                                           aaTCGGGGAACGGCAGCACGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATgg                                                                                                      <  7:306191/132‑1 (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                                                               <  6:245196/149‑1 (MQ=255)
                                                                                 cAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGgttt                                                               <  1:260542/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)
                                                                                              aagaagCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGCGGGCGGCACCGGGCAGGTTGCCCCAGGTGTGTTTGACGTCCG‑t                                                 >  4:267875/1‑149 (MQ=255)
                                                                                               agaagCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTGCCGGCAG‑Tg                                                >  8:314153/1‑149 (MQ=255)
                                                                                                  cgcgcAGCTCCTGCTGCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTGCCGGCCG‑TGGcc                                             <  4:21257/148‑1 (MQ=255)
                                                                                                  agCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTGCCGGCCG‑TGGcc                                             >  8:58261/1‑149 (MQ=255)
                                                                                                      cAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTGCCGGCCG‑TGGCCGACa                                         <  8:89176/149‑1 (MQ=255)
                                                                                                         cTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTgccg                                                       >  2:292359/1‑133 (MQ=255)
                                                                                                         cTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTgccg                                                       <  1:292359/133‑1 (MQ=255)
                                                                                                          tCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTg                                                                                              >  7:99931/1‑93 (MQ=255)
                                                                                                          tCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTg                                                                                              <  8:99931/93‑1 (MQ=255)
                                                                                                          tCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTGCCGGCCG‑TGGCCGACATTTg                                     >  8:55954/1‑149 (MQ=255)
                                                                                                          cccTGCTGCTGACGGCCGACAGCCAGGCCGCCGACCGTCTGGAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTGCCGGCCG‑TGGCCGACATTTg                                     >  6:71013/2‑149 (MQ=25)
                                                                                                          cccTGCTGCTGACGGCCGACAGCCAGGCCGCCGACCGTCTGGAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTGCCGGCCG‑TGGCCGACATTTg                                     >  6:335939/2‑149 (MQ=25)
                                                                                                           ccTGCTGCTGACGGCCGACAGCCAGGCCGCCGACCGTCTGGAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTGCCGGCCG‑TGGCCGACATTTGc                                    >  6:71628/1‑149 (MQ=37)
                                                                                                           ccTGCTGCTGACGGCCGACAGCCAGGCCGCCGACCGTCTGGAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTGCCGGCCG‑TGGCCGACATTTGc                                    >  3:54745/1‑149 (MQ=37)
                                                                                                            ctgctgctgACGGCCGACAGCCAGGCCGCCGACCGTCTGGAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTGCCGGCCG‑TGGCCGACATTTGCg                                   <  5:71628/149‑1 (MQ=37)
                                                                                                               cTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAACCGCCAGCAGCGCTggc                                                                                                                    >  6:75877/1‑66 (MQ=255)
                                                                                                               cTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAACCGCCAGCAGCGCTggc                                                                                                                    <  5:75877/66‑1 (MQ=255)
                                                                                                               cTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGACAGCAGCAGGGTGAAGCCGCCAGCGCTCCTGGCGGCTTCGGCCATGGCAAGCCTGAGGGCCCCACCTGGCAGGTTGCCCCAGCTTAGTTTGCCGGCCC‑TGTCCGACATTAGCCGTa                                >  7:302447/1‑149 (MQ=255)
                                                                                                                tCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTGCCGGCCG‑TGGCCGACATTCGCGGTAg                               >  6:319079/1‑149 (MQ=255)
                                                                                                                  ctgACGGCCGACAGCCAGGCCGCCGACCGTCTGGAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTGCCGGCCG‑TGGCCGACATTTGCGGTAGgc                             <  4:54745/147‑1 (MQ=38)
                                                                                                                  cAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTGCCGGCCG‑TGGCCGACATTTGCGGTAGgc                             <  7:58261/149‑1 (MQ=255)
                                                                                                                        gTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTGCCGGCCG‑TGGCCGACATTCGCGGTAGGCGCAGaa                       <  5:319079/149‑1 (MQ=255)
                                                                                                                         ccgacagcCAGGCCGCCGACCGTCTGGAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTGCCGGCCG‑TGGCCGACATTTGCGGTAGGCGCAGAAc                      <  5:335939/143‑1 (MQ=255)
                                                                                                                          cggcggCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTGCCGGCCG‑TGGCCGACATTTGCGGTAGGCGCAGAACt                     <  5:115727/149‑1 (MQ=255)
                                                                                                                                  tGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTGCCGGCCG‑TGGCCGACATTTGCGGTAGGCGCAGAACTGACACgg              >  2:72280/1‑148 (MQ=255)
                                                                                                                                    gCTGTCGGCCGTCAGCACCACGCTGAAGCGGCCAGCGCTGCTGGCGGCTACGGCGATGGCAAGGCTGAGGGCGGCACCCGGCAGGTTGCCCCACGTTTGTTTGCCGCCCG‑TGGCCGACATTTGCGGTAGCCGCAGAACTGACACgg              >  4:321893/1‑146 (MQ=255)
                                                                                                                                      tGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGA‑GGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTGCCGGCCG‑TGGCCGACATTTGCGGTAGGCGCAGAACTGACACGGGAGAtt        >  4:303931/1‑149 (MQ=255)
                                                                                                                                       gTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTGCCGGCCG‑TGGCCGACATTTGCGGTAGGCGCAGAACTGACACGGGAGAtt        >  4:40361/1‑149 (MQ=255)
                                                                                                                                        tCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTGCCGGCCG‑TGGCCGACATTTGCGGTAGGCGCAGAACTGACACGGGAGATTg       <  7:55954/149‑1 (MQ=255)
                                                                                                                                          ggCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTGCCGGCCGATGGCCGACATTTGCGGTAGGCGCAGaa                       <  1:9887/132‑1 (MQ=255)
                                                                                                                                          ggCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTGCCGGCCG‑TGGCCGACATTTGCGGTAGGCGCAGaa                       >  2:9887/1‑131 (MQ=255)
                                                                                                                                          ggCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTGCCGGCCG‑TGGCCGACATTTGCGGTAGGCGCAGAACTGACACGGGAGATTGc      <  3:267875/148‑1 (MQ=255)
                                                                                                                                          ggCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTGCCGGCCG‑TGGCCGACATTTGCGGTAGGCGCAGAACTGACACGGGAGATTGCa     <  2:151573/149‑1 (MQ=255)
                                                                                                                                           gCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTGCCGGCCG‑TGGCCGACATTTGCGGTAGGCGCAGAACTGACACGGGAGATTGccc    >  8:259752/1‑147 (MQ=255)
                                                                                                                                              gTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCGGGTTTGTTTGCCGGCCG‑TGGCCGACATTTGCGGTAGGCGCAGAACTGACACGGGAGATTGCACTc  >  8:152778/1‑148 (MQ=255)
                                                                                                                                              gTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTGCCGGCCG‑TGGCCGACATTTGCGGTAGGCGCAGAACTGACACGGGAGATTGCACTc  >  1:288127/1‑148 (MQ=255)
                                                                                                                                               tCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTGCCGGCCG‑TGGCCGACATTTGCg                                   >  3:20109/1‑114 (MQ=255)
                                                                                                                                               tCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTGCCGGCCG‑TGGCCGACATTTGCg                                   <  4:20109/114‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
GATGCGCTGGGAGATGATGTCCTGGTGCGGCGAGAACAGGTCGTAGGGCAGGGTTTCCCAATCGGGGAACGGCAGCACGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTGCCGGCCG‑TGGCCGACATTTGCGGTAGGCGCAGAACTGACACGGGAGATTGCACTC  >  NC_002947/2454843‑2455132

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