Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 425,115 G→A 6.1% S414S (AGC→AGT)  U712_01935 ← putative ABC transporter permease yclI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP006881425,1150GA 6.1% 202.1 / 2.4 66S414S (AGC→AGT) U712_01935putative ABC transporter permease yclI
Reads supporting (aligned to +/- strand):  ref base G (35/25);  new base A (2/2);  total (39/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.71e-01

GGATGCCGATCCCGCCAAGTATAAGCATGTCGTTCATAGAAACCGCGACGTTCAGTGAATCGATAACGTCCACGTTGGATGAACTGTGGCCGAACATTTTGCCGCCCATGCCGCCACCGCCGCCGGGCATTTGTCCGCTTGCGGTTTGCGTTGAATCGGTTGATGAACTGATTTGCTGGGATAAAAGCTGATTTCCAAGCTGATTCGCGACGAGATTTCCCGTCACAGATGCAAGGCCGATAGCGATAACGGCCACAATCAGAATCTCAGTCA  >  CP006881/424979‑425251
                                                                                                                                        |                                                                                                                                        
ggATGCCGATCCCGCCAAGTATAAGCATGTCGTTCATAGAAACCGCGACGTTCAGTGAATCGATAACGTCCACGTTGGATGAACTGTGGCCGAACATTTTGCCGCCCATGCCGCCACCGCCGCCGGGCATTTGTCCACt                                                                                                                                        >  6:181552/1‑139 (MQ=255)
     ccGATCCCGCCAAGTATAAGCATGTCGTTCATAGAAACCGCGACGTTCAGTGAATCGATAACGTCCACGTTGGATGAACTGTGGCCGAACATTTTGCCGCCCATGCCGCCACCGCCGCCGGGCATTTGTCCGCTTGCgg                                                                                                                                   >  8:158436/1‑139 (MQ=255)
     ccGATCCCGCCAAGTATAAGCATGTCGTTCATAGAAACCGCGACGTTCAGTGAATCGATAACGTCCACGTTGGATGAACTGTGGCCGAACATTTTGCCGCCCATGCCGCCACCGCCGCCGGGCATTTGTCCGCTTGCgg                                                                                                                                   >  4:404981/1‑139 (MQ=255)
      cGATCCCGCCAAGTATAAGCATGTCGTTCATAGAAACCGCGACGTTCAGTGAATCGATAACGTCCACGTTGGATGAACTGTGGCCGAACATTTTGCCGCCCATGCCGCCACCGCCGCCGGGCATTTGTCCACTTGCGGt                                                                                                                                  <  5:181552/139‑1 (MQ=255)
                  gTATAAGCATGTCGTTCATAGAAACCGCGACGTTCAGTGAATCGATAACGTCCACGTTGGATGAACTGTGGCCGAACATTTTGCCGCCCATGCCGCCACCGCCGCCGGGCATTTGTCCGCTTGCGGTTTGCGTTGAATc                                                                                                                      <  3:404981/139‑1 (MQ=255)
                       aGCATGTCGTTCATAGAAACCGCGACGTTCAGTGAATCGATAACGTCCACGTTGGATGCACTGTGGCCGAACATTTTGCAGCCCATGCCGCCACCGCCGCCGGGAATTTGTCCTCTTGCGGTTTGCGGTGAATCGGTtg                                                                                                                 >  8:185103/1‑139 (MQ=255)
                       aGCATGTCGTTCATAGAAACCGCGACGTTCAGTGAATCGATAACGTCCACGTTGGATGAACTGTGGCCGAACATTTTGCCGCCCATGCCGCCACCGCCGCCGGGCATTTGTCCGCTTGCGGTTTGCGTTGAATCGGtt                                                                                                                  >  2:215485/1‑138 (MQ=255)
                       aGCATGTCGTTCATAGAAACCGCGACGTTCAGTGAATCGATAACGTCCACGTTGGATGAACTGTGGCCGAACATTTTGCCGCCCATGCCGCCACCGCCGCCGGGCATTTGTCCGCTTGCGGTTTGCGTTGAATCGGtt                                                                                                                  >  8:228663/1‑138 (MQ=255)
                       aGCATGTCGTTCATAGAAACCGCGACGTTCAGTGAATCGATAACGTCCACGTTGGATGAACTGTGGCCGAACATTTTGCCGCCCATGCCGCCACCGCCGCCGGGCATTTGTCCGCTTGCGGTTTGCGTTGAATCGGTtg                                                                                                                 >  2:347025/1‑139 (MQ=255)
                             tCGTTCATAGAAACCGCGACGTTCAGTGAATCGTTAACGTCCACGTTGGATGAACTGTGGCCGAACATTTTGCCGCCCATGCCGCCACCGCCGCCGGGCATTTGTCCGCTTGCGGTTTGCGTTGGATCGGTTGATGAAc                                                                                                           >  7:276573/1‑139 (MQ=255)
                             tCGTTCATAGAAACCGCGACGTTCAGTGAATCGATAACGTCCCCGGTGGATGAACTGTGGCAGAACATTTTGCCGCCCATGCCGCCACCGCCGCCGGGCATTTGTCCGCTTGCGGTTTGCGTTGAATCGGTTGATGAAc                                                                                                           >  7:148003/1‑139 (MQ=255)
                             tCGTTCATAGAAACCGCGACGTTCAGTGAATCGATAACGTCCACGTTGGATGCACTGTGGCCGAACATTTTGCCGCCCATGCCGCCACCGCCGCCGGGCATTTGTCCGCTTGCGGTTTGCGTTGAATCGGTTGATGAAc                                                                                                           >  6:305925/1‑139 (MQ=255)
                             tCGTTCATAGAAACCGCGACGTTCAGTGAATCGATAACGTCCACGTTGGATGAACTTTGGCCGAACATTTTGCCGCCCATGCCGCCACCGCCGCCGGGCATTTGTCCGCTTTCGGTTTGCGTTGAATCGGTTGATGAAc                                                                                                           >  3:370539/1‑139 (MQ=255)
                             tCGTTCATAGAAACCGCGACGTTCAGTGAATCGATAACGTCCACGTTGGATGAACTGTGGCCGAACATTTTGCCGCCCATGCCGCCCCCGCCGCCGGGCCTTTGTCCGCGTGCGGGTGGCGTTGAATCGGGTGATGAAc                                                                                                           >  2:188961/1‑139 (MQ=255)
                             tCGTTCATAGAAACCGCGACGTTCAGTGAATCGATAACGTCCACGTTGGATGAACTGTGGCCGAACATTTTGCCGCCCATGCCGCCACCGCCGCCGGGCATTTGTCTGCTTGCGGTTTGCGTTGAATCGGTTGATGAAc                                                                                                           >  1:106908/1‑139 (MQ=255)
                             tCGTTCATAGAAACCGCGACGTTCAGTGAATCGATAACGTCCACGTTGGATGAACTGTGGCCGAACATTTTGCCGCCCATGCCGCCACCGCCGCCGGGCATTTGTCGGCTTGCGGTTTGCGTTGAATCGGTTGATGAAc                                                                                                           >  5:122904/1‑139 (MQ=255)
                             tCGTTCATAGAAACCGCGACGTTCAGTGAATCGATAACGTCCACGTTGGATGAACTGTGGCCGAACATTTTGCCGCCCATGCCGCCACCGCCGCCGGGCATTTGTCCGCTTGCGGTTTGCGTTGAATCGGTTGATGAAc                                                                                                           >  5:379416/1‑139 (MQ=255)
                             tCGTTCATAGAAACCGCGACGTTCAGTGAATCGATAACGTCCACGTTGGATGAACTGTGGCCGAACATTTTGCCGCCCATGCCGCCACCGCCGCCGGGCATTTGTCCGCTTGCGGTTTGCGTTGAATCGGTTGATGAAc                                                                                                           >  8:158850/1‑139 (MQ=255)
                             tCGTTCATAGAAACCGCGACGTTCAGTGAATCGATAACGTCCACGTTGGATGAACTGTGGCCGAACATTTTGCCGCCCATGCCGCCACCGCCGCCGGGCATTTGTCCGCTTGCGGTTTGCGTTGAATCGGTTGATGAAc                                                                                                           >  8:65689/1‑139 (MQ=255)
                             tCGTTCATAGAAACCGCGACGTTCAGTGAATCGATAACGTCCACGTTGGATGAACTGTGGCCGAACATTTTGCCGCCCATGCCGCCACCGCCGCAGGGCATTTGTCCGCTTGCGGTTTGCGTTGAATCGGTTGATGAAc                                                                                                           >  6:292296/1‑139 (MQ=255)
                                    tAGAAACCGCGACGTTCAGTGAATCGATAACGTCCACGTTGGATGAACTGTGGCCGAACATTTTGCCGCCCATGCCGCCACCGCCGCCGGGCATTTGTCCGCTTGCGGTTTGCGTTGAATCGGTTGATGAACTGATTTg                                                                                                    >  2:216157/1‑139 (MQ=255)
                                      gAAACCGCGACGTTCAGTGAATCGATAACGTCCACGTTGGATGAACTGTGGCCGAACATTTTGCCGCCCATGCCGCCACCGCCGCGGGGCATTTGTCCGCTTGCGGTTTGCGTTGAATCGGTTGATGAACTGATTTGCt                                                                                                  <  5:292296/139‑1 (MQ=255)
                                      gAAACCGCGACGTTCAGTGAATCGATAACGTCCACGTTGGATGAACTGTGGCCGAACATTTTGCCGCCCATGCCGCCACCGCCGCCGGGCATTTGTCCGCTTGCGGTTTGCGTTGAATCGGTTGATGAACTGATTTGCt                                                                                                  <  7:228663/139‑1 (MQ=255)
                                      gAAACCGCGACGTTCAGTGAATCGATAACGTCCACGTTGGATGAACTGTGGCCGAACATTTTGCCGCCCATGCCGCCACCGCCGCCGGGCATTTGTCCGCTTGCGGTTTGCGTTGAATCGGTTGATGAACTGATTTGCt                                                                                                  <  7:111543/139‑1 (MQ=255)
                                         aCCGCGACGTTCAGTGAATCGATAACGTCCACGTTGGATGAACTGTGGCCGAACATTTTGCCGCCCATGCCGCCACCGCCGCCGGGCATTTGTCCGCTTGCGGt                                                                                                                                  >  8:314684/1‑104 (MQ=255)
                                         aCCGCGACGTTCAGTGAATCGATAACGTCCACGTTGGATGAACTGTGGCCGAACATTTTGCCGCCCATGCCGCCACCGCCGCCGGGCATTTGTCCGCTTGCGGt                                                                                                                                  <  7:314684/104‑1 (MQ=255)
                                          cctcgACGTTCAGTGAATCGATAACGTCCACGTTGGATGAACTGTGGCCGAACATTTTGCCGCCCATGACGCCACCGCCGCCGGGCATTTGTCCGCTTGCGGTTTGCGTTGAATCGGTTGATGAACTGATTTGCTGGGa                                                                                              >  6:23393/4‑139 (MQ=255)
                                          ccGCGACGTTCAGTGAATCGATAACGTCCACGTTGGATGAACTGTGGCCGAACATTTTGCCGCCCATGCCGCCACCGCCGCCGGGCATTTGTCCGCTTGCGGTTTGCGTTGAATCGGTTGATGAACTGATTTGCt                                                                                                  <  5:286146/135‑1 (MQ=255)
                                          ccGCGACGTTCAGTGAATCGATAACGTCCACGTTGGATGAACTGTGGCCGAACATTTTGCCGCCCATGCCGCCACCGCCGCCGGGCATTTGTCCGCTTGCGGTTTGCGTTGAATCGGTTGATGAACTGATTTGCt                                                                                                  >  6:286146/1‑135 (MQ=255)
                                             cgACGTTCAGTGAATCGATAACGTCCACGGTGGATGAACTGTGGCCGAACATTTTGCCGCCCATGCCGCCACCGCCGCCGGGCATTTGTCCGCTTGCGGTTTGCGTTGAATCGGTTGATGAACTGATTTGCTGGGATaa                                                                                           <  8:152586/139‑1 (MQ=255)
                                                  ttCAGTGAATCGATAACGTCCACGTTGGATGAACTGTGGCCGAACATTTTGCCGCCCATGCCGCCACCGCCGCCGGGCATTTGTCCGCTTGCGGTTTGCGTTGAATCGGTTGATGAACTGATTTGCTGGGATAAAAGCt                                                                                      <  5:314982/139‑1 (MQ=255)
                                                          aTCGATAACGTCCACGTTGGATGAACTGTGGCCGAACATTTTGCCGCCCATGCCGCCACCGCCGCCGGGCATTTGTCCGCTTGCGGTTTGCGTTGAATCGGTTGATGAACTGATTTGCTGGGATAAAAGCTGATTTcc                                                                               >  4:69783/1‑138 (MQ=255)
                                                                         gTTGGATGAACTGTGGCCGAACATTTTGCCGCCCATGCCGCCACCGCCGCCGGGCATTTGTCCGCTTGCGGTTTGCGTTGAATCGGTTGATGAACTGATTTGCTGGGATAAAAGCTGATTTCCAAGCTGATTCGcgac                                                                >  4:168369/1‑138 (MQ=255)
                                                                            ggATGAACTGTGGCCGAACATTTTGCCGCCCATGCCGCCACCGCCGCCGGGCATTTGTCCGCTTGCGGTTTGCGTTGAATCGGTTGATGAACTGATTTGCTGGGATAAAAGCTGATTTCCAAGCTGATTCGcgac                                                                >  7:17789/1‑135 (MQ=255)
                                                                             gATGAACTGTGGCCGAACATTTTGCCGCCCATGCCGCCACCGCCGCCGGGCATTTGTCCGCTTGCGGTTTGCGTTGAATCGGTTGATGAACTGATTTGCTGGGATAAAAGCTGATTTCCAAGCTGATTCGCGACGAGAt                                                           >  3:334922/1‑139 (MQ=255)
                                                                             gATGAACTGTGGCCGAACATTTTGAAGCCAATGCCGCCACCGCCGCCGGGCATTTGTCCGCTTGCGGTTTGCGTTGAATCGGTTGATGAACTGATTTGCTGGGATAAAAGCTGATTTCCAAGCTGATTCGCGACGAGAt                                                           >  3:341645/1‑139 (MQ=255)
                                                                              aTGAACTGTGGCCGAACATTTTGCCGCCCATGCCGCCACCGCCGCCGGGCATTTGTCCGCTTGCGGTTTGCGTTGAATCGGTTGATGAACTGATTTGCTGGGATAAAAGCTGATTTCCAAGCTGATTCGCGACGAGAtt                                                          >  8:310299/1‑139 (MQ=255)
                                                                              aTGAACTGTGGCCGAACATTTTGCCGCCCATGCCGCCACCGCCGCCGGGCATTTGTCCGCATGCGGTTTGCGTGGAATCGGGTGATGAACTGATTTGCTGGGATAAAAGCTGATTTCCAAGCTCATTCGCGACgagctg                                                          >  7:285585/1‑136 (MQ=255)
                                                                                   cTGTGGCCGAACATTTTGCCGCCCATGCCGCCACCGCCGCCGGGCATTTGACCTATTGCTGTTTGCGTTGAATCGGTTGATGAACTGATTTGCTGGGATAAAAGCTGATTTCCAAGCTGATTCGCGACGAGATTTccag                                                     >  3:358950/1‑137 (MQ=255)
                                                                                       ggCCGAACATTTTGCCGCCCATGCCGCCACCGCCGCCGGGCATTTGTCCGCTTGCGGTTTGCGTTGAATCGGTTGATGAACTGATTTGCTGGGATAAAAGCTGATTTCCAAGCTGATTCGCGACGAGATTTCCCGTcac                                                 <  8:17789/139‑1 (MQ=255)
                                                                                       ggCCGAACATTTTGCCGCCCATGCCGCCACCGCCGCCGGGCATTTGTCCGCTTGCGGTTTGCGTTGAATCGGTTGATGAACTGATTTGCTGGGATAAAAGCTGATTTCCAAGCTGATTCGCGACGAGATTTCCCGTcac                                                 <  3:55166/139‑1 (MQ=255)
                                                                                       ggCCGAACATTTTGCCGCCCATGCCGCCACCGCCGCCGGGCATTTGTCCGCTTGCGGTTTGCGTTGAATCGGTTGATGAACTGATTTGCTGGGATAAAAGCTGATTTCCAAGCTGATTCGCGACGAGATTTCCCGTcac                                                 <  7:158436/139‑1 (MQ=255)
                                                                                       ggCCGAACATTTTGCCGCCCATGCCGCCACCGCCGCCGGGCATTTGTCCGCTTGCGGTTTGCGTTGAATCGGTTGATGAACTGATTTGCTGGGATAAAAGCTGATTTCCAAGCTGATTCGCGACGAGATTTCCCGTca                                                  >  3:338133/1‑138 (MQ=255)
                                                                                            aaCATTTTGCCGCCCATGCCGCCACCGCCGCCGGGCATTTGTCCGCtt                                                                                                                                       <  5:159400/48‑1 (MQ=255)
                                                                                            aaCATTTTGCCGCCCATGCCGCCACCGCCGCCGGGCATTTGTCCGCtt                                                                                                                                       >  6:159400/1‑48 (MQ=255)
                                                                                            aaCATTTTGCCGCCCATGCCGCCACCGCCGCCGGGCATTTGTCCGCTTGCGGTTTGCGTTGAATCGGTTGATGAACTGATTTGCTGGGATAAAAGCTGATTTCCAAGCTGATTCGCGACGAGATTTCCCGTCACAGATg                                            >  1:276008/1‑139 (MQ=255)
                                                                                            aaCATTTTGCCGCCCATGCCGCCACCGCCGCCGGGCATTTGTCCGCTTGCGGTTTGCGTTGAATCGGTTGATGAACTGATTTGCTGGGATAAAAGCTGATTTCCAAGCTGATTCGCGACGAGATTTCCCGTCACAGATg                                            >  7:289312/1‑139 (MQ=255)
                                                                                            aaCATTTTGCCGCCCATGCCGCCACCCCCGCCGGCCATTTGTCCGCTAGCGGTTTGCGTTGAATCGGTAGATGAACTCATTTGCAGGGATAAAAGCTGATTTCCAAGCTGATTCGCGACGAGATTTCCCGTCACAGATg                                            >  7:262864/1‑139 (MQ=255)
                                                                                             aCATTTTGCCGCCCATGCCGCCACCGCCGCCGGGCATTTGTCCGCTTGCGGTTTGCGTTGAATCGGTTGATGAACTGATTTGCTGGGATAAAAGCTGATTTCCAAGCTGATTCGCGACGAGATTTCCCGTCACAGATGc                                           <  3:6719/139‑1 (MQ=255)
                                                                                                ttttGCCGCCCATGCCGCCACCGCCGCCGGGCATTTGTCCGCTTGCGGTTTGCGTTGAATCGGTTGATGAACTGATTTGCTGGGATAAAAGCTGATTTCCAAGCTGATTCGCGACGAGATTTCCCGTCACAGATGc                                           >  8:195960/1‑136 (MQ=255)
                                                                                                ttttGCCGCCCATGCCGCCACCGCCGCCGGGCATTTGTCCGCTTGCGGTTTGCGTTGAATCGGTTGATGAACTGATTTGCTGGGATAAAAGCTGATTTCCAAGCTGATTCGCGACGAGATTTCCCGTCACAGATGc                                           <  7:195960/136‑1 (MQ=255)
                                                                                                                   aCCGCCGCCGGGCATTTGTCCACTTGCGGTTTGCGTTGAATCGGTTGATGAACTGATTTGCTGGGATAAAAGCTGATTTCCAAGCTGATTCGCGACGAGATTTCCCGTCACAGATGCAAGGCCGATAGCGATAACGGc                      >  4:271014/1‑138 (MQ=255)
                                                                                                                    ccgccgccgGGCATTTGTCCACTTGCGGTTTGCGTTGAATCGGTTGATGAACTGATTTGCTGGGATAAAAGCTGATTTCCAAGCTGATTCGCGACGAGATTTCCCGTCACAGATGCAAGGCCGATAGCGATAACGGCca                    <  3:271014/139‑1 (MQ=255)
                                                                                                                       ccgccgGGCATTTGTCCGCTTGCGGTTTGCGTTGAATCGGTTGATGAACTGATTTGCTGGGATAAAAGCTGATTTCCAAGCTGATTCGCGCCGAGATTTCCAGTCACAGATGCAAGGCCGATAGCGATATCGGCcacag                 >  1:205310/1‑138 (MQ=255)
                                                                                                                       ccgccgGGCATTTGTCCGCTTGCGGTTTGCGTTGAATCGGTTGATGAACTGATGTGCTGGGATAAAAGCTGATTTCCAAGCTGATTCGCGACGAGATTTCCCGTCACAGATGCAAGGCCGATAGCGATGACGGCCACaa                 >  4:327278/1‑139 (MQ=255)
                                                                                                                                aTTTGTCCGGTTGGGGTTTGCGTTGATTCGGTTGATGAACTGATTTGCTGGGATAAAAGCTGATTTCCAAGCTGATTCGCGACGAGATTTCCCGTCACAGATGCAAGGCCGATAGCGATAACGGCCACAATCAGAAtct        <  7:160501/139‑1 (MQ=255)
                                                                                                                                aTTTGTCCGCTTGCGTTTTGCGTTGTATCGGTTTATTAACTTATTTGCTTCGAAAAAAGCTGATTACCAAGCTGAATCGCGACGAGATTTCCCGTCACAGATGCAAGGCCGATAGCGATAACGGCCACAATCAGAAtat        <  6:287472/139‑3 (MQ=255)
                                                                                                                                aTTTGTCCGCTTGCGGTTTGCGTTGAATCGGTTGATGAACTGATTTTCTGGGATAAAAGCTGATTTCCAAGCTGATTCGCGACGAGATTTCCCGTCACAGATGCAAGGCCGATAGCGATAACGGCCACAATCAGAAtct        <  3:168369/139‑1 (MQ=255)
                                                                                                                                aTTTGTCCGCTTGCGGTTTGCGTTGAATCGGTTGATGAACTGATTTGCTGGGATAAAAGCTGATTTCCAAGCTGATTCGCGACGAGATTTCCCGTCACAGATGCAAGGCCGATAGCGATAACGGCCACAATCAGAAtct        <  7:310299/139‑1 (MQ=255)
                                                                                                                                aTTTGTCCGCTTGCGGTTTGCGTTGAATCGGTTGATGAACTGATTTGCTGGGATAAAAGCTGATTTCCAAGCTGATTCGCGACGAGATTTCCCGTCACAGATGCAAGGCCGATAGCGATAACGGCCACAATCAGAAtct        <  7:65689/139‑1 (MQ=255)
                                                                                                                                aTTTGTCCGCTTGCGGTTTGCGTTGAATCGGTTGATGAACTGATTTGCTGGGATAAAAGCTGATTTCCAAGCTGATTCGCGACGAGATTTCCCGTCACAGATGCAAGGCCGATAGCGATAACGGCCACAATCAGAAtct        <  5:273621/139‑1 (MQ=255)
                                                                                                                                aTTTGTCCGCTTGCGGTTTGCGTTGAATCGGTTGATGAACTGATTTGCTGGGATAAAAGCTGATTTCCAAGCTGATTCGCGACGAGATTTCCCGTCACAGATGCAAGGCCGATAGCGATAACGGCCACAATCAGAAtct        <  2:257585/139‑1 (MQ=255)
                                                                                                                                aTTTGTCCGCTTGCGGTTTGCGTTGAATCGGTTGATGAACTGATTTGCTGGGATAAAAGCTGATTTCCAAGCTGATTCGCGACGAGATTTCCCGTCACAGATGCAAGGCCGATAGCGATAACGGCCACAATCAGAAtct        <  1:216157/139‑1 (MQ=255)
                                                                                                                                 tttGTCCGCTTGCGGTTTGCGTTGAATCGGTTGATGAACTGATTTGCTGGGATAAAAGCTGATTTCCAAGCTGATTCGCGACGAGATTTCCCGTCACAGATGCAAGGCCGATAGCGATAACGGCCACAATCAGAAtct        <  7:158850/138‑1 (MQ=255)
                                                                                                                                      ccGTTTGCGGTTTGCGTTGAATGGGTTGATGAACTGATTTGCTGGCATAAACTCTGATTTCCAAGCTGATTCGCGACGAGATTTCCCGTCACAGATTCTAGGCCGATAGCGATAACGGCCACAATCAGAATCTCAGTc   <  4:358950/138‑1 (MQ=255)
                                                                                                                                      ccGCTTGCGGTTTGCGTTGAATCGGTTGATGAACTGATTTGCTGGGATAAAAGCTGATTTCCAAGCTGATTCGCGACGAGATTTCCCGTCACAGATGCAAGGCCGATAGCGATAACGGCCACAATCAGAATCTCAGTCa  <  5:380889/139‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                        
GGATGCCGATCCCGCCAAGTATAAGCATGTCGTTCATAGAAACCGCGACGTTCAGTGAATCGATAACGTCCACGTTGGATGAACTGTGGCCGAACATTTTGCCGCCCATGCCGCCACCGCCGCCGGGCATTTGTCCGCTTGCGGTTTGCGTTGAATCGGTTGATGAACTGATTTGCTGGGATAAAAGCTGATTTCCAAGCTGATTCGCGACGAGATTTCCCGTCACAGATGCAAGGCCGATAGCGATAACGGCCACAATCAGAATCTCAGTCA  >  CP006881/424979‑425251

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