Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,641,804 C→A 5.6% R23S (CGC→AGC)  PP_2312 → lytic transglycosylase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,641,8040CA 5.6% 230.7 / 2.8 71R23S (CGC→AGC) PP_2312lytic transglycosylase family protein
Reads supporting (aligned to +/- strand):  ref base C (30/37);  new base A (2/2);  total (32/39)
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.08e-01

TTGCCGTGGTCTTTTACCGTCGTTTATCACTTGGCGGTGCAACCCATGGATGCCGCTCACCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCG‑‑‑‑CCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTGGTGAGGTGAAGGGTGAGCCAGTCGGCATCGAGTA  >  NC_002947/2641657‑2641949
                                                                                                                                                   |                                                                                                                                                     
ttGCCGTGGTCTTTTACCGTCGTTTATCACTTGGCGGTGCAACGCATGGATGCCGCTCACCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGgcgc                                                                                                                                                       <  8:284979/148‑1 (MQ=255)
ttGCCGTGGTCTTTTACCGTCGTTTATCACTTGGCGGTGCAACCCATGGATGCCGCTCACCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGgcgc                                                                                                                                                       >  7:284979/1‑148 (MQ=255)
   ccGTGGTCTTTTACCGTCGTTTATCACTTGGCGGTGCAACCCATGGATGCCGCTCACCGACCACGGACACGTTTTCGCATGCAGCGATGCTGACTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTc                                                                                                                                                   >  1:291428/1‑149 (MQ=255)
    cGTGGTCTTTTACCGTCGTTTATCACTTGGCGGTGCAACCCATGGATGCCGCTCACCGACCACGGACACGTGTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCCGGCGAGcgcg                                                                                                                                                  >  1:197243/1‑149 (MQ=255)
       ggTCTTTTACCGTCGCTTATCACTGCGCGGGGCAACCCATGGATCCCGCTCACCGACCACCGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCGGTGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCG‑‑‑‑CCg                                                                                                                                           <  8:187232/149‑1 (MQ=255)
        gTCTTTTACCGTCGTTTATCACTTGGCGGTGCAACCCATGGATGCCGCTCACCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCG‑‑‑‑CCgg                                                                                                                                          <  8:113336/149‑1 (MQ=255)
        gTCGTTTACCGTCGTTTATCACTTGGCGGTGCAACCCATGGATGCCGCTCACCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCG‑‑‑‑CCgg                                                                                                                                          <  3:260643/149‑1 (MQ=255)
         tCTTTTACCGTCGTTTATCACTTGGCGGTGCAACCCATGGATGCCGCTCACCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCG‑‑‑‑CCggg                                                                                                                                         <  6:337653/149‑1 (MQ=255)
             ttACCGTCGTTTATCACTTGGCGGTGCAACCCATGGATGCCGCTCACCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCG‑‑‑‑CCgggg                                                                                                                                        <  1:346579/146‑1 (MQ=255)
             ttACCGTCGTTTATCACTTGGCGGTGCAACCCATGGATGCCGCTCACCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCG‑‑‑‑CCGGGGcc                                                                                                                                      <  3:1890/148‑1 (MQ=255)
             ttACCGTCGTTTATCACTTGGCGGTGCAACCCATGGATGCCGCTCACCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCG‑‑‑‑CCGGGGc                                                                                                                                       <  1:43092/147‑1 (MQ=255)
                ccGTCGTTTATCACTTGGCGGTGCAACCCATGGATGCCGCTCACCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCCGGCGAGCGCG‑‑‑‑CCTGGGcggga                                                                                                                                   <  2:197243/148‑5 (MQ=255)
                ccGTCGTTTATCACTTGGCGGTGCAACCCATGGATGCCGCTCACCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCG‑‑‑‑CCGGGGCCgcag                                                                                                                                  <  3:301939/149‑1 (MQ=255)
                ccGTCGTTTATCACTTGGCGGTGCAACCCATGGATGCCGCTCACCGACCACGGACACGTTTTCGCATGCAGCGATGCTGACTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCG‑‑‑‑CCGGGGCCgca                                                                                                                                   <  2:291428/148‑1 (MQ=255)
                ccGTCGTTTATCACTTGGCGGTGCAACCCATGGATGCCGCTCACCGACCACGGACACGGTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCG‑‑‑‑CCGGGGCCgca                                                                                                                                   <  6:248192/148‑1 (MQ=255)
                 cgacgtTTATCACAACCCGGTGCAACCCATGGATGCCGCTCACCGACCACCGACACCTTTTCGCATGCAGCCATGCACGCACATTTTTCTGCTGACGCTGGGGTTGGCCTATGCAGGTCCCGCCCAGGCGCGCTCG‑‑‑‑CCGGGGCCccag                                                                                                                                  >  6:391583/4‑148 (MQ=255)
                      tttATCACTTGGCGGTGCAACCCATGGATGCCGCTCACCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCG‑‑‑‑CCGGGGCCGCAGCaaaa                                                                                                                             >  2:122380/1‑148 (MQ=255)
                         aTCACTTGGCGGTGCAACCCATGGATGCCGCTCACCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCG‑‑‑‑CCGGGGCCGCAGCAAAACAt                                                                                                                          <  2:350537/148‑1 (MQ=255)
                                 gCGGTGCAACCCATGGATGCCGCTCACCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCG‑‑‑‑CCGGGGCCGCAGCAAAACATACCGCCTGc                                                                                                                 >  1:290248/1‑149 (MQ=255)
                                     tGCAACCCATGGATGCCGCTCCCCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCG‑‑‑‑CCGGGGCCGCTGCAAAACATACCGCCTGCCAAg                                                                                                             >  8:198544/1‑149 (MQ=255)
                                        aaCCCATGGATGCCGCTCACCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCG‑‑‑‑CCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTg                                                                                                          >  2:250253/1‑149 (MQ=255)
                                                  tGCCGCTCACCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCG‑‑‑‑CCGGGGCCGCAGCAAAACAt                                                                                                                          <  5:340981/123‑1 (MQ=255)
                                                  tGCCGCTCACCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCG‑‑‑‑CCGGGGCCGCAGCAAAACAt                                                                                                                          >  6:340981/1‑123 (MQ=255)
                                                    ccGCTCACCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCG‑‑‑‑CCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGgc                                                                                               <  2:290248/148‑1 (MQ=255)
                                                          aCCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCAGGTCCCGCCCAGGCGCGCTCG‑‑‑‑CCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGAtt                                                                                        >  5:62038/1‑149 (MQ=255)
                                                            cGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCG‑‑‑‑CCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCg                                                                                      >  7:152836/1‑149 (MQ=255)
                                                             gACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCG‑‑‑‑CCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGc                                                                                     <  3:100691/149‑1 (MQ=255)
                                                                  cGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCG‑‑‑‑CCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAg                                                                                <  1:250253/149‑1 (MQ=255)
                                                                     acacGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCG‑‑‑‑CCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCa                                                                              >  1:246483/1‑148 (MQ=255)
                                                                      cacGTTTTCGCATGCAGCGATGCTGGGTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCG‑‑‑‑CCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCAc                                                                                   >  1:18555/1‑142 (MQ=255)
                                                                      cacGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCG‑‑‑‑CCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCAc                                                                                   <  2:18555/142‑1 (MQ=255)
                                                                         gTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCG‑‑‑‑CCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTa                                                                         >  5:52864/1‑149 (MQ=255)
                                                                             tagCATGCAGCGATGCTGGCTCATTTTTCTGCTTACGCTTGTGTTGGCCTTTGCCGGTCGCGCCCAGGCGCGCTCG‑‑‑‑CCGGTGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGTCGCAGATTCGCACCAGCAAAGTACTg                                                                      <  3:16515/146‑1 (MQ=255)
                                                                              cGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCG‑‑‑‑CCGGGGCCGCAGCAAAACATACCGCCTGCCaa                                                                                                              <  5:84927/107‑1 (MQ=255)
                                                                              cGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCG‑‑‑‑CCGGGGCCGCAGCAAAACATACCGCCTGCCaa                                                                                                              >  6:84927/1‑107 (MQ=255)
                                                                              cGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCG‑‑‑‑CCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGc                                                                     >  5:158865/1‑148 (MQ=255)
                                                                                cATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCG‑‑‑‑CCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAg                                                                  >  8:161018/1‑149 (MQ=255)
                                                                                  tGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCG‑‑‑‑CCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTg                                                                >  5:39984/1‑149 (MQ=255)
                                                                                  agCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCG‑‑‑‑CCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTTCGAGTg                                                                <  5:391583/148‑1 (MQ=255)
                                                                                     aGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCG‑‑‑‑CCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTc                                                                                       >  3:167356/1‑123 (MQ=255)
                                                                                     aGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCG‑‑‑‑CCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTc                                                                                       <  4:167356/123‑1 (MQ=255)
                                                                                     aGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCG‑‑‑‑CCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTg                                                             <  5:18096/149‑1 (MQ=255)
                                                                                     aGCGATGCTGGCTCATTTTTCTGCTGACGCGGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCG‑‑‑‑CCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGAGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTg                                                             >  1:122984/1‑149 (MQ=255)
                                                                                             tGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCG‑‑‑‑CCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAAc                                                       >  7:200363/1‑147 (MQ=255)
                                                                                               gCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCG‑‑‑‑CCGGGGCCGCAGCTAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCa                                                     >  3:34713/1‑147 (MQ=255)
                                                                                               gCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCG‑‑‑‑CCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCaga                                                   <  6:158865/149‑1 (MQ=255)
                                                                                               gCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCG‑‑‑‑CCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCaga                                                   <  6:39984/149‑1 (MQ=255)
                                                                                                cTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCG‑‑‑‑CCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGAt                                                                                         >  7:252904/1‑110 (MQ=255)
                                                                                                cTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCG‑‑‑‑CCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGAt                                                                                         <  8:252904/110‑1 (MQ=255)
                                                                                                   aTTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCG‑‑‑‑CCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGc                                                 <  8:152836/147‑1 (MQ=255)
                                                                                                         ctgctgACGCTGGGGTTGTCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCCTCCCCGGCCGTCCCCATACTTACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGc                                              <  1:122380/148‑1 (MQ=255)
                                                                                                         ctgctgACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCG‑‑‑‑CCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACa                                          <  1:148203/148‑1 (MQ=255)
                                                                                                            ctgACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCG‑‑‑‑CCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTc                                      <  2:311972/149‑1 (MQ=255)
                                                                                                                cGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCG‑‑‑‑CCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTggt                                  <  3:264539/149‑1 (MQ=255)
                                                                                                                    ggggTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCG‑‑‑‑CCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTggtgagg                              >  8:306354/1‑149 (MQ=255)
                                                                                                                           gCCTTTGCCGGTCCCGCCCAGGCGCGCTCG‑‑‑‑CCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTGGTGAGGTGAAgg                        <  6:52864/148‑1 (MQ=255)
                                                                                                                            ccTTTGCCGGTCCCGCCCAGGCGCGCTCG‑‑‑‑CCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTGGTGAGGTGAAGGGt                      >  4:347431/1‑149 (MQ=255)
                                                                                                                             cTTTGCCGGTCCCGCCCAGGCGCGCTCG‑‑‑‑CCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTGGTGAGGTGAAGGGt                      >  7:234228/1‑148 (MQ=255)
                                                                                                                               ttGCCGGTCCCGCCCAGGCGCGCTCG‑‑‑‑CCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTGGTGAGGTGAAGGGTGAg                   >  2:194779/1‑149 (MQ=255)
                                                                                                                                tGCCGGTCCCGCCCAGGCGCGCTCG‑‑‑‑CCGGGGCCGCAGGAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTGGTGAGGTGAAGGGTGa                    <  2:246483/147‑1 (MQ=255)
                                                                                                                                tGCCGGTCCCGCCCAGGCGCGCTCG‑‑‑‑CCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTGGTGAGGTGAAGGGTGAg                   <  8:200363/148‑1 (MQ=255)
                                                                                                                                tGCCGGTCCCGCCCAGGCGCGCTCG‑‑‑‑CCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTGGTGAGGTGAAGGGTGAg                   <  4:34713/148‑1 (MQ=255)
                                                                                                                                tGCCGGTCCCGCCCAGGCGCGCTCG‑‑‑‑CCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTGGTGAGGTGAAGGGTGAg                   <  7:306354/148‑1 (MQ=255)
                                                                                                                                  ccGGTCCCGCCCAGGCGCGCTCG‑‑‑‑CCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTGGTGAGGTGAAGGGTGAGc                  >  2:363445/1‑147 (MQ=255)
                                                                                                                                    ggTCCCGCCCAGGCGCGCTCG‑‑‑‑CCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTGGTGAGGTGAAGGGTGAGCCAGt              >  2:180814/1‑149 (MQ=255)
                                                                                                                                    ggTCCCGCCCAGGCGCGCTCG‑‑‑‑CCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTGGTGAGGTGAAGGGTGAGCCAGt              >  5:216606/1‑149 (MQ=255)
                                                                                                                                      tCCCGCCCAGGCGCGCTCG‑‑‑‑CCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTGGTGAGGTGAAGGGTGAGCCAGTCg            <  3:347431/149‑1 (MQ=255)
                                                                                                                                        ccGCCCAGGCGCGCTCG‑‑‑‑CCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTGGTGAGGTGAAGGGTGAGCCAGTCGGc          <  6:216606/149‑1 (MQ=255)
                                                                                                                                        ccGCCCAGGCGCGCTCG‑‑‑‑CCGGGCCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGAGTCGCACCAGCAAAGTACTGCGAGTCCTGGTCAACCCGAGCCGCAACAGTTCTGGTGAGGTGAAGGGTGAGCCAGTCgg           >  4:159814/1‑148 (MQ=255)
                                                                                                                                         cGCCCAGGCGCGCTCG‑‑‑‑CCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTGGTGAGGTGAAGGGTGAGCCAGTCGGCa         >  7:279535/1‑149 (MQ=255)
                                                                                                                                          gCCCCGGCGAGCGCG‑‑‑‑CCTGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTGGTGAGGTGAAGGGt                      <  3:214958/135‑1 (MQ=255)
                                                                                                                                          gCCCCGGCGAGCGCG‑‑‑‑CCTGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTGGTGAGGTGAAGGGt                      >  4:214958/1‑135 (MQ=255)
                                                                                                                                               ggCGCGCTCG‑‑‑‑CCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTGGTGAGGTGAAGGGTGAGCCAGTCGGCATCGAGt   <  1:180814/149‑1 (MQ=255)
                                                                                                                                                gcgcgcTCG‑‑‑‑CCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTGGTGAGGTGAAGGGTGAGCCAGTCGGCATCGAGta  <  7:357402/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TTGCCGTGGTCTTTTACCGTCGTTTATCACTTGGCGGTGCAACCCATGGATGCCGCTCACCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCG‑‑‑‑CCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTGGTGAGGTGAAGGGTGAGCCAGTCGGCATCGAGTA  >  NC_002947/2641657‑2641949

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: