Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,073,536 G→T 5.2% H240Q (CAC→CAA PP_5327 ← phosphate ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,073,5360GT 5.2% 254.1 / 3.3 77H240Q (CAC→CAAPP_5327phosphate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base G (35/38);  new base T (2/2);  total (37/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.02e-01

ACTCCCAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGTTCAGTTCGGCACGGTCGGCGTCCATGTCGGCCTGGGCGGCGGCGTCCAGCTTGCCGTCCAGCTCCAGCTTGCGCGCTTGCAGGCGCAGGCGCTCAAGGCCAT  >  NC_002947/6073388‑6073681
                                                                                                                                                    |                                                                                                                                                 
aCTCCCAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTgg                                                                                                                                                   <  4:273761/149‑1 (MQ=255)
 cTCCCAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGt                                                                                                                                                  >  8:92387/1‑149 (MQ=255)
      aGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAgg                                                                                                                                             >  1:131944/1‑149 (MQ=255)
      aGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGACTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAgg                                                                                                                                             >  8:284490/1‑149 (MQ=255)
         ccTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCg                                                                                                                                          <  3:205665/149‑1 (MQ=255)
           ttGCCGAAGAAGGTGCCCAACATGGTGAACGCCGACATGACGATTGGCTGAACGGCATCCACTACCTTGCACAGGTTGATTTCCACTTCACGGCCGATGCCATACCGGGACACAAGGCTGACGCGGGAAAAAGCCTGGTGAAGGCAgg                                                                                                                                         >  4:70205/1‑148 (MQ=255)
               cGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAgg                                                                                                                                    <  1:286943/149‑1 (MQ=255)
               cGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAgg                                                                                                                                    >  5:99595/1‑149 (MQ=255)
               cGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAgg                                                                                                                                    >  2:28670/1‑149 (MQ=255)
               cGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAgg                                                                                                                                    >  2:255598/1‑149 (MQ=255)
                 aaGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCg                                                                                                                                  <  4:64711/149‑1 (MQ=255)
                   gTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGt                                                                                                                                <  4:1646/149‑1 (MQ=255)
                    tAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTc                                                                                                                               <  6:108302/149‑1 (MQ=255)
                     aGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCt                                                                                                                              >  6:82622/1‑149 (MQ=255)
                      ggTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCtt                                                                                                                             <  7:83325/149‑1 (MQ=255)
                      ggTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCtt                                                                                                                             <  5:82622/149‑1 (MQ=255)
                           ccAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCCGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGa                                                                                                                          >  8:44007/1‑147 (MQ=255)
                           ccAACTTGGTGAACGCCGACATGCCGGTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGa                                                                                                                          >  6:40574/1‑147 (MQ=255)
                             aaCTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGc                                                                                                                      >  5:200427/1‑149 (MQ=255)
                               cTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCTCGGCCTTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCt                                                                                                                    >  1:223393/1‑149 (MQ=255)
                                  ggTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGATAAAGCCTGGTGCAGGCCGGTGAGGGTGTCTTCGATAGCCtt                                                                                                                   <  2:25605/147‑1 (MQ=255)
                                  ggTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGa                                                                                                                          <  3:29669/140‑1 (MQ=255)
                                  ggTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGa                                                                                                                          >  4:29669/1‑140 (MQ=255)
                                  ggTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCtt                                                                                                                   <  2:131944/147‑1 (MQ=255)
                                  ggTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTg                                                                                                                  <  6:200427/148‑1 (MQ=255)
                                     gAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGTTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGc                                                                                                              >  7:221394/1‑149 (MQ=255)
                                     gAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGTTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGc                                                                                                              >  7:225674/1‑149 (MQ=255)
                                      aaCGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGGCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGc                                                                                                              <  3:70205/148‑1 (MQ=255)
                                         gCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGCt                                                                                                          >  8:167339/1‑149 (MQ=255)
                                            gACCTGCCGTTTGGCTGAACGGCATACACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGtt                                                                                                       <  8:295028/149‑1 (MQ=255)
                                              cATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGTTCa                                                                                                     <  7:167339/149‑1 (MQ=255)
                                                 gCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGTTCAGtt                                                                                                  >  2:83625/1‑149 (MQ=255)
                                                 gCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGTTCAGtt                                                                                                  >  4:211119/1‑149 (MQ=255)
                                                    gTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGTTCAGTTCgg                                                                                               <  6:99595/149‑1 (MQ=255)
                                                    gTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGTTCAGTTCgg                                                                                               <  3:211119/149‑1 (MQ=255)
                                                    gTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGTTCAGTTCgg                                                                                               >  6:12568/1‑149 (MQ=255)
                                                     tttGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGTTCAGTTCGGc                                                                                              <  5:44451/149‑1 (MQ=255)
                                                         gCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGTTCAGTTCGGCACgg                                                                                          >  7:170499/1‑149 (MQ=255)
                                                             aaCGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGTTCAGTTCGGCACGGTCg                                                                                       <  7:284490/148‑1 (MQ=255)
                                                              aCGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGTTCAGTTCGGCACGGTCGGc                                                                                     <  5:318196/149‑1 (MQ=255)
                                                                  cATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGTTCAGTTCGGCa                                                                                             <  6:203728/137‑1 (MQ=255)
                                                                  cATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCCGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGTTCAGTTCGGCACGGTCGGCGTcc                                                                                 <  7:44007/149‑1 (MQ=255)
                                                                    tGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGTTCAGTTCGGCACGGTCGGCGTCCAt                                                                               <  7:92387/149‑1 (MQ=255)
                                                                       aCTACCTTGCTCAGGTTGATTTCCACTTCGCGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGTTCAGTTCGGCACGGTCg                                                                                       <  2:209008/138‑1 (MQ=255)
                                                                       aCTACCTTGCTCAGGTTGATTTCCACTTCGCGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGTTCAGTTCGGCACGGTCg                                                                                       >  1:209008/1‑138 (MQ=255)
                                                                           ccTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGTTCAGTTCGGCACGGTCGGCGTCCATGTCGGcc                                                                        >  1:235937/1‑149 (MQ=255)
                                                                            cTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGTTCAGGCgg                                                                                                                                 <  4:39107/91‑1 (MQ=255)
                                                                            cTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGTTCAGTTCGGCACGGTCGGCGTCCATGTCGGCCt                                                                       <  1:255598/149‑1 (MQ=255)
                                                                            cTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGATCAGGCgg                                                                                                                                 >  3:39107/1‑91 (MQ=255)
                                                                               gCTCAGGTTGATTTCCACTTCTCGGCCTTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGTTCAGTTCGGCACGGTCGGCGTCCATGTCGGCCTgg                                                                     <  2:223393/148‑1 (MQ=255)
                                                                               gCTCAGGTTGATTTCCACTTCACGGCAGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGTTCAGTTCGGCACGGTCGGCGTCCATGTCGGCCTggg                                                                    >  1:4978/1‑149 (MQ=255)
                                                                               gCTCAGGTTGAATTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGTTCAGTTCGGCACGGTCGGCGTCCATGTCGGCCTg                                                                      <  6:273550/147‑1 (MQ=255)
                                                                                   aGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGTTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGTTCAGTTCGGCACGGTCGGCGTCCATGTCGGCCTGggcggc                                                                <  8:225674/149‑1 (MQ=255)
                                                                                   aGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGTTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGTTCAGTTCGGCACGGTCGGCGTCCATGTCGGCCTGggcggc                                                                <  8:221394/149‑1 (MQ=255)
                                                                                   aGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGTTCAGTTCGGCACGGTCGGCGTCCATGTCGGCCTGggcggc                                                                <  8:170499/149‑1 (MQ=255)
                                                                                    ggTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGTTCAGTTCGGCACGGTCGGCGTCCATGTCGGCCTGggcggcg                                                               <  2:235937/149‑1 (MQ=255)
                                                                                            tacACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGTTCAGTTCGGCACGGTCGGCGTCCATGTCGGCCTGGGCGGCGGCGTCCAg                                                       >  4:228250/3‑149 (MQ=255)
                                                                                            tCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGGTCAGTTCGGCACGGTCGGCGTCCATGTCGGCCTGGGCGGCGGCGTCCAg                                                       >  2:307095/1‑149 (MQ=255)
                                                                                                cTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGTTCAGTTCGGCACGGTCGGCGTCCATGTCGGCCTGGGCGGCGGCGTCCAGCTTg                                                   >  4:93287/1‑149 (MQ=255)
                                                                                                    aCGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGTTCAGTTCGGCACGGTCGGCGTCCATGTCGGCCTGGGCGGCGGCGTCCAGCTTGcc                                                 >  2:295729/1‑147 (MQ=255)
                                                                                                    aCGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGTTCAGTTCGGCACGGTCGGCGTCCATGTCGGCCTGGGCGGCGGCGTCCAGCTTGcc                                                 <  1:295729/147‑1 (MQ=255)
                                                                                                    aCGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGTTCAGTTCGGCACGGTCGGCGTCCATGTCGGCCTGGGCGGCGGCGTCCAGCTTGCCGt                                               >  6:25941/1‑149 (MQ=255)
                                                                                                    aCGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGTTCAGTTCGGCACGGTCGGCGTCCATGTCGGCCTGGGCGGCGGCGTCCAGCTTGCCGt                                               >  6:131579/1‑149 (MQ=255)
                                                                                                     cGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGGGGTCTTCGATAGCCTTGTAGCGGCTGTTCAGTTCGGCACGGTCGGCGTCCATGTCGGCCTGGGCGGCGGCGTCCAGCTTGCCGtc                                              >  5:75604/1‑149 (MQ=255)
                                                                                                     cGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGTTCAGTTCGGCACGGTCGGCGTCCATGTCGGCCTGGGCGGCGGCGTCCAGCTTGCCGtc                                              <  2:4978/149‑1 (MQ=255)
                                                                                                        ccGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAATGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGTTCAGTTCGGCACGGTCGGCGTCCATGTCGGCCTGGGCGGCGGCGTCCAGCTTGCCGtccag                                           <  6:4631/149‑1 (MQ=255)
                                                                                                          gTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGTTCAGTTCGGCACGGTCGGCGTCCATGTCGGCCTGGGCGGCGGCGTCCAGCTTGCCGtccagct                                         >  7:136423/1‑149 (MQ=255)
                                                                                                                 tCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGTTCAGTTCGGCACGGTCGGCGTCCATGTCGGCCTGggcggcg                                                               >  6:205096/1‑120 (MQ=255)
                                                                                                                 tCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGTTCAGTTCGGCACGGTCGGCGTCCATGTCGGCCTGggcggcg                                                               <  5:205096/120‑1 (MQ=255)
                                                                                                                 tCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGTTCAGTTCGGCACGGTCGGCGTCCATGTCGGCCTGGGCGGCGGCGTCCAGCTTGCCGtccagctcc                                       <  1:114130/144‑1 (MQ=255)
                                                                                                                 tCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGTTCAGTTCGGCACGGTCGGCGTCCATGTCGGCCTGGGCGGCGGCGTCCAGCTTGCCGtccagctcc                                       >  2:114130/1‑144 (MQ=255)
                                                                                                                        ccaccaGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGTTCAGTTCGGCACGGTCGGCGTCCATGTCGGCCTGGGCGGCGGCGTCCAGCTTGCCGTCCAGCTCCAGCTTGCGCGCt                           >  6:91842/1‑149 (MQ=255)
                                                                                                                         caccaGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGTTCAGTTCGGCACGGTCGGCGTCCATGTCGGCCTGGGCGGCGGCGTCCAGCTTGCCGTCCAGCTCCAGCTTGCGCGCtt                          <  8:136423/149‑1 (MQ=255)
                                                                                                                          accaGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCTATAGCCTTGTAGCGGCTGTTCAGTTCGGCACGGTCGGCGTCCATGTCGGCCTGGGCGGCGGCGTCCAGCTTGCCGTCCAGCTCCAGCTTGCGCGCTTg                         <  3:93287/149‑1 (MQ=255)
                                                                                                                          accaGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGTTCAGTTCGGCACGGTCGGCGTCCATGTCGGCCTGGGCGGCGGCGTCCAGCTTGCCGTCCAGCTCCAGCTTGCGCGCTTg                         <  1:83625/149‑1 (MQ=255)
                                                                                                                                           aaaaGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGTTCAGTTCGGCACGGTCGGCGTCCATGTCGGCCTGGGCGGCGGCGTCCAGCTTGCCGTCCAGCTCCAGCTTGCGCGCTTGCAGGCGCAGGCGCTc          >  2:304336/1‑147 (MQ=255)
                                                                                                                                               gCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGTTCAGTTCGGCACGGTCGGCGTCCATGTCGGCCTGGGCGGCGGCGTCCAGCTTGCCGTCCAGCTCCAGCTTGCGCGCTTGCAGGCGCAGGCGCTCAAGGc     >  1:295293/1‑148 (MQ=255)
                                                                                                                                                 cTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGTTCAGTTCGGCACGGTCGGCGTCCATGTCGGCCTGGGCGGCGGCGTCCAGCTTGCCGTCCAGCTCCAGCTTGCGCGCTTGCAGGCGCAGGCGCTCAAGGCCAt  <  5:91842/149‑1 (MQ=255)
                                                                                                                                                   ggTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGTTCAGTTCGGCACGGTCGGCGTCCATGTCGGCCTGGGCGGCGGCGTCCAGCTTGCCGTCCAGCTCCAGCTTGCGCGCTTgcaggcgcag                >  6:307051/1‑133 (MQ=255)
                                                                                                                                                   ggTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGTTCAGTTCGGCACGGTCGGCGTCCATGTCGGCCTGGGCGGCGGCGTCCAGCTTGCCGTCCAGCTCCAGCTTGCGCGCTTgcaggcgcag                <  5:307051/133‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
ACTCCCAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGTTCAGTTCGGCACGGTCGGCGTCCATGTCGGCCTGGGCGGCGGCGTCCAGCTTGCCGTCCAGCTCCAGCTTGCGCGCTTGCAGGCGCAGGCGCTCAAGGCCAT  >  NC_002947/6073388‑6073681

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: