Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,826,195 C→T 8.7% R568H (CGC→CAC)  yicI ← putative alpha‑glucosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,826,1950CT 8.7% 138.5 / 3.0 48R568H (CGC→CAC) yicIputative alpha‑glucosidase
Reads supporting (aligned to +/- strand):  ref base C (22/21);  new base T (2/2);  total (25/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01

TTTATGCCAGCGACTACCGTCGAGTTCATCGTTGTGCCACAGGTGTGTCCAGCGACCTTCCGGCAGGTAGAACTGCACATCGCCCGCTTCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAACGCGCACCAGCGTTTGTAAACGTGCGCCGGAGCGGTATTTTCAAAGCCGCCGATATCGTGGCTCCAGAAGCCAAAACCTGAAAGGCCAATAGACAAACCACCGCGCAGGCTTTCCGCCATTGATTCGTAGTTAGCGTAACAATCGCCAC  >  CP009273/3825945‑3826459
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                        
tttATGCCAGTGAATACCGTCGGGTTCATCGTTCTGACACAGGTGTGTCCAGCGACCTTCCGGCAGGTAGAACTGCACATCGCCCGCTTCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGaagaag                                                                                                                                                                                                                                                   <  1:274372/274‑1 (MQ=255)
   aTGCCTGCTACTCCCTTCGTGTTCATCGTTGTTCCACAGTTGTGTCCATCGACCTTCCGGCTGGTAGAACTGCACATCGCCCGCTTCTGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGGCGGTCAAGGTAATCACAATCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGAGACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGca                                                                                                                                                                                                                                               <  1:170580/275‑1 (MQ=255)
     gCCAGCGACTACCGTCGAGTTCATCGTTGTGCCACAGGTGTGTCCAGCGACCTTCCGGCAGGTAGAACTGCACATCGCCCGCTTCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGTGGCATTTCAGTTGCGTGAAGAAGCGcac                                                                                                                                                                                                                                              >  2:171959/1‑274 (MQ=255)
          cGACTACCGTCGAGTTCATCGTTGTGCCACAGGTGTGTCCAGCGACCTTCCGGCAGGTAGAACTGCACATCGCCCGCTTCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATAGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCCCGAAACACATAACGCCGCCAGGGGCATTTCACTTGCGTGacgaagacaaccagcgc                                                                                                                                                                                                                                        >  2:170575/1‑264 (MQ=255)
                                  tgCCACAGGTGTGTCCAGCGACCTTCCGGCAGGTAGAACTGCACATCGCCCGCTTCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGTGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCa                                                                                                                                                                                                                <  1:171965/275‑1 (MQ=255)
                                                ccAGCGACCTTCCTGCAGGTANAACTGCACATCGCCCGCTTCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCgg                                                                                                                                                                                                                                                                        <  1:391405/205‑1 (MQ=255)
                                                                                   ccGCGTCAGTGAACGCCGGAGAACCCATCACGTTGTCGACTAACATGTATTGTCGGTCAAGGTAATCACAAGCAGGATCGTCCGGGAACTCCATCATCATGGCGCGCATCATCGGCGTACCCCGCGCGTTCGAACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGagag                                                                                                                                                                <  1:316802/274‑1 (MQ=255)
                                                                                                 acCGGCGCAACCATCACGTTGTCGCCTACGATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCTGCATTTCAGTTGCGTGACGAAGCGCACCNCATCACAGGa                                                                                                                                                                                                                                  >  1:371846/1‑194 (MQ=255)
                                                                                                 acCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCCCGGCACACGATAAGCATTGCTACCGTGTAAACGGCTATGCATGGAGAGAAAACCAAACGCGCa                                                                                                                                                 >  1:24383/1‑275 (MQ=255)
                                                                                                 acCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAAcgcga                                                                                                                                                  >  1:224077/1‑273 (MQ=255)
                                                                                                                     gTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATAACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAACGCGCACCAGCGGATGTAAACGTgcg                                                                                                                             >  1:223155/1‑275 (MQ=255)
                                                                                                                                         ggTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCacac                                                                                                                                                                                                        >  2:270053/1‑180 (MQ=255)
                                                                                                                                           tCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGAGACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGagag                                                                                                                                                                <  1:193269/218‑1 (MQ=255)
                                                                                                                                             aaGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATa                                                                                                                                                                                                    <  1:270064/180‑1 (MQ=255)
                                                                                                                                                        aaGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGc                                                                                                                                                                                                                                                                          <  1:50119/99‑1 (MQ=255)
                                                                                                                                                                     cGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTgaa                                                                                                                                                                                                                                                       >  1:86552/1‑105 (MQ=255)
                                                                                                                                                                         aaCTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGaagaag                                                                                                                                                                                                                                                   <  2:86551/105‑1 (MQ=255)
                                                                                                                                                                              catcatcatGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTaaa                                                                                                                                                                                 >  1:174820/1‑166 (MQ=255)
                                                                                                                                                                                  atcatGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGc                                                                                                                                                                             <  2:174813/166‑1 (MQ=255)
                                                                                                                                                                                      tGGCACGCATCATCGGCGTAACCCGCGCGTTCGCACGTGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTTAGTTGCGTGAAGAAGCGCACCACATCATAGGACTCACGAACGTAGGCCCACgtc                                                                                                                                                                                                            >  1:390132/1‑129 (MQ=255)
                                                                                                                                                                                         cccGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATcac                                                                                                                                                                                                                                      >  1:78473/1‑102 (MQ=255)
                                                                                                                                                                                          ccGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCacac                                                                                                                                                                                                        <  2:390117/131‑1 (MQ=255)
                                                                                                                                                                                           cGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAgcgc                                                                                                                                                                                                                                                >  1:226642/1‑90 (MQ=255)
                                                                                                                                                                                             catcatCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGa                                                                                                                                                                                                                                  <  2:78472/102‑1 (MQ=255)
                                                                                                                                                                                             catcatCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGGAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAACGCGCACCACCGTACGTAAACGTGCGCCGGAGCGCGATATTCAAAACCGCCGAAACCGTGGCCCCAGCAACCAAAACCCTAAACGACAATAAAc                                                         >  2:203560/1‑271 (MQ=255)
                                                                                                                                                                                               tcatcGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGcacca                                                                                                                                                                                                                                            <  2:226633/90‑1 (MQ=255)
                                                                                                                                                                                                   cGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATcac                                                                                                                                                                                                                                      >  1:108330/1‑92 (MQ=255)
                                                                                                                                                                                                       gTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGa                                                                                                                                                                                                                                  <  2:108329/92‑1 (MQ=255)
                                                                                                                                                                                                        tACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGTGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAAcgc                                                                                                                                                    >  1:124425/1‑169 (MQ=255)
                                                                                                                                                                                                        tACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGa                                                                                                                                                                                                      >  2:347078/1‑119 (MQ=255)
                                                                                                                                                                                                        tACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAACGCGCACCAGCGTTTGTAACCGTGCGCCGGAGCGGTATTTTCAAAGCCGCTGATATCGTGGCTCCAGAAGCCAAAACCTGAAAGGCCAATAGACAAACCACCGCACAgg                                          >  1:263138/1‑275 (MQ=255)
                                                                                                                                                                                                            ccGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGTGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAACGCGCAc                                                                                                                                                <  2:124424/169‑1 (MQ=255)
                                                                                                                                                                                                            ccGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAg                                                                                                                                                                                                  <  1:347090/119‑1 (MQ=255)
                                                                                                                                                                                                             cgcgcgTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTTTAAACGGCTATGGCTGGAGAGCAAACCAAACGTCCACCAGCGTTTGTAAACGNGCGCCGGAGCGGTATTTTCAAAGCCGCCGATATCGTGGCTCCAGAAGCCAAAACCTAACACGCCAATAGAAAAACCACCGCGCAGGCTTTc                                     >  1:393651/1‑275 (MQ=255)
                                                                                                                                                                                                             cgcgcgTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTaaa                                                                                                                                                                                 >  1:34283/1‑135 (MQ=255)
                                                                                                                                                                                                                gcgTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAgcgc                                                                                                                                                                                                                                                >  2:42413/1‑69 (MQ=255)
                                                                                                                                                                                                                 cgTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGc                                                                                                                                                                             <  2:34283/135‑1 (MQ=255)
                                                                                                                                                                                                                    tCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGcacca                                                                                                                                                                                                                                            <  1:42413/69‑1 (MQ=255)
                                                                                                                                                                                                                     cGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAACGCGCACCAGCGTTTGTAAACGTGCGCCGGAGCGGTATTTTCAAAGCCGCCGATATCGTGGCTCCAGAAGCCAAAACCTGAAAGGCCAATAGACAAACCACCGCGCAGGCTTTCCGCCAtt                              <  1:203568/274‑1 (MQ=255)
                                                                                                                                                                                                                                           aGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTaaa                                                                                                                                                                                 >  1:350203/1‑105 (MQ=255)
                                                                                                                                                                                                                                             aTACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTaaa                                                                                                                                                                                 >  1:127338/1‑103 (MQ=255)
                                                                                                                                                                                                                                               aCGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGc                                                                                                                                                                             <  2:350191/105‑1 (MQ=255)
                                                                                                                                                                                                                                                 ggCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGc                                                                                                                                                                             <  2:127337/103‑1 (MQ=255)
                                                                                                                                                                                                                                                  gCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCCCGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAACGCGCACCAGCGTTTGTAAACGTGCGCCGGAGCGGTATTTTCAAAGCCGCCGATATCGTGGCTCCAGAAGCCAAAACCTTAAAGGCCAATAAACAAACCACCGCACAGGCTTTCCGCCATTGATTCGTAGTTAGCGTAACACTCAcccc  >  2:133585/1‑271 (MQ=255)
                                                                                                                                                                                                                                                    atcatGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAACGCGCACCAGCGTTTGTaaa                                                                                                                                   >  2:433593/1‑142 (MQ=255)
                                                                                                                                                                                                                                                     tcatGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTa                                                                                                                                                                                                                      >  2:418778/1‑58 (MQ=255)
                                                                                                                                                                                                                                                      catGCGGCATTTCAGTGGCGTGAAGAAGCGCACCACATCAGAGGACTAAACATCGTAGGCCGACGGCACACGATAAGATTTGCTACCGGGTAAACGGCTATGGCTGGAGAGCAAACCAAACGCGCACCAGCGTTTGTAAACGTGCGCCGGAGCGGTATTTTCAAAGCCGCCGATATCGTGGCTCCAGAAGCCAAAACCTGAAAGGCCAATAGACACACCACCGCGCAGGCTTTCCGCCATTGAt                           <  1:142383/244‑1 (MQ=255)
                                                                                                                                                                                                                                                        tGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAACGCGCACCAGCGTTTCTAAACGTg                                                                                                                               <  1:433608/142‑1 (MQ=255)
                                                                                                                                                                                                                                                         gCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGcc                                                                                                                                                                                                                  <  1:418793/58‑1 (MQ=255)
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                        
TTTATGCCAGCGACTACCGTCGAGTTCATCGTTGTGCCACAGGTGTGTCCAGCGACCTTCCGGCAGGTAGAACTGCACATCGCCCGCTTCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAACGCGCACCAGCGTTTGTAAACGTGCGCCGGAGCGGTATTTTCAAAGCCGCCGATATCGTGGCTCCAGAAGCCAAAACCTGAAAGGCCAATAGACAAACCACCGCGCAGGCTTTCCGCCATTGATTCGTAGTTAGCGTAACAATCGCCAC  >  CP009273/3825945‑3826459

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