Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,454,979 T→C 7.1% D47G (GAC→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,9790TC 7.1% 281.6 / 8.2 98D47G (GAC→GGC) mfdtranscription‑repair coupling factor
Reads supporting (aligned to +/- strand):  ref base T (43/48);  new base C (4/3);  total (47/51)
Fisher's exact test for biased strand distribution p-value = 7.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.45e-01

TACAGGCTGGCGATGCGCTGGGAGATGATGTCCTGGTGCGGCGAGAACAGGTCGTAGGGCAGGGTTTCCCAATCGGGGAACGGCAGCACGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTGCCGGCCGTGGCCGACATTTGCGGTAGGCGCAGAACTGACACGGGAGATTG  >  NC_002947/2454832‑2455127
                                                                                                                                                   |                                                                                                                                                    
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                                                                                                                                     <  3:338841/147‑1 (MQ=255)
                  tGGGAGATGATGTCCTGGTGCGGCGAGAACAGGTCGTAGGGCAGGGTTTCCCAATCGGGGAACGGCAGCACGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAgg                                                                                                                                     >  4:338841/1‑147 (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)
                                         cgAGAACAGGTCGTAGGGCAGGGTTTCCCAATCGGGGAACGGCAGCACGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTcag                                                                                                                                            <  1:326413/117‑1 (MQ=255)
                                         cgAGAACAGGTCGTAGGGCAGGGTTTCCCAATCGGGGAACGGCAGCACGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTcag                                                                                                                                            >  2:326413/1‑117 (MQ=255)
                                          gagaACAGGTCGTAGGGCAGGGTTTCCCAATCGGGGAACGGCAGCACGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTggcggc                                                                                                           >  3:65364/1‑149 (MQ=255)
                                          gagaACAGGTCGTAGGGCAGGGTTTCCCAATCGGGGAACGGCAGCACGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCCGCCAGCAGCGCTggccgc                                                                                                           <  8:100046/149‑1 (MQ=255)
                                              aCAGGTCGTAGGGCAGGGTTTCCCAATCGGGGAACGGCAGCACGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTcag                                                                                                                                            <  2:34626/112‑1 (MQ=255)
                                              aCAGGTCGTAGGGCAGGGTTTCCCAATCGGGGAACGGCAGCACGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTcag                                                                                                                                            >  1:34626/1‑112 (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                                                                                                           >  2:115742/1‑135 (MQ=255)
                                                        gggCAGGGTTTCCCAATCGGGGAACGGCAGCACGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCCGCCAGCAGCGCTggccgc                                                                                                           <  1:115742/135‑1 (MQ=255)
                                                         ggCAGGGTTTCCCAATCGGGGAACGGCAGCACGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCgctgc                                                                                                                  <  1:35115/127‑1 (MQ=255)
                                                         ggCAGGGTTTCCCAATCGGGGAACGGCAGCACGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCgctgc                                                                                                                  >  2:35115/1‑127 (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                                                                           <  5:72738/141‑1 (MQ=255)
                                                                                  gcagcaCGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCgg                                                                           >  6:72738/1‑141 (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                                                      <  4:320122/149‑1 (MQ=255)
                                                                                               gTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGgtttgtt                                                      >  7:89176/1‑149 (MQ=255)
                                                                                                  tGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTgc                                                   >  6:115727/1‑149 (MQ=255)
                                                                                                    gggCAAAGAAGCCCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCTGCACCGGGCAGGTTGCCCCAGGtttgtttg                                                    <  4:65364/146‑1 (MQ=255)
                                                                                                      gCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTgccggc                                               >  3:156354/1‑149 (MQ=255)
                                                                                                      gCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTgccggc                                               >  2:88002/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                                                                                        >  7:99931/1‑93 (MQ=255)
                                                                                                                     tCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTg                                                                                        <  8:99931/93‑1 (MQ=255)
                                                                                                                     tCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTGCCGGCCGTGGCCGACATTTg                                >  8:55954/1‑149 (MQ=255)
                                                                                                                     cccTGCTGCTGACGGCCGACAGCCAGGCCGCCGACCGTCTGGAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTGCCGGCCGTGGCCGACATTTg                                >  6:335939/2‑149 (MQ=25)
                                                                                                                     cccTGCTGCTGACGGCCGACAGCCAGGCCGCCGACCGTCTGGAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTGCCGGCCGTGGCCGACATTTg                                >  6:71013/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)
                                                                                                                                                 tGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGA‑GGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTGCCGGCCGTGGCCGACATTTGCGGTAGGCGCAGAACTGACACGGGAGAtt   >  4:303931/1‑149 (MQ=255)
                                                                                                                                                  gTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTGCCGGCCGTGGCCGACATTTGCGGTAGGCGCAGAACTGACACGGGAGAtt   >  4:40361/1‑149 (MQ=255)
                                                                                                                                                   tCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTGCCGGCCGTGGCCGACATTTGCGGTAGGCGCAGAACTGACACGGGAGATTg  <  7:55954/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TACAGGCTGGCGATGCGCTGGGAGATGATGTCCTGGTGCGGCGAGAACAGGTCGTAGGGCAGGGTTTCCCAATCGGGGAACGGCAGCACGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTGCCGGCCGTGGCCGACATTTGCGGTAGGCGCAGAACTGACACGGGAGATTG  >  NC_002947/2454832‑2455127

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