Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,289,589 A→T 6.7% L246H (CTC→CAC) ‡ PP_4663 ← transmembrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,289,5890AT 6.7% 162.5 / 2.8 60L246H (CTC→CAC) ‡PP_4663transmembrane protein
Reads supporting (aligned to +/- strand):  ref base A (34/21);  new base T (2/2);  total (36/24)
Fisher's exact test for biased strand distribution p-value = 6.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.17e-01

GGCCGCGCTGGACGTGCACTCGGCCCTTCCGCACATGGCCGCCTTCTTCGAAAAGCGCATGGACTACCTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCAGCATCAGCAGCGACGCCCCCAGCGCCTGGCCACCCAAGGCATCCTG  >  NC_002947/5289442‑5289734
                                                                                                                                                   |                                                                                                                                                 
ggCCGCGCTGGACGTGCACTCGGCCCTTCCGCACATGGCCGCCTTCTTCGAAAAGCGCATGGACTACCTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAg                                                                                                                                                  >  8:193226/1‑149 (MQ=255)
ggCCGCGCTGGACGTGCACTCGGCCCTTCCGCACATGGCCGCCTTCTTCGAAAAGCGCATGGACTACCTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAg                                                                                                                                                  >  7:32615/1‑149 (MQ=255)
        tGGTCGTGCACTAGTGCCTTCCTCACATGGCCGCCTTCTTCGAAAAGCGCATGGACTTCCTGGAGCTCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCCTCGGCGGCCa                                                                                                                                          <  8:205341/149‑1 (MQ=255)
              tGCACTCGGCCCTTCCGCACATGGCCGCCTTCTTCGAAAAGCGCATGGACTACCTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCCTCGGCGGCCacctgct                                                                                                                                    <  3:307507/149‑7 (MQ=255)
               gCACTCGGCCCTTCCGCACATGGCCGCCTTCTTCGAAAAGCGCATGGACTACCTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGAt                                                                                                                                   >  6:169425/1‑149 (MQ=255)
                       cccTTCCTCACAGCGCCGCCTTCTTCCAAATCCCCATCGACTCCCTGGAGCGCTATTTCCATATCCCCCGGGTCGACATGCTCCGCGCTTCGTCCGCTAACCGTTGATCACCAGGTGGCCGCCGATGGCGGCCAGGCCGATGAAAAAGc                                                                                                                           <  8:154613/149‑1 (MQ=255)
                        ccTTCCGACCATGGCCGCCTTCTTCCAAAAGCGCATGGACTACCTGGAGCGCGATTTCTATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCTTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCa                                                                                                                          <  4:224104/149‑1 (MQ=255)
                          ttCCGCACATGGCCGCCTTCTTCGAAAAGCGCATGGACTACCTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGATGATCAGCAGGTGGCCGCCGAGCGCGGCCAGGCCGATGAAAAAGCAAc                                                                                                                        >  5:65867/1‑149 (MQ=255)
                           tCCGCACATGGCCGCCTTCTTCGAAAAGCGCATGGACTACCTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGAAGGTGGCCACCGAGGGCGGCCAGGCCGATGAAAAAGCAACg                                                                                                                       >  3:63262/1‑149 (MQ=255)
                            ccGCACATGGCCGCCTTCTTCGAAAAGCGCATGGACTACCTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTTGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGc                                                                                                                      >  4:13179/1‑149 (MQ=255)
                            ccGCACATGGCCGCCTTCTTCGAAAAGCGCATGGACTACCTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGACGGCCAGGCCGATGAAAAAGCAAc                                                                                                                        <  3:149690/147‑1 (MQ=255)
                               cacaTGGCCGCCTTCTTCGAAAAGCGCATGGACTACCTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCtt                                                                                                                    >  3:220262/1‑148 (MQ=255)
                                acaTGGCCGCCTTCTTCGAAAAGCGCATGGACTACCTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTg                                                                                                                   >  4:26732/1‑148 (MQ=255)
                                     gccgccTTCTTCGAAAAGCGCATGGACTACCTGGAGTGCGATTTCGATATCCAGCGTGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAAcag                                                                                                              <  5:291953/148‑1 (MQ=255)
                                      ccgccTTCTTCGAAAAGCGCATGGACTACCTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAAcag                                                                                                              >  5:329305/1‑147 (MQ=255)
                                        gccTTCTTCGAAAAGCGCATGGACTACCTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACagcagc                                                                                                          <  5:163579/149‑1 (MQ=255)
                                          cttcttCGAAAAGCGCATGGACTACCTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAgcgc                                                                                                        >  4:295712/1‑149 (MQ=255)
                                                  aaaaGCGCATGGACTACCTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATcc                                                                                                >  7:235776/1‑149 (MQ=255)
                                                                tACCTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTg                                                                                                                   >  7:66740/1‑116 (MQ=255)
                                                                tACCTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTg                                                                                                                   <  8:66740/116‑1 (MQ=255)
                                                                  ccTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTg                                                                                >  3:187504/1‑149 (MQ=255)
                                                                  ccTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTg                                                                                >  7:56909/1‑149 (MQ=255)
                                                                     ggAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGcc                                                                             >  5:259507/1‑149 (MQ=255)
                                                                     ggAGCGCGATTTCGATATCAAGCGGGTCGACATGCTAAGCGCTTCGTCCGATAACCGTTGATCAGAAGGTGGCCGCCGAGGGAGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATACGCTGACGCAGCCACTGGc                                                                              >  6:314393/1‑148 (MQ=255)
                                                                                        cAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAg                                                          >  2:269339/1‑149 (MQ=255)
                                                                                        cAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAg                                                          <  8:206073/149‑1 (MQ=255)
                                                                                         aGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGc                                                                              >  8:99386/1‑128 (MQ=255)
                                                                                         aGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGc                                                                              <  7:99386/128‑1 (MQ=255)
                                                                                         aGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAgg                                                         >  1:174702/1‑149 (MQ=255)
                                                                                            gggTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCg                                                      <  5:169425/149‑1 (MQ=255)
                                                                                              gTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGgc                                                    <  3:295712/149‑1 (MQ=255)
                                                                                                 gACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGgcgc                                                  >  2:206199/1‑148 (MQ=255)
                                                                                                  aCATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCa                                                <  2:174702/149‑1 (MQ=255)
                                                                                                  aCATGCCCTGCGCTGCGCCCCATACCCGTTCATCCGCTGGGGCCCGCCGCGGGCGCCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCa                                                <  3:2525/149‑1 (MQ=255)
                                                                                                      gCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCAGCAt                                            >  4:97732/1‑149 (MQ=255)
                                                                                                       cTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCAGCAt                                            >  6:48457/1‑148 (MQ=255)
                                                                                                        tCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCAGCAt                                            >  3:302686/1‑147 (MQ=255)
                                                                                                        tCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCAGCAt                                            >  1:295308/1‑147 (MQ=255)
                                                                                                        tCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCAGCATca                                          >  8:116938/1‑149 (MQ=255)
                                                                                                          aGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCAGCAt                                            >  3:152765/1‑145 (MQ=255)
                                                                                                              cTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCAGCATCagcagc                                     >  4:183626/1‑148 (MQ=255)
                                                                                                               ttCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCAGCATcagcag                                      <  2:214608/146‑1 (MQ=255)
                                                                                                               ttCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCAGCATcagcag                                      >  1:214608/1‑146 (MQ=255)
                                                                                                                     cGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCAGCATCAGCAGCGACGccc                              <  5:227414/148‑1 (MQ=255)
                                                                                                                       cTAACCGTTGATCAGCAGGTGGCCACCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCAGCATCAGCAGCGACGCCCCCa                           <  5:60718/149‑1 (MQ=255)
                                                                                                                              ttGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCAGCATCAGCAGCGACGCCCCCAGCGCCTg                    >  1:323752/1‑149 (MQ=255)
                                                                                                                              ttGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCAGCATCAGCAGCGACGCCCCCAGCGCCTg                    <  8:32615/149‑1 (MQ=255)
                                                                                                                              ttGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCAGCATCAGCAGCGACGCCCCCAGCGCCTg                    <  6:329305/149‑1 (MQ=255)
                                                                                                                                gATCCGCACGGGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCAGCATCAGCAGCGACGCCCCCAGCGCCTg                    <  4:220262/147‑1 (MQ=255)
                                                                                                                                  tCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCAGCATCAGCAGCGACGCCCCCAGCGCCTGGCCa                >  4:43797/1‑149 (MQ=255)
                                                                                                                                    agcagGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCAGCATCAGCAGCGACGCCCCCAGCGCCTGGCCAcc              >  4:85287/1‑149 (MQ=255)
                                                                                                                                       gcctGGCCGCCCTCGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGGGCCAGCATCAGCAGCGACGCCCCCAGCGCCTGGCCAccc             >  5:280286/4‑147 (MQ=255)
                                                                                                                                       gcctGGCCGCCCTCGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCAGCATCAGCAGCGACGCCCCCAGCGCCTGGCCAccc             >  5:136850/4‑147 (MQ=255)
                                                                                                                                       agGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGGTCATTGATTTGCTCTATATGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCAGCATCAGCAGCGACGCCCCCAGCGCCTGGCCACCCaa           <  1:269339/149‑1 (MQ=255)
                                                                                                                                        ggTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCAGCATCAGCAGCGACGCCCCCAGCGCCTGGCCACCCAAg          <  8:25137/149‑1 (MQ=255)
                                                                                                                                             ccgccgAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCAGCATCAGCAGCGACGCCCCCAGCGCCTGGCCACCCAAGGCATc     >  8:198398/1‑149 (MQ=255)
                                                                                                                                             ccgccgAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCAGCATCAGCAGCGACGCCCCCAGCGCCTGGCCACCCAAGGCATc     <  7:198398/149‑1 (MQ=255)
                                                                                                                                             ccgccgAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCAGCATCAGCAGCGACGCCCCCAGCGCCTGGCCACCCAAGGCATc     >  4:176245/1‑149 (MQ=255)
                                                                                                                                              cgccgAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCAGCATCAGCAGCGACGCCCCCAGCGCCTGGCCACCCAAGGCATcc    >  3:126564/1‑149 (MQ=255)
                                                                                                                                                ccgAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCAGCATCAGCAGCGACGCCCCCAGCGCCTGGCCACCCAAGGCATCCTg  <  7:116938/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GGCCGCGCTGGACGTGCACTCGGCCCTTCCGCACATGGCCGCCTTCTTCGAAAAGCGCATGGACTACCTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCAGCATCAGCAGCGACGCCCCCAGCGCCTGGCCACCCAAGGCATCCTG  >  NC_002947/5289442‑5289734

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