Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,201,830 C→T 8.9% A46V (GCC→GTC)  PP_4583 → putative Peptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,201,8300CT 8.9% 188.7 / 6.8 67A46V (GCC→GTC) PP_4583putative Peptidase
Reads supporting (aligned to +/- strand):  ref base C (29/32);  new base T (2/4);  total (31/36)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01

CTTCAATGCCAACCAAGCCCCAACCCCCGATTGCCCAAGCCGACCAGGCCACCGACCCGTATGCCTGGCTGCAACAGCGCGACA‑CCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGG‑GGCCCCTACCTGTACTACACCCGTACCACCGCCGGCGACGAATATCCGCGTCAC  >  NC_002947/5201689‑5201969
                                                                                                                                              |                                                                                                                                            
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                                         >  3:142016/1‑130 (MQ=255)
                                                                                                                 tGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGG‑GGCCCCTACCTGtac                                         <  4:142016/130‑1 (MQ=255)
                                                                                                                  gAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGG‑GGCCCCTACCTGTACTACACCCGTa                               <  5:125140/139‑1 (MQ=255)
                                                                                                                  gAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGG‑GGCCCCTACCTGTACTACACCCGTa                               >  6:125140/1‑139 (MQ=255)
                                                                                                                   aaaaCGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGG‑GGCCCCTACCTGTACTACACCCGTacca                            >  2:43835/1‑141 (MQ=255)
                                                                                                                   aaaaCGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGG‑GGCCCCTACCTGTACTACACCCGTacca                            <  1:43835/141‑1 (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             <  7:101483/136‑1 (MQ=255)
                                                                                                                                gcaacggtgcctgGTCGGCCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGG‑GGCCCCTACCTGTACTACACCCGTACCACCGCCGGCGACGAat             >  8:101483/8‑143 (MQ=255)
                                                                                                                                    cggtgcctgGTCGGCCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGG‑GGCCCCTACCTGTACTACACCCGTACCACCGCCGGCGACGAATATCCGCGTCa   >  7:340394/4‑149 (MQ=255)
                                                                                                                                     cctgcctgGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGG‑GGCCCCTACCTGTACTACACCCGTACCACCGCCGGCGACGAATATCCGCGTCAc  >  8:152299/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                            
CTTCAATGCCAACCAAGCCCCAACCCCCGATTGCCCAAGCCGACCAGGCCACCGACCCGTATGCCTGGCTGCAACAGCGCGACA‑CCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGG‑GGCCCCTACCTGTACTACACCCGTACCACCGCCGGCGACGAATATCCGCGTCAC  >  NC_002947/5201689‑5201969

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