Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,517,232 A→C 7.3% intergenic (+69/+17) infA → / ← clpA translation initiation factor IF‑1/ATP‑dependent serine protease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,517,2320AC 7.3% 208.2 / 7.8 69intergenic (+69/+17)infA/clpAtranslation initiation factor IF‑1/ATP‑dependent serine protease
Reads supporting (aligned to +/- strand):  ref base A (36/28);  new base C (2/3);  total (38/31)
Fisher's exact test for biased strand distribution p-value = 6.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.04e-01

CTCACTGGCGACAAGGTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCGGCCAGCTCGCCAAACAGGATCTCCTCGGCCAGCGGCCGCTTGATCT  >  NC_002947/4517086‑4517376
                                                                                                                                                   |                                                                                                                                                
ctcACTGGCGACAAGGTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCaa                                                                                                                                                 <  8:61090/148‑1 (MQ=255)
 tcACTGGCGACAAGGTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAgg                                                                                                                                               >  1:144934/1‑149 (MQ=255)
   aCTGGCGACAACGGCCGCGTCGACCTGACGCCGGACGACCTGAGCAAGGGCCGCCGCACCTTCCGTGCGCGCTAAGCGCAAGCCATTAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCa                                                                                                                                             <  6:258596/149‑1 (MQ=255)
     tGGCGACAAGGTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGc                                                                                                                                           >  7:330276/1‑149 (MQ=255)
         gACAAGGTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCt                                                                                                                                       <  3:224307/149‑1 (MQ=255)
          aCAAGGTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTg                                                                                                                                      <  5:234850/149‑1 (MQ=255)
            aaGGTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCCTTGCAGCCCCTGca                                                                                                                                    <  2:49738/149‑1 (MQ=255)
             aGGTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGca                                                                                                                                    >  2:269911/1‑148 (MQ=255)
              ggTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGcac                                                                                                                                   >  2:169337/1‑148 (MQ=255)
                   gcgTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCCTTGCAGCCCCTGcaatttctt                                                                                                                             <  6:345481/149‑8 (MQ=255)
                    cgTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCAcgc                                                                                                                            <  4:8998/149‑1 (MQ=255)
                        gAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGAcc                                                                                                                        >  5:81890/1‑149 (MQ=255)
                        gAACTGACGCCGTACGACCTGAGCAAGGGCCGCAACCCGAACCTTGCCCGAAAACAACAACCCATAAAAAACCCCGCCCAGTTGCCGGCCTTTTTTGTGCGTAAAGAAATTGAAGGGCATCCAAGGCAGCCCATGAACAACACGCGAc                                                                                                                         >  8:16900/1‑148 (MQ=255)
                            tGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAg                                                                                                                    >  6:133102/1‑149 (MQ=255)
                              aCGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGcc                                                                                                                  >  7:125511/1‑149 (MQ=255)
                                gCCGTACGACCTGAGCAAGGGCCGCATCACCTACCATGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGt                                                                                                                >  7:230655/1‑149 (MQ=255)
                                       gACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGTAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCg                                                                                                         <  5:405564/149‑1 (MQ=255)
                                        aCCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTg                                                                                                               >  1:68077/1‑142 (MQ=255)
                                        aCCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTg                                                                                                               <  2:68077/142‑1 (MQ=255)
                                                aaGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAac                                                                                                <  3:17400/149‑1 (MQ=255)
                                                 aGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAaca                                                                                               >  2:41001/1‑149 (MQ=255)
                                                   ggCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGTGCTGCAAg                                                                                                                                                <  3:364800/98‑1 (MQ=255)
                                                   ggCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAg                                                                                                                                                >  4:364800/1‑98 (MQ=255)
                                                    gCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCa                                                                                            <  8:125511/149‑1 (MQ=255)
                                                     ccGCATCACCTACCGTGCGCGCTAAGCTCAAGCGCATGAAAAAGCCCCGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCg                                                                                                         >  1:260254/1‑136 (MQ=255)
                                                     ccGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAATAGCCCTGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCg                                                                                                         <  2:260254/135‑1 (MQ=255)
                                                     ccGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAg                                                                                           >  4:115149/1‑149 (MQ=255)
                                                         aTCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCg                                                                                       <  3:368443/149‑1 (MQ=255)
                                                           cACCTACCGTGCGCGCTAAGCTCAAGGCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGc                                                                                      >  1:10964/1‑148 (MQ=255)
                                                            aCCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCa                                                                                    <  7:16900/149‑1 (MQ=255)
                                                            aCCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCa                                                                                    <  3:115149/149‑1 (MQ=255)
                                                                 ccGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATcgcg                                                                               >  8:307031/1‑149 (MQ=255)
                                                                  catGCTCGATAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGGCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCgcgc                                                                              <  8:230655/147‑1 (MQ=255)
                                                                   gTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTTTGCCTCAAGAAATTGCAGGGGCTAGAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCa                                                                             <  1:330769/149‑1 (MQ=255)
                                                                    tGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGaa                                                                                                       <  6:357008/122‑1 (MQ=255)
                                                                    tGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGaa                                                                                                       >  5:357008/1‑122 (MQ=255)
                                                                             aGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACAt                                                                    <  7:307031/148‑1 (MQ=255)
                                                                                    gCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCGCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACg                                                            >  3:215567/1‑149 (MQ=255)
                                                                                    agacacaaaaaaGCCCGGCACCTGGCCGGGCTTTTTCATGGCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGcaccac                                                             >  1:385150/8‑149 (MQ=37)
                                                                                        aTGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCCTCGCGCAGGTCGACATGCACCACGcca                                                         >  7:198032/1‑149 (MQ=255)
                                                                                        aTGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGcca                                                         >  8:212537/1‑149 (MQ=255)
                                                                                          aaaaaaGCCCGGCACCTGGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGccac                                                        >  7:207020/2‑148 (MQ=255)
                                                                                                gCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTc                                                 >  1:249454/1‑149 (MQ=255)
                                                                                                       cAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCGGCCaac                                          >  5:19539/1‑147 (MQ=255)
                                                                                                           tGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCGGCCAACTCGc                                      <  6:19539/149‑1 (MQ=255)
                                                                                                              cGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCGGCCAGCTCGCCa                                    >  4:356971/1‑148 (MQ=255)
                                                                                                                ggCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGTCACCATGCTCGGCCAGCTCGCCAAAc                                 <  8:250480/149‑1 (MQ=255)
                                                                                                                ggCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCGGCCAGCTCGCCAAAc                                 <  7:212537/149‑1 (MQ=255)
                                                                                                                   ttttttGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCGGCCAGCTCGCCAAACAg                               >  1:179287/1‑148 (MQ=255)
                                                                                                                       ttGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCGGCCAGCTCGCCAAACAGGAtc                           <  1:41001/148‑1 (MQ=255)
                                                                                                                            ccTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCGGCCAGCTCGCCAAACAGGATctcctc                      <  8:330276/148‑1 (MQ=255)
                                                                                                                            ccTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCGGCCAGCTCGCCAAACAGGATctcctc                      <  1:171392/148‑1 (MQ=255)
                                                                                                                            ccTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCGGCCAGCTCGCCAAACAGGATCTCCTCg                     <  1:169337/149‑1 (MQ=255)
                                                                                                                              tCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCGGCTAGCTCGCCAAACAGGATCTCCTCGGc                   >  1:362619/1‑149 (MQ=255)
                                                                                                                              tCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCGGCCAGCTCGCCAAACAGGATCTCCTCGGc                   >  1:154330/1‑149 (MQ=255)
                                                                                                                              tCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCAAGCCACCATGCTCGGCCAGCTCGCCAAACAGCATCTCCTCgcc                   >  8:185299/1‑147 (MQ=255)
                                                                                                                                aaGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCGGCCAGCTCGCCAAACAGGATCTCCTCGGc                   >  1:43776/1‑147 (MQ=255)
                                                                                                                                 aGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCCCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCg                                                <  4:87703/117‑1 (MQ=255)
                                                                                                                                 aGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCCCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCg                                                >  3:87703/1‑117 (MQ=255)
                                                                                                                                    aaTTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCGGCCAGCTCGCCAAACAGGATCTCCTCGGCCAGCg              <  2:179287/148‑1 (MQ=255)
                                                                                                                                         cAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCGGCCAGCTCGCCAAACAGGATCTCCTCGGCCAGCGGCCGc         <  2:249454/148‑1 (MQ=255)
                                                                                                                                          aGGGGCTGCCTTGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCGGCCAGCTCGCCAAACAGGATCTCCTCGGCCAGCGGCCGCtt       >  3:299109/1‑149 (MQ=255)
                                                                                                                                           ggggCTGCCTTGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCGGCCAGCTCGCCAAAc                                 <  6:129082/122‑1 (MQ=255)
                                                                                                                                           ggggCTGCCTTGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCGGCCAGCTCGCCAAAc                                 >  5:129082/1‑122 (MQ=255)
                                                                                                                                            gggCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCGGCCAGCTCGc                                      >  1:186616/1‑116 (MQ=255)
                                                                                                                                            gggCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCGGCCAGCTCGc                                      <  2:186616/116‑1 (MQ=255)
                                                                                                                                            gggCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCGGCCAGCTCGCCAAACAGGATCTCCTCGGCCAGCGGCCGCTTGa     >  1:230952/1‑149 (MQ=255)
                                                                                                                                            gggCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCGGCCAGCTCGCCAAACAGGATCTCCTCGGCCAGCGGCCGCTTGa     >  6:207676/1‑149 (MQ=255)
                                                                                                                                              gCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCGGCCAGCTCGCCAAACAGGATCTCCTCGGCCAGCGGCCGCTTGATc   >  5:32757/1‑149 (MQ=255)
                                                                                                                                               cTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTAGCCATAGCGCAGGTCGACATGCACCACGCCACCATGCTCGGCCAGCTCGCCAAACAGGATCTCCTCGGCCAGCGGCCGCTTGATCt  >  3:37935/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
CTCACTGGCGACAAGGTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCGGCCAGCTCGCCAAACAGGATCTCCTCGGCCAGCGGCCGCTTGATCT  >  NC_002947/4517086‑4517376

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: