Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,073,444 G→A 8.5% P271L (CCA→CTA) ‡ 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,4440GA 8.5% 266.8 / 10.2 94P271L (CCA→CTA) ‡PP_5327phosphate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base G (46/39);  new base A (6/2);  total (52/42)
Fisher's exact test for biased strand distribution p-value = 4.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

CGGCCATGATCAGGGTCATCATCACGGTACCGAAGATGGCCGGGAAGATACCGCCCTCGGTGTTGGCTTCACGCGGGTCGTCGCTGAGGAACTCCCAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGTTCAGTTC  >  NC_002947/6073298‑6073586
                                                                                                                                                  |                                                                                                                                              
cGGCCATGATCAGGGTCATCATCACGGTACCGAAGATGGCCGGGAAGATACCGCCCTCGGTGTTGGCTTCACGCGGGTCGTCGCTGAGGAACTCCCAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGc                                                                                                                                              <  1:304669/149‑1 (MQ=255)
cGGCCATGATCAGGGTCATCATCACGGTACCGAAGATGGCCGGGAAGATACCGCCCTCGGTGTTGGCTTCACGCGGGTCGTCGCTGAGGAACTCCCAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGc                                                                                                                                              <  5:143009/149‑1 (MQ=255)
         tCAGGGTCATCATCACGGTACCGAAGATGGCCGGGAAGATACCGCCCTCGGTGTTGGCTTCACGCGGGTCGTCGCTGAGGAACTCCCAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCa                                                                                                                                     >  4:256421/1‑149 (MQ=255)
                    atcaCGGTACCGAAGATGGCCGGGAAGATACCGCCCTCGGTGTTGGCTTCACGCGGGTCGTCGCTGAGGAACTCCCAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCtt                                                                                                                          >  8:294306/1‑149 (MQ=255)
                     tcaCGGTACCGAAGATGGCCGGGAAGATACCGCCCTCGGTGTTGGCTTCACGCGGGTCGTCGCTGAGGAACTCCCAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTg                                                                                                                         <  4:176710/149‑1 (MQ=255)
                      caCGGTACCGAAGATGGCCGGGAAGATACCGCCCTCGGTGTTGGCTTCACGCGGGTAGTCCCTGAGGAACTCCCAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGc                                                                                                                        >  2:57613/1‑149 (MQ=255)
                           tACCGAAGATGGCCGGGAAGATACCGCCCTCGGTGTTGGCTTCACGCGGGTCGTCGCTGAGGAACTCCCAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAgg                                                                                                                   <  6:360765/149‑1 (MQ=255)
                             ccGAAGATGGCCGGGAAGATACCGCCCTCGGTGTTGGCTTCACGCGGGTCGTCGATGAGGAACTCCCAGACATTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGATGAACGGCATGCACTACCTTGCTCAGGtt                                                                                                                 >  5:318483/1‑149 (MQ=255)
                               gAAGATGGCCGGGAAGATACCGCCCTCGGTGTTGGCTTCACGCGGGTCGTCGCTGAGGAACTCCCAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGa                                                                                                               >  7:90601/1‑149 (MQ=255)
                                  gATGGCCGGGAAGATACCGCCCTCGGTGTTGGCTTCACGCGGGTCGTCGCTGAGGAACTCCCAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGAttt                                                                                                            >  3:89349/1‑149 (MQ=255)
                                     ggCCGGGAAGATACCGCCCTCGGTGTTGGCTTCACGCGGGTCGTCGCTGAGGAACTCCCAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCa                                                                                                         >  3:38334/1‑149 (MQ=255)
                                      gCCGGGAAGATACCGCCCTCGGTGTTGGCTTCACGCGGGTCGTCGCTGAGGAACTCCCAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCa                                                                                                         <  2:64594/148‑1 (MQ=255)
                                      gCCGGGAAGATACCGCCCTCGGTGTTGGCTTCACGCGGGTCGTCGCTGAGGAACTCCCAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCa                                                                                                         <  7:99427/148‑1 (MQ=255)
                                       ccGGGAAGATACCGCCCTCGGTGTTGGCTTCACGCGGGTCGTCGCTGAGGAACTCCCAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACt                                                                                                       >  1:205681/1‑149 (MQ=255)
                                        cGGGAAGATACCGCCCTCGGTGTTGGCTTCACGCGGGTCGTCGCTGAGGAACTCCCAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACtt                                                                                                      >  6:334747/1‑149 (MQ=255)
                                           gAAGATACCGCCCTCGGTGTTGGCTTCACGCGGGTCGTCGCTGAGGAACTCCCAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCAc                                                                                                   <  3:256421/149‑1 (MQ=255)
                                               aTACCGCCCTCGGTGTTGGCTTCACGCGGGTCGTCGCTGAGGAACTCCCAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGcc                                                                                               >  8:278189/1‑149 (MQ=255)
                                                tACCGCCCTCGGTTGTGGCTTCACGCGGGTCCTCGCTGAGGAACTGCCAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCg                                                                                              <  8:34012/149‑1 (MQ=255)
                                                 aCCGCCCTCGGTGTTGGCTTCACGCGGGTCGTCGCTGAGGAACTCCCAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGt                                                                                             >  2:95223/1‑149 (MQ=255)
                                                   cGCCCTCGGTGTTGGCTTCACGCGGGTCGTCGCTGAGGAACTCCCAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTg                                                                                           >  5:165536/1‑149 (MQ=255)
                                                   cGCCCTCGGTGTTGGCTTCACGCGGGTCGTCGCTGAGGAACTCCCAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTg                                                                                           >  6:172239/1‑149 (MQ=255)
                                                    gCCCTCGGTGTTGGCTTCACGCGGGTCGTCGCTGAGGAACTCCCAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGc                                                                                          <  4:155579/149‑1 (MQ=255)
                                                            tgtTGGCTTCACGCGGGTCGTCGCTGAGGAACTCCCAGACCTTGCCGAAGTAGGTGCCCAAATTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCAGTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCgg                                                                                  >  1:106904/1‑149 (MQ=255)
                                                              ttGGCTTCACGCGGGTCGTCGCTGAGGAACTCCCAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCggg                                                                                 <  1:95223/148‑1 (MQ=255)
                                                                 gCTTCACGCGGGTCGTCGCTGAGGAACTCCCAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGccac                                                                             >  3:238089/1‑149 (MQ=255)
                                                                   ttCACGCGGGTCGTCGCTGAGGAACTCCCAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGccac                                                                             >  6:32203/1‑147 (MQ=255)
                                                                   ttCACGCGGGTCGTCGCTGAGGAACTCCCAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGccac                                                                             >  1:180986/1‑147 (MQ=255)
                                                                      aCGCGGGTCGTCGCTGAGGAACTCCCAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGc                                                                        >  7:353432/1‑149 (MQ=255)
                                                                      aCGCGGGTCGTCGCTGAGGAACTCCCAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGc                                                                        >  7:224448/1‑149 (MQ=255)
                                                                      aCGCGGGTCGTCGCTGAGGAACTCCCAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGc                                                                        >  3:136717/1‑149 (MQ=255)
                                                                      aCGCGGGTCGTCGCTGAGGAACTCCCAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGc                                                                        >  6:213723/1‑149 (MQ=255)
                                                                            gtcgtcGCTGAGGAACTCCCAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTcgc                                                                  >  7:214310/1‑149 (MQ=255)
                                                                               gtcgCTGAGGAACTCCCAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCgg                                                                <  1:57613/148‑1 (MQ=255)
                                                                               gtcgCTGAGGAACTCCCAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCgg                                                                <  4:38334/148‑1 (MQ=255)
                                                                                  gCTGAGGAACTCCCAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGaa                                                             <  2:180986/148‑1 (MQ=255)
                                                                                     gAGGAACTCCCAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCgg                                                                                  <  6:321031/124‑1 (MQ=255)
                                                                                     gAGGAACTCCCAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCgg                                                                                  >  5:321031/1‑124 (MQ=255)
                                                                                      aGGAACTCCCAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGaaaaa                                                          <  6:165536/147‑1 (MQ=255)
                                                                                      aGGAACTCCCAGACCTTGCCGAAGTAGGTGCCCAACTTGGGGAACGCCGACCTGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCTGGaaaaacc                                                        <  5:116443/149‑3 (MQ=255)
                                                                                        gAACTCCCAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCt                                                      >  4:171671/1‑149 (MQ=255)
                                                                                        gAACTCCCAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCt                                                      >  3:115320/1‑149 (MQ=255)
                                                                                          aCTCCCAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTgg                                                    <  5:172239/149‑1 (MQ=255)
                                                                                            tCCCAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGt                                                   >  3:226437/1‑148 (MQ=255)
                                                                                            tCCCAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTg                                                  >  5:280587/1‑149 (MQ=255)
                                                                                            tCCCAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTg                                                  >  2:328418/1‑149 (MQ=255)
                                                                                              ccAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCgg                                                                                  <  5:182967/115‑1 (MQ=255)
                                                                                              ccAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCgg                                                                                  >  6:182967/1‑115 (MQ=255)
                                                                                              ccAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGccac                                                                             <  3:255054/120‑1 (MQ=255)
                                                                                              ccAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGccac                                                                             >  4:255054/1‑120 (MQ=255)
                                                                                              ccAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGaaaaa                                                          >  7:166237/1‑139 (MQ=255)
                                                                                              ccAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGaaaaa                                                          <  8:166237/139‑1 (MQ=255)
                                                                                              ccAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGt                                                   <  4:226437/146‑1 (MQ=255)
                                                                                              ccAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTg                                                  >  2:91844/1‑147 (MQ=255)
                                                                                               cAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAg                                               <  5:32203/149‑1 (MQ=255)
                                                                                               cAGACCTTGCCGAAGTAGGTGCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAgg                                              >  6:281053/1‑149 (MQ=255)
                                                                                                aGACCTTGCCGAAGTAGGTGCCCAACTCGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAgg                                              >  3:352538/1‑149 (MQ=255)
                                                                                                     ttGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGt                                         >  2:101044/1‑149 (MQ=255)
                                                                                                        ccGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCgg                                                                                  >  7:291749/1‑105 (MQ=255)
                                                                                                        ccGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTGTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCgg                                                                                  <  8:291749/105‑1 (MQ=255)
                                                                                                        acGAAGTAGGTGACCAAATTGGTGAACGCCGACATGCCGTTTTGCTGAACGGCATGCACTACCTTTCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCa                                       <  6:318483/147‑1 (MQ=255)
                                                                                                         cGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAgg                                     >  2:63493/1‑149 (MQ=255)
                                                                                                          gAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGc                                    <  5:334747/149‑1 (MQ=255)
                                                                                                           aaGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCg                                   >  2:292630/1‑149 (MQ=255)
                                                                                                                   gCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGa                           >  2:189381/1‑149 (MQ=255)
                                                                                                                    cccAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGAt                          <  7:294306/149‑1 (MQ=255)
                                                                                                                     ccAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATa                         <  1:91844/149‑1 (MQ=255)
                                                                                                                     ccAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATa                         <  8:353432/149‑1 (MQ=255)
                                                                                                                       aaCTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGc                       <  8:214310/149‑1 (MQ=255)
                                                                                                                        aCTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGcc                      <  8:224448/149‑1 (MQ=255)
                                                                                                                        aCTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGcc                      <  4:136717/149‑1 (MQ=255)
                                                                                                                          ttGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCtt                    <  4:89349/149‑1 (MQ=255)
                                                                                                                            ggTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCtt                    <  1:189381/147‑1 (MQ=255)
                                                                                                                            ggTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTg                   <  4:238089/148‑1 (MQ=255)
                                                                                                                              tcaaCTCCGACATGACTTTTGGCTGAACGGCATGCACAAACTTCCTAAGGTTGATTTCCACTTCACGGACGTTGACATCCCCGGACACAAGGATGTATCGCGAAAAAGCCTGGTGCAGGCCGCTAAGGCGGTCTTCGATAGCCttc                   >  6:222348/3‑145 (MQ=255)
                                                                                                                                aaCGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGTTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCg              >  1:60068/1‑149 (MQ=255)
                                                                                                                                 aCGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCgg             <  6:82848/149‑1 (MQ=255)
                                                                                                                                  cGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGc            <  5:281053/149‑1 (MQ=255)
                                                                                                                                   gCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGCt           >  6:41795/1‑149 (MQ=255)
                                                                                                                                      tgcATGCCGTTCAGCCAAACGGCATGTACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCt                                                      <  6:262799/101‑1 (MQ=38)
                                                                                                                                      tgcATGCCGTTCAGCCAAACGGCATGTACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCt                                                      >  5:262799/3‑103 (MQ=38)
                                                                                                                                      tgcATGCCGTTCAGCCAAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGtt        >  5:328767/3‑149 (MQ=255)
                                                                                                                                      gACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGtt        <  4:115320/149‑1 (MQ=255)
                                                                                                                                        cATGCCGTTCAGCCAAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGtt        >  6:179796/1‑147 (MQ=255)
                                                                                                                                        cATGCCGTTCAGCCAAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGTTCa      >  3:21349/1‑149 (MQ=255)
                                                                                                                                         aTGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAGAGCCTGGTGCAGGCCGGTCAGGc                                    >  8:97013/1‑118 (MQ=255)
                                                                                                                                         aTGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAACAGCCTGGTGCAGGCCGGTCAGGc                                    <  7:97013/118‑1 (MQ=255)
                                                                                                                                         aTGCCGTTCAGCCAAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGTTCAg     >  4:184491/1‑149 (MQ=255)
                                                                                                                                         aTGCCGTTCAGCCAAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGTTCAg     <  3:184491/149‑1 (MQ=255)
                                                                                                                                          tGCCGTTCAGCCAAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGTTCAGt    >  2:252861/1‑149 (MQ=255)
                                                                                                                                           gCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTa                 <  4:158922/135‑1 (MQ=255)
                                                                                                                                           gCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTa                 >  3:158922/1‑135 (MQ=255)
                                                                                                                                           gCCGTTTGGCTGAACGGCATGCACTACCGTGCTCAGGGTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAACGCCTGGTGCAGGCCGGTCAGGCGGACTACGATAGCCTTGTAGCGGCTGtt        >  5:282544/1‑144 (MQ=255)
                                                                                                                                           gCCGGTTGGCTGATCGGCATGCACTACCTTGCTCAGGTTGTTTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGtt        <  6:282544/144‑1 (MQ=255)
                                                                                                                                            ccGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGTTCAGTTc  >  1:98909/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                              
CGGCCATGATCAGGGTCATCATCACGGTACCGAAGATGGCCGGGAAGATACCGCCCTCGGTGTTGGCTTCACGCGGGTCGTCGCTGAGGAACTCCCAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGTTCAGTTC  >  NC_002947/6073298‑6073586

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: