Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,541,781 C→A 8.3% Y166* (TAC→TAA) ‡ nhaB → Na(+)/H(+) antiporter NhaB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,541,7810CA 8.3% 186.8 / 6.8 60Y166* (TAC→TAA) ‡nhaBNa(+)/H(+) antiporter NhaB
Reads supporting (aligned to +/- strand):  ref base C (31/24);  new base A (2/3);  total (33/27)
Fisher's exact test for biased strand distribution p-value = 6.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.98e-01

TCACGCATTTTGCTCGGGGTGCGTTCCAAGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCTGTGGGTACCGCCCTG  >  NC_002947/4541638‑4541928
                                                                                                                                               |                                                                                                                                                   
tcACGCATTTTGCTCGGGGTGCGTTCCAAGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACcg                                                                                                                                                >  2:222229/1‑149 (MQ=255)
tcACGCATTTTGCTCGGGGTGCGTTCCAAGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACcg                                                                                                                                                >  8:217517/1‑149 (MQ=255)
  aCGCATTTTGCTCGGGGTGCGTTCCAAGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGGTACACcgcg                                                                                                                                              <  6:30349/149‑1 (MQ=255)
   cGCATTTTGCTCGGGGTGCGTTCCAAGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGGTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGGTACACcgcgt                                                                                                                                             <  5:358671/149‑1 (MQ=255)
       ttttGCTCGGGGTGCGTTCCAAGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACcgcgtcgcg                                                                                                                                         >  8:149202/1‑149 (MQ=255)
            cTCGGGGTGCGTTCCAAGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCgg                                                                                                                                    >  5:67859/1‑149 (MQ=255)
               ggggTGCGTTCCAAGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCgggg                                                                                                                                  <  6:106887/148‑1 (MQ=255)
               ggggTGCGTTCCAAGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCgggg                                                                                                                                  <  4:321792/148‑1 (MQ=255)
                      gTTCCAAGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCcg                                                                                                                          >  5:268649/1‑149 (MQ=255)
                           aaGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGa                                                                                                                     <  3:32439/149‑1 (MQ=255)
                           aaGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGa                                                                                                                     <  7:149202/149‑1 (MQ=255)
                            aGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAg                                                                                                                    <  1:222229/149‑1 (MQ=255)
                            aGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAg                                                                                                                    <  5:323435/149‑1 (MQ=255)
                              gCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCgcgtcgcgtc                                                                                                                                       <  5:287334/128‑1 (MQ=255)
                              gCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCgcgtcgcgtc                                                                                                                                       >  6:287334/1‑128 (MQ=255)
                                  ttCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAgc                                                                                                              >  5:212126/1‑149 (MQ=255)
                                        cgcTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTg                                                                                                        <  8:140526/149‑1 (MQ=255)
                                           tgctgTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTcc                                                                                                                                      >  3:89187/1‑116 (MQ=255)
                                           tgctgTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTcc                                                                                                                                      <  4:89187/116‑1 (MQ=255)
                                           tgctgTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCGACGCGGTGGTACACcgcgtagaa                                                                                                                                         >  2:290064/1‑109 (MQ=255)
                                           tgctgTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCGACGCGGTGGTACACcgcgtagaa                                                                                                                                         <  1:290064/113‑5 (MQ=255)
                                                ttttGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGa                                                                                                <  6:268649/149‑1 (MQ=255)
                                                  ttGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGAc                                                                                               >  3:323500/1‑148 (MQ=255)
                                                   tGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCa                                                                                             >  8:320108/1‑149 (MQ=255)
                                                     cGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGc                                                                                           <  7:34185/149‑1 (MQ=255)
                                                       tACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAgc                                                                                                              >  2:123585/1‑128 (MQ=255)
                                                       tACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCaa                                                                                         <  4:323500/149‑1 (MQ=255)
                                                            tCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTc                                                                                      <  7:217517/147‑1 (MQ=255)
                                                            tCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGc                                                                                    >  5:121592/1‑149 (MQ=255)
                                                                      tATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCAc                                                                          >  7:244299/1‑149 (MQ=255)
                                                                          cgcgTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCgg                                                                       >  4:164041/1‑148 (MQ=255)
                                                                           gcgTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCggg                                                                      >  4:272959/1‑148 (MQ=255)
                                                                            cgTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGaa                                                                    <  7:320108/149‑1 (MQ=255)
                                                                              tttCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGa                                                                  >  8:324472/1‑149 (MQ=255)
                                                                                 ctcGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCt                                                               >  2:308002/1‑149 (MQ=255)
                                                                                  tcGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGTGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTc                                                              >  4:313720/1‑149 (MQ=255)
                                                                                     aTGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTCCCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCCCCGGGAAGACCTCGAc                                                           >  3:90911/1‑149 (MQ=255)
                                                                                        cccTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCa                                                         <  6:121592/148‑1 (MQ=255)
                                                                                          cTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGtt                                                      >  3:62528/1‑149 (MQ=255)
                                                                                             aCCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTGcc                                                    <  7:289285/148‑1 (MQ=255)
                                                                                              ccGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCcgc                                                  <  4:62528/149‑1 (MQ=255)
                                                                                               cGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCcgcg                                                 <  5:159147/149‑1 (MQ=255)
                                                                                                gTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCgcgc                                                >  4:68063/1‑149 (MQ=255)
                                                                                                     ggcggTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTcc                                           <  7:324472/149‑1 (MQ=255)
                                                                                                        ggTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGTGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGc                                        <  3:313720/149‑1 (MQ=255)
                                                                                                         gTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCAACGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCGCACCTGCTCCGGGCAGACCTCGACCAGTTCCGCGCCTTCCTGCg                                       >  3:52731/1‑149 (MQ=255)
                                                                                                               atcaGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCct                                 >  6:114301/1‑149 (MQ=255)
                                                                                                                   gcgcTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCctgct                              >  5:299681/1‑148 (MQ=255)
                                                                                                                    cgcTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCTCACTGGAGAGCGAGCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGa                            <  6:212126/149‑1 (MQ=255)
                                                                                                                    cgcTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGa                            <  3:164041/149‑1 (MQ=255)
                                                                                                                     gcTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGAt                           <  5:114301/149‑1 (MQ=255)
                                                                                                                         cGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTcc                                                    >  4:305916/1‑120 (MQ=255)
                                                                                                                            tCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGgc                    >  4:226387/1‑149 (MQ=255)
                                                                                                                               gCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCt                 >  8:309053/1‑149 (MQ=255)
                                                                                                                                 ttCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGCTGCATGGCGCtgt               >  5:331707/1‑149 (MQ=255)
                                                                                                                                   cgaCGCGGTGGTACACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCTGTgg             >  1:147569/3‑149 (MQ=255)
                                                                                                                                    tACGCTGTGTACCTCCGCGTCGCGTCCGTGGCCAATCCTCGCGAGGATATTGCACTGGAGAGCGACCAGCAAGTGGCTCAGGTGGAGCGGGAAGAGCTTGACCAGTTTCTAGGCTTGGTGGGTAGGGTGCTGATGCATGTCGGTGTggg            <  4:40785/149‑1 (MQ=255)
                                                                                                                                     aCGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCTGTggg            >  5:152871/1‑148 (MQ=255)
                                                                                                                                     aCGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCTGTggg            >  5:216641/1‑148 (MQ=255)
                                                                                                                                     aCGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCTGTggg            >  5:221749/1‑148 (MQ=255)
                                                                                                                                              accaccGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCTGTGGGTACCGCCCTg  >  1:286510/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
TCACGCATTTTGCTCGGGGTGCGTTCCAAGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCTGTGGGTACCGCCCTG  >  NC_002947/4541638‑4541928

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