Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,305,623 C→A 5.2% V93L (GTG→TTG)  phnN ← ribose 1,5‑bisphosphokinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,305,6230CA 5.2% 1390.3 / 30.7 347V93L (GTG→TTG) phnNribose 1,5‑bisphosphokinase
Reads supporting (aligned to +/- strand):  ref base C (153/176);  new base A (10/8);  total (163/184)
Fisher's exact test for biased strand distribution p-value = 4.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 350 total aligned reads displayed.

GGCGTTAATTTCACTGGCGTTTTCACGGCCACGGTTTTCCAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTGCTCGCTCAGGGCGATATGGTTTTCACTTCCGGC  >  CP009273/4305474‑4305770
                                                                                                                                                     |                                                                                                                                                   
 gcgTTAATTTCACTGGCGTTTTCACGGCCACGGTTTTCCAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACg                                                                                                                                                    >  2:4515291/1‑150 (MQ=255)
  cgTTAATTTCACTGGCGTTTTCACGGCCACGGTTTTCCAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGt                                                                                                                                                   >  2:964882/1‑150 (MQ=255)
   gTTAATTTCACTGGCGTTTTCACGGCCACGGTTTTCCAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTc                                                                                                                                                  >  2:5107647/1‑150 (MQ=255)
     tAATTTCACTGGCGTTTTCACGGCCACGGTTTTCCAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGa                                                                                                                                                >  1:657391/1‑150 (MQ=255)
     tAATTTCACTGGCGTTTTCACGGCCACGGTTTTCCAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGa                                                                                                                                                >  1:2841896/1‑150 (MQ=255)
            aCTGGCGTTTTCACGGCCACGGTTTTCCAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGc                                                                                                                                         <  2:2222243/150‑1 (MQ=255)
             cTGGCGTTTTCACGGCCACGGTTTTCCAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCg                                                                                                                                        >  1:1063590/1‑150 (MQ=255)
               ggCGTTTTCACGGCCACGGTTTTCCAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTg                                                                                                                                      <  2:3841011/150‑1 (MQ=255)
               ggCGTTTTCACGGCCACGGTTTTCCAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTg                                                                                                                                      <  1:15777/150‑1 (MQ=255)
               ggCGTTTTCACGGCCACGGTTTTCCAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCAAGTCGAATCCGGCGTg                                                                                                                                      >  1:3230688/1‑150 (MQ=255)
                gCGTTTTCACGGCCACGGTTTTCCAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAGTCCGGCGTGc                                                                                                                                     <  2:2857002/150‑1 (MQ=255)
                  gTTTTCACGGCCACGGTTTTCCAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAg                                                                                                                                   <  2:2064103/150‑1 (MQ=255)
                    tttCACGGCCACGGTTTTCCAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGcc                                                                                                                                 <  1:1731844/150‑1 (MQ=255)
                     ttaaCGGCCACGGTTTTCAAGGCGCTGTCGGAGGATCTCCGGCGAAAACTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCTCGTGAGCCGTAGACCAGCACGTCGAATCCGGCGTGCAGCCa                                                                                                                                >  2:2719681/4‑150 (MQ=255)
                     ttCACGGCCACGGTTTTCCAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCa                                                                                                                                <  2:13491/150‑1 (MQ=255)
                      tCACGGCCACGGTTTTCCAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCag                                                                                                                               <  2:3183973/150‑1 (MQ=255)
                      tCACGGCCACGGTTTTCCAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCag                                                                                                                               <  1:5071504/150‑1 (MQ=255)
                      tCACGGCCACGGTTGTCCAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCag                                                                                                                               <  2:4865566/150‑1 (MQ=255)
                       cACGGCCACGGTTTTCCAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCaga                                                                                                                              >  2:1841024/1‑150 (MQ=255)
                        aCGGCCACGGTTTTCCAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCagag                                                                                                                             <  2:4983806/150‑1 (MQ=255)
                          ggCCACGGTTTTCCAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGAt                                                                                                                           <  2:1008897/150‑1 (MQ=255)
                           gCCACGGTTTTCCAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATc                                                                                                                          >  2:525699/1‑150 (MQ=255)
                             cACGGTTTTCCAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCaa                                                                                                                        >  2:1693160/1‑150 (MQ=255)
                             cACGGTTTTCCAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCaa                                                                                                                        >  1:3658975/1‑150 (MQ=255)
                                 gTTTTCCAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAAtctc                                                                                                                    <  2:1342232/150‑1 (MQ=255)
                                  ttttCCAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTcg                                                                                                                   >  1:4383045/1‑150 (MQ=255)
                                  ttttCCAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTcg                                                                                                                   >  1:4574182/1‑150 (MQ=255)
                                   tttCCAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGACTCCGGCGTGCAGCCAGAGATCAATCTcga                                                                                                                  >  2:4384207/1‑150 (MQ=255)
                                   tttCCAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTcga                                                                                                                  <  2:1111162/150‑1 (MQ=255)
                                   tttCCAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTcga                                                                                                                  >  1:3637194/1‑150 (MQ=255)
                                       cAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTcgacgcc                                                                                                              <  2:363473/150‑1 (MQ=255)
                                       cAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTcgacgcc                                                                                                              <  2:307183/150‑1 (MQ=255)
                                        aGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTcgacgccg                                                                                                             >  1:879246/1‑150 (MQ=255)
                                        aGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTcgacgccg                                                                                                             >  2:1563064/1‑150 (MQ=255)
                                         ggCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTcgacgccga                                                                                                            >  2:2516770/1‑150 (MQ=255)
                                                 cGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCataa                                                                                                    >  1:4706967/1‑150 (MQ=255)
                                                 cGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCataa                                                                                                    <  1:2435865/150‑1 (MQ=255)
                                                  ggaggaTCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCataat                                                                                                   <  2:3014507/150‑1 (MQ=255)
                                                  ggaggaTCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCataat                                                                                                   >  2:4308966/1‑150 (MQ=255)
                                                    aggaTCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATAc                                                                                                 >  1:1911416/1‑150 (MQ=255)
                                                      gaTCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAg                                                                                               <  1:3061199/150‑1 (MQ=255)
                                                       aTCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGa                                                                                              >  2:2961653/1‑150 (MQ=255)
                                                         ctcCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGAcc                                                                                            <  1:303328/150‑1 (MQ=255)
                                                            cGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGtt                                                                                         >  1:3712119/1‑150 (MQ=255)
                                                            cGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGtt                                                                                         >  1:2494124/1‑150 (MQ=255)
                                                             ggCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCCGAGATCAATCTCGACGCCGACGCCATCATACAGACCGTTa                                                                                        >  2:372755/1‑150 (MQ=255)
                                                               cGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGc                                                                                      >  1:3430316/1‑150 (MQ=255)
                                                                    ccTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGcc                                                                                 <  2:2288175/150‑1 (MQ=255)
                                                                     cTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCa                                                                                <  1:3941770/150‑1 (MQ=255)
                                                                      tGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAg                                                                               >  2:3895006/1‑150 (MQ=255)
                                                                           aCAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCAAGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTaa                                                                          <  1:291343/150‑1 (MQ=255)
                                                                           aCAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCAAGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTaa                                                                          <  1:2887030/150‑1 (MQ=255)
                                                                              gACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGc                                                                       <  2:1945814/150‑1 (MQ=255)
                                                                               aCGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGcc                                                                      <  1:621947/150‑1 (MQ=255)
                                                                                  ggCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAAc                                                                   >  2:1147376/1‑150 (MQ=255)
                                                                                  ggCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAAc                                                                   >  1:4783828/1‑150 (MQ=255)
                                                                                   tcagcagCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACa                                                                  >  2:1086012/2‑150 (MQ=255)
                                                                                    cagcagCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAg                                                                 >  2:3537404/1‑150 (MQ=255)
                                                                                     agcagcGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGa                                                                >  2:2455305/1‑150 (MQ=255)
                                                                                       cagcGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGAtt                                                              <  2:1063318/150‑1 (MQ=255)
                                                                                            ccGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGcc                                                         <  1:2113769/150‑1 (MQ=255)
                                                                                                 tGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTACGGCCAACAGATTTTGCCCcgcg                                                    >  2:314550/1‑150 (MQ=255)
                                                                                                 tGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCcgcg                                                    <  1:1944951/150‑1 (MQ=255)
                                                                                                 tGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCcgcg                                                    <  2:2599689/150‑1 (MQ=255)
                                                                                                 tGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCcgcg                                                    <  1:2214034/150‑1 (MQ=255)
                                                                                                 tGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCcgcg                                                    >  1:5130075/1‑150 (MQ=255)
                                                                                                  gATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCcgcgc                                                   >  1:3369264/1‑150 (MQ=255)
                                                                                                   aTAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCcgcgcg                                                  <  2:832477/150‑1 (MQ=255)
                                                                                                   aTAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCcgcgcg                                                  <  1:313454/150‑1 (MQ=255)
                                                                                                   aTAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCcgcgcg                                                  <  1:989320/150‑1 (MQ=255)
                                                                                                   aTAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCAAGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCcgcgcg                                                  >  1:3559984/1‑150 (MQ=255)
                                                                                                    tAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCgg                                                 >  1:1112729/1‑150 (MQ=255)
                                                                                                    tAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCgg                                                 >  2:3155260/1‑150 (MQ=255)
                                                                                                    tAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGACGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTCAGGCCAACAGATGATGCCCCGCGCgg                                                 >  2:2924737/1‑150 (MQ=255)
                                                                                                     aGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCggg                                                <  2:1639605/150‑1 (MQ=255)
                                                                                                      gcgcgcCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGAtttt                                                            <  1:1032646/137‑1 (MQ=255)
                                                                                                         cgcccgcgccTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTaaa                                            <  1:454625/150‑1 (MQ=255)
                                                                                                         cgcccgcgccTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTaaa                                            <  2:151546/150‑1 (MQ=255)
                                                                                                          gcccgcgccTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTaaaa                                           >  1:2005838/1‑150 (MQ=255)
                                                                                                           cccgcgccTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTaaaaa                                          <  1:4062574/150‑1 (MQ=255)
                                                                                                             cgcgccTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAAc                                        >  2:1624015/1‑150 (MQ=255)
                                                                                                              gcgccTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACt                                       <  2:1494044/150‑1 (MQ=255)
                                                                                                                 ccTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCt                                    <  2:4814672/150‑1 (MQ=255)
                                                                                                                   tGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTgc                                  >  1:1864913/1‑150 (MQ=255)
                                                                                                                      ggCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTgctcg                               >  1:2298403/1‑150 (MQ=255)
                                                                                                                      ggCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTgctcg                               <  1:2618913/150‑1 (MQ=255)
                                                                                                                        cAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTgctcgct                             <  1:4267906/150‑1 (MQ=255)
                                                                                                                         aGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTgctcgctc                            >  1:745741/1‑150 (MQ=255)
                                                                                                                         aGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTgctcgctc                            <  1:1162773/150‑1 (MQ=255)
                                                                                                                           aTGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTGCTCGCTCAg                          <  1:3189523/150‑1 (MQ=255)
                                                                                                                            tGGGCGCGTGAGCCGTTGACCAGCAAGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTGCTCGCTCAgg                         <  2:649839/150‑1 (MQ=255)
                                                                                                                                cgcgTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTGCTCGCTCAGGGCGa                     <  1:2690494/150‑1 (MQ=255)
                                                                                                                                cgcgTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTGCTCGCTCAGGGCGa                     <  1:4598485/150‑1 (MQ=255)
                                                                                                                                    tGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTGCTCGCTCAGGGCGATATg                 >  2:478126/1‑150 (MQ=255)
                                                                                                                                         cGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTGCTCGCTCAGGGCGATATGGtttt            <  2:1872844/150‑1 (MQ=255)
                                                                                                                                           ttGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTGCTCGCTCAGGGCGATATGGTTTTCa          <  2:4167327/150‑1 (MQ=255)
                                                                                                                                            tGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTGCTCGCTCAGGGCGATATGGTTTTCAc         >  2:3302697/1‑150 (MQ=255)
                                                                                                                                                cAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTGCTCGCTCAGGGCGATATGGTTTTCACTTcc     <  2:4011853/150‑1 (MQ=255)
                                                                                                                                                cAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTGCTCGCTCAGGGCGATATGGTTTTCACTTcc     <  1:4224458/150‑1 (MQ=255)
                                                                                                                                                cAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTGTGCCCCGCGCGGGTAAAAAACTCCTGCTCGCTCAGGGCGATAGGGTTTTCACTTcc     >  1:1382496/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GGCGTTAATTTCACTGGCGTTTTCACGGCCACGGTTTTCCAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTGCTCGCTCAGGGCGATATGGTTTTCACTTCCGGC  >  CP009273/4305474‑4305770

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