Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,267,765:1 +C 100% coding (93/1254 nt) galP‑IV ← porin‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,267,7651.C80.0% 147.2 / 37.8 60coding (93/1254 nt)galP‑IVporin‑like protein
Reads supporting (aligned to +/- strand):  ref base . (8/4);  new base C (23/23);  total (32/28)
Fisher's exact test for biased strand distribution p-value = 3.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.

CGGGGTGTAGCCAGAGCTGAACTTGAAGATGAAACCCTGGGCCCACTCTTCAACCTTGCCCTTGGCGGCCCCCGGGTCGCGGAAATCGCGGTTGAAGTAGTAGTTGCGCAATGTCAGCCCACC‑‑CTTGGCATCTTCCAAGAAG‑‑CCCCCCTCTTGCGCCCAGGCCGGCAAACCAAGGCCAAACGCCAGCGAAACAGCCAGGACCGAACGGGCCAGCGTGCAAGTTGTGTCCATCATTATTATTGTTCTCCACGTTCAATAGGCGCAGCAGACCGCCCCAGGCCAGC  >  NC_002947/4267624‑4267907
                                                                                                                                                  |                                                                                                                                                  
cGGGGTGTAGCCAGAGCTGAACTTGAAGATGAAACCCTGGGCCCACTCTTCAACCTTGCCCTTGGCGGCCCCCGGGTCGCGGAAATCGCGGTTGAAGTAGTAGTTGCGCAATGTCAGCCCACC‑‑CTTGGCTTCTTGGAAGATG‑‑ccaaggg                                                                                                                                            <  7:18327/149‑6 (MQ=255)
cGGGGTGTAGCCAGAGCTGAACTTGAAGATGAAACCCTGGGCCCACTCTTCAACCTTGCCCTTGGCGGCCCCCGGGTCGCGGAAATCGCGGTTGAAGTAGTAGTTGCGCAATGTCAGCCCACC‑‑CTTGGCATCTTCCAAGAAG‑‑ccccccc                                                                                                                                           <  5:322670/149‑2 (MQ=255)
       tAGCCAGAGCTGAACTTGAAGATGAAACCCTGGGCCCACTCTTCAACCTTGCCCTTGGCGGCCCCCGGGTCGCGGAAATCGCGGTTGAAGTAGTAGTTGCGCAATGTCAGCCCACC‑‑CTTGGCATCTTCCAAGAAGCCCCCCCTCTTgcg                                                                                                                                   >  5:286893/1‑149 (MQ=255)
             gagCTGAACTTGAAGATGAAACCCTGGGCCCACTCTTCAACCTTGCCGTTGGCGGCCCCCGGGTCGCGGAAATCGCGGTTGAAGTAGTAGTTGCGCAATGTCAGCCCACC‑‑CTTGGCATCTTCCAAGAAG‑‑CCCCCCTCTTGCTCCCAGGc                                                                                                                            <  8:127561/149‑1 (MQ=255)
               gCTGAACTTGAAGATGAAACCCTGGGCCCACTCTTCAACCTTGCCCTTGGCGGCCCCCGGGTCGCGGAAATCGCGGTTGAAGTAGTAGTTGCGCAATGTCAGCCCACC‑‑CTTGGCATCTTCCAAGAAGCCCCCCCTCTTGCGCCCAGGcc                                                                                                                           >  7:136325/1‑149 (MQ=255)
                cTGAACTTGAAGATGAAACCCTGGGCCCACTCTTCAACCTTGCCCTTGGCGGCCCCCGGGTCGCGGAAATCGCGGTTGAAGTAGTAGTTGCGCAATGTCAGCCCACC‑‑CTTGGCATCTTCCAAGAAGG‑‑CCCCCCTCGTGCGCCCCGGCCg                                                                                                                          >  8:213236/1‑149 (MQ=255)
                cTGAACTTGAAGATGAAACCCTGGGCCCACTCTTCAACCTTGCCCTTGGCGGCCACCGGGTCGCGGAAATCGCGGTTGAAGTAGTAGTTGGGCAATGTCAGCGCACC‑‑CTTGGCCTCGTCCAAGAAG‑‑CCCCCCGTCGTGCGCCCGGGCCg                                                                                                                          >  8:207797/1‑149 (MQ=255)
                     cTTGAAGATGAAACCCTGGGCCCACTCTTCAACCTTGCCCTTGGCGGCCCCCGGGTCGCGGAAATCGCGGTTGAAGTAGTAGTTGCGCAATGTCAGCCCACC‑‑CTTGGCATCTTCCAAGAAGCCCCCCCTCTTGCGCCCAGGCCGGCaaa                                                                                                                     <  2:98786/149‑1 (MQ=255)
                      ttGAAGATGAAACCCTGGGCCCACTCTTCAACCTTGCCCTTGGCGGCCCCCGGGTCGCGGAAATCGCGGTTGAAGTAGTAGTTGCGCAATGTCAGCCCACC‑‑CTTGGCATCTTCCAAGAAGCCCCCCCTCTTGCGCCCAGGCCGGCAAAc                                                                                                                    >  4:84120/1‑149 (MQ=255)
                        gAAGATGAAACCCTGGGCCCACTCTTCAACCTTGCCCTTGGCGGCCCCCGGTTCGCGGAAATCGCGGTTGAAGTAGTAGTTGCGCAATGTCAGCCCAAC‑‑CTTGGCATCTTCCAAGAAGCCCCCCCTCTTCCGCCCAGGCAGCCCAAcc                                                                                                                   >  4:141191/1‑148 (MQ=255)
                          aGATGAAACCCTGGGCCCACTCTTCAACCTTGCCCTTGGCGGCCCCCGGGTCGCGGAAATCGCGGTTGAAGTAGTAGTTGCGCAATGTCAGCCCACC‑‑CTTGGCATCTTCCAAGAAGCCCCCCCTCTTGCGCCCAGGCCGGCAAACCAAg                                                                                                                >  8:166975/1‑149 (MQ=255)
                          aGATGAAACCCTGGGCCCACTCTTCAACCTTGCCCTTGGCGGCCCCCGGGTCGCGGAAATCGCGGTTGAAGTAGTAGTTGCGCAATGTCAGCCCACC‑‑CTTGGCATCTTCCAAGAAGCCCCCCCTCTTGCGCCCAGGCCGGCAAACCAAg                                                                                                                <  7:213236/149‑1 (MQ=255)
                          aGATGAAACCCTGGGCCCACTCTTCAACCTTGCCCTTGGCGGCCCCCGGGTCGCGGAAATCGCGGTTGAAGTAGTAGTTGCGCAATGTCAGCCCACC‑‑CTTGGCATCTTCCAAGAAGCCCCCCCTCTTGCGCCCAGGCCGGCAAACCAAg                                                                                                                <  7:207797/149‑1 (MQ=255)
                          aGATGAAACCCTGGGCCCACTCTTCAACCTTGCCCTTGGCGGCCCCCGGGTCGCGGAAATCGCGGTTGAAGTAGTAGTTGCGCAATGTCAGCCCACC‑‑CTTGGCATCTTCCAAGAAGCCCCCCCTCTTGCGCCCAGGCCGGCAAACCAAg                                                                                                                <  6:286893/149‑1 (MQ=255)
                              gAAACCCTGGGCCCACTCTTCAACCTTGCCCTTGGCGGCCCCCGGGTCGCGGAAATCGCGGTTGAACTAGTAGTTTCGCAATTTCATCCCACC‑‑CTTGGCATCTTCCAAGAAT‑‑CCCCCCTATTGCGCCCAGGCCGGCAAACCAAGGACa                                                                                                            >  2:167282/1‑148 (MQ=255)
                                aaCCCTGGGCCCACTCTTCAACCTTGCCCTTGGCGGCCCCCGGGTCGCGGAAATCGCGGTTGAAGTAGTAGTTGCGCAATGTCAGCCCACC‑‑CTTGGCATCTTCCAAGAAG‑‑CCACCCCTCTTGCGCCCCGGCCGGCAAAACAACGCCaa                                                                                                           >  4:288195/1‑148 (MQ=255)
                                         cccACTCTTCAACCTTGCCCTTGGCGGCCCCCGGGTCGCGGAAATCGCGGTTGAAGTAGTAGTTGCGCAATGTCAGCCCACC‑‑CTTGGCATCTTCCAAGAAGCCCCCCCTCTTGCGCCCAGGCCGGCAAACCAAGGCCAAACGCCAGCg                                                                                                  >  7:218863/1‑148 (MQ=255)
                                          ccACTCTTCAACCTTGCCCTTGGCGGCCCCCGGGTCGCGGAAATCGCGGTTGAAGTAGTAGTTGCGCAATGTCAGCCCACC‑‑CTTGGCATCTTCCAAGAAG‑‑CCCCCCTCTTGCGCCCAGGCCGGCAAACCAAGGCCAAACGCCAGCGaaa                                                                                               >  4:36485/1‑149 (MQ=255)
                                           cACTCTTCAACCTTGCCCTTGGCGGCCCCCGGGTCGCGGAAATCGCTGTTGAAGTAGTAGTTGCGCAATGTCAGCCCACC‑‑CTTGGCATCTTCCAAGAAGCCCCCCCTCTTGCGCCCAGGCCGGCAAACCAAGGCCAAACGCCAGCGaaa                                                                                               <  3:288195/149‑1 (MQ=255)
                                                  aaaCATTGCCCTGGCGGGCCCCCGGGGGCGGCAAACGCGGGTGAACTATTATTTGCGCAATTTTACCCCCCCCTCTTCGGCATCTTCCAAGAAG‑‑CCCCCCTCTTGCGCCCAGGCCGGCAAACCAAGGCCAAACGCCAGCGAAACAGcc                                                                                          <  1:293512/147‑1 (MQ=255)
                                                     ccTTGCCCTTGGCGGCCCCCGGGTCGCGGAAATCGCGGTTGAAGTAGTAGTTTCGCAATGTCAGCCCACC‑‑CTTGGCATCTTCCAAGAAGCCCCCCCTCTTGCGCCCAGGCCGGCAAACCAAGGCCAAACGCCAGCGAAACAGCCAGGAc                                                                                     <  5:206858/149‑1 (MQ=255)
                                                     ccTTGCCCTTGGCGGCCCCCGGGTCGCGGAAATCGCGGTTGAAGTAGTAGTTGCGCAATGTCAGCCCACC‑‑CTTGGCATCTTCCAAGAAGCCCCCCCTCTTGCGCCCAGGCCGGAAAACCAAGGCCAAACGCCAGACAAACAGCCAGCAc                                                                                     >  6:206858/1‑149 (MQ=255)
                                                       ttGCCCTTGGCGGCCCCCGGGTCGCGGAAATCGCGGTTGAAGTAGTAGTTGCGCAATGTCAGCCCACC‑‑CTTGGCATCTTCCAAGAAGCCCCCCCTCTTGCGCCCAGGCCGGCAAACCAAGGCCAAACGCCAGCGAAACAGCCAGGACCg                                                                                   >  4:61022/1‑149 (MQ=255)
                                                          cccTTGGCGGCCCCCGGGTCGCGGAAATCGCGGTTGAAGTAGTAGTTGCGCAATGTCAGCCCACC‑‑CTTGGCATCTTCCAAGAACC‑‑CCCCCCTATTGCGCCCAGGCCGCCCAACCAAGGCCAACCGCCAGCGAAACAGCCAGGACCGa                                                                                  >  6:102898/1‑147 (MQ=255)
                                                           ccTTGGCGGCCCCCGGTTCGGGGAAATCGCTGTTGAAGTAGTAGTTGCGCAATGTCAGCCCACC‑‑CTTGGCATCTTCCAAGAAG‑‑CCCCCCTCTTGCGCCCAGGCCGGCAAACCAAGGCCAAACGCCAGCGAAACAGCCAGGACCGAACgg                                                                              <  1:167282/149‑1 (MQ=255)
                                                           ccTTGGCGGCCCCCGGGTCGCGGAAATCGCGGTTGAAGTAGTAGTTGCGCAATGTCAGCCCACC‑‑CTTGGCATCTTCCAAGAAGTCCCCCCTCTTGCGCCCAGGCCGGCAAACCAAGGCCAAACGCCAGCGAAACAGCCAGGACCGAACg                                                                               >  3:334308/1‑149 (MQ=255)
                                                           ccTTGGCGGCCCCCGGGTCGCGGAAATCGCGGTTGAAGTAGTAGTTGCGCAATGTCAGCCCACC‑‑CTTGGCATCTTCCAAGAAGCCCCCCCTCTTGCGCCCAGGCCGGCAAACCAAGGCCAAACGCCAGCGAAACAGCCAGGACCGAACg                                                                               >  8:205840/1‑149 (MQ=255)
                                                           ccTTGGCGGCCCCCGGGTCGCGGAAATCGCGGTTGAAGTAGTAGTTGCGCAATGTCAGCCCACC‑‑CTTGGCATCTTCCAAGAAGCCCCCCCTCTTGCGCCCAGGCCGGCAAACCAAGGCCAAACGCCAGCGAAACAGCCAGGACCGAACg                                                                               >  1:272527/1‑149 (MQ=255)
                                                           ccTTGGCGGCCCCCGGGTCGCGGAAATCGCGGTTGAAGTAGTAGTTGCGCAATGTCAGCCCACC‑‑CTTGGCATCGTCCAAGAAG‑‑CCCCCCTCTTGCGCCCAGGCCGGCAAACCAAGGCCAACCGCCAGCGAAAGAGCCAGGACCGAACgg                                                                              >  6:365852/1‑149 (MQ=255)
                                                                    cccccGGGTCGCGGAAATCGCGGTTGAAGTAGTAGTTGCGCAATGTCAGCCCACC‑‑CTTGGCATCTTCCAAGAAGCCCCCCCTCTTGCGCCCAGGCCGGCAAACCAAGGCCAAACGCCAGCGAAACAGCCAGGACCGAACGGGCCAGCGt                                                                      <  3:326082/149‑1 (MQ=255)
                                                                      cccGGGTCGCGGAAATCGCGGTTGAAGTAGTAGTTGCGCAATGTCAGCCCACC‑‑CTTGGCATCTTCCAAGAAGCCCCCCCTCTTGCGCCCAGGCCGGCAAACCAAGGCCAAACGCCAGCGAAACAGCCAGGACCGAACGGGCCAGCGTg                                                                     >  4:304131/1‑148 (MQ=255)
                                                                        cGGGTCGCGGAAATCGCGGTTGAAGTAGTAGTTGCGCAATGTCAGCCCACC‑‑CTTGGCATCTTCCAAGAAGCCCCCCCTCTTGCGCCCAGGCCGGCAAACCAAGGCCAAACGCCAGCGAAACAGCCAGGACCGAACGGGCCAGCGTGc                                                                    >  1:163888/1‑147 (MQ=255)
                                                                              gcgGAAATCGCGGTTGAAGTAGTAGTTGCGCAATGTCAGCCCACC‑‑CTTGGCATCTTCCAAGAAGCCCCCCCTCTTGCGCCCAGGCCGGCAAACCAAGGCCAAACGCCCGCGAAACAGCCAGGACCGAACGGGCCAGCGTGCAAGTtgtg                                                            >  2:227911/1‑149 (MQ=255)
                                                                                       gcgGTTGAAGTAGTAGTTGCGCAATGTCAGCCCACC‑‑CTTGGCATCTTCCAAGAAGCCCCCCCTCTTGCGCCCAGGCCGGCAAACCAAGGCCAAACGCCAGCGAAACAGCCAGGACCGAACGGGCCAGCGTGCAAGTTGTGTCCATCatt                                                   <  3:160371/149‑1 (MQ=255)
                                                                                          gTTGAAGTAGTAGTTGCGCAATGTCAGCCCACC‑‑CTTGGCATCTTCCAAGAAGCCCCCCCTCTTGCGCCCAGGCCGGCAAACCAAGGCCAAACGCCAGCGAAACAGCCAGGACCGAACGGGCCAGCGTGCAAGTTGTGTCCATCattatt                                                >  7:139578/1‑149 (MQ=255)
                                                                                             gAAGTAGTAGTTGCGCAATGTCAGCCCACC‑‑CTTGGCATCTTCCAAGAAGCCCCCCCTCTTGCGCCCAGGCCGGCAAACCAAGGCCAAACGCCAGCGAAACAGCCa                                                                                         >  6:286269/1‑105 (MQ=255)
                                                                                             gAAGTAGTAGTTGCGCAATGTCAGCCCACC‑‑CTTGGCATCTTCCAAGAAGCCCCCCCTCTTGCGCCCAGGCCGGCAAACCAAGGCCAAACGCCAGCGAAACAGCCa                                                                                         <  5:286269/105‑1 (MQ=255)
                                                                                              aaGTAGTAGTTGCGCAATGTCAGCCCACC‑‑CTTGGCATCTTCCAAGAAGCCCCCCCTCTTGCGCCCAGGCCGGCAAACCAAGGCCAAACGCCAGCGAAACAGCCAGGACCGAACGGGCCAGCGTGCAAGTTGTGTCCATCattattatt                                             <  3:61022/148‑1 (MQ=255)
                                                                                                 tagtagtTGCGCAATGTCAGCCCACC‑‑CTTGGCATCTTCCAAGAAGCCCCCCCCTCTTGCGCCCAGGCCGGCAAACCAAGGCCAAACGCCAGCGAAACAGCCAGGACCGAACGGGCCAGCGTGCAAGTTGTGTCCATCattattatt                                             <  7:1873/146‑1 (MQ=255)
                                                                                                 tagtagtTGCGCAATGTCAGCCCACC‑‑CTTGGCATCTTCCAAGAAGCCCCCCCCTCTTGCGCCCAGGCCGGCAAACCAAGGCCAAACCCCACCGAAACACCCAGGACCGAACGGGCCAGCGTGCAAGATGAGTGCAAAattattattta                                           >  8:1873/1‑146 (MQ=255)
                                                                                                 tagtagtTGCGCAATGTCAGCCCACC‑‑CTTGGCATCTTCCAAGAAGCCCCCCCTCTTAGCGCCCAGGCCGGCAAACCAAGGCCAAACGCCAGCGAAACAGCCAGGACCGAACGGGCCAGCGTGCAAGTTGTGTCCATCATTATTATTGt                                           >  1:280714/1‑148 (MQ=255)
                                                                                                   gtagtTGCGCAATGTCAGCCCACC‑‑CTTGGCATCTTCCAAGAAGTCCCCCCTCTTGCGCCCAGGCCGGCAAACCAAGGCCAAACGCCAGCGAAACAGCCAGGACCGAACGGGCCAGCGTGCAAGTTGTGTCCATCATTATTATTGTtctc                                       <  4:334308/149‑1 (MQ=255)
                                                                                                   gtagtTGCGCAATGTCAGCCCACC‑‑CTTGGCATCTTCCAAGAAGCCCCCCCTCTTGCGCCCAGGCCGTCAAACCAAGGCCAAACGCCAGCGAAACAGCCAGGACCGAACGGGCCAGCGTGCAAGTTGTGTCCATCATTATTATTGTtctc                                       <  3:84120/149‑1 (MQ=255)
                                                                                                          cgcAATGTCAGCCCACC‑‑CTTGGCATCTTCCAAGAAGCCCCCCCTCTTGCGCCCAGGCCGGCAAACCAAGGCCAAACGCCAGCGAAACAGCCAGGACCGAACGGGCCAGCGTGCAAGTTGTGTCCATCATTATTATTGTTCTCCACGtt                                 <  1:227911/148‑1 (MQ=255)
                                                                                                            cAATGTCAGCCCACA‑‑CTTGGCATCTTCCAAGAAGCCCCCCCTCTTGCGCCCAGGCCCGCAAACCAAGGCCAAACGCCAGCGAAACAGCCAGGACCGAACGGGCCAGAGTGAAAGTTGTGTCCATCATTATTATTGTTATCAACGTTCa                               >  4:216766/1‑148 (MQ=255)
                                                                                                               tGTCAGCCCACC‑‑CTTGGCATCTTCCAAGAAGCCCCCCCTCTTGCGCCCAGGCCGGCAAACCAAGGCCAAACGCCAGCGAAACAGCCAGGACCGAACGGGCCAGCGTGCAAGTTGTGTCCATCATTATTATTGTTCTCCACGTTCAATAg                           <  8:292394/149‑1 (MQ=255)
                                                                                                               tGTCAGCCCACC‑‑CTTGGCATCTTCCAAGAAGCCCCCCCTCTTGCGCCCAGGCCGGCAAACCAAGGCCAAACGCCAGCGAAACAGCCAGGACCGAACGGGCCAGCGTGCAAGTTGTGTCCATCATTATTATTGTTCTCCACGTTCAATAg                           <  2:163888/149‑1 (MQ=255)
                                                                                                               tGTCAGCCCACC‑‑CTTGGCATCTTCCAAGAAGCCCCCCCTCTTGCGCCCAGGCCGGCAAACCAAGGCCAAACGCCAGCGAAACAGCCAGGACCGAACGGGCCAGCGTGCAAGTTGTGTCCATCATTATTATTGTTCTCCACGTTCAATAg                           <  3:304131/149‑1 (MQ=255)
                                                                                                               tGTCAGCCCACC‑‑CTTGGCATCTTCCAAGAAGCCCCCCCTCTTGCGCCCAGGCCGGCAAACCAAGGCCAAACGCCAGCGAAACAGCCAGGACCGAACGGGCCAGCGTGCAAGTTGTGTCCATCATTATTATTGTTCTCCACGTTCAATAg                           <  2:280714/149‑1 (MQ=255)
                                                                                                                gTCAGCCCACC‑‑CTTGGCATCTTCCAAGAAGCCCCCCCTCTTGCGCCCAGGCCGGCAAACCAAGGCCAAACGCCAGCGAAACAGCCAGGACCGAACGGGCCAGCGTGCAAGTTGTGTCCATCATTATTATTGTTCTCCACGTTCAATAgg                          >  8:32180/1‑149 (MQ=255)
                                                                                                                   aGCCCACC‑‑CTTGGCATCTTCCAAGAAGCCCCCCCTCTTGCGCCCAGGCCGGCAAACCAAGGCCAAACGCCAGCGAAACAGCCAGGACCGAACGGGCCAGCGTGCAAGTTGTGTCCATCATTATTATTGTTCTCCACGTTCAATa                            <  6:25189/144‑1 (MQ=255)
                                                                                                                   aGCCCACC‑‑CTTGGCATCTTCCAAGAAGCCCCCCCTCTTGCGCCCAGGCCGGCAAACCAAGGCCAAACGCCAGCGAAACAGCCAGGACCGAACGGGCCAGCGTGCAAGTTGTGTCCATCATTATTATTGTTCTCCACGTTCAATa                            >  5:25189/1‑144 (MQ=255)
                                                                                                                   aGCCCACC‑‑CTTGGCATCTTCCAAGAAGCCCCCCCTCTTGCGCCCAGGCCGGCAAACCAAGGCCAAACGCCAGCGAAACAGCCAGGACCGAACGGGCCAGCGTGCAAGTTGTGTCCATCATTATTATTGTTCTCCACGTTCAATAGgcg                        >  5:6090/1‑148 (MQ=255)
                                                                                                                         cc‑‑cTTGGCATCTTCCAAGAAGCCCCCCCTCTTGCGCCCAGGCCGGCAAACCAAGGCCAAACGCCAGCGAAACAGCCAGGACCGAACGGGCCAGCGTGCAAGTTGTGTCCATCATTATTATTGTTCTCCACGTTCAATAGGCGCAGCAGa                 <  8:218863/149‑1 (MQ=255)
                                                                                                                         cc‑‑cTTGGCATCTTCCAAGAAG‑‑CCACCCCTCTTGCGCCCAGGCCTGCAAAACAAGGCCAAACGCCAGCGAAACAGCCACGACCGAACGGGCCAGCGTGCAAGTTGTGTCCATCATTATTATTGTTCTCCACGTTCAATAGGCGCAGCAGa                 <  3:20615/149‑1 (MQ=255)
                                                                                                                              ttGGCATCTTCCAAGAAGCCCCCCCTCTTGCGCCCAGGCCGGCAAACCAAGGCCAAACGCCAGCGAAACAGCCAGGACCGAACGGGCCAGCGTGCAAGTTGTGTCCATCATTATTATTGTTCTCCACGTTCAATAGGCGCAGCAGACCg              >  7:131543/1‑149 (MQ=255)
                                                                                                                                 ggCATCTTCCAAGAAGCCCCCCCTCTTGCGCCCAGGCCGGCAAACCAAGGCCAAACGCCAGCGAAACAGCCAGGACCGAACGGGCCAGCGTGCAAGTTGTGTCCATCATTATTATTGTTCTCCACGTTCAATAGGCGCAGCAGACCg              <  7:166975/147‑1 (MQ=255)
                                                                                                                                 ggCATCTTCCAAGAAGCCCCCCCTCTTGCGCCCAGGCCGGCAAACCAAGGCCAAACGCCAGCGAAACAGCCAGGACCGAACGGGCCAGCGTGCAAGTTGTGTCCATCATTATTATTGTTCTCCACGTTCAATAGGCGCAGCAGACCGc             <  8:136325/148‑1 (MQ=255)
                                                                                                                                 ggCATCTTCCAAGAAGCCCCCCCTCTTGCGCCCAGGCCGGCAAACCAAGGCCAAACGCCAGCGAAACAGACAGGACCGAACGGGCCAGCGTGCAAGTTGTGTCCATCATTATTATTGTTCTCCACGTTCAATAGGCGCAGCAGACCgac            >  3:254935/1‑147 (MQ=255)
                                                                                                                                   cATCTTCCAAGAAGCCCCCCCTCTTGCGCCCAGGCCGGCAAACCAAGGCCAAACGCCAGCGAAACAGCCAGGACCGAACGGGCCAGCGTGCAAGTTGTGTCCATCATTATTATTGTTCTCCACGTTCAATAGGCGCAGCAGACCGcccc          <  5:102898/149‑1 (MQ=255)
                                                                                                                                         ccAAGAAGCCCCCCCTCTTGCGCCCAGGCCGGCAAACCAAGGCCAAACGCCAGCGAAACAGCCAGGACCGAACGGGCCAGCGTGCAAGTTGTGTCCATCATTATTATTGTTCTCCACGTTCAATAGGCGCAGCAGACCGCCCCAg        <  7:205840/145‑1 (MQ=255)
                                                                                                                                         ccAAGAAG‑‑CCCCCCTCTTGCGCCCAGGCCGGCAAACCAAGGCCAAACGCCAGCGAAACAGCCAGGACCGAACGGGCCAGCGTGCAAGTTGTGTCCATCATTATTATTGTTCTCCACGTGCAATAGGCGAAGCAGACGCCCCCAGGCCAg   >  8:141678/1‑149 (MQ=255)
                                                                                                                                           aagaagcCCCCCCTCTTGCGCCCAGGCCGGCAAACCAAGGCCAAACGCCAGCGAAACAGCCAGGACCGAACGGGCCAGCGTGCAAGTTGTGTCCATCATTATTATTGTTCTCCACGTTCAATAGGCGCAGCAGACCGCCCCAGGCCAg   >  1:48990/1‑148 (MQ=255)
                                                                                                                                           aagaagcCCCCCCTCTTGCGCCCAGGCCGGCAAACCAAGGCCAAACGCCAGCGAAACAGCCAGGACCGAACGGGCCAGCGTGCAAGTTGTGTCCATCATTATTATTGTTCTCCACGTTCAATAGGCGCAGCAGACCGCCCCAGGCCAGc  <  8:139578/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
CGGGGTGTAGCCAGAGCTGAACTTGAAGATGAAACCCTGGGCCCACTCTTCAACCTTGCCCTTGGCGGCCCCCGGGTCGCGGAAATCGCGGTTGAAGTAGTAGTTGCGCAATGTCAGCCCACC‑‑CTTGGCATCTTCCAAGAAG‑‑CCCCCCTCTTGCGCCCAGGCCGGCAAACCAAGGCCAAACGCCAGCGAAACAGCCAGGACCGAACGGGCCAGCGTGCAAGTTGTGTCCATCATTATTATTGTTCTCCACGTTCAATAGGCGCAGCAGACCGCCCCAGGCCAGC  >  NC_002947/4267624‑4267907

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: