Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 973,313 A→G 6.2% E535E (GAA→GAG secD → protein translocase subunit SecD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947973,3130AG 6.2% 244.4 / 4.7 81E535E (GAA→GAGsecDprotein translocase subunit SecD
Reads supporting (aligned to +/- strand):  ref base A (35/41);  new base G (2/3);  total (37/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.56e-01

CCTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGCCTACTCGGCAATCGTCGACGCCAACCTGACCACCCTGTTGGTCGGCGGCATCCTGTTTGCCATGGGTACCGGCCCGGTCAAAGGCTTTGC  >  NC_002947/973172‑973459
                                                                                                                                             |                                                                                                                                                  
ccTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTCTTCGCGAATACGCGAGAAGATca                                                                                                                                             <  3:343804/149‑2 (MQ=255)
    gCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGc                                                                                                                                         <  5:313688/149‑1 (MQ=255)
      gggCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTg                                                                                                                                       <  5:219486/149‑1 (MQ=255)
      gggCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTg                                                                                                                                       <  2:359380/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                                                                                                                                           <  3:215316/115‑4 (MQ=38)
                                    gTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTCTTCGCGAATACGCGAGAAGATcagc                                                                                                                                           >  4:215316/1‑112 (MQ=38)
                                     tCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAgg                                                                                                        >  8:40066/1‑149 (MQ=255)
                                       gctgctGGGGGCGAGCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAaga                                                                                                                                                  >  7:328988/1‑105 (MQ=255)
                                       gctgctGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAaga                                                                                                                                                  <  8:328988/105‑1 (MQ=255)
                                       gctgctGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCt                                                                                                      <  4:236520/149‑1 (MQ=255)
                                       gctgctGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCt                                                                                                      >  3:246497/1‑149 (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                                                                                           <  7:40066/149‑1 (MQ=255)
                                                  cGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGcc                                                                                           <  2:264424/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                                                                                             >  5:172385/1‑118 (MQ=255)
                                                                               aTCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATcgcg                                                                                             <  6:172385/118‑1 (MQ=255)
                                                                                tCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGCCTACTCGGCAATCGTCGACGCCAACCTGac                                                              >  4:173514/1‑148 (MQ=255)
                                                                                   tgtTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGCCTACTCGGCAATCGTCGACGCCAACCTGaccacc                                                          >  5:263954/1‑149 (MQ=255)
                                                                                   tgtTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGCCTACTCGGCAATCGTCGACGCCAACCTGaccacc                                                          <  3:332532/149‑1 (MQ=255)
                                                                                   tgtTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGCCTACTCGGCAATCGTCGACGCCAACCTGaccacc                                                          <  4:246497/149‑1 (MQ=255)
                                                                                        aCCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTg                                                                                                                                       >  5:186647/1‑67 (MQ=255)
                                                                                        aCCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTg                                                                                                                                       <  6:186647/67‑1 (MQ=255)
                                                                                        aCCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATcg                                                                                               >  4:293901/1‑107 (MQ=255)
                                                                                        aCCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATcg                                                                                               <  3:293901/107‑1 (MQ=255)
                                                                                        aCCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGCCTACTCGGCAATCGTCGACGCCAACCTGACCACCCTGtt                                                     <  3:173514/149‑1 (MQ=255)
                                                                                        aCCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGCCTACTCGGCAATCGTCGACGCCAACCTGACCACCCTGtt                                                     <  8:27974/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                                                        >  3:177297/6‑112 (MQ=255)
                                                                                                                          tctcttcGCGAATACGCGAGAAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGCCTACTCGGCAATCGTCGACGCCAACCTGACCACCCt                                                        <  4:177297/107‑1 (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)
                                                                                                                                      ttCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGCCTACTCGGCAATCGTCGACGCCAACCTGACCACCCTGTTGGTCGGCGGCATCCTGTTTGCCATGGGTACCGGCCCGGTCAAAgg        <  8:673/148‑1 (MQ=255)
                                                                                                                                        cgcgAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGCCTACTCGGCAATCGTCGACGCCAACCTGACCACCCTGTTGGTCGGCGGCATCCTGTTTGCCATGGGTACCGGCCCGGTCAAAGGCt      <  8:132105/148‑1 (MQ=255)
                                                                                                                                         gcgAAGAGATCGCGGCTGGGATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCCATCGCGCCTACTCGGCAATCGTCGCCGCCAACCTGACCACCCTGTTGGTCGGCGGCATCCGGTTGGCCATGGGTCCCGGCCCGGTCAAAGGCttt    >  5:281424/1‑149 (MQ=255)
                                                                                                                                          cgAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGCCTACTCGGCAATCGTCGACGCCAACCTGACCACCCTGTTGGTCGGCGGCATCCTGTTTGCCATGGGTACCGGCCCGGTCAAAGGCTTTg   >  3:132751/1‑149 (MQ=255)
                                                                                                                                           gAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGCCTACTCGGCAATCGTCGACGCCAACCTGACCACCCTGTTGGTCGGCGGCATCCTGTTTGCCATGGGTACCGGCCCGGTCAAAGGCTTTGc  >  1:161873/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                                  
CCTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGCCTACTCGGCAATCGTCGACGCCAACCTGACCACCCTGTTGGTCGGCGGCATCCTGTTTGCCATGGGTACCGGCCCGGTCAAAGGCTTTGC  >  NC_002947/973172‑973459

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: