Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 225,366 C→T 7.5% A99T (GCC→ACC)  PP_0173 ← transcriptional factor‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947225,3660CT 7.5% 191.4 / 3.9 66A99T (GCC→ACC) PP_0173transcriptional factor‑like protein
Reads supporting (aligned to +/- strand):  ref base C (27/34);  new base T (2/3);  total (29/37)
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.08e-01

CAATGGCGCATCCCGCAACTGCCGCGCCAGTGCATCGCTGATTCGCGCCACCAGCTCCACGCCCTCACCCTGCGATGGCTGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAGGTCACGCTCGACAGTGCGTTTGCTGACCGTGAAGCCAACGTCACGCAGGCGCCAGACCAGCTCGG  >  NC_002947/225220‑225503
                                                                                                                                                  |                                                                                                                                           
cAATGGCGCATCCCGCAACTGCCGCGCCAGTGCATCGCTGATTCGCGCCACCAGCTCCACGCCCTCACCCTGCGATGGCTGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGTCg                                                                                                                                           >  2:149935/1‑149 (MQ=255)
 aaTGGCGCATCCCGCAACTGCCGCGCCAGTGCATCGCTGATTCGCGCCACCAGCTCCACGCCCTCACCCTGCGATGGCTGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCg                                                                                                                                           <  6:262367/148‑1 (MQ=255)
 aaTGGCGCATCCCGCAACTGCCGCGCCAGTGCATCGCTGATTCGCGCCACCAGCTCCACGCCCTCACCCTGCGATGGCTGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGa                                                                                                                                          >  4:160537/1‑149 (MQ=255)
        cATCCCGCAACTGCCGCGCCAGTGCATCGCTGATTCGCGCCACCAGCTCCACGCCCTCACCCTGCGATGGCTGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGc                                                                                                                                   >  7:15458/1‑149 (MQ=255)
            ccGCAACTGCCGCGCCAGTGCATCGCTGATTCGCGCCACCAGCTCCACGCCCTCACCCTGCGATGGCTGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGTCGGCCAGTGccgt                                                                                                                                <  1:149935/148‑3 (MQ=255)
            ccGCAACTGCCGCGCCAGTGCATCGCTGATTCGCGCCACCAGCTCCACGCCCTCACCCTGCGATGGCTGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGTCGGCCAGTGccgt                                                                                                                                <  3:149458/148‑3 (MQ=255)
            ccGCAACTGCCGCGCCAGTGCATCGCTGATTCGCGCCACCAGCTCCACGCCCTCACCCTGCGATGGCTGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAAc                                                                                                                              <  4:167388/149‑1 (MQ=255)
             cGCAACTGCCGCGCCAGTGCATCGCTGATTCGCGCCACCAGCTCCACGCCCTCACCCTGCGATGGCTGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAAcc                                                                                                                             >  4:208619/1‑149 (MQ=255)
                aaCTGCCGCGCCAGTGCATCGCTGATTCGCGCCACCAGCTCCACGCCCTCACCCTGCGATGGCTGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCgaa                                                                                                                          <  1:280452/149‑1 (MQ=255)
                    gCCGCGCCAGTGCATCGCTGATTCGCGCCACCAGCTCCACGCCCTCACCCTGCGATGGCTGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCgaatgg                                                                                                                       >  5:147832/1‑148 (MQ=255)
                       gcgcCAGTGCATCGCTGATTCGCGCCACCAGCTCCACGCCCTCACCCTGCGATGGCTGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCgaatggaatg                                                                                                                   >  2:51788/1‑149 (MQ=255)
                         gcCAGTGCATCGCTGATTCGCGCCACCAGCTCCACGCCCTCACCCTGCGATGGCTGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGTCACTGGCCGACCAGTGCCAACCGAATGGAATGct                                                                                                                 <  6:43428/149‑1 (MQ=255)
                                 aTCGCTGATTCGCGCCACCAGCTCCACGCCCTCACCCTGCGATGGCTGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGAc                                                                                                                                         <  7:218164/118‑1 (MQ=255)
                                 aTCGCTGATTCGCGCCACCAGCTCCACGCCCTCACCCTGCGATGGCTGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGAc                                                                                                                                         >  8:218164/1‑118 (MQ=255)
                                       gATTCGCGCCACCAGCTCCACGCCCTCACCCTGCGATGGCTGCAACGGATCCCCTCGCAGATGTCTCAGCAAGTCGAAATTACCGCGCAATTAGCCCAAGCCACTGGCCGACCAGTTCCAACCGAAAGGAATGCTCTTGTAATTTCGtt                                                                                                   >  3:57956/1‑149 (MQ=255)
                                           cgcgCCACCAGCTCCACGCCCTCACCCTGCGATGGCTGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCa                                                                                                <  3:160537/148‑1 (MQ=255)
                                           cgcgCCACCAGCTCCACGCCCTCACCCTGCGATGGCTGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCa                                                                                                <  3:208619/148‑1 (MQ=255)
                                             cgcCACCAGCTCCACGCCCTCACCCTGCGATGGCTGCGACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGc                                                                                              <  1:200926/149‑1 (MQ=255)
                                             cgcCACCAGCTCCACGCCCTCACCCTGCGATGGCTGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCg                                                                                             <  5:132188/149‑1 (MQ=255)
                                               ccaccaGCTCCACGCCGTGACCCTGCGATGGCTGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCggg                                                                                           <  4:57956/149‑1 (MQ=255)
                                               ccaccaGCTCCACGCCCTCACCCTGCGATGGCTGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCgg                                                                                            >  1:271997/1‑148 (MQ=255)
                                                   caGCTCCACGCCCTCACCCTGCGATGGCTGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAata                                                                                       >  6:239213/1‑149 (MQ=255)
                                                         cACGCCCTCACCCTGCGATGGCTGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATc                                                                                     <  2:78303/145‑1 (MQ=255)
                                                         cACGCCCTCACCCTGCGATGGCTGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATc                                                                                     >  1:78303/1‑145 (MQ=255)
                                                           cGCCCTCACCCTGCGATGGCTGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGa                                                                                >  8:27211/1‑148 (MQ=255)
                                                            gCCCTCACCCTGCGATGGCTGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGaaa                                                                              <  2:9069/149‑1 (MQ=255)
                                                                   cccTGCGATGGCTGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCggg                                                                                           >  5:57310/1‑129 (MQ=255)
                                                                   cccTGCGATGGCTGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCggg                                                                                           <  6:57310/129‑1 (MQ=255)
                                                                    ccTGCGATGGCTGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGtt                                                                       >  4:112313/1‑148 (MQ=255)
                                                                     cTGCGATGGCTGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGa                                                                     <  7:118278/149‑1 (MQ=255)
                                                                        cGATGGCTGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAg                                                                    <  5:239213/147‑1 (MQ=255)
                                                                              cTGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAGGTCACGCt                                                            >  6:221234/1‑149 (MQ=255)
                                                                               tGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAGGTCACGCTc                                                           >  1:154983/1‑149 (MQ=255)
                                                                                  aaCGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAGGTCACGCTCGa                                                        >  3:244925/1‑148 (MQ=255)
                                                                                  aaCGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAGGTCACGCTCGa                                                        >  6:107034/1‑148 (MQ=255)
                                                                                  aaCGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAGGTCACGCTCGAc                                                       <  4:244925/149‑1 (MQ=255)
                                                                                     gCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAAAGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAGGTCACGCTCGACAGt                                                    >  7:121653/1‑149 (MQ=255)
                                                                                      cATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAGGTCACGCTCGACAGTg                                                   <  8:25922/149‑1 (MQ=255)
                                                                                           ccTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAGGTCACGCTCGACAGTGCGttt                                              <  5:221234/149‑1 (MQ=255)
                                                                                               gCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAGGTCACGCTCGACAGTGCGTTTGc                                            <  5:107034/147‑1 (MQ=255)
                                                                                               gCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAGGTCACGCTCGACAGTGCGTTTGCt                                           >  7:285606/1‑148 (MQ=255)
                                                                                                cAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAGGTCACGCTCGACAGTGCGTTTGCTGa                                         <  3:167852/149‑1 (MQ=255)
                                                                                                 aGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAGGTCACGCTCGACAGTGCGTTTGCTGAc                                        <  1:51788/149‑1 (MQ=255)
                                                                                                      tgtCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAGGTCACGCTCGACAGTGCGTTTGCTGACCGTg                                    <  3:112313/148‑1 (MQ=255)
                                                                                                      tgtCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAGGTCACGCTCGACAGTGCGTTTGCTGACCGTg                                    <  2:154983/148‑1 (MQ=255)
                                                                                                       gtCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAGGTCACGCTCGACAGTGCGTTTg                                             <  7:112043/138‑1 (MQ=255)
                                                                                                       gtCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAGGTCACGCTCGACAGTGCGTTTg                                             >  8:112043/1‑138 (MQ=255)
                                                                                                        tCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAGGTCACGCTCGACAGTGCGTTTGCTGACCGTGaa                                  <  7:27211/148‑1 (MQ=255)
                                                                                                        tCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAGGTCACGCTCGACAGTGCGTTTGCTGACCGTGAAg                                 >  6:38019/1‑149 (MQ=255)
                                                                                                        tCAGCAACTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACACAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAGGTCACGCTCTACAGTGCGTTTGCTGACCGTGaa                                  <  2:282686/148‑1 (MQ=255)
                                                                                                           gCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAGGTCACGCTCGACAGTGCGTTTGCTGACCGTGAAGcc                               >  3:103727/1‑148 (MQ=255)
                                                                                                                  gAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAGGTCACGCTCGACAGTGCGTTTGCTGACCGTGAAGCCAACGTCAc                       >  6:47747/1‑149 (MQ=255)
                                                                                                                    aaTTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAGGTCACGCTCGACAGTGCGTTTGCTGACCGTGAAGCCAACGTCACGc                     <  4:226800/149‑1 (MQ=255)
                                                                                                                            cgcAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAGGTCACGCTCGACAGTGCGTTTGCTGACCGTGAAGCCAACGTCACGCAGgcgc               >  6:104054/1‑147 (MQ=255)
                                                                                                                            cgcAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAGGTCACGCTCGACAGTGCGTTTGCTGACCGTGAAGCCAACGTCACGCAGgcgc               >  2:238260/1‑147 (MQ=255)
                                                                                                                             gcAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAgg                                                                   >  2:268652/1‑95 (MQ=255)
                                                                                                                             gcAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAgg                                                                   >  2:272932/1‑95 (MQ=255)
                                                                                                                             gcAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAgg                                                                   <  1:272932/95‑1 (MQ=255)
                                                                                                                             gcAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAgg                                                                   <  1:268652/95‑1 (MQ=255)
                                                                                                                               aaTTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAGGTCACGCTCGACAGTGCGTTTGCTGACCGTGAAGCCAACGTCACGCAGGCGCCAGAc          <  5:104054/149‑1 (MQ=255)
                                                                                                                                   cGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAGGTCACGCTCGACAGTGCGTTTGCTGACCGTGAAGCCAACGTCACGCAGGCGCCAGACCAg       <  4:103727/148‑1 (MQ=255)
                                                                                                                                       cACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAGGTCACGCTCGACAGTGCGTTTGCTGACCGTGAAGCCAACGTCACGCAGGCGCCAGACCAGCTCgg  <  2:271997/149‑1 (MQ=255)
                                                                                                                                       cACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAGGTCACGCTCGACAGTGCGTTTGCTGACCGTGAAGCCAACGTCACGCAGGCGCCAGACCAGCTCgg  <  8:15458/149‑1 (MQ=255)
                                                                                                                                        aCGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAGGTCACGCTCTTACAGTGCGTTTGCTGACCGTGAAGCCAACGTCACGCAGGCGCCAGACCAGCTCgg  <  1:238260/149‑1 (MQ=255)
                                                                                                                                              cTGGTCGGCCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCgggg                                                                                          >  1:79697/1‑54 (MQ=255)
                                                                                                                                              cTGGTCGGCCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCggg                                                                                           <  2:79697/54‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                           
CAATGGCGCATCCCGCAACTGCCGCGCCAGTGCATCGCTGATTCGCGCCACCAGCTCCACGCCCTCACCCTGCGATGGCTGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAGGTCACGCTCGACAGTGCGTTTGCTGACCGTGAAGCCAACGTCACGCAGGCGCCAGACCAGCTCGG  >  NC_002947/225220‑225503

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: