Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,455,012 C→G 6.7% S36T (AGC→ACC) ‡ mfd ← transcription‑repair coupling factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,455,0120CG 6.7% 283.9 / 8.5 89S36T (AGC→ACC) ‡mfdtranscription‑repair coupling factor
Reads supporting (aligned to +/- strand):  ref base C (40/43);  new base G (3/3);  total (43/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.27e-01

GTGCGGCGAGAACAGGTCGTAGGGCAGGGTTTCCCAATCGGGGAACGGCAGCACGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTGCCGGCCGTGGCCGACATTTGCGGTAGGCGCAGAACTGACACGGGAGATTGCACTCCAAGCGTTGCGGCAAAGAGAC  >  NC_002947/2454867‑2455153
                                                                                                                                                 |                                                                                                                                              
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                                                                                                                                             <  1:35115/127‑1 (MQ=255)
                      ggCAGGGTTTCCCAATCGGGGAACGGCAGCACGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCgctgc                                                                                                                                             >  2:35115/1‑127 (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                                                                                                                            <  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                                                                                                                           >  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                                                                                    <  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)
                                                                      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                                                                            >  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)
                                                                                  tCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTGCCGGCCGTGGCCGACATTTg                                                          >  8:55954/1‑149 (MQ=255)
                                                                                  cccTGCTGCTGACGGCCGACAGCCAGGCCGCCGACCGTCTGGAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTGCCGGCCGTGGCCGACATTTg                                                          >  6:71013/2‑149 (MQ=25)
                                                                                  cccTGCTGCTGACGGCCGACAGCCAGGCCGCCGACCGTCTGGAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTGCCGGCCGTGGCCGACATTTg                                                          >  6:335939/2‑149 (MQ=25)
                                                                                   ccTGCTGCTGACGGCCGACAGCCAGGCCGCCGACCGTCTGGAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTGCCGGCCGTGGCCGACATTTGc                                                         >  6:71628/1‑149 (MQ=37)
                                                                                   ccTGCTGCTGACGGCCGACAGCCAGGCCGCCGACCGTCTGGAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTGCCGGCCGTGGCCGACATTTGc                                                         >  3:54745/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)
                                                                                                               gTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTGCCGGCCGTGGCCGACATTTGCGGTAGGCGCAGAACTGACACGGGAGAtt                             >  4:40361/1‑149 (MQ=255)
                                                                                                                tCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTGCCGGCCGTGGCCGACATTTGCGGTAGGCGCAGAACTGACACGGGAGATTg                            <  7:55954/149‑1 (MQ=255)
                                                                                                                  ggCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTGCCGGCCGATGGCCGACATTTGCGGTAGGCGCAGaa                                            <  1:9887/132‑1 (MQ=255)
                                                                                                                  ggCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTGCCGGCCGTGGCCGACATTTGCGGTAGGCGCAGaa                                            >  2:9887/1‑131 (MQ=255)
                                                                                                                  ggCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTGCCGGCCGTGGCCGACATTTGCGGTAGGCGCAGAACTGACACGGGAGATTGc                           <  3:267875/148‑1 (MQ=255)
                                                                                                                  ggCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTGCCGGCCGTGGCCGACATTTGCGGTAGGCGCAGAACTGACACGGGAGATTGCa                          <  2:151573/149‑1 (MQ=255)
                                                                                                                   gCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTGCCGGCCGTGGCCGACATTTGCGGTAGGCGCAGAACTGACACGGGAGATTGccc                         >  8:259752/1‑147 (MQ=255)
                                                                                                                      gTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCGGGTTTGTTTGCCGGCCGTGGCCGACATTTGCGGTAGGCGCAGAACTGACACGGGAGATTGCACTc                       >  8:152778/1‑148 (MQ=255)
                                                                                                                      gTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTGCCGGCCGTGGCCGACATTTGCGGTAGGCGCAGAACTGACACGGGAGATTGCACTc                       >  1:288127/1‑148 (MQ=255)
                                                                                                                       tCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTGCCGGCCGTGGCCGACATTTGCg                                                        >  3:20109/1‑114 (MQ=255)
                                                                                                                       tCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTGCCGGCCGTGGCCGACATTTGCg                                                        <  4:20109/114‑1 (MQ=255)
                                                                                                                            agcagGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTGCCGGCCGTGGCCGACATTTGCGGTAGGCGCAGAACTGACACGGGAGATTGCACTCCAAGCGtt               <  3:303931/149‑1 (MQ=255)
                                                                                                                              cagGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTGCCGGCCGTGGCCGACATTTGGGGTAGGCGCAGAACTGACACGGGAGATTGCACTCCAAGCGTTg              <  7:3936/149‑1 (MQ=255)
                                                                                                                              cagGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTGCCGGCCGTGGCCGACATTTGCGGTAGGCGCAGAACTGACACGGGAGATTGCACTCCAAGCGTTg              >  1:220760/1‑149 (MQ=255)
                                                                                                                               agGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTGCCGGCCGTGGCCGACATTTGCGGTAGGCGCAGAACTGACACGGGAGATTGCACTCCAAGCGTTGc             <  1:72280/149‑1 (MQ=255)
                                                                                                                                gggTGAAGCGGGCAGGGCTGCTGGCGGCTTCGGCGATGGCAAGGATGAGGGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTGCCGGCCGTGGCCGACATTTGCGGTAGGCGCAGAACTGACACGGGAGATTGCACTCCAAGCGTTGc             <  8:302447/148‑1 (MQ=255)
                                                                                                                                 gccgAAGCCGCCAGCAGCGCTGGCCGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTGCCGGCCGTGGCCGACATTTGCGGTAGGCGCAGAACTGACACGGGAGATTGCACTCCAAGCGTTGCgg           >  5:96873/4‑149 (MQ=255)
                                                                                                                                     aaGCGGCCAGCGCTGCTGGCGGCTTCGGCTATGGCAAGGCTGAGGGCGGCACCGGTCAGGTTGCCCCAGGTTTGTTTGCCGGCCGTGGCCGACATTTGCGGTAGGCGGAGAACTGACACGGGAGATTGCACTCCAAGCGTTGCGGCaaa       <  3:321893/149‑1 (MQ=255)
                                                                                                                                      agtggCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGGCCCAGGTTTGTTTGCCGGCCGTGGCCGACATTTGCGGTAGGCGCAGAACTGACACGGGAGATTGCACTCCAAGCGTTGCGGCa         >  4:10937/4‑146 (MQ=255)
                                                                                                                                       gCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTGCCGGCCGTGGCCGACATTTGCGGTAGGCGCAGAACTGACACGGGAGATTGCACTCCAAGCGTTGCGGCAAaga     <  1:88002/149‑1 (MQ=255)
                                                                                                                                          gCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTGCCGGCCGTGGCCGACATTTGCGGTAGGCGCAGAACTGACACGGGAGATTGCACTCCAAGCGTTGCGGCAAaga     >  6:1552/1‑146 (MQ=255)
                                                                                                                                          gCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTGCCGGCCGTGGCCGACATTTGCGGTAGGCGCAGAACTGACACGGGAGATTGCACTCCAAGCGTTGCGGCAAaga     <  5:1552/146‑1 (MQ=255)
                                                                                                                                          gCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTGCCGGCCGTGGCCGACATTTGCGGTAGGCGCAGAACTGACACGGGAGATTGCACTCCAAGCGTTGCGGCAAAGAGAc  <  2:288127/149‑1 (MQ=255)
                                                                                                                                           ccAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTGCCGGCCGTGGCCGACATTTGCGGTAGGCGCAGAACTGACACGGGAGATTGCACTCCAAGCGTTGCGGCa         <  3:343407/141‑1 (MQ=255)
                                                                                                                                           ccAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTGCCGGCCGTGGCCGACATTTGCGGTAGGCGCAGAACTGACACGGGAGATTGCACTCCAAGCGTTGCGGCa         >  4:343407/1‑141 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
GTGCGGCGAGAACAGGTCGTAGGGCAGGGTTTCCCAATCGGGGAACGGCAGCACGGGCAGGTCTGGGGCAAAGAAGCGCAGCTCCTGCTCCAGACGGTCGGCGGCCTGGCTGTCGGCCGTCAGCAGCAGGGTGAAGCGGCCAGCGCTGCTGGCGGCTTCGGCGATGGCAAGGCTGAGGGCGGCACCGGGCAGGTTGCCCCAGGTTTGTTTGCCGGCCGTGGCCGACATTTGCGGTAGGCGCAGAACTGACACGGGAGATTGCACTCCAAGCGTTGCGGCAAAGAGAC  >  NC_002947/2454867‑2455153

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: