Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 805,138 G→C 6.7% S731C (TCC→TGC) ‡ PP_0693 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947805,1380GC 6.7% 241.6 / 5.2 74S731C (TCC→TGC) ‡PP_0693hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (34/35);  new base C (2/3);  total (36/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.77e-01

TTGTCCTGATCGGCACTCTGGGTCAGGTCGAAG‑ACCGACAGGGTTTCGCCTTCTTCGTACTGGCCGGCACTGTGGATAATCTGGTCGATGTGCTTGAGCGCTTCCTTGTTCGGCGCCAGTACCACGCGGAAGCTTTCCAGGTTGGAAACCTGGCGCTTGTTGATCTCGCGGTTGAACAGGGCCAGCTGGTGCTCGAGGCTGTCGTAGTCGCTGCGGATGTTGCGCAGGGTTCGGGCGATGTCGGTGACCGCAGCGCGGCGTGCCTTGGCCAAGGTCAGGGCTTCGTCGGTACG  >  NC_002947/804994‑805286
                                                                                                                                                 |                                                                                                                                                    
ttgtCCTGATCGGCACTCTGGGTCAGGTCGAAG‑ACCGACAGGGTTTCGCCTTCTTCGTACTGGCCGGCACTGTGGATAATCTGGTCGATGTGCTTGAGCGCTTCCTTGTTCGGCGCCAGTACCACGCGGAAGCTTTCCAGGTTGGAAAc                                                                                                                                                  >  4:131579/1‑149 (MQ=255)
   tCCTGATCGGGACTCTGGGTCAGGTCGAAG‑ACCGACAGGGTTTCGCCTTCTTCGTACTGGCCGGCACTGTGGATAATCTGGTCGATGTGCTTGAGCGCTTCCTTGTTCGGCGCCAGTACCACGCGGAAGCTTTCCAGGTTGGAAACCt                                                                                                                                                <  6:215141/148‑1 (MQ=255)
   tCCTGATCGGCACTCTGGGTCAGGTCGAAG‑ACCGACAGGGTTTCGCCTTCTTCGTACTGGCCGGCACTGTGGATAATCTGGTCGATGTGCTTGAGCGCTTCCTTGTTCGGCGCCAGTACCACGCGGAAGCTTTCCAGGTTGGAAACCTg                                                                                                                                               >  6:33725/1‑149 (MQ=255)
     cTGATCGGCACTCTGGGTCAGGTCGAAG‑ACCGACAGGGTTTCGCCTTCTTCGTACTGGCCGGCACTGTGGATAATCTGGTCGATGTGCTTGAGCGCTTCCTTGTTCGGCGCCAGTACCACGCGGAAGCTTTCCAGGTTGGAAACCTgg                                                                                                                                              <  2:215068/148‑1 (MQ=255)
     cTGATCGGCACTCTGGGTCAGGTCGAAG‑ACCGACAGGGTTTCGCCTTCTTCGTACTGGCCGGCACTGTGGATAATCTGGTCGATGTGCTTGAGCGCTTCCTTGTTCGGCGCCAGTACCACGCGGAAGCTTTCCAGGTTGGAAACCTgg                                                                                                                                              <  3:131579/148‑1 (MQ=255)
     cTGATCGGCACTCTGGGTCAGGTCGAAG‑ACCGACAGGGTTTCGCCTTCTTCGTACTGGCCGGCACTGTGGATAATCTGGTCGATGTGCTTGAGCGCTTCCTTGTTCGGCGCCAGTACCACGCGGAAGCTTTCCAGGTTGGAAACCTgg                                                                                                                                              <  6:176514/148‑1 (MQ=255)
     cTGATCGGCACTCTGGGTCAGGTCGAAG‑ACCGACAGGGTTTCGCCTTCTTCGTACTGGCCGGCACTGTGGATAATCTGGTCGATGTGCTTGAGCGCTTCCTTGTTCGGCGCCAGTACCACGCGGAAGCTTTCCAGGTTGGAAACCTGgc                                                                                                                                             <  3:65791/149‑1 (MQ=255)
      tGATCGGCACTCTGGGTCAGGTCGAAG‑ACCGACAGGGTTTCGCCTTCTTCGTACTGGCCGGCACTGTGGATAATCTGGTCGATGTGCTTGAGCGCTTCCTTGTTCGGCGCCAGTACCACGCGGAAGCTTTCCAGGTTTCCAACCTgga                                                                                                                                             >  1:313494/1‑147 (MQ=255)
      tGATCGGCACTCTGGGTCAGGTCGAAG‑ACCGACAGGGTTTCGCCTTCTTCGTACTGGCCGGCACTGTGGATAATCTGGTCGATGTGCTTGAGCGCTTCCTTGTTCGGCGCCAGTACCACGCGGAAGCTTTCCAGGTTGGAAACCTGgc                                                                                                                                             >  2:314174/1‑148 (MQ=255)
        aTCGGCACTCTGGGTCAGGTCGAAGAACCGACAGGGTTTCGCCTTCTTCGTACTGGCCGGCACTGTGGATAATCTGGTCGATGTGCTTGAGCGCTTCCTTGTTCGGCGCCAGTACCACGCGGAAGCTTTCCAGGTTGGATACCTGgcgc                                                                                                                                           >  1:3568/1‑149 (MQ=255)
        aTCGGCACTCTGGGTCAGGTCGAAG‑ACCGACAGGGTTTCGCCTTCTTCGTACTGGCCGGCACTGTGGATAATCTGGTCGATGTGCTTGAGCGCTTCCTTGTTCGGCGCCAGTACCACGCGGAAGCTTTCCAGGTTTCCAACCTggaaa                                                                                                                                           <  2:313494/148‑4 (MQ=255)
        aTCGGCACTCTGGGTCAGGTCGAAG‑ACCGACAGGGTTTCGCCTTCTTCGTACTGGCCGGCACTGTGGATAATCTGGTCGATGTGCTTGAGCGCTTCCTTGTTCGGCGCCAGTACCACGCGGAAGCTTTCCAGGTTTCCAACCTggaaa                                                                                                                                           <  2:31949/148‑4 (MQ=255)
        aTCGGCACTCTGGGTCAGGTCGAAG‑ACCGACAGGGTTTCGCCTTCTTCGTACTGGCCGGCACTGTGGATAATCTGGTCGATGTGCTTGAGCGCTTCCTTGTTCGGCGCCAGTACCACGCGGAAGCTTTCCAGGTTTCCAACCTggaaa                                                                                                                                           <  5:222242/148‑4 (MQ=255)
        aTCGGCACTCTGGGTCAGGTCGAAG‑ACCGACAGGGTTTCGCCTTCTTCGTACTGGCCGGCACTGTGGATAATCTGGTCGATGTGCTTGAGCGCTTCCTTGTTCGGCGCCAGTACCACGCGGAAGCTTTCCAGGTTGGAAACCTGgcgc                                                                                                                                           <  6:90061/148‑1 (MQ=255)
           ggCACTCTGGGTCAGGTCGAAG‑ACCGACAGGGTTTCGCCTTCTTCGTACTGGCCGGCACTGTGGATAATCTGGTCGATGTGCTTGAGCGCTTCCTTGTTCGGCGCCAGTACCACGCGGAAGCTTTCCAGGTTGGAAACCTGGCGCttgt                                                                                                                                       <  1:51851/149‑1 (MQ=255)
            gCACTCTGGGTCAGGTCGAAG‑CCCGACAGGGTTTCGCCTTCTTCGTACTGGCCGGCACTGTGGATAATCTGGTCGATGTGCTTGAGCGCTTCCTTGTTCGGCGCCAGTACCACGCGGAAGCTTTCCAGGTTGGAAACCTGGCGCttgtt                                                                                                                                      >  6:85245/1‑149 (MQ=255)
             cACTCTGGGTCAGGTCGAAG‑ACCGACAGGGTTTCGCCTTCTTCGTACTGGCCGGCACTGTGGATAATCTGGTCGATGTGCTTGAGCGCTTCCTTGTTCGGCGCCAGTACCACGCGGAAGCTTTCCAGGTTGGAAACCTGGCGCttgttg                                                                                                                                     >  8:161401/1‑149 (MQ=255)
               ctctGGGTGAGGTCTAAG‑ACCGACAGGGGTTCGCCTTCTTCGTACTGGCCGGCACTGGGGATAATCTGGTCGATGTGCTTGAGCGCTTCCTTGTTCGGCTCCAGTACCACGCGGAAGCTTTCCAGGTTGGAAACCTGGCGCTTGTTGAt                                                                                                                                   <  2:289036/149‑1 (MQ=255)
                tctGGGTCAGGTCGAAG‑ACCGACAGGGTTTCGCCTTCTTCGTACTGGCCGGCACTGTGGATAATCTGGTCGATGTGCTTGAGCGCTTCCTTGTTCGGCGCCAGTACCACGCGGAAGCTTTCCAGGTTGGAAACCTGGCGCTTGTTGAtc                                                                                                                                  >  6:341237/1‑149 (MQ=255)
                    ggTCAGGTCGAAG‑ACCGACAGGGTTTCGCCTTCTTCGTACTGGCCGGCACTGTGGATAATCTGGTCGATGTGCTTGAGCGCTTCCTTGTTCGGCGCCAGTACCACGCGGAAGCTTTCCAGGTTGGAAACCTGGCGCTTGTTGATCTcgc                                                                                                                              >  8:227281/1‑149 (MQ=255)
                           tCGAAG‑ACCGACAGGGTTTCGCCTTCTTCGTACTGGCCGGCACTGTGGATAATCTGGTCGATGTGCTTGAGCGCTTCCTTGTTCGGCGCCAGTACCACGCGGAAGCTTTCCAGGTTGGAAACCTGGCGCTTGTTGATCTCGCGGTTGaa                                                                                                                       <  4:106417/149‑1 (MQ=255)
                             gAAG‑ACCGACAGGGTTTCGCCTTCTTCGTACTGGCCGGCACGGTGGATAATCTGGTCGATGTGCTTGAGCGCTTCCTTGTTCGGCGCCAGTACCACGCGGAAGCTTTCCAGGTTGGAAACCTGGCGCTTGTTGATCTCGCgg                                                                                                                            >  2:44609/1‑142 (MQ=255)
                             gAAG‑ACCGACAGGGTTTCGCCTTCTTCGTACTGGCCGGCACGGTGGATAATCTGGTCGATGTGCTTGAGCGCTTCCTTGTTCGGCGCCAGTACCACGCGGAAGCTTTCCAGGTTGGAAACCTGGCGCTTGTTGATCTCGCgg                                                                                                                            <  1:44609/142‑1 (MQ=255)
                                  aCCGACAGGGTTTCGCCTTCTTCGTACTGGCCGGCACTGTGGATAATCTGGTCGATGTGCTTGAGCGCTTCCTTGTTCGGCGCCAGTACCACGCGGAAGCTTTCCAGGTTGGAAACCTGGCGCTTGTTGATCTCGCGGTTGAACAGGGc                                                                                                                 <  5:341237/149‑1 (MQ=255)
                                       cAGGGTTTCGCCTTCTTCGTACTGGCCGGCACTGTGGATAATCTGGTCGATGTGCTTGAGCGCTTCCTTGTTCGGCGCCAGTACCACGCGGAAGCTTTCCAGGTTGGAAACCTGGCGCTTGTTGATCTCGCGGTTGAACAGGGCCAGCt                                                                                                            <  3:137494/149‑1 (MQ=255)
                                           gTTTCGCCTTCTTCGTACTGGCCGGCACTGTGGATAATCTGGTCGATGTGCTTGAGCGCTTCCTTGTTCGGCGCCAGTACCACGCGGAAGCTTTCCAGGTTGGAAACCTGGCGCTTGTTGATCTCGCGGTTGAACAGGGCCAGCTGGTg                                                                                                        <  5:85245/149‑1 (MQ=255)
                                                gCCTTCTTCGTACTGGCCGGCACTGTGGATAATCTGGTCGATGTGCTTGAGCGCTTCCTTGTTCGGCGCCAGTACCACGCGGAAGCTTTCCAGGTTGGAAACCTGGCGCTTGTTGATCTCGCGGTTGAACAGGGCCAGCTGGTGCTc                                                                                                     >  4:51522/1‑147 (MQ=255)
                                                 ccTTCTTCGTACTGGCCGGCACTGTGGATAATCTGGTCGATGTGCTTGAGCGCTTCCTTGTTCGGCGCCAGTACCACGCGGAAGCTTTCCAGGTTGGAAACCTGGCGCTTGTTGATCTCGCGGTTGAACAGGGCCAGCTGGTGCTCGAg                                                                                                  >  4:100008/1‑149 (MQ=255)
                                                  cttcttCGTACTGGCCGGCACTGTGGATAATCTGGTCGATGTGCTTGAGCGCTTCCTTGTTCGGCGCCAGTACCACGCGGAAGCTTTCCAGGTTGGAAACCTGGCGCTTGTTGATCTCGCGGTTGAACAGGGCCAGCTGGTGCTCGAgg                                                                                                 <  3:100008/149‑1 (MQ=255)
                                                      ttcGTACTGGCCGGCACTGTGGATAATCTGGTCGATGTGCTTGAGCGCTTCCTTGTTCGGCGCCAGTACCACGCGGAAGCTTTCCAGGTTGGAAACCTGGCGCTTGTTGATCTCGCGGTTGAACAGGGCCAGCTGGTGCTCGAGGCTGt                                                                                             >  6:100324/1‑149 (MQ=255)
                                                       tcGTACTGGCCGGCACTGTGGATAATCTGGTCGATGTGCTTGAGCGCTTCCTTGTTCGGCGCCAGTACCACGCGGAAGCTTTCCAGGTTGGAAACCTGGCGCTTGTTGATCTCGCGGTTGAAc                                                                                                                      <  8:153601/123‑1 (MQ=255)
                                                        cGTACTGGCCGGCACTGTGGATAATCTGGTCGATGTGCTTGAGCGCTTCCTTGTTCGGCGCCAGTACCACGCGGAAGCTTTCCAGGTTGGAAACCTGGCGCTTGTTGATCTCGCGGTTGAACAGGGCCAGCTGGTGCTCGAGGCTGTc                                                                                            <  3:176933/148‑1 (MQ=255)
                                                         gTACTGGCCGGCACTGTGGATAATGTGGTCGATGTGCTTGAGCGCTTCCTTGTTCGGCGCCAGTACCACGCGGAAGCTTTCCAGGTTGGAAACCTGGCGCTTGTTGATCTCGCGGTTGAACAGGGCCAGCTGGTGCTCGAGGCTGTCg                                                                                           <  5:304229/148‑1 (MQ=255)
                                                          tACTGGCCGGCACTGTGGATAATCTGGTCGATGTGCTTGAGCGCTTCCTTGTTCGGCGCCAGTACCACGCGGAAGCTTTCCAGGTTGGAAACCTGGCGCTTGTTGATCTCGCGGTTGAACAGGGCCAGCTGGTGCTCGAGGCTGTCGTa                                                                                         <  5:33725/149‑1 (MQ=255)
                                                               gCCGGCACTGTGGATAATCTGGTCGATGTGCTTGAGCGCTTCCTTGTTCGGCGCCAGTACCACGCGGAAGCTTTCCAGGTTGGAAACCTGGCGCTTGTTGATCTCGCGGTTGAACAGGGCCAGCTGGTGCTCGAGGCTGTCGTAGTCGc                                                                                    >  8:77140/1‑149 (MQ=255)
                                                                      cTGTGGATAATCTGGTCGATGTGCTTGAGCGCTTCCTTGTTCGGCGCCAGTACCACGCGGAAGCTTTCCAGGTTGGAAACCTGGCGCTTGTTGATCTCGCGGTTGAACAGGGCCAGCTGGTGCTCGAGGCTGTCGTAGTCGCTGCGGAt                                                                             >  8:22920/1‑149 (MQ=255)
                                                                             tAATCTGGTCGATGTGCTTGAGCGCTTCCTTGTTCGGCGCCAGTACCACGCGGAAGCTTTCCAGGTTGGAAACCTGGCGCTTGTTGATCTCGCGGTTGAACAGGGCCAGCTGGTGCTCGAGGCTGTCGTAGTCGCTGCGGATGTTgcgc                                                                      <  7:161401/149‑1 (MQ=255)
                                                                                           tgCTTGAGCGCTTCCTTGTTCGGCGCCAGTACCACGCGGAAGCTTTCCAGGTTGGAAACCTGGCGCTTGTTGATCTCGCGGTTGAACAGGGCCAGCTGGTGCTCGAGGCTGTCGTAGTCGCTGCGGATGTTGCGCAGGGTTCGGGCGa                                                         <  5:189736/148‑1 (MQ=255)
                                                                                           tgCTTGAGCGCTTCCTTGTTCGGCGCCAGTACCACGCGGAAGCTTTCCAGGTTGGAAACCTGGCGCTTGTTGATCTCGCGGTTGAACAGGGCCAGCTGGTGCTCGAGGCTGTCGTAGTCGCTGCGGATGTTGCGCAGGGTTCGGGCGAt                                                        <  7:22920/149‑1 (MQ=255)
                                                                                           tgCTTGAGCGCTTCCTTGTTCGGCGCCAGTACCACGCGGAAGCTTTCCAGGTTGGAAACCTGGCGCTTGTTGATCTCGCGGTTGAACAGGGCCAGCTGGTGCTCGAGGCTGTCGTAGTCGCTGCGGATGTTGCGCAGGGTTCGGGCGAt                                                        <  7:227281/149‑1 (MQ=255)
                                                                                           tgCTTGAGCGCTTCCTTGTTCGGCGCCAGTACCACGCGGAAGCTTTCCAGGTTGGAAACCTGGCGCTTGTTGATCTCGCGGTTGAACAGGGCCAGCTGGTGCTCGAGGCTGTCGTAGTCGCTGCGGATGTTGCGCAGGGTTCGGGCGAt                                                        <  1:314174/149‑1 (MQ=255)
                                                                                            gCTTGAGCGCTTCCTTGTTCGGCGCCAGTACCACGCGGAAGCTTTCCAGGTTGGAAACCTGGCGCTTGTTGATCTCGCGGTTGAACAGGGCCAGCTGGTGCt                                                                                                      <  8:101751/102‑1 (MQ=255)
                                                                                            gCTTGAGCGCTTCCTTGTTCGGCGCCAGTACCACGCGGAAGCTTTCCAGGTTGGAAACCTGGCGCTTGTTGATCTCGCGGTTGAACAGGGCCAGCTGGTGCt                                                                                                      >  7:101751/1‑102 (MQ=255)
                                                                                                   cgcTTCCTTGTTCGGCGCCAGTACCACGCGGAAGCTTTCCAGGTTGGAAACCTGGCGCTTGTTGATCTCGCGGTTGAACAGGGCCAGCTGGTGCTCGAGGCTGTCGTAGTCGCTGCGGATGTTGCGCAGGGTTCGGGCGATGTCGGTGa                                                >  8:193769/1‑149 (MQ=255)
                                                                                                      ttCCTTGTTCGGCGCCAGTACCACGCGGAAGCTTTCCAGGTTGGAAACCTGGCGCTTGTTGATCTCGCGGTTGAACAggg                                                                                                                  >  8:189237/1‑80 (MQ=255)
                                                                                                      ttCCTTGTTCGGCGCCAGTACCACGCGGAAGCTTTCCAGGTTGGAAACCTGGCGCTTGTTGATCTCGCGGTTGAACAggg                                                                                                                  <  7:189237/80‑1 (MQ=255)
                                                                                                      ttCCTTGTTCGGCGCCAGTACCACGCGGAAGCTTTCCAGGTTGGAAACCTGGCGCTTGTTGATCTCGCGGTTGAACAGGGCCAGCTGGTGCTCGAGGCTGTCGTAGTCGCTGCGGATGtt                                                                          <  3:135898/120‑1 (MQ=255)
                                                                                                      ttCCTTGTTCGGCGCCAGTACCACGCGGAAGCTTTCCAGGTTGGAAACCTGGCGCTTGTTGATCTCGCGGTTGAACAGGGCCAGCTGGTGCTCGAGGCTGTCGTAGTCGCTGCGGATGtt                                                                          >  4:135898/1‑120 (MQ=255)
                                                                                                      ttCCTTGTTCGGCGCCAGTACCACGCGGAAGCTTTCCAGGTTGGAAACCTGGCGCTTGTTGATCTCGCGGTTGAACAGGGCCAGCTGGTGCTCGAGGCTGTCGTAGTCGCTGCGGATGTTGCGCa                                                                     <  2:125561/125‑1 (MQ=255)
                                                                                                      ttCCTTGTTCGGCGCCAGTACCACGCGGAAGCTTTCCAGGTTGGAAACCTGGCGCTTGTTGATCTCGCGGTTGAACAGGGCCAGCTGGTGCTCGAGGCTGTCGTAGTCGCTGCGGATGTTGCGCa                                                                     >  1:125561/1‑125 (MQ=255)
                                                                                                      ttCCTTGTTCGGCGCCAGTACCACGCGGAAGCTTTCCAGGTTGGAAACCTGGCGCTTGTTGATCTCGCGGTTGAACAGGGCCAGCTGGTGCTCGAGGCTGTCGTAGTCGCTGCGGATGTTGCGCAGGGTTCGGGCGATGTCGGt                                                  <  1:181475/144‑1 (MQ=255)
                                                                                                      ttCCTTGTTCGGCGCCAGTACCACGCGGAAGCTTTCCAGGTTGGAAACCTGGCGCTTGTTGATCTCGCGGTTGAACAGGGCCAGCTGGTGCTCGAGGCTGTCGTAGTCGCTGCGGATGTTGCGCAGGGTTCGGGCGATGTCGGt                                                  >  2:181475/1‑144 (MQ=255)
                                                                                                      ttCCTTGTTCGGCGCCAGTACCACGCGGAAGCTTTCCAGGTTGGAAACCTGGCGCTTGTTGATCTCGCGGTTGAACAGGGCCAGCTGGTGCTCGAGGCTGTCGTAGTCGCTGCGGATGTTGCGCAGGGTTCGGGCGATGTCGGTGACCg                                             >  1:246545/1‑149 (MQ=255)
                                                                                                      ttCCTTGTTCGGCGCCAGTACCACGCGGAAGCTTTCCAGGTTGGAAACCTGGCGCTTGTTGATCTCGCGGTTGAACAGGGCCAGCTGGTGCTCGAGGCTGTCGTAGTCGCTGCGGATGTTGCGCAGGGTTCGGGCGATGTCGGTGACCg                                             >  2:289159/1‑149 (MQ=255)
                                                                                                      ttCCTTGTTCGGCGCCAGTACCACGCGGAAGCTTTCCAGGTTGGAAACCTGGCGCTTGTTGATCTCGCGGTTGAACAGGGCCAGCTGGTGCTCGAGGCTGTCGTAGTCGCTGCGGATGTTGCGCAGGGTTCGGGCGATGTCGGTGACCg                                             >  2:292490/1‑149 (MQ=255)
                                                                                                         cTTGTTCGGCGCCAGTACCACGCGGAAGCTTTCGAGGTTGGAAACCTGGCGCTTGTTGATCTCGCGGTTGAACAGGGCCAGCTGGTGCTCGAGGCTGTCGTAGTCGCTGCGGATGTTGCGCAGGGTTCGGGCGATGTCGGTGACCGCAg                                          <  3:51522/149‑1 (MQ=255)
                                                                                                           tGTTCGGCGCCAGTACCACGCGGAAGCTTTCCAGGTTGGAAACCTGGCGCTTGTTGATCTCGCGGTTGAACAGGGCCAGCTGGTGCTCGAGGCTGTCGTAGTCGCTGCGGATGTTGCGCAGGGTTCGGGCGATGTCGGTGACCGCAgtg                                        >  4:353998/1‑147 (MQ=255)
                                                                                                                ggCGCCAGTACCACGCGGAAGCTTTCCAGGTTGGAAACCTGGCGCTTGTTGATCTCGCGGTTGAACAGGGCCAGCTGGTGCTCGAGGCTGTCGTAGTCGCTGCGGATGTTGCGCAGGGTTCGGGCGATGTCGGTGACCGCAGTgcggc                                    <  3:353998/148‑1 (MQ=255)
                                                                                                                ggCGCCAGTACCACGCGGAAGCTTTCCAGGTTGGAAACCTGGCGCTTGTTGATCTCGCGGTTGAACAGGGCCAGCTGGTGCTCGAGGCTGTCGTAGTCGCTGCGGATGTTGCGCAGGGTTCGGGCGATGTCGGTGACCGCAGCgcggcg                                   >  2:257233/1‑149 (MQ=255)
                                                                                                                       gTACCACGCGGAAGCTTTCCAGGTTGGAAACCTGGCGCTTGTTGATCTCGCGGTTGAACAGGGCCAGCTGGTGCTCGAGGCTGTCGTAGTCGCTGCGGATGTTGCGCAGGGTTCGGGCGATGTCGGTGACCGCAGCGCGGCGTGCCTTg                            <  7:77140/149‑1 (MQ=255)
                                                                                                                          ccACGCGGAAGCTTTCCAGGTTGGAAACCTGGCGCTTGTTGATCTCGCGGTTGAACAGGGCCAGCTGGTGCTCGAGGCTGTCGTAGTCGCTGCGGATGTTGCGCAGGGTTCGGGCGATGTCGGTGACCGCAGCGCGGCGTGCCTTGGc                          >  2:259990/1‑148 (MQ=255)
                                                                                                                           cACGCGGAAGCTTTCCAGGTTGGAAACCTGGCGCTTGTTGATCTCGCGGTTGAACAGGGCCAGCTGGTGCTCGAGGCTGTCGTAGTCGCTGCGGATGTTGCGCAGGGTTCGGGCGATGTCGGTGACCGCAGCGCGGCGTGCCTTGGcc                         >  2:2263/1‑148 (MQ=255)
                                                                                                                           cACGCGGAAGCTTTCCAGGTTGGAAACCTGGCGCTTGTTGATCTCGCGGTTGAACAGGGCCAGCTGGTGCTCGAGGCTGTCGTAGTCGCTGCGGATGTTGCGCAGGGTTCGGGCGATGTCGATGACCGCAGCGCGGCGTGCCTTGGcc                         >  2:18348/1‑148 (MQ=255)
                                                                                                                             cgcgGAAGCTTTCCAGGTTGGAAACCTGGCGCTTGTTGATCTCGCGGTTGAACAGGGCCAGCTGGTGCTCGAGGCTGTCGTAGTCGCTGCGGATGTTGCGCAGGGTTCGGGCGATGTCGGTGACCGCAGCGCGGCGTGCCTTGGCCAAg                      <  8:252894/149‑1 (MQ=255)
                                                                                                                                 gAAGCTTTCCAGGTTGGAAACCTGGCGCTTGTTGATCTCGCGGTTGAACAGGGCCAGCTGGTGCTCGAGGCTGTCGTAGTCGCTGCTGATGTTGCGCAGGGTTCGGGCGATGTCGGTGACCGCAGCGCGGCGTGCCTTGGCCAAGGTCa                  <  1:292490/149‑1 (MQ=255)
                                                                                                                                     agcgccAGGTTTCCAACCTGGCGCTTGTTGATCTCGCGGTTGAACAGGGCCAGCTGGTGCTCGAGGCTGTCGTAGTCGCTGCGGATGTTGCGCAGGGTTCGGGCGATGTCGGTGACCGCAGCGCGGCGTGCCTTGGCCAAGGTCAGGGc              >  2:79177/5‑149 (MQ=255)
                                                                                                                                       ttCCAGGTTGGAAACCTGGCGCTTGTTGATCTCGCGGTTGAACAGGGCCAGCTGGTGCTCGAGGCTGTCGTAGTCGCTGCGGATGTTGCGCAGGGTTCGGGCGATGTCGGTGACCGCAGCGCGGCGTGCCTTGGCCAAGGTCAGGGCtt            >  1:324352/1‑149 (MQ=255)
                                                                                                                                         ccAGGTTGGAAACCTGGCGCTTGTTGATCTCGCGGTTGAACAGGGCCAGCTGGTGCTCGAGGCTGTCGTAGTCGCTGCGGATGTTGCGCAGGGTTCGGGCGATGTCGGTGACCGCAGCGCGGCGTGCCTTGGCCAAGGTCAGGGCTtcg          >  3:62013/1‑149 (MQ=255)
                                                                                                                                          cAGGTTGGAAACCTGGCGCTTGTTGATCTCGCGGTTGAACAGGGCCAGCTGGTGCTCGAGGCTGTCGTAGTCGCTGCGGATGTTGCGCAGGGTTCGGGCGATGTCGGTGACCGCAGCGCGGCGTGCCTTGGCCAAGGTCAGGGCTtcgt         >  5:269857/1‑149 (MQ=255)
                                                                                                                                          cAGGATGGAAACCTGGCGCTTGTTGATATCGCGGTTGAACAGGGCCAGCTGGTGCTCGAGGCTGTCGAAGTCGCTGCGGATGAAGCGCAGGGTTCGCGCGATGTCGGTGACCGCAGCGCGGCGTGCCTTGGCCAAGGTCAGGCCTtcgt         >  8:253799/1‑149 (MQ=255)
                                                                                                                                           aGGTTGGAAACCTGGCGCTTGTTGATCTCGCGGTTGAACAGGGCCAGCTGGTGCTCGAGGCTGTCGTAGTCGCTGCGGATGTTGCGCAGGGTTCGGGCGATGTCGGTGACCGCAGCGCGGCGTGCCTTGGCCAAGGTCAGGGCTtcgtc        <  4:62013/149‑1 (MQ=255)
                                                                                                                                            ggTTGGAAACCTGGCGCTTGTTGATCTCGCGGTTGAACAGGGCCAGCTGGTGCTCGAGGCTGTCGTAGTCGCTGCGGATGTTGCGCAGGGTTCGGGCGATGTCGGAGACCGCAGCGCGGCGTGCCTTGGCCAAGGTCAGGGCTtcgtcg       >  4:276082/1‑149 (MQ=255)
                                                                                                                                             gTTGGAAACCTGGCGCTTGTTGATCTCGCGGTTGAACAGGGCCAGCTGGTGCTCGAGGCTGTCGTAGTCGCTGCGGATGTTGCGCAGGGTTCGGGCGATGTCGGTGACCGCAGCGCGGCGTGCCTTGGCCAAGGTCAGGGCTTCGTCgg      >  4:155653/1‑149 (MQ=255)
                                                                                                                                              ttGGAAACCTGGCGCTTGTTGATCTCGCGGTTGAACAGGGCCAGCTGGTGCTCGAGGCTGTCGTAGTCGCTGCGGATGTTGCGCAGGGGTCGGGCGATGTCGGTGACCGCAGCGCGGCGTGCCTTGGCCAAGGTCAGGGCTTCGTCGGt     <  2:3568/149‑1 (MQ=255)
                                                                                                                                                 gAAACCTGGCGCTTGTTGATCTCGCGGTTGAACAGGGCCAGCTGGTGCTCGAGGCTGTCGTAGTCGCTGCGGATGTTGCGCAGGGTTCGGGCGATGTCGGTGACCGCAGCGCGGCGTGCCTTGGCCAAGGTCAGGGCTTCGTCGGTACg  <  1:289159/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
TTGTCCTGATCGGCACTCTGGGTCAGGTCGAAG‑ACCGACAGGGTTTCGCCTTCTTCGTACTGGCCGGCACTGTGGATAATCTGGTCGATGTGCTTGAGCGCTTCCTTGTTCGGCGCCAGTACCACGCGGAAGCTTTCCAGGTTGGAAACCTGGCGCTTGTTGATCTCGCGGTTGAACAGGGCCAGCTGGTGCTCGAGGCTGTCGTAGTCGCTGCGGATGTTGCGCAGGGTTCGGGCGATGTCGGTGACCGCAGCGCGGCGTGCCTTGGCCAAGGTCAGGGCTTCGTCGGTACG  >  NC_002947/804994‑805286

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