Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,517,233 A→T 7.1% intergenic (+70/+16) 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,2330AT 7.1% 208.3 / 5.3 70intergenic (+70/+16)infA/clpAtranslation initiation factor IF‑1/ATP‑dependent serine protease
Reads supporting (aligned to +/- strand):  ref base A (34/31);  new base T (3/2);  total (37/33)
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.37e-01

TCACTGGCGACAAGGTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCGGCCAGCTCGCCAAACAGGATCTCCTCGGCCAGCGGCCGCTTGATCTTGTC  >  NC_002947/4517087‑4517380
                                                                                                                                                  |                                                                                                                                                   
tcACTGGCGACAAGGTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAgg                                                                                                                                                   >  6:393476/1‑149 (MQ=255)
      gCGACAAGGTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGccc                                                                                                                                             >  6:336771/1‑149 (MQ=255)
      gCGACAAGGTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGc                                                                                                                                               <  5:414954/147‑1 (MQ=255)
        gACAAGGTCCGCGTCGAACTGACGCCGTACGCCCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCt                                                                                                                                           >  5:362797/1‑149 (MQ=255)
        gACAAGGTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGccc                                                                                                                                             <  1:63456/147‑1 (MQ=255)
        gACAAGGTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGccc                                                                                                                                             >  2:63456/1‑147 (MQ=255)
            aGGTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGcac                                                                                                                                       <  5:322643/149‑1 (MQ=255)
             ggTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCCTTGCAGCCCCTGcaat                                                                                                                                      <  6:78291/149‑3 (MQ=255)
             ggTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCCTTGCAGCCCCTGcaat                                                                                                                                      >  5:78291/1‑147 (MQ=255)
                       gAACTGACGCCGTACGACCTGAGCAAGGACCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGAc                                                                                                                             <  6:19164/148‑1 (MQ=255)
                          cTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCa                                                                                                                         <  2:33040/149‑1 (MQ=255)
                          cTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCa                                                                                                                         <  2:251910/149‑1 (MQ=255)
                            gACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTGAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGc                                                                                                                       >  7:150293/1‑149 (MQ=255)
                               gCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGt                                                                                                                    >  5:13009/1‑149 (MQ=255)
                               gCCGAACGACCTGAGCAAGGGCCGCATCACCTACCGAGCGCGCTAAGCTCAAGCAATGAAAAACCCCCGCCAGGTCCCGGGCTATTATGTGACTCAAGAAATACCAAGGGCTGCAAGGCAGCCCCTGCACATCACGCTACATCACCCGt                                                                                                                    >  2:103665/1‑149 (MQ=255)
                                     cGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTc                                                                                                              <  6:13009/149‑1 (MQ=255)
                                      gACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCg                                                                                                             <  6:209866/149‑1 (MQ=255)
                                        ccTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGaa                                                                                                           <  7:55337/149‑1 (MQ=255)
                                         cTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAg                                                                                                          >  1:394432/1‑149 (MQ=255)
                                              cAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGaa                                                                                                     <  5:13310/149‑1 (MQ=255)
                                               aaGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAgg                                                                                                                                                   <  3:8060/102‑1 (MQ=255)
                                               aaGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAgg                                                                                                                                                   >  4:8060/1‑102 (MQ=255)
                                               aaGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAac                                                                                                    >  5:110157/1‑149 (MQ=255)
                                                 gggCCGTATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCCCATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAacac                                                                                                  >  2:327318/1‑149 (MQ=255)
                                                 gggCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAacac                                                                                                  >  3:176263/1‑149 (MQ=255)
                                                 gggCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAacac                                                                                                  >  5:237861/1‑149 (MQ=255)
                                                   gCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAacac                                                                                                  >  3:423675/1‑147 (MQ=255)
                                                   gCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAacac                                                                                                  >  3:416961/1‑147 (MQ=255)
                                                      gCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGtt                                                                                             <  5:393476/149‑1 (MQ=255)
                                                       cATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTc                                                                                            <  4:176263/149‑1 (MQ=255)
                                                         tCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGc                                                                                          <  2:65731/149‑1 (MQ=255)
                                                            ccTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAg                                                                                                          <  2:368349/130‑1 (MQ=255)
                                                               aCCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATcgc                                                                                    >  2:306177/1‑149 (MQ=255)
                                                                ccGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGTGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATcg                                                                                     <  1:134717/147‑1 (MQ=255)
                                                                         cTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGa                                                                           >  1:82975/1‑148 (MQ=255)
                                                                            aGCTCAAGCCATCAAAAAGCCCGGCCATGTGCAGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGAACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCAATCGCGCAGGTCGACAt                                                                        <  6:334075/148‑1 (MQ=255)
                                                                             gCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGc                                                                      <  2:82975/149‑1 (MQ=255)
                                                                              cTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGca                                                                     >  4:384549/1‑149 (MQ=255)
                                                                              cTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGca                                                                     >  5:251256/1‑149 (MQ=255)
                                                                                 aaGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGcacca                                                                  >  1:407423/1‑149 (MQ=255)
                                                                                      aTGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGcca                                                             >  2:48103/1‑149 (MQ=255)
                                                                                      aTGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGcca                                                             <  6:362797/149‑1 (MQ=255)
                                                                                        gAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGccac                                                            >  4:201404/1‑148 (MQ=255)
                                                                                          aaaaGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCAt                                                         >  2:288620/1‑149 (MQ=255)
                                                                                            aaGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGc                                                       >  6:28047/1‑149 (MQ=255)
                                                                                            aaGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGc                                                       >  8:409215/1‑149 (MQ=255)
                                                                                              gCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTc                                                     >  7:13324/1‑149 (MQ=255)
                                                                                                  ggCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCGGc                                                  >  7:142757/1‑148 (MQ=255)
                                                                                                   gCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCGGCCa                                                <  5:398430/149‑1 (MQ=255)
                                                                                                         tGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCGGCCAGCTCGc                                          <  8:13324/149‑1 (MQ=255)
                                                                                                         tGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCGGCCAGCTCGc                                          <  3:384549/149‑1 (MQ=255)
                                                                                                           ccGGGCTTTTTTGTGCCTCAAGAAATTGCATGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCGGCCAGCTCGcc                                         <  5:336771/148‑1 (MQ=255)
                                                                                                           ccGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCACCATAGCGCAGGTCGACATGCACCACGCCACCATGCTCGGCCAGCTCGcc                                         >  5:24265/1‑148 (MQ=255)
                                                                                                              ggCTTTTTTGTGCCTCAAGGAATTGCAGGGGCTGCAAGGCAGGCCCTGCAGATAACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCGGCCAGATCGCCaaa                                      >  7:197489/1‑148 (MQ=255)
                                                                                                                 ttttttGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCGGCCAGCTCGCCAAACAg                                   >  3:72233/1‑148 (MQ=255)
                                                                                                                   ttttGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCGGCCAGCTCGCCAAACAgg                                  <  8:142757/147‑1 (MQ=255)
                                                                                                                    tttGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCACCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCGGCCAGCTCGCCAAACAGGAt                                <  6:24265/148‑1 (MQ=255)
                                                                                                                          ccTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCGGCCAGCTCGCCAAACAGGATCTCCTCg                         >  7:216776/1‑149 (MQ=255)
                                                                                                                              aaGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCGGCCAGCTCGCCAAACAGGATCTCCTCGGc                       >  3:345148/1‑147 (MQ=255)
                                                                                                                              aaGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCGGCCAGCTCGCCAAACAGGATCTCCTCGGc                       <  1:288620/147‑1 (MQ=255)
                                                                                                                              aaGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCGGCCAGCTCGCCAAACAGGATCTCCTCGGc                       >  6:197846/1‑147 (MQ=255)
                                                                                                                                     tGCAGGGGCTGCCTTGTAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCGGCCAGCTCGCCAAACAGGATCTCCTCGGCCAGCGGCCg              >  4:332392/1‑149 (MQ=255)
                                                                                                                                       cAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCGGCCAGCTCGCCAAACAGGATCTCCTCGGCCAGCGGCCGc             <  4:72233/148‑1 (MQ=255)
                                                                                                                                         ggggCTGCCTTGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCCAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCGGCCAGCTCGCCAAACAGGATCTCCTCGGCCAGCGGCCGCt            <  3:161782/147‑1 (MQ=255)
                                                                                                                                         ggggCTGCCTTGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCCAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCGGCCAGCTCGCCAAACAGGATCTCCTCGGCCAGCGGCCGCt            >  4:161782/1‑147 (MQ=255)
                                                                                                                                         ggggCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCGGCCAGCTCGCCAAACAGGATCTCCTCGGCCAGCGGCCGCTTg          <  8:150293/149‑1 (MQ=255)
                                                                                                                                         ggggCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCGGCCAGCTCGCCAAACAGGATCTCCTCGGCCAGCGGCCGCTTg          <  5:28047/149‑1 (MQ=255)
                                                                                                                                          gggCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCACGGCCAGCTCGCCAAACAGGATCTCCTCGGCCAGCGGCCGCTTGa         <  1:306177/149‑1 (MQ=255)
                                                                                                                                               gCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCGGCCAGCTCGCCAAACAGGATCTCCTCGGCCAGCGGCCGCTTGATCTTg    >  3:189947/1‑149 (MQ=255)
                                                                                                                                                cAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCGGCCAGCTCGCCAAACAGGATCTCCTCGGCCAGCGGCCGCTTGATCtt     <  2:407423/147‑1 (MQ=255)
                                                                                                                                                 aaGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACATCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCGGCCAGCTCGCCAAACAGGATCTCCTCGGCCAGCg                  >  3:58314/1‑133 (MQ=255)
                                                                                                                                                 aaGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACATCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCGGCCAGCTCGCCAAACAGGATCTCCTCGGCCAGCg                  <  4:58314/133‑1 (MQ=255)
                                                                                                                                                 aaGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCGGCCAGCTCGCCAAACAGGATCTCCTCGGCCAGCGGCCGCTTGATCTTGTc  <  8:216776/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
TCACTGGCGACAAGGTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCGGCCAGCTCGCCAAACAGGATCTCCTCGGCCAGCGGCCGCTTGATCTTGTC  >  NC_002947/4517087‑4517380

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: