Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 973,304 T→A 5.8% R532R (CGT→CGA secD → protein translocase subunit SecD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947973,3040TA 5.8% 271.5 / 4.7 86R532R (CGT→CGAsecDprotein translocase subunit SecD
Reads supporting (aligned to +/- strand):  ref base T (38/43);  new base A (2/3);  total (40/46)
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.77e-01

TGGCCACCATCGCCCTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGCCTACTCGGCAATCGTCGACGCCAACCTGACCACCCTGTTGGTCGGCGGCATCCTGTTTGCCATGGGTACCGGCCCGGTC  >  NC_002947/973159‑973448
                                                                                                                                                 |                                                                                                                                                
tGGCCACCATCGCCCTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTAtt                                                                                                                                               >  1:261268/1‑149 (MQ=255)
tGGCCACCATCGCCCTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTAtt                                                                                                                                               >  4:82209/1‑149 (MQ=255)
  gCCACCATCGCCCTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTcg                                                                                                                                             >  6:313688/1‑149 (MQ=255)
  gCCACCATCGCCCTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTc                                                                                                                                              <  5:145317/148‑1 (MQ=255)
    caccaTCGCCCTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTcgcg                                                                                                                                           <  1:367333/149‑1 (MQ=255)
    caccaTCGCCCTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTcgc                                                                                                                                            >  1:216719/1‑148 (MQ=255)
    caccaTCGCCCTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCAGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTcgc                                                                                                                                            >  8:39600/1‑148 (MQ=255)
     accaTCGCCCTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGa                                                                                                                                          <  8:345199/149‑1 (MQ=255)
             ccTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTCTTCGCGAATACGCGAGAAGATca                                                                                                                                  <  3:343804/149‑2 (MQ=255)
                 gCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGc                                                                                                                              <  5:313688/149‑1 (MQ=255)
                   gggCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTg                                                                                                                            <  2:359380/149‑1 (MQ=255)
                   gggCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTg                                                                                                                            <  5:219486/149‑1 (MQ=255)
                     gCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGt                                                                                                                          <  1:183514/149‑1 (MQ=255)
                       aaCATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTAt                                                                                                                        <  2:191369/149‑1 (MQ=255)
                           tGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCg                                                                                                                    <  8:172093/149‑1 (MQ=255)
                              tgctgctgcTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTg                                                                                                                 <  6:152012/149‑1 (MQ=255)
                                ctgctgctgGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTg                                                                                                                 >  4:67344/1‑147 (MQ=255)
                                             tGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAgcgc                                                                                                           <  7:175449/140‑1 (MQ=255)
                                             tGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAgcgc                                                                                                           >  8:175449/1‑140 (MQ=255)
                                              gATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCAc                                                                                                  >  4:216186/1‑148 (MQ=255)
                                                tGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGaa                                                                                               <  3:82209/149‑1 (MQ=255)
                                                 gTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTCTTCGCGAATACGCGAGAAGATcagc                                                                                                                                >  4:215316/1‑112 (MQ=38)
                                                 gTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTCTTCGCGAATACGCGAGAAGATcagc                                                                                                                                <  3:215316/115‑4 (MQ=38)
                                                  tCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAgg                                                                                             >  8:40066/1‑149 (MQ=255)
                                                    gctgctGGGGGCGAGCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAaga                                                                                                                                       >  7:328988/1‑105 (MQ=255)
                                                    gctgctGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAaga                                                                                                                                       <  8:328988/105‑1 (MQ=255)
                                                    gctgctGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCt                                                                                           >  3:246497/1‑149 (MQ=255)
                                                    gctgctGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCt                                                                                           <  4:236520/149‑1 (MQ=255)
                                                    gctgctGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCt                                                                                           >  6:344835/1‑149 (MQ=255)
                                                     ctgctgGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCtt                                                                                          <  2:331504/149‑1 (MQ=255)
                                                      tgctgGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCtt                                                                                          >  3:96150/1‑148 (MQ=255)
                                                        ctgGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCaa                                                                                       <  5:344835/149‑1 (MQ=255)
                                                            gggCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATcgc                                                                                   <  3:67344/149‑1 (MQ=255)
                                                             ggCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATcgcg                                                                                  <  2:179256/149‑1 (MQ=255)
                                                              gCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCgcgc                                                                                 <  4:305257/149‑1 (MQ=255)
                                                              gCGACCCTGACCCTGCCGGGTATTGCCAGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCgcgc                                                                                 <  7:39600/149‑1 (MQ=255)
                                                               cGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGcc                                                                                <  2:264424/149‑1 (MQ=255)
                                                               cGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGcc                                                                                <  7:40066/149‑1 (MQ=255)
                                                                         cctgCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGCCTACTCGGCaa                                                                      <  3:234241/149‑1 (MQ=255)
                                                                                ggTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGCCTACTCGGCAATCGTCGa                                                               <  5:486/149‑1 (MQ=255)
                                                                                ggTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGGGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGCCTACTCGGCAATCGTCGa                                                               <  5:221348/149‑1 (MQ=255)
                                                                                 gTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGCCTACTCGGCAATCGTCGAc                                                              <  2:261268/149‑1 (MQ=255)
                                                                                    ttGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGCCTACTCGGCAATCGTCGACGcc                                                           >  4:202918/1‑149 (MQ=255)
                                                                                        cGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGCCTACTCGGCAAtcgtcg                                                                <  1:257190/140‑1 (MQ=255)
                                                                                        cGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGCCTACTCGGCAAtcgtcg                                                                >  2:257190/1‑140 (MQ=255)
                                                                                         ggCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGCCTACTCGGCAATCGTCGACGCCAACCt                                                      <  3:216186/149‑1 (MQ=255)
                                                                                          gCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTACTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGCCTACTCGGCAATCGTCGACGCCAACCt                                                      >  2:255581/1‑148 (MQ=255)
                                                                                            aTCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATcgcg                                                                                  <  6:172385/118‑1 (MQ=255)
                                                                                            aTCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATcgcg                                                                                  >  5:172385/1‑118 (MQ=255)
                                                                                             tCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGCCTACTCGGCAATCGTCGACGCCAACCTGac                                                   >  4:173514/1‑148 (MQ=255)
                                                                                                tgtTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGCCTACTCGGCAATCGTCGACGCCAACCTGaccacc                                               <  4:246497/149‑1 (MQ=255)
                                                                                                tgtTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGCCTACTCGGCAATCGTCGACGCCAACCTGaccacc                                               <  3:332532/149‑1 (MQ=255)
                                                                                                tgtTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGCCTACTCGGCAATCGTCGACGCCAACCTGaccacc                                               >  5:263954/1‑149 (MQ=255)
                                                                                                     aCCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTg                                                                                                                            <  6:186647/67‑1 (MQ=255)
                                                                                                     aCCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTg                                                                                                                            >  5:186647/1‑67 (MQ=255)
                                                                                                     aCCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATcg                                                                                    <  3:293901/107‑1 (MQ=255)
                                                                                                     aCCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATcg                                                                                    >  4:293901/1‑107 (MQ=255)
                                                                                                     aCCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGCCTACTCGGCAATCGTCGACGCCAACCTGACCACCCTGtt                                          <  8:27974/149‑1 (MQ=255)
                                                                                                     aCCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGCCTACTCGGCAATCGTCGACGCCAACCTGACCACCCTGtt                                          <  3:173514/149‑1 (MQ=255)
                                                                                                      ccATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGCCTACTCGGCAATCGTCGACGCCAACCTGACCACCCTGTTg                                         <  3:88572/149‑1 (MQ=255)
                                                                                                       cATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGGGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGCCTACTCGGCAATCGTCGACGCCAACCTGACCACCCTGTTgg                                        >  3:48586/1‑149 (MQ=255)
                                                                                                        aTGGGTATGGCGGTGGACGCCAACGTTCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGCCTACTCGGCAATCGTCGACGCCAACCTGACCACCCTGTTGGt                                       <  1:255581/149‑1 (MQ=255)
                                                                                                          gggTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGCCTACTCGGCAATCGTCGACGCCAACCTGACCACCCTGTTGGTcg                                     >  4:197790/1‑149 (MQ=255)
                                                                                                             tATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGCCTACTCGGCAATCGTCGACGCCAACCTGACCACCCTGTTGGTcggcg                                  >  4:122207/1‑149 (MQ=255)
                                                                                                               tGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGCCTACTCGGCAATCGTCGACGCCAACCTGACCACCCTGTTGGTcggcgg                                 >  7:64555/1‑148 (MQ=255)
                                                                                                                       gACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGCCTACTCGGCAATCGTCGACGCCAACCTGACCACCCTGTTGGTCGGCGGCATCCTGtt                        >  7:85206/1‑149 (MQ=255)
                                                                                                                            cAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGCCTACTCGGCAATCGTCGACGCCAACCTGACCACCCTGTTGGTCGGCGGCATCCTGTTTGCCa                   >  5:83952/1‑149 (MQ=255)
                                                                                                                             aaCGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTGCAATCGCGCCTACTCGGCAATCGTCGACGCCAACCTGACCACCCTGTTGGTCGGCGGCATCCTGTTTGCCAt                  <  2:216719/149‑1 (MQ=255)
                                                                                                                              aCGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGCCTACTCGGCAATCGTCGACGCCAACCTGACCACCCTGTTGGTCGGCGGCATCCTGTTTGCCATg                 <  7:142884/149‑1 (MQ=255)
                                                                                                                               cGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGCCTACTCGGCAATCGTCGACGCCAACCTGACCACCCTGTTGGTCGGCGGCATCCTGTTTGCCATgg                <  3:202918/149‑1 (MQ=255)
                                                                                                                                 tGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGCCTACTCGGCAATCGTCGACGCCAACCTGACCACCCTGTTGGTCGGCGGCATCCTGTTTGCCATggg               >  7:132105/1‑148 (MQ=255)
                                                                                                                                  gCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGCCTACTCGGCAATCGTCGACGCCAACCTGACCACCCTGTTGGTCGGCGGCATCCTGTTTGCCATGGGTa             >  2:353256/1‑149 (MQ=255)
                                                                                                                                  gCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGCCTACTCGGCAATCGTCGACGCCAACCTGACCACCCTGTTGGTCGGCGGCATCCTGTTTGCCATGGGTa             >  5:199987/1‑149 (MQ=255)
                                                                                                                                    tGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGCCTACTCGGCAATCGTCGACGCCAACCTGACCACCCTGTTGGTCGGCGGCATCCTGTTTGCCATGGGTAcc           >  1:165582/1‑149 (MQ=255)
                                                                                                                                    tGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGCCTACTCGGCAATCGTCGACGCCAACCTGACCACCCTGTTGGTCGGCGGCATCCTGTTTGCCATGGGTAcc           >  3:375312/1‑149 (MQ=255)
                                                                                                                                     gATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGCCTACTCGGCAATCGTCGACGCCAAc                                                        >  6:152164/1‑103 (MQ=255)
                                                                                                                                     gATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGCCTACTCGGCAATCGTCGACGCCAAc                                                        <  5:152164/103‑1 (MQ=255)
                                                                                                                                     gATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGCCTACTCGGCAATCGTCGACGCCAACCTGACCACCCTGTTGGTCGGCGGCATCCTGTTTGCCATGGGTACCg          >  5:254418/1‑149 (MQ=255)
                                                                                                                                      aTCTTCTCGCGTATTCGCGAAGAGATCGGGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGCCTACTCGGCAATCGTCGACGCCAACCTGACCACCCTGTTGGTCGGCGGCATCCTGTTTGCCATGGGTAcc           <  4:48586/147‑1 (MQ=255)
                                                                                                                                       tctcttcGCGAATACGCGAGAAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGCCTACTCGGCAATCGTCGACGCCAACCTGACCACCCt                                             <  4:177297/107‑1 (MQ=255)
                                                                                                                                       tctcttcGCGAATACGCGAGAAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGCCTACTCGGCAATCGTCGACGCCAACCTGACCACCCt                                             >  3:177297/6‑112 (MQ=255)
                                                                                                                                         ttctCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGCCTACTCGGCAATCGTCGACGCCAACCTGACCACCCTGTTGGTCGGCGGCATCCTGTTTGCCATGGGTACCGGccc      >  7:673/1‑149 (MQ=255)
                                                                                                                                           ctcGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGCCTACTCGGCAATCGTCGACGCCAACCTGACCACCCTGTTGGTCGGCGGCATCCTGTTTGCCATGGGTACCGGc        <  6:263954/145‑1 (MQ=255)
                                                                                                                                            tcGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGCCTACTCGGCAATCGTCGACGCCAACCTGACCACCCTGTTGGTCGGCGGCATCCTGTTTGCCATGGGTACCGGCCCGGt   >  5:102546/1‑149 (MQ=255)
                                                                                                                                            tcGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGCCTACTCGGCAATCGTCGACGCCAACCTGACCACCCTGTTGGTCGGCGGCATCCTGTTTGCCATGGGTACCGGCCCGGt   >  3:353271/1‑149 (MQ=255)
                                                                                                                                              gcgTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGCCTACTCGGCAATCGTCGACGCCAACCTGACCACCCTGTTGGTCGGCGGCATCCTGTTTGCCATGGGTACCGGCCCGGTc  >  3:166624/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
TGGCCACCATCGCCCTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGCCTACTCGGCAATCGTCGACGCCAACCTGACCACCCTGTTGGTCGGCGGCATCCTGTTTGCCATGGGTACCGGCCCGGTC  >  NC_002947/973159‑973448

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: