Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 941,939 C→T 5.1% G1228D (GGC→GAC)  PP_0806 ← surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947941,9390CT 5.1% 299.0 / 3.3 97G1228D (GGC→GAC) PP_0806surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base C (46/46);  new base T (2/3);  total (48/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.47e-01

CGCCTCGGGCGGTGTGATATCCGGGGCTGTGAACTCCAGCGGTATCGAACTGTTGCCGGCGGCATCCACCAGACGCACGTCGAGTGCTTCACCATTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGGCGGCAGCAAGTCGGCCGAGGTGACCGAGCCCGGTTG  >  NC_002947/941791‑942087
                                                                                                                                                    |                                                                                                                                                    
cGCCTCGGGCGGTGTGATATCCGGGGCTGTGAACTCCAGCGGTATCGAACTGTTGCCGGCGGCATCCACCAGACGCACGTCGAGTGCTTCACCATTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGc                                                                                                                                                      >  2:83060/1‑149 (MQ=38)
   cTCGGGCGGTGTGATATCCGGGGCTGTGAACTCCAGCGGTATCGAACTGTTGCCGGCGGCATCCACCAGACGCACGTCGAGTGCTTCACCATTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCgg                                                                                                                                                   <  5:142671/149‑1 (MQ=39)
      gggCGGTGTGATATCCGGGGCTGTGAACTCCAGCGGTATCGAACTGTTGCCGGCGGCATCCACCAGACGCTCGTCGAGTGCTTCACCATTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGc                                                                                                                                                <  1:102469/149‑1 (MQ=255)
      gggCGGTGTGATATCCGGGGCTGTGAACTCCAGCGGTATCGAACTGTTGCCGGCGGCATCCACCAGACGCACGTCGAGTGCTTCACCATTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGGCGTCGTCGGTGc                                                                                                                                                <  7:285673/149‑1 (MQ=255)
      gggCGGTGTGATATCCGGGGCTGTGAACTCCAGCGGTATCGAACTGTTGCCGGCGGCATCCACCAGACGCACGTCGAGTGCTTCACCATTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGt                                                                                                                                                  <  8:216702/147‑1 (MQ=39)
      gggCGGTGTGATATCCGGGGCTGTGAACTCCAGCGGTATCGAACTGTTGCCGGCGGCATCCACCAGACGCACGTCGAGTGCTTCACCATTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGt                                                                                                                                                  <  4:321164/147‑1 (MQ=39)
      gggCGGTGTGATATCCGGGGCTGTGAACTCCAGCGGTATCGAACTGTTGCCGGCGGCATCCACCAGACGCACGTCGAGTGCTTCACCATTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTg                                                                                                                                                 <  7:35007/148‑1 (MQ=40)
         cGGTGTGATATCCGGGGCTGTGAACTCCAGCGGTATCGAACTGTTGCCGGCGGCATCCACCAGACGCACGTCGAGTGCTTCACCATTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGGCGTCGTCGGTGccga                                                                                                                                             <  7:84286/149‑3 (MQ=255)
          ggTGTGATATCCGGGGCTGTGAACTCCAGCGGTATCGAACTGTTGCCGGCGGCATCCACCAGACGCACGTCGAGTGCTTCACCATTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTTGGCACCGGCGTCGTCGGTGccgac                                                                                                                                            <  7:262416/149‑4 (MQ=11)
                atatCCGGGGCTGTGAACTCCAGCGGTATCGAACTGTTGCCGGCGGCATCCACCAGACGCACGTCGAGTGCTTCACCATTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGtt                                                                                                                                      >  5:344056/1‑149 (MQ=255)
                  atCCGGGGCTGTGAACTCCAGCGGAATCGAACTGTTGCCGGCGGCATCCACCAGACGCACGTCGAGTGCTTCACCATTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGc                                                                                                                                    >  3:60685/1‑149 (MQ=255)
                   tCCGGGGCTGTGAACTCCAGCGGTATCGAACTGTTGCCGGCGGCATCCACCAGACGCACGTCGAGTGCTTCACCATTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGcc                                                                                                                                   >  4:368314/1‑149 (MQ=255)
                      ggggCTGTGAACTCCAGCGGTATCGAACTGTTGCCGGCGGCATCCACCAGACGCACGTCGAGTGCTTCACCATTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCgg                                                                                                                                 >  7:16351/1‑148 (MQ=255)
                         gCTGTGAACTCCAGCGGTATCGAACTGTTGCCGGCGGCATCCACCAGACGCACGTCGAGTGCTTCACCATTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGgcag                                                                                                                              >  4:181242/1‑148 (MQ=255)
                           tgtgAACTCCAGCGGTATCGAACTGTTGCCGGCGGCATCCACCAGACGCACGTCGAGTGCTTCACCATTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCAt                                                                                                                           >  1:125461/1‑149 (MQ=255)
                           agtgAACTCCAGCGGTATCGAACTGTTGCCGGCGGCATCCACCAGACGCACGTCGAGTGCTTCACCAATGGCCTGGGACGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGCCGACGCCCGTGCCAAATAGGATGCCGGCAGCAt                                                                                                                           >  6:256517/2‑149 (MQ=17)
                             tgAACTCCAGCGGTATCGAACTGTTGCCGGCGGCATCCACCAGACGCACGTCGAGTGCTTCACCATTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATcg                                                                                                                         >  8:80177/1‑149 (MQ=255)
                                 cTCCAGCGGTATCGAACTGTTGCCGGCGGCATCCACCAGACGCACGTCGAGTGCTTCACCATTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCa                                                                                                                     <  1:148450/149‑1 (MQ=255)
                                 cTCCAGCGGTATCGAACTGTTGCCGGCGGCATCCACCAGACGCACGTCGAGTGCTTCACCATTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCa                                                                                                                     >  3:74278/1‑149 (MQ=255)
                                  tCCAGCGGTATCGAACTGTTGCCGGCGGCATCCACCAGACGCACGTCGAGTGCTTCACCATTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCAc                                                                                                                    >  3:294728/1‑149 (MQ=255)
                                   ccAGCGGTATCGACCTGTTGCCGGCGGCATCCACCAGACGCACGTCGAGTGCTTCACCATTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACt                                                                                                                   <  4:38891/149‑1 (MQ=255)
                                       cGGTATCGAACTGTTGCCGGCGGCATCCACCAGACGCACGTCGAGTGCTTCACCATTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGTCGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCa                                                                                                               <  1:83060/149‑1 (MQ=255)
                                           aTCGAACTGTTGCCGGCGGCATCCACCAGACGCACGTCGAGTGCTTCACCATTGACCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACgg                                                                                                            >  1:288630/1‑148 (MQ=255)
                                              gAACTGTTGCCGGCGGCATCCACCAGACGCACGTCGAGTGCTTCACCATTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAgc                                                                                                        <  3:368314/149‑1 (MQ=255)
                                               aaCTGTTGCCGGCGGCATCCACCAGACGCACGTCGAGTGCTTCACCATTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAgcg                                                                                                       >  8:135395/1‑149 (MQ=255)
                                                  tGTTGCCGGCGGCATCCACCAGACGCACGTCGAGTGCTTCACCATTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCg                                                                                                    >  7:314587/1‑149 (MQ=255)
                                                    ttGCCGGCGGCATCCACCAGACGCACGTCGAGTGCTTCACCATTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCgg                                                                                                   <  7:80177/148‑1 (MQ=255)
                                                    ttGCCGGCGGCATCCACCAGACGCACGGCGAGTGCTTCACCATTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGGTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCgg                                                                                                   <  2:120997/148‑1 (MQ=255)
                                                       ccGGCGGCATCCAACAGACGCACGTCGAGTGCTTCCCCATTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTc                                                                                               >  6:187142/1‑149 (MQ=255)
                                                          gcggcATCCACCAGACGCACGTCGAGTGCTTCACCATTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTcgcc                                                                                            <  8:16351/149‑1 (MQ=255)
                                                          gcggcATCCACCAGACGCACGTCGAGTGCTTCACCATTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTcgcc                                                                                            >  7:107962/1‑149 (MQ=255)
                                                            ggcATCCACCAGACGCACGTCGAGTGCTTCACCAATGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGACCGTCGGCACCGACGACGCCGGTGCCAATTAGGATGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTcgccgc                                                                                          <  6:192816/149‑1 (MQ=255)
                                                              cATCCACCAGACGCACGTCGAGTGCTTCACCATTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGa                                                                                        >  2:353366/1‑149 (MQ=255)
                                                               aTCCACCAGACGCACGTCGAGTGCTTCACCTTTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATc                                                                                                                          >  1:33583/1‑114 (MQ=255)
                                                               aTCCACCAGACGCACGTCGAGTGCTTCACCATTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATc                                                                                                                          <  2:33583/114‑1 (MQ=255)
                                                                tcaaccaGACGCACGTCGAGTGCTTCACCATTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCg                                                                                                                                                    <  1:337972/84‑1 (MQ=255)
                                                                tcaaccaGACGCACGTCGAGTGCTTCACCATTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCg                                                                                                                                                    >  2:337972/4‑87 (MQ=255)
                                                                tCCACCAGACGCACGTCGAGTGCTTCACCATTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGAcc                                                                                      <  5:132799/149‑1 (MQ=255)
                                                                  caccaGACGCACGTCGAGTGCTTCACCATTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCgg                                                                                    >  5:238109/1‑149 (MQ=255)
                                                                  caccaGACGCACGTCGAGTGCTTCACCATTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCg                                                                                     <  7:135395/148‑1 (MQ=255)
                                                                   accaGACGCACGTCGAGTGCTTCACCATTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGt                                                                                   >  2:378678/1‑149 (MQ=255)
                                                                    ccaGACGCACGTCGAGTGCTTCACCATTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGt                                                                                   >  2:95132/1‑148 (MQ=255)
                                                                              ggcGAGTGCTTCTCCATTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCAGCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGc                                                                        <  5:256517/147‑1 (MQ=255)
                                                                              gTCGAGTGCTTCACCATTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGaa                                                                                >  3:39258/1‑141 (MQ=255)
                                                                              gTCGAGTGCTTCACCATTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGaa                                                                                <  4:39258/141‑1 (MQ=255)
                                                                                   gTGCTTCACCATTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCg                                                                   <  5:4216/149‑1 (MQ=255)
                                                                                   gTGCTTCACCATTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCg                                                                   >  1:128751/1‑149 (MQ=255)
                                                                                    tGCTTCACCATTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCg                                                                   <  6:344056/148‑1 (MQ=255)
                                                                                      cTTCACCATTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGCCGCCGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGc                                                                 >  2:163921/1‑148 (MQ=255)
                                                                                       ttCACCATTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACg                                                                         >  1:144413/1‑139 (MQ=255)
                                                                                       ttCACCATTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACg                                                                         <  2:144413/139‑1 (MQ=255)
                                                                                       ttCACCATTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGcc                                                                >  7:146444/1‑148 (MQ=255)
                                                                                       ttCACCATTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGcc                                                                >  6:148556/1‑148 (MQ=255)
                                                                                       ttCACCATTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGcc                                                                >  3:6900/1‑148 (MQ=255)
                                                                                       ttCACCATTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGcc                                                                >  7:150827/1‑148 (MQ=255)
                                                                                        tCACCATTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGGCGCCGGTGCCAATTTGGTTGGCGGGAGCATAGCGCACTTGCACGGTAGCGCCGGGCTCGCCGGGACCGGGGGAAGTGACGCCCTCGGcc                                                                >  6:237005/1‑147 (MQ=17)
                                                                                        tCACCATTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGACGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGcc                                                                >  1:268103/1‑147 (MQ=255)
                                                                                         cACCATTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAg                                                             >  4:37521/1‑149 (MQ=255)
                                                                                          aCCATTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAg                                                             <  8:314587/148‑1 (MQ=255)
                                                                                            cATTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGAc                                                                          >  3:234685/1‑133 (MQ=255)
                                                                                            cATTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGAc                                                                          <  4:234685/133‑1 (MQ=255)
                                                                                            cATTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGAAAGGTGAAAGTGACGCCATCGGCCAAGGc                                                           >  6:70088/1‑148 (MQ=255)
                                                                                              ttGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAg                                                        <  1:353366/149‑1 (MQ=255)
                                                                                                 gCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGt                                                      <  3:181242/148‑1 (MQ=255)
                                                                                                 gCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGt                                                      <  5:237005/148‑1 (MQ=255)
                                                                                                 gCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGt                                                      >  4:319822/1‑148 (MQ=255)
                                                                                                 gCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGt                                                      <  4:294728/148‑1 (MQ=255)
                                                                                                   cTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAg                                                   >  3:290111/1‑149 (MQ=255)
                                                                                                        ccGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGt                                              >  6:189456/1‑149 (MQ=255)
                                                                                                         cGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCaa                                                              <  5:348709/132‑1 (MQ=255)
                                                                                                         cGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCaa                                                              >  6:348709/1‑132 (MQ=255)
                                                                                                         cGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGtt                                             <  2:288630/149‑1 (MQ=255)
                                                                                                          ggCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTtg                                            <  3:319822/149‑1 (MQ=255)
                                                                                                           gCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTtgt                                           <  8:150827/149‑1 (MQ=255)
                                                                                                             gACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGc                                         <  2:216490/149‑1 (MQ=255)
                                                                                                                 gggTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAggc                                     <  2:253822/149‑1 (MQ=255)
                                                                                                                      aGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGGCAATTATGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGGCGgca                                 <  3:26017/148‑1 (MQ=255)
                                                                                                                      aGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGGCGgca                                 <  5:260928/148‑1 (MQ=255)
                                                                                                                         cTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAg                                                             >  8:316617/1‑117 (MQ=255)
                                                                                                                         cTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAg                                                             <  7:316617/117‑1 (MQ=255)
                                                                                                                         cTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGGCGGCAGCaa                             <  4:74278/149‑1 (MQ=255)
                                                                                                                          tGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGGCGGCAGCAAg                            >  8:334023/1‑149 (MQ=255)
                                                                                                                            aaCAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCTCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGGCGGCAGCAAGt                           <  6:351403/148‑1 (MQ=255)
                                                                                                                            aaCAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGGCGGCATCAAGTc                          >  8:179126/1‑149 (MQ=255)
                                                                                                                             aCAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGGCGGCAGCAAGTCg                         <  5:189456/149‑1 (MQ=255)
                                                                                                                              cAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGGCGGCAGCAAGTCgg                        >  1:218428/1‑149 (MQ=255)
                                                                                                                                cctaattggCACCGGCGTCGTCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGGCGGCAGCAAGTCGGcc                      >  6:14506/8‑149 (MQ=21)
                                                                                                                                 gCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGGCGGCAGCAAGTCGGCCg                     <  4:6900/149‑1 (MQ=255)
                                                                                                                                   cGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGGCGGCAGCAAGTCGGCCGAg                   <  3:37521/149‑1 (MQ=255)
                                                                                                                                        gCACCGGCGTCGTCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGGCGGCAGCAAGTCGGCCGAGGTGAc              >  5:288099/1‑149 (MQ=255)
                                                                                                                                         cACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCg                                                                                     <  8:108919/77‑1 (MQ=255)
                                                                                                                                         cACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCg                                                                                     >  7:108919/1‑77 (MQ=255)
                                                                                                                                           ccGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGGCGGCAGCAAGTCGGCCGAGGTGACCGa           <  2:218428/149‑1 (MQ=255)
                                                                                                                                            cgacgaCGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGGCGGCAGCAAGTCGGCCGAGGTGACCGAg          >  4:171611/1‑149 (MQ=255)
                                                                                                                                             gacgacGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGGCGGCAGCAAGTCGGCCGAGGTGACCGAGc         >  2:95571/1‑149 (MQ=255)
                                                                                                                                             gacgacGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGGCGGCAGCAAGTCGGCCGAGGTGACCGAGc         >  3:21973/1‑149 (MQ=255)
                                                                                                                                                   gCCGGTGCCAATTAGTTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGTCGGCAGCAAGTCGGCCGAGGTGACCGAGCCCGGtt   <  8:162903/149‑1 (MQ=255)
                                                                                                                                                   gCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGGCGGCAGCAAGTCGGCCGAGGTGACCGAGCCCGGtt   <  1:378678/149‑1 (MQ=255)
                                                                                                                                                    ccGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGGCGGCAGCAAGTCGGCCGAGGTGACCGAGCCCGGTTg  <  2:125461/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CGCCTCGGGCGGTGTGATATCCGGGGCTGTGAACTCCAGCGGTATCGAACTGTTGCCGGCGGCATCCACCAGACGCACGTCGAGTGCTTCACCATTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGGCGGCAGCAAGTCGGCCGAGGTGACCGAGCCCGGTTG  >  NC_002947/941791‑942087

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