Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,856,454 T→C 6.6% G127G (GGT→GGC)  PP_5133 → membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,856,4540TC 6.6% 168.5 / 2.3 60G127G (GGT→GGC) PP_5133membrane protein
Reads supporting (aligned to +/- strand):  ref base T (29/27);  new base C (2/2);  total (31/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.88e-01

TGCTGGCCCTGCTGCTGATGGCCGTGCTCAGCGGCACGGGCCTGGGGCTGGCCGCCCTCGGTGACGGGCAGCGCCCGGTGAATGTGTTCTGGGCGCTGGGCAGCCTGCTCGGGCTGAACCTGCTGATGCTACTTGGCTGGGGCATGGGTTTTGCCCTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTT‑TGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGTGCTGCTGCAACGTCAACGCCTCAACCGGTGGCTGCTC  >  NC_002947/5856306‑5856601
                                                                                                                                                    |                                                                                                                                                    
tgctgGCCCTGCTGCTGATGGCCGTGCTCAGCGGCACGGGCCTGGGGCTGGCCGCCCTCGGTGACGGGCAGCGCCCGGTGAATGTGTTCTGGGCGCTGGGCAGCCTGCTCGGGCTGAACCTGCTGATGCTACTTGGCTGGGGCATGGGt                                                                                                                                                      <  4:21448/149‑1 (MQ=255)
     gCCCTGCTGCTGATGGCCGTGCTCAGCGGCACGGGCCTGGGGCTGGCCGCCCTCGGTGACGGGCAGCGCCCGGTGAATGTGTTCTGGGCGCTGGGCAGCCTGCTCGGGCTGAACCTGCTGATGCTACTTGGCTGGGGCATGGGtttt                                                                                                                                                   <  1:16644/147‑1 (MQ=255)
     gCCCTGCTGCTGATGGCCGTGCTCAGCGGCACGGGCCTGGGGCTGGCCGCCCTCGGTGACGGGCAGCGCCCGGTGAATGTGTTCTGGGCGCTGGGCAGCCTGCTCGGGCTGAACCTGCTGATGCTACTTGGCTGGGGCATGGGTTTTGc                                                                                                                                                 >  6:102543/1‑149 (MQ=255)
     gCCCTGCTGCTGATGGCCGTGCTCAGCGGCACGGGCCTGGGGCTGGCCGCCCTCGGTGACGGGCAGCGCCCGGTGAATGTGTTCTGGGCGCTGGGCAGCCTGCTCGGGCTGAACCTGCTGATGCTACTTGGCTGGGGCATGGGTTTTGc                                                                                                                                                 <  2:14424/149‑1 (MQ=255)
      cccTGCTGCTGATGGCCGTGCTCAGCGGCACGGGCCTGGGGCTGGCCGCCCTCGGTGACGGGCAGCGCCCGGTGAATGTGTTCTGGGCGCTGGGCAGCCTGCTCGGGCTGAACCTGCTGATGCTACTTGGCTGGGGCATGGGTTTTGcc                                                                                                                                                >  4:163810/1‑149 (MQ=255)
       ccTGCTGCTGATGGCCGTGCTCAGCGGCACGGGCCTGGGGCTGGCCGCCCTCGGTGACGGGCAGCGCCCGGTGAATGTGTTCTGGGCGCTGGGCAGCCTGCTCGGGCTGAACCTGCTGATGCTACTTGGCTGGGGCATGGGTTTTGccc                                                                                                                                               >  6:57733/1‑149 (MQ=255)
            tgctgATGGCCGTGCTCAGCGGCACGGGCCTGGGGCTGGCCGCCCTCGGTGACGGGCAGCGCCCGGTGAATGTGTTCTGGGCGCTGGGCAGCCTGCTCGGGCTGAACCTGCTGATGCTACTTGGCTGGGGCATGGGTTTTGCCCTGAgc                                                                                                                                          >  7:7499/1‑149 (MQ=255)
              ctgATGGCCGTGCTCAGCGGCACGGGCCTGGGGCTGGCCGCCCTCGGTGACGGGCAGCGCCCGGTGAATGTGTTCTGGGCGCTGGGCAGCCTGCTCGGGCTGAACCTGCTGATGCTACTTGGCTGGGGCATGGGTTTTGCCCTGAgcgg                                                                                                                                        >  6:350598/1‑149 (MQ=255)
                gATGGCCGTGCTCAGCGGCACGGGCCTGGGGCTGGCCGCCCTCGGTGACGGGCAGCGCCCGGTGAATGTGTTCTGGGCGCTGGGCAGCCTGCTCGGGCTGAACCTGCTGATGCTACTTGGCTGGGGCATGGGTTTTGCCCTGAgcggc                                                                                                                                       >  1:81415/1‑148 (MQ=255)
                  tGGCCGTGCTCAGCGGCACGGGCCTGGGGCTGGCCGCCCTCGGTGACGGGCAGCGCCCGGTGAATGTGTTCTGGGCGCTGGGCAGCCTGCTCGGGCTGAACCTGCTGATGCTACTTGGCTGGGGCATGGGTTTTGCCCTGAGCGGCGAg                                                                                                                                    <  2:3574/149‑1 (MQ=255)
                     ccGTGCTCAGCGGCACGGGCGTGGGGCTGGCCGCCCTCGGTGACGGGCAGCGCCCGGTGAATGTGTTCTGGGCGCTGGGCAGCCTGCTCGGGCTGAACCTGCTGATGCTACTTGGCTGGGGCATGGGTTTTGCCCTGAGCGGCGa                                                                                                                                     <  4:271467/145‑1 (MQ=255)
                     ccGTGCTCAGCGGCACGGGCCTGGGGCTGGCCGCCCTCGGTGACGGGCAGCGCCCGGTGAATGTGTTCTGGGCGCTGGGCAGCCTGCTCGGGCTGCACCTGCTGATGCTACTTGGCTGGGGCATGGGTTTTGCCCTGAGCGGCGa                                                                                                                                     >  3:271467/1‑145 (MQ=255)
                     ccGTGCTCAGCGGCACGGGCCTGGGGCTGGCCGCCCTCGGTGACGGGCAGCGCCCGGTGAATGTGTTCTGGGCGCTGGGCAGCCTGCTCGGGCTGAACCTGCTGATGCTACTTGGCTGGGGCATGGGTTTTGCCCTGAGCGGCGAGCAt                                                                                                                                 >  1:1658/1‑149 (MQ=255)
                      cGTGCTCAGCGGCACGGGCCTGGGGCTGGCCGCCCTCGGTGACGGGCAGCGCCCGGTGAATGTGTTCTGGGCGCTGGGCAGCCTGCTCGGGCTGAACCTGCTGATGCTACTTGGCTGGGGCATGGGTTTTGCCCTGAGCGGCGAGCATg                                                                                                                                <  8:7499/149‑1 (MQ=255)
                      cGTGCTCAGCGGCACGGGCCTGGGGCTGGCCGCCCTCGGTGACGGGCAGCGCCCGGTGAATGTGTTCTGGGCGCTGGGCAGCCTGCTCGGGCTGAACCTGCTGATGCTACTTGGCTGGGGCATGGGCTTTGCCCTGAGCGGCGAGCATg                                                                                                                                >  8:272076/1‑149 (MQ=255)
                      cGTGCTCAGCGGCACGGGCCTGGGGCTGGCCGCCCTCGGTGACGGGCAGCGCCCGGTGAATGTGTTCTGGGCGCTGGGCAGCCTGCTCGGGCTGAACCTGCTGATGCTACTTGGCTGGGGCATGGGCTTTGCCCTGAGCGGCGAGCATg                                                                                                                                >  8:276905/1‑149 (MQ=255)
                          cTCAGCGGCACGGGCCTGGGGCTGGCCGCCCTCGGTGACGGGCAGCGCCCGGTGAATGTGTTCTGGGCGCTGGGCAGCCTGCTCGGGCTGAACCTGCTGATGCTACTTGGCTGGGGCATGGGTTTTGCCCTGAGCGGCGAGCATGgcgc                                                                                                                            >  6:195078/1‑149 (MQ=255)
                             aGCGGCACGGGCCTGGGGCTGGCCGCCCTCGGTGACGGGCAGCGCCCGGTGAATGTGTTCTGGGCGCTGGGCAGCCTGCTCGGGCTGAACCTGCTGATGCTACTTGGCTGGGGCATGGGTTTTGCCCTGAGCGGCGAGCATGGCGCAgg                                                                                                                         >  8:64147/1‑149 (MQ=255)
                                       gCCTGGGGCTGGCCGCCCTCGGTGACGGGCAGCGCCCGGTGAATGTGTTCTGGGCGCTGGGCAGCCTGCTCGGGCTGAACCTGCTGATGCTACTTGGCTGGGGCATGGGTTTTGCCCTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGc                                                                                                               >  8:214075/1‑149 (MQ=255)
                                         cTGGGGCTGGCCGCCCTCGGTGACGGGCAGCGCCCGGTGAATGTGTTCTGGGCGCTGGGCAGCCTGCTCGGGCTGAACCTGCTGATGCTACTTGGCTGGGGCATGGGTTTTGCCCTGAGCGGCGAGCATGGCGCAGGGc                                                                                                                       >  3:176030/1‑139 (MQ=255)
                                         cTGGGGCTGGCCGCCCTCGGTGACGGGCAGCGCCCGGTGAATGTGTTCTGGGCGCTGGGCAGCCTGCTCGGGCTGAACCTGCTGATGCTACTTGGCTGGGGCATGGGTTTTGCCCTGAGCGGCGAGCATGGCGCAGGGc                                                                                                                       <  4:176030/139‑1 (MQ=255)
                                          tGGGGCTGGCCGCCCTCGGTGACGGGCAGCGCCCGGTGAATGTGTTCTGGGCGCTGGGCAGCCTGCTCGGGCTGAACCTGCTGATGCTACTTGGCTGGGGCATGGGTTTTGCCCTGAGCGGCGAGCATGGCGCAg                                                                                                                          <  8:192569/135‑1 (MQ=255)
                                          tGGGGCTGGCCGCCCTCGGTGACGGGCAGCGCCCGGTGAATGTGTTCTGGGCGCTGGGCAGCCTGCTCGGGCTGAACCTGCTGATGCTACTTGGCTGGGGCATGGGTTTTGCCCTGAGCGGCGAGCATGGCGCAg                                                                                                                          >  7:192569/1‑135 (MQ=255)
                                             ggCTGGCCGCCCTCGGTGACGGGCAGCGCCCGGTGAATGTGTTCTGGGCGCTGGGCAGCCTGCTCGGGCTGAACCTGCTGATGCTACTTGGCTGGGGCATGGGTTTTGCCCTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTtgtgg                                                                                                         <  8:316532/149‑1 (MQ=255)
                                               cTGGCCGCCCTCGGTGACGGGCAGCGCCCGGTGAATGTGTTCTGGGCGCTGGGCAGCCTGCTCGGGCTGAACCTGCTGATGCTACTTGGCTGGGGCATGGGTTTTGCCCTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTtgtggct                                                                                                       <  5:195078/149‑1 (MQ=255)
                                                  gccgccCTCGGTGACGGGCAGCGCCCGGTGAATGTGTTCTGGGCGCTGGGCAGCCTGCTCGGGCTGAACCTGCTGATGCTACTTGGCTGGGGCATGGGTTTTGCCCTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTtgtggctgtg                                                                                                    <  5:57733/149‑1 (MQ=255)
                                                  gccgccCTCGGTGACGGGCAGCGCCCGGTGAATGTGTTCTGGGCGCTGGGCAGCCTGCTCGGGCTGAACCTGCTGATGCTACTTGGCTGGGGCATGGGCTTTGCCCTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTtgtggctg                                                                                                      <  7:272076/147‑1 (MQ=255)
                                                  gccgccCTCGGTGACGGGCAGCGCCCGGTGAATGTGTTCTGGGCGCTGGGCAGCCTGCTCGGGCTGAACCTGCTGATGCTACTTGGCTGGGGCATGGGCTTTGCCCTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTtgtggctg                                                                                                      <  7:276905/147‑1 (MQ=255)
                                                   ccgccCTCGGTGACGGGCAGCGCCCGGTGAATGTGTTCTGGGCGCTGGGCAGCCTGCTCGGGCTGAACCTGCTGATGCTACTTGGCTGGGGCATGGGTTTTGCCCTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTtgtggctgtgg                                                                                                   <  7:64147/149‑1 (MQ=255)
                                                     gccCTCGGTGACGGGCAGCGCCCGGTGAATGTGTTCTGGGCGCTGGGCAGCCTGCTCGGGCTGAACCTGCTGATGCTACTTGGCTGGGGCATGGGTTTTGCCCTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTtgtggctgtggc                                                                                                  >  1:354981/1‑148 (MQ=255)
                                                     gccCTCGGTGACGGGCAGCGCCCGGTGAATGTGTTCTGGGCGCTGGGCAGCCTGCTCGGGCTGAACCTGCTGATGCTACTTGGCTGGGGCATGGGTTTTGCCCTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTgtggctgtggct                                                                                                 <  2:354981/149‑1 (MQ=255)
                                                        cTCGGTGACGGGCAGCGCCCGGTGAATGTGTTCTGGGCGCTGGGCAGCCTGCTCGGGCTGAACCTGCTGATGCTACTTGGCTGGGGCATGGGTTTTGCCCTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCAAAg                                                                                              >  8:90359/1‑149 (MQ=255)
                                                            gTGACGGGCAGCGCCCGGTGAATGTGTTCTGGGCGCTGGGCAGCCTGCTCGGGCTGAACCTGCTGATGCTACTTGGCTGGGGCATGGGTTTTGCCCTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGaa                                                                                          >  8:88772/1‑149 (MQ=255)
                                                            gTGACGGGCAGCGCCCGGTGAATGTGTTCTGGGCGCTGGGCAGCCTGCTCGGGCTGAACCTGCTGATGCTACTTGGCTGGGGCATGGGTTTTGCCCTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGaa                                                                                          <  1:206881/149‑1 (MQ=255)
                                                                cGGGCAGCGCCCGGTGAATGTGTTCTGGGCGCTGGGCAGCCTGCTCGGGCTGAACCTGCTGATGCTACTTGGCTGGGGCATGGGTTTTGCCCTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCt                                                                                      <  7:88772/149‑1 (MQ=255)
                                                                 gggCAGCGCCCGGTGAATGTGTTCTGGGCGCTGGGCAGCCTGCTCGGGCTGAACCTGCTGATGCTACTTGGCTGGGGCATGGGTTTTGCCCTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCtt                                                                                     >  6:71107/1‑149 (MQ=255)
                                                                      gcgcCCGTGGAATGTGTTCTGGGCGCTGGGCAGCCGGCTCGGGCTGAACCTGCTGATGCTACTTGGCTGGGGCATGGGTTTTGCCCTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTT‑TGccc                                                                               <  2:81415/149‑1 (MQ=255)
                                                                         cccGGTGAATGTGTTCTGGGCGCTGGGCAGCCTGCTCGGGCTGAACCTGCTGATGCTACTTGGCTGGGGCATGGGTTTTGCCCTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCAAAGCGTTCGCTT‑AGccacag                                                                            <  2:303039/149‑5 (MQ=255)
                                                                             gTGAATGTGTTCTGGGCGCTGGGCAGCCTGCTCGGGCTGAACCTGCTGATGCTACTTGGCTGGGGCATGGGTTTTGCCCTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTT‑TGCCCGCGATg                                                                         >  8:115683/1‑148 (MQ=255)
                                                                                  tgtgTTCTGGGCGCTGGGCAGCCTGCTCGGGCTGAACCTGCTGATGCTACTTGGCTGGGGCATGGGTTTTGCCCTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTT‑TGCCCGCGATGCCAAgg                                                                   <  5:71107/149‑1 (MQ=255)
                                                                                     gtTCTGGGCGCTGGGCAGCCTGCTCGGGCTGAACCTGCTGATGCTACTTGGCTGGGGCATGGGTTTTGCCCTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTT‑TGCCCGCGATGCCAAGGCcg                                                                >  2:96285/1‑149 (MQ=255)
                                                                                          gggCGCTGGGCAGCCTGCTCGGGCTGAACCTGCTGATGCTACTTGGCTGGGGCATGGGTTTTGCCCTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTT‑TGCCCGCGATGCCAAGGCCGCGCAc                                                           >  3:270730/1‑149 (MQ=255)
                                                                                           ggCGCTGGGCAGCCTGCTCGGGCTGAACCTGCTGATGCTACTTGGCTGGGGCATGGGTTTTGCCCTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTT‑TGCCCGCGATGCCAAGGCCGCGCAcc                                                          >  4:336580/1‑149 (MQ=255)
                                                                                              gcTGGGCAGCCTGCTCGGGCTGAACCTGCTGATGCTACTTGGCTGGGGCATGGGTTTTGCCCTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTT‑TGCCCGCGATGCCAAGGCCGCGCACCTgg                                                       >  4:23568/1‑149 (MQ=255)
                                                                                               cTGGGCAGCCTGCTCGGGCTGAACCTGCTGATGCTACTTGGCTGGGGCATGGGTTTTGCCCTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTT‑TGCCCGCGATGCCAAGGCCGCGCACCTGgc                                                      <  7:115683/149‑1 (MQ=255)
                                                                                                  ggCAGGCTGCTCGGGCTGAACCTGCTGATGCTACTTGGCTGGGGCATGGGTTTTGCCCTGAGCGGCGAGCATGGCGCAGGGCTg                                                                                                                     >  7:326148/1‑84 (MQ=255)
                                                                                                  ggCAGCCTGCTCGGGCTGAACCTGCTGATGCTACTTGGCTGGGGCATGGGTTTTGCCCTGAGCGGCGAGCATGGCGCAGGGCTg                                                                                                                     <  8:326148/84‑1 (MQ=255)
                                                                                                       ccTGCTCGGGCTGAACCTGCTGATGCTACTTGGCTGGGGCATGGGTTTTGCCCTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTT‑TGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAgcgc                                              >  3:58503/1‑149 (MQ=255)
                                                                                                       ccTGCTCGGGCTGAACCTGCTGATGCTACTTGGCTGGGGCATGGGTTTTGCCCTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTT‑TGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAgcgc                                              >  8:279673/1‑149 (MQ=255)
                                                                                                           cTCGGGCTGAACCTGCTGATGCTACTTGGCTGGGGCATGGGTTTTGCCCTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTT‑TGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCgctgct                                          <  5:102543/149‑1 (MQ=255)
                                                                                                              gggCTGAACCTGCTGATGCTACTTGGCTGGGGCATGGGTTTTGCCCTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTT‑TGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGctgctg                                         <  8:46632/147‑1 (MQ=255)
                                                                                                              gggCTGAACCTGCTGATGCTACTTGGCTGGGGCATGGGTTTTGCCCTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTT‑TGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTgg                                        <  1:349811/148‑1 (MQ=255)
                                                                                                                   gAACCTGCTGATGCTACTTGGCTGGGGCATGGGTTTTGCCCTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTT‑TGCCCGCGATACCAAGGCCGCGCACCTGgcgc                                                    <  1:138864/131‑1 (MQ=255)
                                                                                                                    aaCCTGCTGATGCTACTTGGCTGGGGCATGGGTTTTGCCCTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTT‑TGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGtgctgc                                  <  1:96285/148‑1 (MQ=255)
                                                                                                                                cTACTTGGCTGGGGCATGGGTTTTGCCCTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTTATGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGTGCTGCTGCAACGTCAAc                      <  1:287492/149‑1 (MQ=255)
                                                                                                                                      ggCTGGGGCATGGGTTTTGCCCTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTT‑TGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGTGCTGCTGCAACGTCAACGCCTCaa               >  1:146364/1‑149 (MQ=255)
                                                                                                                                           gggCATGGGTTTTGCCCTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTT‑TGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGTGCTGCTGCAACGTCAACGCCTCAACCgg           <  7:214075/148‑1 (MQ=255)
                                                                                                                                              cATGGGTTTTGCCCTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTT‑TGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGTGCTGCTGCAACGTCAACGCCTCAACCGGTGgc       <  4:270730/149‑1 (MQ=255)
                                                                                                                                               aTGGGTTTTGCCCTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTT‑TGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGTGCTGCTGCAACGTCAACGCCTCAACCGGTGgct      >  7:117517/1‑149 (MQ=255)
                                                                                                                                                tGGGTTTTGCCCTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTT‑TGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGTGCTGCTGCAACGTCAACGCCTCAACCGGTGgctg     >  5:185053/1‑149 (MQ=255)
                                                                                                                                                 gggTTTTGCCCTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTT‑TGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGTGCTGCTGCAACGTCAACGCCTCAACCGGTGgctgc    >  1:34552/1‑149 (MQ=255)
                                                                                                                                                   gTTTTGCCCTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTT‑TGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGTGCTGCTGCAACGTCAACGCCTCAACCGGTGGCTGCTc  <  3:163810/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TGCTGGCCCTGCTGCTGATGGCCGTGCTCAGCGGCACGGGCCTGGGGCTGGCCGCCCTCGGTGACGGGCAGCGCCCGGTGAATGTGTTCTGGGCGCTGGGCAGCCTGCTCGGGCTGAACCTGCTGATGCTACTTGGCTGGGGCATGGGTTTTGCCCTGAGCGGCGAGCATGGCGCAGGGCTGGGCCGCTTGTGGCTGTGGCTAAGCGAACGCTT‑TGCCCGCGATGCCAAGGCCGCGCACCTGGCGCCAGCGCTGCTGGTGCTGCTGCAACGTCAACGCCTCAACCGGTGGCTGCTC  >  NC_002947/5856306‑5856601

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