Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,099,695 G→A A277A (GCC→GCT)  rhaT ← L‑rhamnose:proton symporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,099,6950GA98.4% 240.2 / ‑3.7 63A277A (GCC→GCT) rhaTL‑rhamnose:proton symporter
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (25/37);  total (25/38)
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.77e-01

AGGCTCAACACCGTTACCGGACGGCGTCCTGCATTGTTCCACTCTTTCAGCACCAGCCCGACGATACCGCCGCACAATACATAGAAACTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGTGGCCCCAGGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCGAGAGTAACACATTGTGAATGATCAGCGATTTTGCCAGCGAGAAGTCGGCTTTTAGCGACAAATCCTTCACTTTTGCCAGACGAATAAAACAGAAACCGAG  >  NC_000913/4099549‑4099841
                                                                                                                                                  |                                                                                                                                                  
aGGCTCAACACCGTTACCGGACGGCGTCCTGCATTGTTCCACTCTTTCAGCACCAGCCCGACGATACCGCCGCACAATACATAGAAACTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGTGGCCCCAAGCAt                                                                                                                                                <  2:883427/151‑1 (MQ=255)
 ggCTCAACACCGTTACCGGACGGCGTCCTGCATTGTTCCACTCTTTCAGCACCAGCCCGACGATACCGCCGCACAATACATAGAAACTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGTGGCCCCAAGCATa                                                                                                                                               <  2:418820/151‑1 (MQ=255)
  gCTCAACACCGTTACCGGACGGCGTCCTGCATTGTTCCACTCTTTCAGCACCAGCCCGACGATACCGCCGCACAATACATAGAAACTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGTGGCCCCAAGCATAg                                                                                                                                              <  1:232570/151‑1 (MQ=255)
    tCAACACCGTTACCGGACGGCGTCCTGCATTGTTCCACTCTTTCAGCACCAGCCCGACGATACCGCCGCACAATACATAGAAACTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGTGGCCCCAAGCATAGa                                                                                                                                             >  1:205908/1‑150 (MQ=255)
        cacCGTTACCGGACGGCGTCCTGCATTGTTCCACTCTTTCAGCACCAGCCCGACGATACCGCCGCACAATACATAGAAACTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGTGGCCCCAAGCATAGaaaa                                                                                                                                          >  1:844523/1‑149 (MQ=255)
         acCGTTACCGGACGGCGTCCTGCATTGTTCCACTCTTTCAGCACCAGCCCGACGATACCGCCGCACAATACATAGAAACTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGTGGCCCCAAGCATAGaaaa                                                                                                                                          >  1:50692/1‑148 (MQ=255)
              tACCGGACGGCGTCCTGCATTGTTCCACTCTTTCAGCACCAGCCCGACGATACCGCCGCACAATACATAGAAACTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGTGGCCCCAAGCATAGAAAAAGAATTg                                                                                                                                   <  1:794425/150‑1 (MQ=255)
               aCCGGACGGCGTCCTGCATTGTTCCACTCTATCAGCACCAGCCCGACGATACCGCCGCACAATACATAGAAACTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGTGGCCCCAAGCATAGAAAAAGAATTGc                                                                                                                                  <  1:290503/150‑1 (MQ=255)
                  ggACGGCGTCCTGCATTGTTCCACTCTTTCAGCACCAGCCCGACGATACCGCCGCACAATACATAGAAACTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGTGGCCCCAAGCATAGAAAAAGAATTGc                                                                                                                                  >  2:120312/1‑147 (MQ=255)
                  ggACGGCGTCCTGCATTGTTCCACTCTTTCAGCACCAGCCCGACGATACCGCCGCACAATACATAGAAACTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGTGGCCCCAAGCATAGAAAAAGAATTGCAGAt                                                                                                                              <  1:830706/151‑1 (MQ=255)
                   gACGGCGTCCTGCATTGTTCCACTCTTTCAGCACCAGCCCGACGATACCGCCGCACAATACATAGAAACTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGTGGCCCCAAGCATAGAAAAAGAATTGCAGATa                                                                                                                             <  1:326770/151‑1 (MQ=255)
                   gACGGCGTCCTGCATTGTTCCACTCTTTCAGCACCAGCCCGACGATACCGCCGCACAATACATAGAAACTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGTGGCCCCAAGCATAGAAAAAGAATTGCAGATa                                                                                                                             <  2:937174/151‑1 (MQ=255)
                       gCGTCCTGCATTGTTCCACTCTTTCAGCACCAGCCCGACGATACCGCCGCACAATACATAGAAACTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGTGGCCCCAAGCATAGAAAAAGAATTGCAGATACcac                                                                                                                         <  1:309147/151‑1 (MQ=255)
                        cGTCCTGCATTGTTCCACTCTTTCAGCACCAGCCCGACGATACCGCCGCACAATACATAGAAACTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGTGGCCCCAAGCATAGAAAAAGAATTGCAGATACcac                                                                                                                         >  2:376291/1‑150 (MQ=255)
                         gTCCTGCATTGTTCCACTCTTTCAGCACCAGCCCGACGATACCGCCGCACAATACATAGAAACTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGTGGCCCCAGGCATAGAAAAAGAATTGCAGATACCACAt                                                                                                                       <  1:811785/151‑1 (MQ=255)
                           ccTGCATTGTTCCACTCTTTCAGCACCAGCCCGACGATACCGCCGCACAATACATAGAAACTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGTGGCCCCAAGCATAGAAAAAGAATTGCAGATACCACATCa                                                                                                                     <  1:70967/151‑1 (MQ=255)
                           ccTGCATTGTTCCACTCTTTCAGCACCAGCCCGACGATACCGCCGCACAATACATAGAAACTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGTGGCCCCAAGCATAGAAAAAGAATTGCAGATACCACATCa                                                                                                                     <  1:634893/151‑1 (MQ=255)
                              gCATTGTTCCACTCTTTCAGCACCAGCCCGACGATACCGCCGCACAATACATAGAAACTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGTAATGCGGGCGTGGCCCCAAGCATAGAAAAAGAATTGCAGATACCACATCaa                                                                                                                    <  1:464677/149‑1 (MQ=255)
                              gCATTGTTCCACTCTTTCAGCACCAGCCCGACGATACCGCCGCACAATACATAGAAACTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGTGGCCCCAAGCATAGAAAAAGAATTGCAGATACCACATCAAcc                                                                                                                  <  2:612122/151‑1 (MQ=255)
                              gCATTGTTCCACTCTTTCAGCACCAGCCCGACGATACCGCCGCACAATACATAGAAACTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGTGGCCCCAAGCATAGAAAAAGAATTGCAGATACCACATCAAcc                                                                                                                  <  2:285263/151‑1 (MQ=255)
                               cATTGTTCCACTCTTTCAGCACCAGCCCGACGATACCGCCGCACAATACATAGAAACTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGTGGCCCCAAGCATAGAAAAAGAATTGCAGATACCACATCAAccc                                                                                                                 <  2:322520/151‑1 (MQ=255)
                                 ttGTTCCACTCTTTCAGCACCAGCCCGACGATACCGCCGCACAATACATAGAAACTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGTGGCCCCAAGCATAGAAAAAGAATTGCAGATACCACATCAACCCGc                                                                                                               >  2:985046/1‑151 (MQ=255)
                                    ttCCACTCTTTCAGCACCAGCCCGACGATACCGCCGCACAATACATAGAAACTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGTGGCCCCAAGCATAGAAAAAGAATTGCAGATACCACATCAACCCGccc                                                                                                             >  2:968077/1‑150 (MQ=255)
                                       cACTCTTTCAGCACCAGCCCGACGATACCGCCGCACAATACATAGAAACTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGTGGCCCCAAGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAgtg                                                                                                         <  2:204555/151‑1 (MQ=255)
                                         ctctTTCAGCACCAGCCCGACGATACCGCCGCACAATACATAGAAACTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGTGGCCCCAAGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAgtgt                                                                                                        <  2:253764/150‑1 (MQ=255)
                                         ctctTTCAGCACCAGCCCGACGATACCGCCGCACAATACATAGAAACTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGTGGCCCCAAGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTc                                                                                                       >  2:644888/1‑151 (MQ=255)
                                           ctTTCAGCACCAGCCCGACGATACCGCCGCACAATACATAGAAACTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGTGGCCCCAAGAATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCg                                                                                                      >  1:143051/1‑150 (MQ=255)
                                             ttCAGCACCAGCCCGACGATACCGCCGCACAATACATAGAAACTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGTGGCCCCAAGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCga                                                                                                     <  2:471019/149‑1 (MQ=255)
                                             ttCAGCACCAGCCCGACGATACCGCCGCACAATACATAGAAACTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGTGGCCCCAAGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCg                                                                                                      <  2:191517/148‑1 (MQ=255)
                                                 gCACCAGCCCGACGATACCGCCGCACAATACATAGAAACTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGTGGCCCCAAGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCGAGAGTaa                                                                                               <  2:611496/151‑1 (MQ=255)
                                                         ccGACGATACCGCCGCACAATACATAGAAACTCATAGGCCGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGTGGCCCCAAGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCGAGAGTAACACATtgt                                                                                       <  1:492702/151‑1 (MQ=255)
                                                              gaTACCGCCGCACAATACATAGAAACTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGTGGCCCCAAGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCGAGAGTAACACATTGTGAATg                                                                                  >  1:470304/1‑151 (MQ=255)
                                                               aTACCGCCGCACAATACATAGAAACTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGTGGCCCCAAGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCGAGAGTAACACATTGTGAATg                                                                                  >  1:789930/1‑150 (MQ=255)
                                                                  ccgccgCACAATACATAGAAACTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGTGGCCCCAAGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCGAGAGTAACACATTGTGAATGATCa                                                                              <  2:334572/151‑1 (MQ=255)
                                                                           aaTACATAGAAACTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGTGGCCCCAAGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCGAGAGTAACACATTGTGAATGATCAGCGATTTTg                                                                     <  1:432342/151‑1 (MQ=255)
                                                                           aaTACATAGAAACGCATATGCAGCATCCAACTGATGTAGTCATATTGCGCCGGAATGCGTGCGTGGCCCCAAGCATAGAAAAAGAATTGCAAATACCACATCAACCCGCCCAGTGTCGAGAGTAACACATTGTGAATGATCAGCGATTTTg                                                                     >  1:570467/1‑151 (MQ=255)
                                                                             tACATAGAAACTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGTGGCCCCAAGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCGAGAGTAACACATTGTGAATGATCAGCGATTTTGcc                                                                   >  2:309683/1‑151 (MQ=255)
                                                                                  aGAAACTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGTGGCCCCAAGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCGAGAGTAACACATTGTGAATGATCAGCGATTTTGCCAGCg                                                               >  2:20073/1‑150 (MQ=255)
                                                                                     aaCTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGTGGCCCCAAGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCGAGAGTAACACATTGTGAATGATCAGCGATTTTGCCAGCgaga                                                            >  1:199383/1‑150 (MQ=255)
                                                                                      aCTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGTGGCCCCAAGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCGAGAGTAACACATTGTGAATGATCAGCGATTTTGCCAGCGAGAAg                                                          <  2:11255/151‑1 (MQ=255)
                                                                                      aCTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGTGGCCCCAAGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCGAGAGTAACACATTGTGAATGATCAGCGATTTTGCCAGCGAGAAg                                                          <  2:757625/151‑1 (MQ=255)
                                                                                               cagcaTCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGTGGCCCCAAGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCGAGAGTAACACATTGTGAATGATCAGCGATTTTGCCAGCGAGAAGTCGGCtttt                                                 >  2:635552/1‑151 (MQ=255)
                                                                                                 gcaTCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGTGGCCCCAAGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCGAGAGTAACACATTGTGAATGATCAGCGATTTTGCCAGCGAGAAGTCGGCTTTTAg                                               <  1:819654/151‑1 (MQ=255)
                                                                                                     ccAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGTGGCCCCAAGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCGAGAGTAACACATTGTGAATGATCAGCGATTTTGCCAGCGAGAAGTCGGCTTTTAGCGa                                            >  1:691403/1‑150 (MQ=255)
                                                                                                      cAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGTGGCCCCAAGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCGAGAGTAACACATTGTGAATGATCAGCGATTTTGCCAGCGAGAAGTCGGCTTTTAGCGACa                                          <  2:669011/151‑1 (MQ=255)
                                                                                                           gATGTAGTCATACTGCGCCGGAATGCGGGCGTGGCCCCAAGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCGAGAGTAACACATTGTGAATGATCAGCGATTTTGCCAGCGAGAAGTCGGCTTTTAGCGACAAATcc                                     >  2:306575/1‑151 (MQ=255)
                                                                                                            aTGTAGTCATACTGCGCCGGAATGCGGGCGTGGCCCCAAGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCGAGAGTAACACATTGTGAATGATCAGCGATTTTGCCAGCGAGAAGTCGGCTTTTAGCGACAAATCCt                                    <  1:497302/151‑1 (MQ=255)
                                                                                                             tGTAGTCATACTGCGCCGGAATGCGGGCGTGGCCCCAAGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCGAGAGTAACACATTGTGAATGATCAGCGATTTTGCCAGCGAGAAGTCGGCTTTTAGCGACAAATCCtt                                   >  2:389170/1‑151 (MQ=255)
                                                                                                              gTAGTCATACTGCGCCGGAATGCGGGCGTGGCCCCAAGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCGAGAGTAACACATTGTGAATGATCAGCGATTTTGCCAGCGAGAAGTCGGCTTTTAGCGACAAATCCTTc                                  >  2:60382/1‑151 (MQ=255)
                                                                                                               tAGTCATACTGCGCCGGAATGCGGGCGTGGCCCCAAGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCGAGAGTAACACATTGTGAATGATCAGCGATTTTGCCAGCGAGAAGTCGGCTTTTAGCGACAAATCCtgc                                  <  1:570402/150‑3 (MQ=255)
                                                                                                                aGTCATACTGCGCCGGAATGCGGGCGTGGCCCCAAGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCGAGAGTAACACATTGTGAATGATCAGCGATTTTGCCAGCGAGAAGTCGGCTTTTAGCGACAAATCCTTCAc                                <  2:16582/151‑1 (MQ=255)
                                                                                                                   cATACTGCGCCGGAATGCGGGCGTGGCCCCAAGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCGAGAGTAACACATTGTGAATGATCAGCGATTTTGCCAGCGAGAAGTCGGCTTTTAGCGACAAATCCTTCACtt                              <  2:461074/150‑1 (MQ=255)
                                                                                                                         gcgcCGGAATGCGGGCGTGGCCCCAAGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCGAGAGTAACACATTGTGAATGATCAGCGATTTTGCCAGCGAGAAGTCGGCTTTTAGCGACAAATCCTTCTCTTTTGCCa                        <  1:659520/150‑1 (MQ=255)
                                                                                                                          cgcCGGAATGCGGGCGTGGCCCCAAGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCGAGAGTAACACATTGTGAATGATCAGCGATTTTGCCAGCGAGAAGTCGGCTTTTAGCGACAAATCCTTCACTTTTGCCAGa                      <  2:385573/151‑1 (MQ=255)
                                                                                                                              ggAATGCGGGCGTGGCCCCAAGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCGAGAGTAACACATTGTGAATGATCAGCGATTTTGCCAGCGAGAAGTCGGCTTTTAGCGACAAATCCTTCACTTTTGCCAGACGa                   >  1:868473/1‑150 (MQ=255)
                                                                                                                                aaTGCGGGCGTGGCCCCAAGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCGAGAGTAACACATTGTGAATGATCAGCGATTTTGCCAGCGAGAAGTCGGCTTTTAGCGACAAATCCTTCACTTTTGCCAGACGAAt                 >  1:207406/1‑150 (MQ=255)
                                                                                                                                aaTGCGGGCGTGGCCCCAAGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCGAGAGTAACACATTGTGAATGATCAGCGATTTTGCCAGCGAGAAGTCGGCTTTTAGCGACAAATCCTTCACTTTTGCCAGACGAATa                <  1:667422/151‑1 (MQ=255)
                                                                                                                                   gCGGGCGTGGCCCCAAGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCGAGAGTAACACATTGTGAATGATCAGCGATTTTGCCAGCGAGAAGTCGGCTTTTAGCGACAAATCCTTCACTTTTGCCAGACGAATaaaa             <  2:941132/151‑1 (MQ=255)
                                                                                                                                   gCGGGCGTGGCCCCAAGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCGAGAGTAACACATTGTGAATGATCAGCGATTTTGCCAGCGAGAAGTCGGCTTTTAGCGACAAATCCTTCACTTTTGCCAGACGAATaaaa             <  1:221636/151‑1 (MQ=255)
                                                                                                                                       gCGTGGCCCCAAGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCGAGAGTAACACATTGTGAATGATCAGCGATTTTGCCAGCGAGAAGTCGGCTTTTAGCGACAAATCCTTCACTTTTGCCAGACGAATAAAAc            >  2:620433/1‑148 (MQ=255)
                                                                                                                                        cGTGGCCCCAAGCATAGAAAAAGAATTGCAGGTACCACATCAACCCGCCCAGTGTCGAGAGTAACACATTGTGAATGATCAGCGATTTTGCCAGCGAGAAGTCGGCTTTTAGCGACAAATCCTTCACTTTTGCCAGACGAATAAAACAGa         >  1:247291/1‑150 (MQ=255)
                                                                                                                                          tGGCCCCAAGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCGAGAGTAACACATTGTGAATGATCAGCGATTTTGCCAGCGAGAAGTCGGCTTTTAGCGACAAATCCTTCACTTTTGCCAGACGAATAAAACAGAAAc      >  1:688690/1‑151 (MQ=255)
                                                                                                                                              cccAAGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCGAGAGTAACACATTGTGAATGATCAGCGATTTTGCCAGCGAGAAGTCGGCTTTTAGCGACAAATCCTTCACTTTTGCCAGACGAATAAAACAGAAACCGAg  <  2:962518/151‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
AGGCTCAACACCGTTACCGGACGGCGTCCTGCATTGTTCCACTCTTTCAGCACCAGCCCGACGATACCGCCGCACAATACATAGAAACTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGTGGCCCCAGGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCGAGAGTAACACATTGTGAATGATCAGCGATTTTGCCAGCGAGAAGTCGGCTTTTAGCGACAAATCCTTCACTTTTGCCAGACGAATAAAACAGAAACCGAG  >  NC_000913/4099549‑4099841

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑