Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 126,115 A→G 6.7% intergenic (+25/+25) PP_0120 → / ← thrB zinc ABC transporter/homoserine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947126,1150AG 6.7% 213.5 / 4.9 74intergenic (+25/+25)PP_0120/thrBzinc ABC transporter/homoserine kinase
Reads supporting (aligned to +/- strand):  ref base A (35/34);  new base G (2/3);  total (37/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.30e-01

CGCCTTGCCGAGCTGGATGCCTTGGGTGGCACCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCA‑‑GGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGATCATCACATCCATCCCGGCGAACGACTCGGCCGCAATCAGG  >  NC_002947/125977‑126257
                                                                                                                                             |                                                                                                                                              
cGCCTTGCCGAGCTGGATGCCTTGGGTGGCACCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCA‑‑GGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCGTGCCCGCgaa                                                                                                                                      <  5:258550/149‑1 (MQ=255)
 gCCTTGCCGAGCTGGATGCCTTGGGTGGCACCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCA‑‑GGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCGTGCCCGCgaag                                                                                                                                     <  2:285949/149‑1 (MQ=255)
 gCCTTGCCGAGCTGGATGCCTTGGGTGGCACCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCA‑‑GGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCcgcg                                                                                                                                        >  6:127072/1‑146 (MQ=255)
 gCCTTGCCGAGCTGGATGCCTTGGGTGGCACCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCA‑‑GGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCcgcg                                                                                                                                        >  5:306194/1‑146 (MQ=255)
 gCCTTGCCGAGCTGGATGCCTTGGGTGGCACCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCA‑‑GGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCgaag                                                                                                                                     <  8:48484/149‑1 (MQ=255)
 gCCTTGCCGAGCTGGATGCCTTGGGTGGCACCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCA‑‑GGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCgaa                                                                                                                                      <  5:127072/148‑1 (MQ=255)
   cTTGCCGAGCTGGATGCCTTGGGTGGCACCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCA‑‑GGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCgaagaa                                                                                                                                   <  7:92425/149‑1 (MQ=255)
   cTTGCCGAGCTGGATGCCTTGGGTGGCACCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCA‑‑GGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCgaagaa                                                                                                                                   <  2:212266/149‑1 (MQ=255)
   attGCCGAGCTGGATGCCTTGGGTGGCACCGATCCGGTGGATGGCAAAGGGTATGAGAGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCA‑‑GGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCgaagaa                                                                                                                                   >  5:73923/2‑149 (MQ=255)
        cGAGCTGCATGCCTTGGGTGGCACCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCA‑‑GGCTGCCTGGAGAAGTTGTAATCCTGTGCCGGCACCTTCGCGGGCACGACAGAGAAGAAGCAGa                                                                                                                              >  6:105888/1‑149 (MQ=255)
            cTGGATGCCTTGGGTGGCACCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCA‑‑GGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGAcgcg                                                                                                                          >  5:294013/1‑149 (MQ=255)
                 tGCCTTGGGTGGCACCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCA‑‑GGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTa                                                                                                                     >  1:11450/1‑149 (MQ=255)
                    cTTGGGTGGCACCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGGGCTGCCTGGAGCAATTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTa                                                                                                                     >  1:107011/1‑148 (MQ=255)
                    cTTGGGTGGCACCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCA‑‑GGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACaa                                                                                                                  >  3:5838/1‑148 (MQ=255)
                             cACCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCA‑‑GGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACgg                                                                                                         >  7:18791/1‑149 (MQ=255)
                              aCCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCATGGAGAAGCTGGGCGGTGACCTGGCA‑‑GGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACgg                                                                                                         >  1:62075/1‑149 (MQ=255)
                                cGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCA‑‑GGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTCCAGCGCGAACGGCAg                                                                                                      >  8:281771/1‑149 (MQ=255)
                                  aTCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCA‑‑GGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGt                                                                                                    <  6:294013/149‑1 (MQ=255)
                                     cGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCA‑‑GGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGc                                                                                                  >  5:194554/1‑148 (MQ=255)
                                      ggTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCA‑‑GGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCaa                                                                                                <  2:223289/149‑1 (MQ=255)
                                      ggTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCA‑‑GGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCaa                                                                                                <  4:307712/149‑1 (MQ=255)
                                       gTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCA‑‑GGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCgg                                                                                                                         >  4:9509/1‑123 (MQ=255)
                                       gTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCA‑‑GGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCgg                                                                                                                         <  3:9509/123‑1 (MQ=255)
                                             ggCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCA‑‑GGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGACCGGCAGGTGCAATGCCAcc                                                                                         >  4:213919/1‑149 (MQ=255)
                                                 aaGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCA‑‑GGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACctgct                                                                                     <  2:1212/149‑1 (MQ=255)
                                                 aaGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCA‑‑GGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACctgct                                                                                     >  6:126863/1‑149 (MQ=255)
                                                  aGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCA‑‑GGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGTTGCAATGCCACctgctg                                                                                    <  6:306194/149‑1 (MQ=255)
                                                      tATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCA‑‑GGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGgcgc                                                                                <  4:163324/149‑1 (MQ=255)
                                                           gcgcTTGCTGGAGAAGCTGGGCGGTGACCTGGCA‑‑GGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGTAGACGCGGTCTACAGCGCGAACGGCAGCTGCAATGCCACCTGCTCGCGCTgagc                                                                           >  4:59183/1‑149 (MQ=255)
                                                           gcgcTTGCTGGAGAAGCTGGGCGGTGACCTGGCA‑‑GGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTgcgc                                                                           <  6:73923/149‑1 (MQ=255)
                                                              cTTGCTGGAGAAGCTGGGCGGTGACCTGGCA‑‑GGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAg                                                                        <  1:299561/149‑1 (MQ=255)
                                                                      gagaAGCTGGGCGGTGACCTGGCA‑‑GGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCATAGCGGGCACGCCCGCGAACAAGCACACGCCGACTACACCGCCAAAGGCACATGCAATCCCACCTGCTGGACCTGCGCCCTACGCACt                                                                 >  2:157213/1‑149 (MQ=255)
                                                                       agaAGCTGGGCGGTGACCTGGCA‑‑GGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACtt                                                                <  1:149568/149‑1 (MQ=255)
                                                                        gaAGCTGGGCGGTGACCTGGCA‑‑GGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTc                                                               >  7:266192/1‑149 (MQ=255)
                                                                        gaAGCTGGGCGGTGACCTGGCA‑‑GGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTc                                                               >  3:115447/1‑149 (MQ=255)
                                                                          aGCTGGGCGGTGACCTGGCA‑‑GGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGTCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCg                                                              <  4:308285/148‑1 (MQ=255)
                                                                           gCTGGGCGGTGACCTGGCA‑‑GGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGaa                                                            <  2:302910/149‑1 (MQ=255)
                                                                            cTGGGCGGTGACCTGGCA‑‑GGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAAc                                                           <  6:194554/149‑1 (MQ=255)
                                                                            cTGGGCGGTGACCTGGCA‑‑GGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAAc                                                           <  3:213919/149‑1 (MQ=255)
                                                                             tGGGCGGTGACCTGGCA‑‑GGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAAc                                                           >  8:263405/1‑148 (MQ=255)
                                                                              gggCGGTGACCTGGCA‑‑GGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAActc                                                         >  8:96716/1‑149 (MQ=255)
                                                                                        cTGGCA‑‑GGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTAGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAAAGCCACCTGCAGGCGCTGCGCCCGACGCACTTAGAACTAGCTCGGATCg                                               >  4:189173/1‑149 (MQ=255)
                                                                                          ggCA‑‑GGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTgcgc                                                                           >  3:128603/1‑119 (MQ=255)
                                                                                          ggCA‑‑GGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTgcgc                                                                           <  4:128603/119‑1 (MQ=255)
                                                                                          ggCA‑‑GGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCAGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTg                                             >  3:130081/1‑149 (MQ=255)
                                                                                           gCA‑‑GGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGt                                              <  5:126863/147‑1 (MQ=255)
                                                                                           gCA‑‑GGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGt                                              <  2:11450/147‑1 (MQ=255)
                                                                                            cA‑‑GGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGa                                           <  3:189173/149‑1 (MQ=255)
                                                                                                ggCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGatc                                         <  5:105888/149‑1 (MQ=255)
                                                                                                 gCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGTGCCAGACGCACTTCGAACTCGCTCGGATCGTGGatca                                        <  1:201240/149‑1 (MQ=255)
                                                                                                 gCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGatca                                        >  2:13999/1‑149 (MQ=255)
                                                                                                       tGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGCGTCGTGGATCATCAcat                                  >  7:253009/1‑149 (MQ=255)
                                                                                                        ggAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGaa                                                            <  8:66734/122‑1 (MQ=255)
                                                                                                        ggAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGaa                                                            >  7:66734/1‑122 (MQ=255)
                                                                                                         gAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGATCATCAcatcc                                <  7:96716/149‑1 (MQ=255)
                                                                                                         gAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCAGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGATCATCAcat                                  <  4:130081/147‑1 (MQ=255)
                                                                                                           gCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACctgctg                                                                                    >  8:154646/1‑95 (MQ=255)
                                                                                                               ttGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCgg                                                   >  1:5652/1‑124 (MQ=255)
                                                                                                               ttGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCgg                                                   <  2:5652/124‑1 (MQ=255)
                                                                                                                 gTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGATCATCACATCCATCCCgg                         <  6:244037/148‑1 (MQ=255)
                                                                                                                 gTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGATCATCACATCCATCCCgg                         >  5:244037/1‑148 (MQ=255)
                                                                                                                  tAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGATCATCACATCCATCCCGGc                        <  2:107011/148‑1 (MQ=255)
                                                                                                                    aGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGATCATCACATCCATCCCGGCGaa                     <  7:281771/149‑1 (MQ=255)
                                                                                                                    aGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGATCATCACATCCATCCCGGCGaa                     <  8:18791/149‑1 (MQ=255)
                                                                                                                      ccTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGATCATCACATCCATCCCGGCGAAc                    >  2:103293/1‑148 (MQ=255)
                                                                                                                       cTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGATCATCACATCCATCCCGGCGAACGa                  >  8:6614/1‑149 (MQ=255)
                                                                                                                          tgCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCACTCGGATCGTGGATCATCACATCCATCCCGGCGAACGAcgc               >  6:75130/1‑147 (MQ=255)
                                                                                                                              ggCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGATCATCACATCCATCCCGGCGAACGACTCGGcc           >  6:9229/1‑149 (MQ=255)
                                                                                                                                 ttcTTCGCGGGCGTGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGa                                                  <  1:80013/105‑1 (MQ=255)
                                                                                                                                 ttcTTCGCGGGCGTGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGa                                                  >  2:80013/3‑107 (MQ=255)
                                                                                                                                 ttcTTCGCGGGCGTGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGATCATCACATCCATCCCGGCGAACGACTCGGCCGc         >  3:189762/3‑148 (MQ=255)
                                                                                                                                  ccTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGATCATCACATCCATCCCGGCGAACGACTCGGCCGCaa       >  8:138521/1‑149 (MQ=255)
                                                                                                                                      cgcgGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGATCATCACATCCATCCCGGCGAACGACTCGGCCGCAATc     <  7:6614/147‑1 (MQ=255)
                                                                                                                                       gcgGGCACGCCCGCGAAGAAGTAGACGCGGTCTACAGCGGGAACGGCAGGTGCAATGCCAGCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGATCATCACATCCATCCCGGCGAACGACTCGGCCGCAATCAgg  <  3:59183/149‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                              
CGCCTTGCCGAGCTGGATGCCTTGGGTGGCACCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCA‑‑GGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGATCATCACATCCATCCCGGCGAACGACTCGGCCGCAATCAGG  >  NC_002947/125977‑126257

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: