Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,641,807 T→G 5.7% S24A (TCG→GCG)  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,8070TG 5.7% 221.1 / 3.1 70S24A (TCG→GCG) PP_2312lytic transglycosylase family protein
Reads supporting (aligned to +/- strand):  ref base T (30/36);  new base G (2/2);  total (32/38)
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.29e-01

CCGTGGTCTTTTACCGTCGTTTATCACTTGGCGGTGCAACCCATGGATGCCGCTCA‑‑‑‑‑‑CCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTGGTGAGGTGAAGGGTGAGCCAGTCGGCATCGA  >  NC_002947/2641660‑2641946
                                                                                                                                                         |                                                                                                                                           
ccGTGGTCTTTTACCGTCGTTTATCACTTGGCGGTGCAACCCATGGATGCCGCTCA‑‑‑‑‑‑CCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTc                                                                                                                                            >  4:88090/1‑149 (MQ=255)
        ttttACCGTCGTTTATCACTTGGCGGTGCAACCCATGGATGCCGCTCA‑‑‑‑‑‑CCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGc                                                                                                                                    >  3:84651/1‑149 (MQ=255)
         tttACCGTCGTTTATCACTTGGCGGTGCAACCCATGGATGCCGCTCA‑‑‑‑‑‑CCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGcc                                                                                                                                   >  3:118707/1‑149 (MQ=255)
          ttACCGTCGTTTATCACTTGGCGGTGCAACCCATGGATGCCGCTCA‑‑‑‑‑‑CCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGc                                                                                                                                    <  4:34038/147‑1 (MQ=255)
           tACCGTCGTTTATCACTTGGCGGTGCAACCCATGGATGCCGCTCA‑‑‑‑‑‑CCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCCGGCGAGCGCGCCTGGGc                                                                                                                                    <  6:362413/146‑1 (MQ=255)
           tACCGTCGTTTATCACTTGGCGGTGCAACCCATGGATGCCGCTCA‑‑‑‑‑‑CCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCCGGCGAGCGCGCCTGGGc                                                                                                                                    >  5:362413/1‑146 (MQ=255)
             ccGTCGTTTATGACTTGGCGGTGCAACCCATGGATGCCGCTCACAGCCCCCTACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCgggg                                                                                                                                     <  1:119594/149‑1 (MQ=255)
             ccGTCGTTTATCACTTGGCGGTGCAACCCATGGATGCCGCTCA‑‑‑‑‑‑CCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCgcag                                                                                                                               >  3:22047/1‑149 (MQ=255)
              cgtcgtTTATCACTTGGCGGTGCAACCCATGGATGCCGCTCA‑‑‑‑‑‑CCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCgcagc                                                                                                                              <  3:88090/149‑1 (MQ=255)
                  gtTTATCACTTGGCGGTGCAACCCATGGATGCCGCTCA‑‑‑‑‑‑CCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCaaa                                                                                                                           >  7:279758/1‑148 (MQ=255)
                   tttATCACTTGGCGGTGCAACCCATGGATGCCGCTCA‑‑‑‑‑‑CCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAAc                                                                                                                         <  4:100639/149‑1 (MQ=255)
                     aaTCACTTGGCGGTGCAACCCATGGATGCCGCTCA‑‑‑‑‑‑CCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACa                                                                                                                        >  8:1385/2‑148 (MQ=255)
                       tCACTTGGCGGTGCAAGCCATGGATGCCGATCA‑‑‑‑‑‑CCGACCAAGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGAGCTCGCCGGGGCCGCAGCAAAACATAc                                                                                                                     <  4:302283/149‑1 (MQ=255)
                       tCACTTGGCGGTGCAACCCATGGATGCCGCTCA‑‑‑‑‑‑CCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATAc                                                                                                                     <  7:337338/149‑1 (MQ=255)
                         aCTTGGCGGTGCAACCCATGGATGCCGCTCA‑‑‑‑‑‑CCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCg                                                                                                                   >  3:166988/1‑149 (MQ=255)
                         aCTTGGCGGTGCAACCCAGGGATGCCGCTCA‑‑‑‑‑‑CCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATAc                                                                                                                     <  3:238314/147‑1 (MQ=255)
                         aCTTGGCGGTGCAACCAATGGATGCCGCTCA‑‑‑‑‑‑CCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCGCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTACAGCCCAGGCGCGCTCGCCGGGCCCGCAGCAAAACATAc                                                                                                                     >  4:238314/1‑147 (MQ=255)
                          cTTGGCGGTGCAACCCATGGATGCCGCTCA‑‑‑‑‑‑CCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCCGGCGAGCGCGCCTGGGc                                                                                                                                    <  1:142446/131‑1 (MQ=255)
                          cTTGGCGGTGCAACCCATGGATGCCGCTCA‑‑‑‑‑‑CCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAAc                                                                                                                         <  7:275950/142‑1 (MQ=255)
                          cTTGGCGGTGCAACCCATGGATGCCGCTCA‑‑‑‑‑‑CCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAAc                                                                                                                         >  8:275950/1‑142 (MQ=255)
                                ggTGCAACCCATGGATGCCGCTCA‑‑‑‑‑‑CCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGc                                                                                                                                    >  7:257708/1‑125 (MQ=255)
                                ggTGCAACCCATGGATGCCGCTCA‑‑‑‑‑‑CCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGc                                                                                                                                    <  8:257708/125‑1 (MQ=255)
                                  tGCAACCCATGGATGCCGCTCA‑‑‑‑‑‑CCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAg                                                                                                          <  8:279758/149‑1 (MQ=255)
                                            ggATGCCGCTCA‑‑‑‑‑‑CCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGAcc                                                                                                >  8:239105/1‑149 (MQ=255)
                                              aTGCCGCTCA‑‑‑‑‑‑CCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTg                                                                                              <  6:238087/149‑1 (MQ=255)
                                                 ccGCTCA‑‑‑‑‑‑CCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCTGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGgcg                                                                                           >  3:177029/1‑149 (MQ=255)
                                                 ccGCTCA‑‑‑‑‑‑CCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGgc                                                                                            <  5:64434/148‑1 (MQ=255)
                                                     tCA‑‑‑‑‑‑CCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTg                                                                                              <  4:4906/142‑1 (MQ=255)
                                                     tCA‑‑‑‑‑‑CCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTg                                                                                              >  3:4906/1‑142 (MQ=255)
                                                       a‑‑‑‑‑‑CCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGAtt                                                                                     >  5:13031/1‑149 (MQ=255)
                                                       a‑‑‑‑‑‑CCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGAt                                                                                      <  2:79281/148‑1 (MQ=255)
                                                                gACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCg                                                                                   >  4:281350/1‑148 (MQ=255)
                                                                   cACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCAc                                                                                >  6:305139/1‑148 (MQ=255)
                                                                      ggACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAg                                                                             >  1:319809/1‑148 (MQ=255)
                                                                          acGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAg                                                                        <  6:13031/149‑1 (MQ=255)
                                                                            gTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCaaaa                                                                                                                          >  7:170278/1‑97 (MQ=255)
                                                                            gTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCaaaa                                                                                                                          <  8:170278/97‑1 (MQ=255)
                                                                            gTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTc                                                                                    >  5:312855/1‑135 (MQ=255)
                                                                            gTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTc                                                                                    <  6:312855/135‑1 (MQ=255)
                                                                             ttttCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTAc                                                                     >  4:340250/1‑149 (MQ=255)
                                                                                tCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTg                                                                   <  2:289068/148‑1 (MQ=255)
                                                                                tCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGc                                                                  >  6:72058/1‑149 (MQ=255)
                                                                                    aTGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGt                                                              <  4:166988/149‑1 (MQ=255)
                                                                                    aTGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCATGACCTGGCGCAGATTCGCGCCAGCAAAGTACTGCGAGt                                                              >  6:336002/1‑149 (MQ=255)
                                                                                         gCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTgg                                                         <  4:8510/149‑1 (MQ=255)
                                                                                         gCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTgg                                                         <  4:84651/149‑1 (MQ=255)
                                                                                             tGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCaa                                                     <  7:1385/149‑1 (MQ=255)
                                                                                               cTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAAc                                                    >  4:340326/1‑148 (MQ=255)
                                                                                               cTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAAc                                                    <  4:118707/148‑1 (MQ=255)
                                                                                                tGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCTGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCa                                                  <  4:177029/149‑1 (MQ=255)
                                                                                                  gCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCaga                                                <  3:281350/149‑1 (MQ=255)
                                                                                                  gCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCaga                                                <  8:155715/149‑1 (MQ=255)
                                                                                                    tCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGc                                              >  6:157398/1‑149 (MQ=255)
                                                                                                    tCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGc                                              >  4:99462/1‑149 (MQ=255)
                                                                                                     cATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGcc                                             >  5:261353/1‑149 (MQ=255)
                                                                                                      aTTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCg                                            <  5:305139/149‑1 (MQ=255)
                                                                                                        ttttCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGa                                                                                       >  1:165293/1‑104 (MQ=255)
                                                                                                        ttttCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGa                                                                                       <  2:165293/104‑1 (MQ=255)
                                                                                                         tttCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCa                                          <  2:227729/148‑1 (MQ=255)
                                                                                                          ttCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCaa                                         <  4:22047/148‑1 (MQ=255)
                                                                                                                  aCGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCt                                  <  7:239105/147‑1 (MQ=255)
                                                                                                                  aCGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCt                                  <  3:340326/147‑1 (MQ=255)
                                                                                                                  aCGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTgg                                <  5:72058/149‑1 (MQ=255)
                                                                                                                    gCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTggt                               >  5:112795/1‑148 (MQ=255)
                                                                                                                               ccTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTGGTGAGGTGAAGGGt                   >  5:274878/1‑149 (MQ=255)
                                                                                                                                cTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTGGTGAGGTGAAGGGTg                  <  3:340250/149‑1 (MQ=255)
                                                                                                                                    gCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATTCCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTGGTGAGGTGAAGGGTGAGc               <  6:274878/148‑1 (MQ=255)
                                                                                                                                       ggTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTGGTGAGGTGAAGGGTGAGCCAGt           >  6:336748/1‑149 (MQ=255)
                                                                                                                                             gCCCCGGCGAGCGCGCCTGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTGGTGAGGTGAAGGGTGAGCCAGTCGGCAt     >  7:91478/1‑149 (MQ=255)
                                                                                                                                                cAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCTCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTGGTGAGGTGAAGGGTGAGCCAGTCGGCATCGa  <  5:157398/149‑1 (MQ=255)
                                                                                                                                                         |                                                                                                                                           
CCGTGGTCTTTTACCGTCGTTTATCACTTGGCGGTGCAACCCATGGATGCCGCTCA‑‑‑‑‑‑CCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTGGTGAGGTGAAGGGTGAGCCAGTCGGCATCGA  >  NC_002947/2641660‑2641946

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: