Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,641,804 C→A 6.7% 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.7% 248.6 / 5.8 85R23S (CGC→AGC) PP_2312lytic transglycosylase family protein
Reads supporting (aligned to +/- strand):  ref base C (39/40);  new base A (2/4);  total (41/44)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.09e-01

GTTGCCGTGGTCTTTTACCGTCGTTTATCACTTGGCGGTGCAACCCATGGATGCCGCTCACCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTGGTGAGGTGAAGGGTGAGCCAGTCGGCATCGAGTATTA  >  NC_002947/2641656‑2641952
                                                                                                                                                    |                                                                                                                                                    
gTTGCCGTTGTCTTTTACCGTCGTTTATCACTTGGCGGTGCAACCCATGGATGCCGCTCACCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTTGCCTTTGCCGGTCCCGCCCAGgcgc                                                                                                                                                      <  5:314247/149‑1 (MQ=255)
 ttGCCGTGGTCTTTTACCGTCGTTTATCACTTGGCGGTGCAACCCATGGATGCCGCTCACCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGgcgc                                                                                                                                                      <  6:322467/148‑1 (MQ=255)
   gCCGTGGTCTTTTACCGTCGTTTATCACTTGGCGGTGCAACCCATGGATGCCGCTCACCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCt                                                                                                                                                   >  1:149397/1‑149 (MQ=255)
    ccGTGGTCTTTTACCGTCGTTTATCACTTGGCGGTGCAACCCATGGATGCCGCTCACCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTc                                                                                                                                                  >  3:208998/1‑149 (MQ=255)
          tCTTTTACCGTCGTTTATCACTTGGCGGTGCAACCCATGGATGCCGCTCACCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCCGGCGAGCGCGCCtgg                                                                                                                                            <  2:154751/149‑1 (MQ=255)
            ttttACCGTCGTTTATCACTTGGCGGTGCAACCCATGGATGCCGCTCACCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGgcgc                                                                                                                                          >  3:276344/1‑149 (MQ=255)
            ttttACCGTCGTTTATCACTTGGCGGTGCAACCCATGGATGCCGCTCACCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGc                                                                                                                                          >  7:197043/1‑149 (MQ=255)
              ttACCGTCGTTTATCACTTGGCGGTGCAACCCATGGATGCCGCTCACCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCCGGCGAGCGCGCCTGGGcg                                                                                                                                         <  6:331336/148‑2 (MQ=255)
              ttACCGTCGTTTATCACTTGGCGGTGCAACCCATGGATGCCGCTCACCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTATTCTGATGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCg                                                                                                                                        >  7:245905/1‑149 (MQ=255)
                 ccGTCGTTTATCACTTGGCGGTGCAACCCATGGATGCCGCTCACCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCgcag                                                                                                                                     >  2:256806/1‑149 (MQ=255)
                 ccGTCGTTTATCACTTGGCGGTGCAACCCATGGATGCCGCTCACCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCgcag                                                                                                                                     >  1:252373/1‑149 (MQ=255)
                       tttATCACTTGGCGGTGCAACCCATGGATGCCGCTCACCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCaaaa                                                                                                                                >  1:68201/1‑148 (MQ=255)
                          aTCACTTGGCGGTGCAACCCATGGATGCCGCTCACCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATa                                                                                                                            >  4:241394/1‑149 (MQ=255)
                          aTCACTTGGCGGTGCAACCCATGGATGCCGCTCACCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGCGGTTGGCCTTTGCCGGTCCCGCCCATGCGCGCTCGCCGGGGCCGCAGCAAACCAt                                                                                                                             >  4:236802/1‑148 (MQ=255)
                                       gCAACCCATGGATGCCGCTCACCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAAc                                                                                                                               >  2:349896/1‑133 (MQ=255)
                                       gCAACCCATGGATGCCGCTCACCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAAc                                                                                                                               <  1:349896/133‑1 (MQ=255)
                                          aCCCATGGATGCCGCTCACCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCCGGCGAGCGCGCCTGGGcg                                                                                                                                         >  6:47928/1‑119 (MQ=255)
                                          aCCCATGGATGCCGCTCACCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCCGGCGAGCGCGCCTGGGcg                                                                                                                                         <  5:47928/120‑2 (MQ=255)
                                             cATGGATGCCGCTCACCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGt                                                                                                          >  6:206129/1‑148 (MQ=255)
                                             cATGGATGCCGCTCACCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGt                                                                                                          >  5:97141/1‑148 (MQ=255)
                                                 gATGCCGCTCACCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCt                                                                                                     >  5:365243/1‑149 (MQ=255)
                                                 gATGCCGCTCACCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCt                                                                                                     >  7:247597/1‑149 (MQ=255)
                                                  aTGCCGCTCACCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCggg                                                                                                                                            >  6:175899/1‑109 (MQ=255)
                                                  aTGCCGCTCACCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCggg                                                                                                                                            <  5:175899/109‑1 (MQ=255)
                                                     ccGCTCACCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGgc                                                                                                  <  7:83152/148‑1 (MQ=255)
                                                      cGCTCACCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGgcgc                                                                                                <  3:341130/149‑1 (MQ=255)
                                                       gCTCACCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGgcg                                                                                                 >  3:336081/1‑147 (MQ=255)
                                                       gCTCACCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGgcg                                                                                                 >  7:268165/1‑147 (MQ=255)
                                                              gACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAg                                                                                                                >  8:328710/1‑125 (MQ=255)
                                                              gACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAg                                                                                                                <  7:328710/125‑1 (MQ=255)
                                                              gACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCg                                                                                         >  6:297041/1‑148 (MQ=255)
                                                                 cACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCAc                                                                                      >  8:15015/1‑148 (MQ=255)
                                                                  aCGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCa                                                                                    <  1:13124/149‑1 (MQ=255)
                                                                    ggACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAg                                                                                   >  1:212225/1‑148 (MQ=255)
                                                                    ggACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAg                                                                                   >  4:293193/1‑148 (MQ=255)
                                                                      acacGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCCGGCGAGCGCGcc                                                                                                                                               <  5:366616/86‑1 (MQ=255)
                                                                      acacGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCCGGCGAGCGCGcc                                                                                                                                               >  6:366616/1‑86 (MQ=255)
                                                                          gTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGATCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTa                                                                            >  2:163685/1‑149 (MQ=255)
                                                                              tCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTg                                                                         <  1:256806/148‑1 (MQ=255)
                                                                              tCGCATGCAGCGATGATGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGAAAAACATACCGCCTGCCAAGGTGCGTGACATGGCGCAGATTCGCACCAGCAAAGTACTGc                                                                        >  6:321295/1‑149 (MQ=255)
                                                                               cGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGc                                                                        <  2:149397/148‑1 (MQ=255)
                                                                               cGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGc                                                                        <  4:208998/148‑1 (MQ=255)
                                                                               cGCATGCAGAGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTTGCCTTTGCCGGTACCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGTCTCAGATTCGCACCAGCAAAGTACTGc                                                                        <  5:321295/148‑1 (MQ=255)
                                                                                gCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGa                                                                      <  3:293193/149‑1 (MQ=255)
                                                                                gCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGa                                                                      <  6:150723/149‑1 (MQ=255)
                                                                                gCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGa                                                                      <  1:178532/149‑1 (MQ=255)
                                                                                 cATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAg                                                                     >  5:346123/1‑149 (MQ=255)
                                                                                        cGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAg                                                                                              <  7:284617/117‑1 (MQ=255)
                                                                                        cGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAg                                                                                              >  8:284617/1‑117 (MQ=255)
                                                                                        cGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGt                                                              <  2:68201/149‑1 (MQ=255)
                                                                                         gATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTc                                                             <  8:197043/149‑1 (MQ=255)
                                                                                         gATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTc                                                             <  7:191161/149‑1 (MQ=255)
                                                                                             cTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAAcc                                                         <  5:38734/149‑1 (MQ=255)
                                                                                             cTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAAc                                                          <  6:97141/148‑1 (MQ=255)
                                                                                                gCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCa                                                        >  8:339262/1‑147 (MQ=255)
                                                                                                 cTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCagag                                                     <  8:268165/149‑1 (MQ=255)
                                                                                                   cATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAAc                                                          >  2:330430/1‑142 (MQ=255)
                                                                                                   cATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAAc                                                          <  1:325105/142‑1 (MQ=255)
                                                                                                   cATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAAc                                                          >  2:325105/1‑142 (MQ=255)
                                                                                                   cATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAAc                                                          <  1:330430/142‑1 (MQ=255)
                                                                                                            gctgACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGtt                                          >  6:75499/1‑149 (MQ=255)
                                                                                                             ctgACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTg                                                                   >  1:184094/1‑123 (MQ=255)
                                                                                                             ctgACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTg                                                                   <  2:184094/123‑1 (MQ=255)
                                                                                                             ctgACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTc                                         <  6:346123/149‑1 (MQ=255)
                                                                                                              tgACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCt                                        >  6:201806/1‑149 (MQ=255)
                                                                                                                aCGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCt                                        <  5:75499/147‑1 (MQ=255)
                                                                                                                aCGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTg                                       <  7:15015/148‑1 (MQ=255)
                                                                                                                  gCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTggtg                                    <  5:165437/149‑1 (MQ=255)
                                                                                                                    ttgggTTGGCCTTTTCCGGTCCCGCCCAGGCGCGGTCGCCGGGGCCGCATCAAAACATACCGCCTGCCAAGGTTCGTGACCTGGCGCAGTTTCCCAACAGCAAAGTACTGCGAGTGCTGGTCAACTAGAGCCGCAACAGTTCTggtgag                                  <  7:227424/147‑1 (MQ=255)
                                                                                                                      gggTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAAGATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTggtgag                                  <  3:236802/147‑1 (MQ=255)
                                                                                                                      gggTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTggtgag                                  <  7:339262/147‑1 (MQ=255)
                                                                                                                      gggTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTggtgag                                  <  5:297041/147‑1 (MQ=255)
                                                                                                                      gggTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTggtgag                                  <  3:241394/147‑1 (MQ=255)
                                                                                                                       ggTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTggtgaggtg                               <  2:212225/149‑1 (MQ=255)
                                                                                                                       ggTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTggtgaggtg                               <  5:206129/149‑1 (MQ=255)
                                                                                                                        gTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTCGCGCAGCTTCGCAAGAGCAAAGTACAGCGAGTGCTGGTCAAGCAGAGCCGCAACAg                                            >  7:237145/1‑135 (MQ=255)
                                                                                                                             ccTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTc                                                             >  3:69578/1‑113 (MQ=255)
                                                                                                                             ccTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTc                                                             <  4:69578/113‑1 (MQ=255)
                                                                                                                                   ccGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTGGTGAGGTGAAGGGTGAGc                     <  3:344399/147‑1 (MQ=255)
                                                                                                                                   ccGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTGGTGAGGTGAAGGGTGAGc                     >  4:344399/1‑147 (MQ=255)
                                                                                                                                      gTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTGGTGAGGTGAAGGGTGAGCCAGTc                >  4:346128/1‑149 (MQ=255)
                                                                                                                                       tCCCGGCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCAGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTGGTGAGGGGAAGGGTGAGCCAGTCg               <  8:245905/149‑1 (MQ=255)
                                                                                                                                         ccGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTGGTGAGGTGAAGGGTGAGCCAGTCGGc             >  5:197437/1‑149 (MQ=255)
                                                                                                                                            cccAGACGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTGGTGAGGTGAAGGGTGAGCCAGTCGGCATc          <  8:147160/149‑1 (MQ=255)
                                                                                                                                             ccAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTGGTGAGGTGAAGGGTGAGCCAGTCGGCATc          <  1:162487/148‑1 (MQ=255)
                                                                                                                                             ccAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTGGTGAGGTGAAGGGTGAGCCAGTCGGCATc          >  2:162487/1‑148 (MQ=255)
                                                                                                                                              cAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTGGTGAGGTGAAGGGTGAGCCAGTCGGCATCGa        <  5:201806/149‑1 (MQ=255)
                                                                                                                                                 gcgcgcTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTGGTGAGGTGAAGGGTGAGCCAGTCGGCATc          <  8:135438/144‑1 (MQ=255)
                                                                                                                                                 gcgcgcTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTGGTGAGGTGAAGGGTGAGCCAGTCGGCATc          >  7:135438/1‑144 (MQ=255)
                                                                                                                                                    cgcTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTGGTGAGGTGAAGGGTGAGCCAGTCGGCATCGAGtatta  >  5:127766/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GTTGCCGTGGTCTTTTACCGTCGTTTATCACTTGGCGGTGCAACCCATGGATGCCGCTCACCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTGGTGAGGTGAAGGGTGAGCCAGTCGGCATCGAGTATTA  >  NC_002947/2641656‑2641952

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: