Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009273 4,568,602 C→A E185* (GAG→TAG)  mcrB ← 5‑methylcytosine‑specific restriction enzyme McrBC, subunit McrB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,568,6020CA90.0% 154.3 / 4.9 60E185* (GAG→TAG) mcrB5‑methylcytosine‑specific restriction enzyme McrBC, subunit McrB
Reads supporting (aligned to +/- strand):  ref base C (4/2);  new base A (27/27);  total (31/29)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.75e-01

GGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCGTGCAACAAAGGTTTTTCCAACGCCGGGCGGCCCCTGGAGGATAATATTTTTTTTGATGGTTAATCGTTTGAGTATCGTCTCTATTGTGGTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCAGTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGATTATGTTGTCCAGCATAGAGGCAAA  >  CP009273/4568458‑4568746
                                                                                                                                                |                                                                                                                                                
ggAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCGTGCAACAAAGGTTTTTCCAACGCCGGGCGGCCCCTGGAGGATAATATTTTTTTTGATGGTTAATCGTTTGAGTATCGTCTATAtt                                                                                                                                              >  2:902209/1‑149 (MQ=255)
   aCTGTACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCGTGCAACAAAGGTTTTTCCAACGCCGGGCGGCCCCTGGAGGATAATATTTTTTTTGATGGTTAATCGTTTGAGTATCGTCTATATtgtg                                                                                                                                           <  2:776749/149‑1 (MQ=255)
   aCTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCGTGCAACAAAGGTTTTTCCAACGCCGGGCGGCCCCTGGAGGATACTATTTTTTTTGATGGTTAATCGTTTGAGTATCGTCTATATtgtg                                                                                                                                           >  2:54927/1‑149 (MQ=255)
     tGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCGTGCAACAAAGGTTTTTCCAACGCCGGGCGGCCCCTGGAGGATAATATTTTTTTTGATGGTTAATCGTTTGAGTATCGTCTATATTGTGGt                                                                                                                                         <  2:187954/149‑1 (MQ=255)
         ccATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCATATCAATAATGTTTTTCCAACGCCGGGCGGCCCCTGGAGGATAATATTTTTTTTGATGGTTAATCGTTTGAGTATCGTCTATATTGTGGTTTCa                                                                                                                                     <  2:1017681/149‑1 (MQ=255)
          cATATTGACGCGTTGCGTAGCTTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCTGGCAACAAAGGTTTTTCCAACGCCGGGCGGCCCCTGGAGGATAATATTTTTTTTGATGGTTAATCGTTTGAGTATCGTCTATATTGTGGt                                                                                                                                         <  1:412012/144‑1 (MQ=255)
              ttGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCGTGCAACAAAGGTTTTTCCAACGCCGGGCGGCCCCTGTAGGATAATATTTTTTTTGATGGTTAATCGTTTGAGTATCGTCTCTATTGTGGTTTCAGGGAt                                                                                                                                >  2:701599/1‑149 (MQ=255)
              ttGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCGTGCAACAAAGGTTTTTCCAACGCCGGGCGGCCCCTGGAGGATAATATTTTTTTTGATGGTTAATCGCTGGAGCATCGTCTATAGTGTGGTTTCa                                                                                                                                     >  1:366830/1‑144 (MQ=255)
              ttGACGCGTTGCGGAGCCTTTTCTCCTGGCAGCACGTAAGCCAGACGGCGTCCAACAAAGGGTTTTCCACCGCCCGGCGGCCCCTGGAGGAGCATATTTTTTTTGATGGTTAATCGTTTGAGTAGCGGGTATATTGTGtttttt                                                                                                                                     >  1:366483/1‑142 (MQ=255)
               tGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCGTGCAACAAAGGTTTTTCCAACGCCGGGCGGCCCCTGGAGGATAATATTTTTTTTGATGGTTAATCGTTTTAGTATCGTCTATATTGTGGTTTCAGGGATa                                                                                                                               >  2:83476/1‑149 (MQ=255)
                gACGCGTTGCGGAGCCTTTTCGCCTGTCAGCAAGTAAGCCAGACGGCGTGCAACAAAGGTTTTTCCAACGCCGGGCGGCCCCTGGAGGATAATATTTTTTTTGATGGTTAATCGTTTGAGTATCGTCTCTATTGTGGtttgcgg                                                                                                                                   <  1:352636/144‑5 (MQ=255)
                    cgTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCGTGCAACAAAGGTTTTCCAACGCCGGGCGGCCCCTGGAGGATAATATTTTTTTTGATGGTTAATCGTTTGAGTATCGTCTATATTGTGGTTTCAGGGATGAACaca                                                                                                                         >  2:438026/1‑147 (MQ=255)
                          ggAGCCTTTTCCCCTGTCAGCAAGTAAGCCAGACGGCGTGCAACAAAGGTTTTTCCAACGCCGGGCGGCCCCTGGAGGATAATATTTTTTTTGATGGTTAATCGTTTGAGTATCGTCTATATTGTGGTTTCAGGGATAAACAAATCAtt                                                                                                                    <  2:853862/149‑1 (MQ=255)
                          ggAGCCTTTGCTCCTGTCAGCAAGTAAGCCAGACGGCGTGCAACAAAGGTTTTTCCAACGCCGGGCGGCCCCTGGAGGATAATATTTTTTTTGATGGTTAATCGTTTGAGTATCGTCTATATTGTGGTTTCAGGGATAAACAAATCAtt                                                                                                                    <  2:443048/149‑1 (MQ=255)
                           gAGCCTTTGCTCCTGTCAGCAAGTAAGCCAGACGGCGTGCAACCAAGGTTTTTCCAACCCCGGGCGGCCCCTGGAGGATAATATTTTTTTTGATGGTTAATCGTTTGAGTATCGTCTATATTGTGGTTTCAGGGATAAACAAAt                                                                                                                        <  1:715779/144‑1 (MQ=255)
                               cTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCGTGCAACAAAGGTTTTTCCAACGCCGGGCGGCACCTGGATGATAATATTTTTTTTGATGGTTAATCGTTTGAGTATCGTCTATATTGTGGTTTCAGGGATAAACAAATCAtt                                                                                                                    <  1:818672/144‑1 (MQ=255)
                                        tGTCAGCAAGTAAGCCAGACGGCGTGCAACAAAGGTTTTTCCAACGCCGGGCGGCCCCTGGAGGATAATATTTTTTTTGATGGTTAATCGTTTGAGTATCGTCTATATTGTGGTTTCAGGGATAAACAACTCATTTAACGCATCTTCCa                                                                                                      >  2:208991/1‑149 (MQ=255)
                                        tGTCAGCAAGTAAGCCAGACGGCGTGCAACAAAGGTTTTTCCAACGCCGGGCGGCCCCTGGAGGATAATATTTTTTTTGATGGTTAATCGTTTGAGTATCGTCTATATTGTGGTTTCAGGGATAAACAAATCATTTAACGCATCTTCCa                                                                                                      >  2:208984/1‑149 (MQ=255)
                                         gTCAGCAAGTAAGCCAGACGGCGTGCAACAAAGGTTTTTCCAACGCCGGGCGGCCCCTGGAGGATAATATTTTTTTTGATGGTTAATCGTTTGAGTATCGTCTATATTGTGGTTTCAGGGATAAACAAATCATTTAACgactct                                                                                                          <  1:831317/144‑6 (MQ=255)
                                          tCAGCAAGTAAGCCAGACGGCGTGCAACAAAGGTTTTTCCAACGCCGGGCGGCCCCGGGAGGATAATATTTTTTTTGCTGGTTACTCGTTTGGGTATCGGCTATATTGTGGTTTCAGGGATAAACAAATCATTTAACGCATCtt                                                                                                         >  1:174962/1‑144 (MQ=255)
                                                aGTAAGCCATACGGCGTGCAACAAAGGTTTTTCAAACGGCGGTCTGCCCCTGGAGGATAATATTTTTTTTGATGGTTAATCGTTTGAGTATCGTCTATATTGTGGTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTAt                                                                                              <  2:453830/149‑1 (MQ=255)
                                                aGTAAGCCAGACGGCGTGGCACAAAGGTTTTTCCAACGCCGGGCGGCCCCTGGAGGATAATATTTTTTTTGATGGTTAATCGTTTGAGTATCGTCTATATTGTGGTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTAt                                                                                              <  2:685661/149‑1 (MQ=255)
                                                     gCCAGACGGCGTGCAACAAAGGTTTTTCCAACGCCGGGCGGCCCCTGGAGGATACTATTTTTTTTGATGGTTAATCGTTTGAGTATCGTCTATATTGTGGTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTag                                                                                              >  1:257953/1‑143 (MQ=255)
                                                         gACGGCGTGCAACAAAGGTTTTTCCAACGCCGGGCGGCCCCTGGAGGATAATATTTTTTTTGATGGTTAATCGTTTGAGTATCGTCTATATTGTGGTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCAGtt                                                                                     >  2:44081/1‑149 (MQ=255)
                                                           ctgCGTGCAACAAAGGTTTTTCCAACGCCGGGCGGCCCCTGGAGGATAATATTTTTTTTGATGGTTAATCGTTTGAGTATCGTCTATATTGTGGTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCa                                                                                        <  1:711236/142‑1 (MQ=255)
                                                           cGGCGTGCAACAAAGGTTTTTCCAACGCCGGGCGGCCCCTGGAGGATAATATTTTTTTTGACGGCCAATGCTGACAGTACCCTCGATAGGGTTTTTTCAGGGATAAACAAATCATTTAACGCATCTGCCTGACAGTATGATTCa                                                                                        <  1:433091/144‑1 (MQ=255)
                                                            tgCGTGCAAACAGGATATTACCGACGCCGGGCGGCCCCTGGAGGATAATATTTTTTTTGATGGGTAATCGTTTGAGTATCGTCTATATTGTGGTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACTGTATGATTCAGttttt                                                                                  <  2:370491/148‑1 (MQ=255)
                                                              cGTGCAACAAATGGTTTTCCAACGCGGGGCGGCCCCTGGAGGATAATATTTTTTTTGATGGTTAATCGTTTGAGTATCGTCTCTATTGTGGTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCAGTTTTTGa                                                                                <  2:699453/149‑1 (MQ=255)
                                                               gTGCAACAAAGGTTTTTCCAACGCCGGGCGGCCCCTGGAGGATAATATTTTTTTTGATGATTAATCGTTTGAGTATCGTCTCTATTGTGGTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCAGttt                                                                                    >  1:704609/1‑144 (MQ=255)
                                                                 gCAACAAAGGTTTTTCCAACGCCGGGCGGCCCGTGGAGGATAATATTTTTTTTGATGGTTAATCGTTTGAGTATCGTCTATATTGTGGTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCAGttttt                                                                                  <  1:49168/144‑1 (MQ=255)
                                                                  caacaaAGGTTTTTCCAACGCCGGGCGGCCCCTGGAGGATAATATTTTTTTTGATGGTTAATCGTTTGAGTATCGTCTATATTGTGGGTTCAGGGATAAACAAACCATTTAACGCATCTTCCAGACAGTATGATTCAGTTTTTGACATa                                                                            >  2:857943/1‑149 (MQ=255)
                                                                   aacaaAGGTTTTTCCAACGCCGGGCGGCCCCTGGAGGATAATATTTTTTTTGATGGTTAATCGTTTGAGTATCGTCTATATTGTGGTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCAGTTTTTGa                                                                                <  1:421083/144‑1 (MQ=255)
                                                                      aaaGGTTTTTCCAACGCCGGGCGGCCCCTGGAGGATAATATTTTTTTTGATGGTTAATCGTTTGAGTATTGTCTATATTGTGGTTTCAGGGATCAACACATCATTTAACGCGTCTTCCAGACCGTCTGGTTCAGTTTTTGatat                                                                             >  1:295000/1‑144 (MQ=255)
                                                                      aaaGGTTTTTCCAACGCCGGGCGGCCCCTGGAGGATAATATTTTTTTTGATGGTTAATCGTTTGAGTATCGTCTATATTGTGGTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCAGTTTTTGACCTCGGTg                                                                        >  2:57338/1‑149 (MQ=255)
                                                                           ttttttCAACGCTTGGCGGCCCCTGGAGGATAATATTTTTTTTGATGGTTAATCGTTTGAGTATCGTCTATATTGTGGTTTCAGGGTTAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCAGTTTTTGACATAGGTg                                                                        >  1:356018/1‑144 (MQ=255)
                                                                             tttCCAACGCCGGGCGGCCCCTGGAGGATAATATTTTTTTTGATGGTTAATCGTTTGAGTATCGTCTATATTGTGGTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCAGTTTTTGACATAGGTGGa                                                                      <  1:947759/144‑1 (MQ=255)
                                                                             tttCCAACGCCGGGCGGCCCCTGGAGGATAATATTTTTTTTGATGGTTAATCGTTTGAGTATCGTCTATATTGTGGTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCAGTTTTTGACATAGGTGGAATAac                                                                 >  2:422479/1‑149 (MQ=255)
                                                                                ccAACGCCGGGCGGCCCCTGGAGGATAATATTTTTTTTGATGGTTAATCGTTTGAGTATCGTCTATATTGTGGTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCAGTTTTTGACATAGGTGGAATAACAct                                                              <  2:325720/149‑1 (MQ=255)
                                                                                ccAACGCCGGGCGGCCCCTGGAGGATAATATTTTTTTTGATGGTTAATCGTTTGAGTATCGTCTATATTGTGGTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCAGTTTTTGACATAAGTAGAATAACAct                                                              >  2:850755/1‑149 (MQ=255)
                                                                                     gCCGGGCGGCCCCTGGAGGATAATATTTTTTTTGATGGTTAATCGTTTGAGTATCGTCTATATTGTGGTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCAGTTTTTGACATAGGTGGAATAACAct                                                              <  1:794944/144‑1 (MQ=255)
                                                                                      ccGGGCGGCCCCTGGAGGATAATATTTTTTTTGATGGTTAATCGTTTGAGTATCGTCTATATTGTGGTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCCGTTTTTGACATAGGTGGAATAACActc                                                             >  1:34909/1‑144 (MQ=255)
                                                                                      ccGGGCGGCCCCTGGAGGATAATATTTTTTTTGATGGTTAATCGTTTGAGTATCGTCTATATTGTGGTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCAGTTTTTGACATAGGTGGAATAACACTCTTGcc                                                        >  2:838795/1‑149 (MQ=255)
                                                                                                    gaggaTAATATTTTTTTTGATGGTTAATCGTTTGAGTATCGTCTATATTGTGGTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCAGTTTTTGACATAGGTGGAACAACACTCTTGCCAGAATTaaa                                               <  1:788625/144‑1 (MQ=255)
                                                                                                       gaTAATATTTTTTTTGATGGTTAATCGTTTGAGTATCGTCTATATTGTGGTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCAGTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAAtt                                       >  2:1002302/1‑149 (MQ=255)
                                                                                                          aattTTTTTTTTGATGGTTAATCGTTTGAGTATCGTCTATATTGTGGTTTCAGGGATTAACAAATCATTTAACGTATCTTCCAGACAGTATGATTCAGTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTaa                                         <  1:432133/144‑1 (MQ=255)
                                                                                                          aataTTTTTTTTGATGGTTAATCGTTTGAGTATCGTCTCTATTGTGGTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCAGTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTtat                                    >  2:964993/1‑149 (MQ=255)
                                                                                                           attgtttttttGATGGTTAAGCGTTTGAGTATCGTCTATATTGTTGTTTCAGGGATAAACAAATCATTTAACGCATCTTCCATACAGTATGATTCAGTTTTTGACATAGGTGGAATAACACTCTTGCCATAGTTAAATATTAAt                                        <  1:829802/140‑1 (MQ=255)
                                                                                                            tatTTTTTTTGTTGTTTTTTCTATTGAGTATCGTCTCTATTGTGGTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTAGGATTCAGTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTaagt                                       >  1:555477/1‑142 (MQ=255)
                                                                                                              ttttttttGATGGTTAATCGTGTGAGTATCGTCTATATTGTGGTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCAGTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATAGTAATTTATAGTc                                >  2:666145/1‑149 (MQ=255)
                                                                                                                  ttttttttttttttCGTTTGAGTATCGTCTATATTGTGTTTTCAGGGCTACACAAATCATTTAACTGATCTTCCAGACAGTATGATTCTGTTTTTGACATTGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGt                                 >  1:14869/14‑144 (MQ=255)
                                                                                                                    ttGAGGGTTAATCGTTTGAGTATCGTCGATATTGTGGTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCAGTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTGtcgtcg                               <  1:137504/144‑1 (MQ=255)
                                                                                                                               tCGTTTGAGTATCGTCTATATTGTGGTGTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCAGTTTTTGACATAGGTGGAATAACACTCTTGCCCGAATTAAATATTAATTTTTAGTCGTGGATTATGTTGTCCa               >  2:478392/1‑149 (MQ=255)
                                                                                                                                cGTTTGAGTATCGTCTATATTGTGGTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCAGTTTTTGACATAGGTGGAATATCACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGATTAttttg                   <  1:777162/144‑1 (MQ=255)
                                                                                                                                 gTTTGAGTATCGTCTATATTGTGGTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCAGTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGATTAtgtggt                  >  1:948883/1‑144 (MQ=255)
                                                                                                                                 gTTTGAGTATCGTCTATATTGTGGTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCAGTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTAGTTGATTAtgttgc                  >  1:151847/1‑143 (MQ=255)
                                                                                                                                   ttGAGTATCGTCTATATTGTGGTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCAGTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGATTATGTTGTCCAGCAt           >  2:69233/1‑149 (MQ=255)
                                                                                                                                      aGTATCGTTTATATTGTGGTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCAGTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGATTATGGTGTCCAGCATaga        <  2:709461/149‑1 (MQ=255)
                                                                                                                                         aTCGTCTATATTGTGGTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCAGTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGATTATGTTGTCCAGCATa          <  1:75877/144‑1 (MQ=255)
                                                                                                                                         aTCGTCTATATTGTGGTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCAGTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGATTATGTTGTCCAGCATAGAGGc     <  2:317366/149‑1 (MQ=255)
                                                                                                                                            gTCTATATTGTGGTTTCAGGGATAACAAACGCAATAAACGCATCTTCCAGACAGTATGATTCAGTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGATTATGTTGTCCAGCATAGAGGCaaa  >  2:1028543/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                
GGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCGTGCAACAAAGGTTTTTCCAACGCCGGGCGGCCCCTGGAGGATAATATTTTTTTTGATGGTTAATCGTTTGAGTATCGTCTCTATTGTGGTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCAGTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGATTATGTTGTCCAGCATAGAGGCAAA  >  CP009273/4568458‑4568746

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: