Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,641,804 C→A 6.9% 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 6.9% 216.3 / 3.4 70R23S (CGC→AGC) PP_2312lytic transglycosylase family protein
Reads supporting (aligned to +/- strand):  ref base C (30/35);  new base A (2/3);  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.47e-01

TTGCCGTGGTCTTTTACCGTCGTTTATCACTTGGCGGTGCAACCCATGGATGCCGCTCA‑‑‑‑‑‑CCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTGGTGAGGTGAAGGGTGAGCCAGTCGGCATCGA  >  NC_002947/2641657‑2641946
                                                                                                                                                         |                                                                                                                                              
ttGCCGTGGTCTTTTACCGTCGTTTATCACTTGGCGGTGCAACCCATGGATGCCGCTCA‑‑‑‑‑‑CCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGgcgcg                                                                                                                                               >  3:100639/1‑149 (MQ=255)
ttGCCGTGGTCTTTTACCGTCGTTTATCACTTGGCGGTGCAACCCATGGATGCCGCTCA‑‑‑‑‑‑CCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGgcgc                                                                                                                                                <  4:309345/148‑1 (MQ=255)
   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                                                                                              >  3:4906/1‑142 (MQ=255)
                                                        tCA‑‑‑‑‑‑CCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTg                                                                                              <  4:4906/142‑1 (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                                                                                    <  6:312855/135‑1 (MQ=255)
                                                                               gTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTc                                                                                    >  5:312855/1‑135 (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:118707/148‑1 (MQ=255)
                                                                                                  cTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAAc                                                    >  4:340326/1‑148 (MQ=255)
                                                                                                   tGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCTGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCa                                                  <  4:177029/149‑1 (MQ=255)
                                                                                                     gCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCaga                                                <  8:155715/149‑1 (MQ=255)
                                                                                                     gCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCaga                                                <  3:281350/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)
                                                                                                                                                         |                                                                                                                                              
TTGCCGTGGTCTTTTACCGTCGTTTATCACTTGGCGGTGCAACCCATGGATGCCGCTCA‑‑‑‑‑‑CCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTGGTGAGGTGAAGGGTGAGCCAGTCGGCATCGA  >  NC_002947/2641657‑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: