Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,641,799 A→C 6.9% Q21P (CAG→CCG)  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,7990AC 6.9% 264.6 / 8.8 87Q21P (CAG→CCG) PP_2312lytic transglycosylase family protein
Reads supporting (aligned to +/- strand):  ref base A (39/42);  new base C (2/4);  total (41/46)
Fisher's exact test for biased strand distribution p-value = 6.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.82e-01

GTTGCCGTGGTCTTTTACCGTCGTTTATCACTTGGCGGTGCAACCCATGGATGCCGCTCACCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTGGTGAGGTGAAGGGTGAGCCAGTCGGCATCGA  >  NC_002947/2641656‑2641946
                                                                                                                                               |                                                                                                                                                   
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                                                                                                                         <  1:349896/133‑1 (MQ=255)
                                       gCAACCCATGGATGCCGCTCACCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAAc                                                                                                                         >  2:349896/1‑133 (MQ=255)
                                          aCCCATGGATGCCGCTCACCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCCGGCGAGCGCGCCTGGGcg                                                                                                                                   <  5:47928/120‑2 (MQ=255)
                                          aCCCATGGATGCCGCTCACCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCCGGCGAGCGCGCCTGGGcg                                                                                                                                   >  6:47928/1‑119 (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                                                                                                                                      <  5:175899/109‑1 (MQ=255)
                                                  aTGCCGCTCACCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCggg                                                                                                                                      >  6:175899/1‑109 (MQ=255)
                                                     ccGCTCACCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGgc                                                                                            <  7:83152/148‑1 (MQ=255)
                                                      cGCTCACCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGgcgc                                                                                          <  3:341130/149‑1 (MQ=255)
                                                       gCTCACCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGgcg                                                                                           >  7:268165/1‑147 (MQ=255)
                                                       gCTCACCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGgcg                                                                                           >  3:336081/1‑147 (MQ=255)
                                                              gACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAg                                                                                                          <  7:328710/125‑1 (MQ=255)
                                                              gACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAg                                                                                                          >  8:328710/1‑125 (MQ=255)
                                                              gACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCg                                                                                   >  6:297041/1‑148 (MQ=255)
                                                                 cACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCAc                                                                                >  8:15015/1‑148 (MQ=255)
                                                                  aCGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCa                                                                              <  1:13124/149‑1 (MQ=255)
                                                                    ggACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAg                                                                             >  4:293193/1‑148 (MQ=255)
                                                                    ggACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAg                                                                             >  1:212225/1‑148 (MQ=255)
                                                                      acacGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCCGGCGAGCGCGcc                                                                                                                                         >  6:366616/1‑86 (MQ=255)
                                                                      acacGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCCGGCGAGCGCGcc                                                                                                                                         <  5:366616/86‑1 (MQ=255)
                                                                          gTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGATCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTa                                                                      >  2:163685/1‑149 (MQ=255)
                                                                              tCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTg                                                                   <  1:256806/148‑1 (MQ=255)
                                                                              tCGCATGCAGCGATGATGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGAAAAACATACCGCCTGCCAAGGTGCGTGACATGGCGCAGATTCGCACCAGCAAAGTACTGc                                                                  >  6:321295/1‑149 (MQ=255)
                                                                               cGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGc                                                                  <  4:208998/148‑1 (MQ=255)
                                                                               cGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGc                                                                  <  2:149397/148‑1 (MQ=255)
                                                                               cGCATGCAGAGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTTGCCTTTGCCGGTACCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGTCTCAGATTCGCACCAGCAAAGTACTGc                                                                  <  5:321295/148‑1 (MQ=255)
                                                                                gCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGa                                                                <  6:150723/149‑1 (MQ=255)
                                                                                gCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGa                                                                <  3:293193/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                                                       <  7:191161/149‑1 (MQ=255)
                                                                                         gATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTc                                                       <  8:197043/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:325105/1‑142 (MQ=255)
                                                                                                   cATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAAc                                                    >  2:330430/1‑142 (MQ=255)
                                                                                                   cATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAAc                                                    <  1:330430/142‑1 (MQ=255)
                                                                                                   cATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAAc                                                    <  1:325105/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                            <  3:241394/147‑1 (MQ=255)
                                                                                                                      gggTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTggtgag                            <  5:297041/147‑1 (MQ=255)
                                                                                                                       ggTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTggtgaggtg                         <  5:206129/149‑1 (MQ=255)
                                                                                                                       ggTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTggtgaggtg                         <  2:212225/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    >  2:162487/1‑148 (MQ=255)
                                                                                                                                             ccAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTGGTGAGGTGAAGGGTGAGCCAGTCGGCATc    <  1:162487/148‑1 (MQ=255)
                                                                                                                                              cAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTGGTGAGGTGAAGGGTGAGCCAGTCGGCATCGa  <  5:201806/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
GTTGCCGTGGTCTTTTACCGTCGTTTATCACTTGGCGGTGCAACCCATGGATGCCGCTCACCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTGGTGAGGTGAAGGGTGAGCCAGTCGGCATCGA  >  NC_002947/2641656‑2641946

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: