Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,201,833 A→G 8.8% D47G (GAC→GGC)  PP_4583 → putative Peptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,201,8330AG 8.8% 188.6 / 7.1 68D47G (GAC→GGC) PP_4583putative Peptidase
Reads supporting (aligned to +/- strand):  ref base A (29/33);  new base G (2/4);  total (31/37)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTCAATGCCAACCAAGCCCCAACCCCCGATTGCCCAAGCCGACCAGGCCACCGACCCGTATGCCTGGCTGCAACAGCGCGACA‑CCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGG‑GGCCCCTACCTGTACTACACCCGTACCACCGCCGGCGACGAATATCCGCGTCACTACCGCTGCC  >  NC_002947/5201689‑5201979
                                                                                                                                                 |                                                                                                                                                   
cTTCAATGCCAACCAAGCCCCAACCCCCGATTGCCCAAGCCGACCAGGCCACCGACCCGTATGCCTGGCTGCAACAGCGCGACA‑CCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGTCGGCCAg                                                                                                                                                 <  2:21831/149‑1 (MQ=255)
cTTCAATGCCAACCAAGCCCCAACCCCCGATTGCCCAAGCCGACCAGGCCACCGACCCGTATGCCTGGCTGCAACAGCGCGACA‑CCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAg                                                                                                                                                 <  8:67784/149‑1 (MQ=255)
 ttCAATGCCAACCAAGCCCCAACCCCCGATTGCCCAAGCCGACCAGGCCACCGACCCGTATGCCTGGCTGCAACAGCGCGACA‑CCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAgg                                                                                                                                                >  4:241908/1‑149 (MQ=255)
  tCAATGCCAACCAAGCCCCAACCCCCGATTGCCCAAGCCGACCAGGCCACCGACCCGTATGCCTGGCTGCAACAGCGCGACA‑CCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAgg                                                                                                                                                >  3:292296/1‑148 (MQ=255)
    aaTGCCAACCAAGCCCCAACCCCCGATTGCCCAAGCCGACCCGGCCACCGTCCCGTATGCCTCGCTGCAACAGCGCGACA‑CCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGTCGGCCAGGCag                                                                                                                                             <  3:5486/149‑2 (MQ=255)
    aaTGCCAACCAAGCCCCAACCCCCGATTGCCCAAGCCGACCAGGCCACCGACCCGTATGCCTGGCTGCAACAGCGCGACA‑CCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGTCGGCCAGGCag                                                                                                                                             <  6:103788/149‑2 (MQ=255)
            ccaaGCCCCAACCCCCGATTGCCCAAGCCGACCAGGCCACCGACCCGTATGCCTGGCTGCAACAGCGCGACA‑CCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGc                                                                                                                                       <  7:172612/147‑1 (MQ=255)
            ccaaGCCCCAACCCCCGATTGCCCAAGCCGACCAGGCCACCGACCCGTATGCCTGGCTGCAACAGCGCGACA‑CCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGc                                                                                                                                       >  8:172612/1‑147 (MQ=255)
             caaGCCCCAACCCCCGATTGCCCAAGCCGACCAGGCCACCGACCCGTATGCCTGGCTGCAACAGCGCGACA‑CCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTg                                                                                                                                    >  5:77241/1‑149 (MQ=255)
                tccccATCCCCCGATTGCCCAAGCCGACCAGGCCACCGACCCGTATGCCTGGCTGCAACAGCGCGACA‑CCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGACGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAgc                                                                                                                                 <  3:221503/148‑1 (MQ=255)
                gCCACAACCACCTATTGCCCAAGCCGACCAGGCCACCGACCCGTATGACTGGCTGCAACAGCGCGACA‑CCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAgc                                                                                                                                 <  3:188734/149‑1 (MQ=255)
                      aCCCCCGATTGCCCAAGCCGACCATGCCACCTACCCGTATGCCTGGCTGCAACAGCGCGACA‑CCCCCGAGGTACTCGCTTAACTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTTAGCGCCTGt                                                                                                                           <  4:330952/149‑1 (MQ=255)
                      aCCCCCGATTGCCCAAGCCGACCAGGCCACCGACCCGTATGCCTGGCTGCAACAGCGCGACA‑CCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGt                                                                                                                           <  4:292296/149‑1 (MQ=255)
                         cccGATTGCCCAAGCCGACCAGGCCACCGACCCGTATGCCTGGCTGCAACAGCGCGACA‑CCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCg                                                                                                                        <  7:208328/149‑1 (MQ=255)
                           cGATTGCCCAAGCCGACCAGGTCACCGACCCGTATGCCTGGCTGCAACAGCGCGACA‑CCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCgaa                                                                                                                      <  2:214686/149‑1 (MQ=255)
                               tGCCCAAGCCGACCAGGCCACCGACCCGTATGCCTGGCTGCAACAGCGCGACA‑CCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCgaagaa                                                                                                                   >  4:12118/1‑148 (MQ=255)
                                 cccAAGCCGACCAGGCCACCGACCCGTATGCCTGGCTGCAACAGCGCGACA‑CCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAAt                                                                                                                 >  7:40114/1‑148 (MQ=255)
                                  ccAAGCCGACCAGGCCACCGACCCGTATGCCTGGCTGCAACAGCGCGACA‑CCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATc                                                                                                                >  7:225377/1‑148 (MQ=255)
                                   cAAGCCGACCAGGCCACCGACCCGTATGCCTGGCTGCAACAGCGCGACA‑CCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCaa                                                                                                              <  8:286650/149‑1 (MQ=255)
                                     aGCCGACCAGGCCACCGACCCGTATGCCTGGCTGCAACAGCGCGACA‑CCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAgg                                                                                                            <  2:65860/149‑1 (MQ=255)
                                               gCCACCGACCCGTATGCCTGGCTGCAACAGCGCGACACCCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATcc                                                                                                   >  5:178093/1‑149 (MQ=255)
                                                 cACCGACCCGTATGCCTGGCTGCAACAGCGCGACA‑CCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCt                                                                                                  <  7:56864/147‑1 (MQ=255)
                                                 cACCGACCCGTATGCCTGGCTGCAACAGCGCGACA‑CCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCt                                                                                                  >  8:56864/1‑147 (MQ=255)
                                                  aCCGACCCGTATGCCTGGCTGCAACAGCGCGACA‑CCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGga                                                                                               <  3:12118/149‑1 (MQ=255)
                                                    cGACCCGTATGCCTGGCTGCAACAGCGCGACA‑CCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGgaga                                                                                             >  8:180431/1‑149 (MQ=255)
                                                         cGTATGCCTGGCTGCAACAGCGCGACA‑CCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAgaccgac                                                                                        <  8:40114/149‑1 (MQ=255)
                                                              gCCTGGCTGCAACAGCGCGACA‑CCCCCGAGGTACTAGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTAGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGt                                                                                    >  5:338140/1‑148 (MQ=255)
                                                               ccTGGCTGCAACAGCGCGACA‑CCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCa                                                                                  <  8:225377/149‑1 (MQ=255)
                                                               ccTGGCTGCAACAGCGCGACA‑CCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCa                                                                                  >  5:234434/1‑149 (MQ=255)
                                                                cTGGCTGCAACAGCGCGACA‑CCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCAt                                                                                 >  3:149937/1‑149 (MQ=255)
                                                                cTGGCTGCAACAGCGCGACA‑CCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCAGCATCCTGGAGACCGACCTGTCAt                                                                                 >  4:148711/1‑149 (MQ=255)
                                                                       cAACAGCGCGACA‑CCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGgc                                                                          >  5:263200/1‑149 (MQ=255)
                                                                       cAACAGCGCGACA‑CCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGgc                                                                          <  6:234434/149‑1 (MQ=255)
                                                                        aaCAGCGCGACA‑CCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCg                                                                                                                        <  1:108850/102‑1 (MQ=255)
                                                                        aaCAGCGCGACA‑CCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCg                                                                                                                        >  2:108850/1‑102 (MQ=255)
                                                                                 aca‑cCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGG‑GGc                                                               <  6:275461/149‑1 (MQ=255)
                                                                                   a‑cCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGG‑GGccc                                                             <  6:77241/149‑1 (MQ=255)
                                                                                       cccGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGG‑GGCCCCTa                                                          >  3:249561/1‑149 (MQ=255)
                                                                                        ccGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGc                                                                                                          <  2:305892/101‑1 (MQ=255)
                                                                                        ccGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGG‑GGCCCCTAc                                                         <  3:241908/149‑1 (MQ=255)
                                                                                        ccGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGG‑GGCCCCTAc                                                         <  6:178093/149‑1 (MQ=255)
                                                                                             gTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCGGGAGACCGACCTGTCATTGCCGGCGCCCTGG‑GGCCCCTACCTGta                                                    >  7:319467/1‑149 (MQ=255)
                                                                                              tACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGG‑GGCCCCTACCTGtac                                                   <  6:263200/149‑1 (MQ=255)
                                                                                               aCTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGG‑GGCCCCTACCTGtac                                                   >  6:188967/1‑148 (MQ=255)
                                                                                                cTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGG‑GGCCCCTACCTGtact                                                  <  4:149937/148‑1 (MQ=255)
                                                                                                        ccTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGG‑GGCCCCTACCTGTACTACACCCGt                                          <  6:338140/148‑1 (MQ=255)
                                                                                                               gCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGG‑GGCCCCTACCTGTACTACACCCGTACCAccgc                                  >  2:262280/1‑149 (MQ=255)
                                                                                                                 tGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGG‑GGCCCCTCCCTGTACTACACCCGTACCAccgccg                                >  8:189402/1‑149 (MQ=255)
                                                                                                                 tGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGG‑GGCCCCTACCTGtac                                                   <  4:142016/130‑1 (MQ=255)
                                                                                                                 tGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGG‑GGCCCCTACCTGtac                                                   >  3:142016/1‑130 (MQ=255)
                                                                                                                  gAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGG‑GGCCCCTACCTGTACTACACCCGTa                                         >  6:125140/1‑139 (MQ=255)
                                                                                                                  gAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGG‑GGCCCCTACCTGTACTACACCCGTa                                         <  5:125140/139‑1 (MQ=255)
                                                                                                                   aaaaCGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGG‑GGCCCCTACCTGTACTACACCCGTacca                                      <  1:43835/141‑1 (MQ=255)
                                                                                                                   aaaaCGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGG‑GGCCCCTACCTGTACTACACCCGTacca                                      >  2:43835/1‑141 (MQ=255)
                                                                                                                   aaaaCGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGG‑GGCCCCCACCTGTACTACACCCGTACCAccgccg                                <  3:148711/147‑1 (MQ=255)
                                                                                                                      aCGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGG‑GGCCCCTACCTGTACTACACCCGTACCACCGCCGGcgac                           <  4:249561/149‑1 (MQ=255)
                                                                                                                        gCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCg                                                                            <  6:365376/99‑1 (MQ=255)
                                                                                                                        gCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCg                                                                            >  5:365376/1‑99 (MQ=255)
                                                                                                                        gCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGG‑GGCCCCTACCTGTACTACACCCGTACCACCGCCGGcgacg                          >  4:151936/1‑148 (MQ=255)
                                                                                                                           tACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGG‑GGCCCCTACCTGTACTACACCCGTACCACCGCCGGCGACGAata                      >  5:122439/1‑149 (MQ=255)
                                                                                                                              cAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGG‑GGCCCCTACCTGTACTACACCCGTACCACCGCCGGCGACGAATATcc                   >  6:115573/1‑149 (MQ=255)
                                                                                                                               aGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGGCGGCCCCTACCTGTACTACACCCGTACCACCGCCGGCGACGAATATcc                   <  1:262280/149‑1 (MQ=255)
                                                                                                                                ggAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGG‑GGCCCCTACCTGTACTACACCCGTACCACCGCCGGCGACGAATATCcgc                 <  6:122439/149‑1 (MQ=255)
                                                                                                                                gcaacggtgcctgGTCGGCCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGG‑GGCCCCTACCTGTACTACACCCGTACCACCGCCGGCGACGAat                       >  8:101483/8‑143 (MQ=255)
                                                                                                                                gcaacggtgcctgGTCGGCCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGG‑GGCCCCTACCTGTACTACACCCGTACCACCGCCGGCGACGAat                       <  7:101483/136‑1 (MQ=255)
                                                                                                                                    cggtgcctgGTCGGCCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGG‑GGCCCCTACCTGTACTACACCCGTACCACCGCCGGCGACGAATATCCGCGTCa             >  7:340394/4‑149 (MQ=255)
                                                                                                                                     cctgcctgGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGG‑GGCCCCTACCTGTACTACACCCGTACCACCGCCGGCGACGAATATCCGCGTCAc            >  8:152299/1‑149 (MQ=255)
                                                                                                                                               cGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGG‑GGCCCCTACCTGTACTACACCCGTACCACCGCCGGCGACGAATATCCGCGTCACTACCGCTGcc  <  5:188967/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
CTTCAATGCCAACCAAGCCCCAACCCCCGATTGCCCAAGCCGACCAGGCCACCGACCCGTATGCCTGGCTGCAACAGCGCGACA‑CCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGG‑GGCCCCTACCTGTACTACACCCGTACCACCGCCGGCGACGAATATCCGCGTCACTACCGCTGCC  >  NC_002947/5201689‑5201979

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: ‑