Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,208,127 A→T 22.7% D364V (GAC→GTC)  aceA → isocitrate lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,208,1270AT22.7% 808.1 / 261.7 352D364V (GAC→GTC) aceAisocitrate lyase
Reads supporting (aligned to +/- strand):  ref base A (142/130);  new base T (49/31);  total (191/161)
Fisher's exact test for biased strand distribution p-value = 1.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

Only 100 of 352 total aligned reads displayed.

GCTTATAACTGCTCGCCGTCGTTCAACTGGCAGAAAAACCTCGACGACAAAACTATTGCCAGCTTCCAGCAGCAGCTGTCGGATATGGGCTACAAGTTCCAGTTCATCACCCTGGCAGGTATCCACAGCATGTGGTTCAACATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGAC  >  CP009273/4207979‑4208274
                                                                                                                                                    |                                                                                                                                                   
gCTTATAACTGCTCGCCGTCGTTCAACTGGCAGAAAAACCTCGACGACAAAACTATTGCCAGCTTCCAGCAGCAGCTGTCGGATATGGGCTACAAGTTCCAGTTCATCACCCTGGCAGGTATCCACAGCATGTGGTTCAACATGTTTGAc                                                                                                                                                    <  2:1546321/150‑1 (MQ=255)
  ttATAACTGCTCGCCGTCGTTCAACTGGCAGAAAAACCTCGACGACAAAACTATTGCCAGCTTCCAGCAGCAGCTGTCGGATATGGGCTACAAGTTCCAGTTCATCACCCTGGCAGGTATCCACAGCATGTGGTTCAACATGTTTGTCCt                                                                                                                                                  >  1:265601/1‑150 (MQ=255)
  ttATAACTGCTCGCCGTCGTTCAACTGGCAGAAAAACCTCGACGACAAAACTATTGCCAGCTTCCAGCAGCAGCTGTCGGATATGGGCTACAAGTTCCAGTTCATCACCCTGGCAGGTATCCACAGCATGTGGTTCAACATGTTTGACCt                                                                                                                                                  >  1:3839637/1‑150 (MQ=255)
   tataACTGCTCGCCGTCGTTCAACTGGCAGAAAAACCTCGACGACAAAACTATTGCCAGCTTCCAGCAGCAGCTGTCGGATATGGGCTACAAGTTCCAGTTCATCACCCTGGCAGGTATCCACAGCATGTGGTTCAACATGTTTGTCCTg                                                                                                                                                 >  1:852256/1‑150 (MQ=255)
     taACTGCTCGCCGTCGTTCAACTGGCAGAAAAACCTCGACGACAAAACTATTGCCAGCTTCCAGCAGCAGCTGTCGGATATGGGCTACAAGTTCCAGTTCATCACCCTGGCAGGTATCCACAGCATGTGGTTCAACATGTTTGACCTGGc                                                                                                                                               >  1:3487269/1‑150 (MQ=255)
      aaCTGCTCGCCGTCGTTCAACTGGCAGAAAAACCTCGACGACAAAACTATTGCCAGCTTCCAGCAGCAGCTGTCGGATATGGGCTACAAGTTCCAGTTCATCACCCTGGCAGGTATCCACAGCATGTGGTTCAACATGTTTGACCTGGCa                                                                                                                                              >  2:667200/1‑150 (MQ=255)
       aCTGCTCGCCGTCGTTCAACTGGCAGAAAAACCTCGACGACAAAACTATTGCCAGCTTCCAGCAGCAGCTGTCGGATATGGGCTACAAGTTCCAGTTCATCACCCTGGCAGGTATCCACAGCATGTGGTTCAACATGTTTGACCTGGCaa                                                                                                                                             >  1:1088971/1‑150 (MQ=255)
        cTGCTCGCCGTCGTTCAACTGGCAGAAAAACCTCGACGACAAAACTATTGCCAGCTTCCAGCAGCAGCTGTCGGATATGGGCTACAAGTTCCAGTTCATCACCCTGGCAGGTATCCACAGCATGTGGTTCAACATGTTTGTCCTGGCaaa                                                                                                                                            >  2:3404586/1‑150 (MQ=255)
         tGCTCGCCGTCGTTCAACTGGCAGAAAAACCTCGACGACAAAACTATTGCCAGCTTCCAGCAGCAGCTGTCGGATATGGGCTACAAGTTCCAGTTCATCACCCTGGCAGGTATCCACAGCATGTGGTTCAACATGTTTGACCTGGCAAAc                                                                                                                                           >  2:3362643/1‑150 (MQ=255)
              gCCGTCGTTCAACTGGCAGAAAAACCTCGACGACAAAACTATTGCCAGCTTCCAGCAGCAGCTGTCGGATATGGGCTACAAGTTCCAGTTCATCACCCTGGCAGGTATCCACAGCATGTGGTTCAACATGTTTGACCTGGCAAACGCCTa                                                                                                                                      >  1:1459179/1‑150 (MQ=255)
               ccGTCGTTCAACTGGCAGAAAAACCTCGACGACAAAACTATTGCCAGCTTCCAGCAGCAGCTGTCGGATATGGGCTACAAGTTCCAGTTCATCACCCTGGCAGGTATCCACAGCATGTGGTTCAACATGTTTGACCTGGCAAACGCCTAt                                                                                                                                     <  2:2322662/150‑1 (MQ=255)
                          cTGGCAGAAAAACCTCGACGACAAAACTATTGCCAGCTTCCAGCAGCAGCTGTCGGATATGGGCTACAAGTTCCAGTTCATCACCCTGGCAGGTATCCACAGCATGTGGTTCAACATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGa                                                                                                                          >  2:4159134/1‑150 (MQ=255)
                               aGAAAAACCTCGACGACAAAACTATTGCCAGCTTCCAGCAGCAGCTGTCGGATATGGGCTACAAGTTCCAGTTCATCACCCTGGCAGGTATCCACAGCATGTGGTTCAACATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTa                                                                                                                     >  2:2537951/1‑150 (MQ=255)
                                  aaaaCCTCGACGACAAAACTATTGCCAGCTTCCAGCAGCAGCTGTCGGATATGGGCTACAAGTTCCAGTTCATCACCCTGGCAGGTATCCACAGCATGTGGTTCAACATGTTTGTCCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGa                                                                                                                  >  2:2499554/1‑150 (MQ=255)
                                    aaCCTCGACGACAAAACTATTGCCAGCTTCCAGCAGCAGCTGTCGGATATGGGCTACAAGTTCCAGTTCATCACCCTGGCAGGTATCCACAGCATGTGGTTCAACATGTTTGTCCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAg                                                                                                                <  2:547689/150‑1 (MQ=255)
                                     aCCTCGACGACAAAACTATTGCCAGCTTCCAGCAGCAGCTGTCGGATATGGGCTACAAGTTCCAGTTCATCACCCTGGCAGGTATCCACAGCATGTGGTTCAACATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGc                                                                                                               <  2:3826158/150‑1 (MQ=255)
                                      ccTCGACGACAAAACTATTGCCAGCTTCCAGCAGCAGCTGTCGGATATGGGCTACAAGTTCCAGTTCATCACCCTGGCAGGTATCCACAGCATGTGGTTCAACATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCa                                                                                                              >  2:4574857/1‑150 (MQ=255)
                                         cgacgaCAAAACTATTGCCAGCTTCCAGCAGCAGCTGTCGGATATGGGCTACAAGTTCCAGTTCATCACCCTGGCAGGTATCCACAGCATGTGGTTCAACATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTa                                                                                                           <  1:3053947/150‑1 (MQ=255)
                                           acgacAAAACTATTGCCAGCTTCCAGCAGCAGCTGTCGGATATGGGCTACAAGTTCCAGTTCATCACCCTGGCAGGTATCCACAGCATGTGGTTCAACATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACg                                                                                                         <  1:615788/150‑1 (MQ=255)
                                             gacAACACTATTGCCAGCTTCCAGCAGCAGCTGTCGGATATGGGCTACAAGTTCCAGTTCATCACCCTGGCAGGTATCCACAGCATGTGGTTCAACATGTTTGTCCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGtt                                                                                                       <  2:1241876/150‑1 (MQ=255)
                                             gacAAAACTATTGCCAGCTTCCAGCAGCAGCTGTCGGATATGGGCTACAAGTTCCAGTTCATCACCCTGGCAGGTATCCACAGCATGTGGTTCAACATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGtt                                                                                                       >  1:2390220/1‑150 (MQ=255)
                                               cAAAACTATTGCCAGCTTCCAGCAGCAGCTGTCGGATATGGGCTACAAGTTCCAGTTCATCACCCTGGCAGGTATCCACAGCATGTGGTTCAACATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTga                                                                                                     <  2:3881643/150‑1 (MQ=255)
                                                     tagtGCCAGCTTCCAGCAGCAGCTGTCGGATATGGGCTACAAGTTCCAGTTCATCACCCTGGCAGGTATCCACAGCATGTGGTTCAACATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGt                                                                                               <  1:3903748/147‑1 (MQ=255)
                                                     tATTGCCAGCTTCCAGCAGCAGCTGTCGGATATGGGCTACAAGTTCCAGTTCATCACCCTGGCAGGTATCCACAGCATGTGGTTCAACATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGt                                                                                               >  1:568961/1‑150 (MQ=255)
                                                        tGCCAGCTTCCAGCAGCAGCTGTCGGATATGGGCTACAAGTTCCAGTTCATCACCCTGGCAGGTATCCACAGCATGTGGTTCAACATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTgca                                                                                            <  2:4267694/150‑1 (MQ=255)
                                                         gCCAGCTTCCAGCAGCAGCTGTCGGATATGGGCTACAAGTTCCAGTTCATCACCCTGGCAGGTATCCACAGCATGTGGTTCAACATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTgcag                                                                                           >  2:794569/1‑150 (MQ=255)
                                                           cAGCTTCCAGCAGCAGCTGTCGGATATGGGCTACAAGTTCCAGTTCATCACCCTGGCAGGTATCCACAGCATGTGGTTCAACATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTgcagca                                                                                         <  2:1569333/150‑1 (MQ=255)
                                                            aGCTTCCAGCAGCAGCTGTCGGATATGGGCTACAAGTTCCAGTTCATCACCCTGGCAGGTATCCACAGCATGTGGTTCAACATGTTTGTCCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGcagcag                                                                                        >  2:2427312/1‑150 (MQ=255)
                                                            aGCTTCCAGCAGCAGCTGTCGGATATGGGCTACAAGTTCCAGTTCATCACCCTGGCAGGTATCCACAGCATGTGGTTCAACATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGcagcag                                                                                        <  1:3011595/150‑1 (MQ=255)
                                                             gCTTCCAGCAGCAGCTGTCGGATATGGGCTACAAGTTCCAGTTCATCACCCTGGCAGGTATCCACAGCATGTGGTTCAACATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCagcagc                                                                                       >  1:4423745/1‑150 (MQ=255)
                                                               ttCCAGCAGCAGCTGTCGGATATGGGCTACAAGTTCCAGTTCATCACCCTGGCAGGTATCCACAGCATGTGGTTCAACATGTTTGTCCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCg                                                                                     <  1:963810/150‑1 (MQ=255)
                                                               ttCCAGCAGCAGCTGTCGGATATGGGCTACAAGTTCCAGTTCATCACCCTGGCAGGTATCCACAGCATGTGGTTCAACATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCg                                                                                     <  1:1800344/150‑1 (MQ=255)
                                                                  cagcagcagCTGTCGGATATGGGCTACAAGTTCCAGTTCATCACCCTGGCAGGTATCCACAGCATGTGGTTCAACATGTTTGTCCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGaa                                                                                  >  1:3667233/1‑150 (MQ=255)
                                                                   agcagcagcTGTCGGATATGGGCTACAAGTTCCAGTTCATCACCCTGGCAGGTATCCACAGCATGTGGTTCAACATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAAt                                                                                 >  1:2574473/1‑150 (MQ=255)
                                                                    gcagcagcTGTCGGATATGGGCTACAAGTTCCAGTTCATCACCCTGGCAGGTATCCACAGCATGTGGTTCAACATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAAtt                                                                                <  1:4612351/150‑1 (MQ=255)
                                                                     cagcagcTGTCGGATATGGGCTACAAGTTCCAGTTCATCACCCTGGCAGGTATCCACAGCATGTGGTTCAACATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAAttt                                                                               >  1:1249783/1‑150 (MQ=255)
                                                                     cagcagcTGTCGGATATGGGCTACAAGTTCCAGTTCATCACCCTGGCAGGTATCCACAGCATGTGGTTCAACATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAAttt                                                                               <  2:1638319/150‑1 (MQ=255)
                                                                       gcagcTGTCGGATATGGGCTACAAGTTCCAGTTCATCACCCTGGCAGGTATCCACAGCATGTGGTTCAACATGTTTGTCCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTgc                                                                             >  1:1127685/1‑150 (MQ=255)
                                                                         agcTGTCGGATATGGGCTACAAGTTCCAGTTCATCACCCTGGCAGGTATCCACAGCATGTGGTTCAACATGTTTGTCCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTgccg                                                                           >  2:4082866/1‑150 (MQ=255)
                                                                         agcTGTCGGATATGGGCTACAAGTTCCAGTTCATCACCCTGGCAGGTATCCACAGCATGTGGTTCAACATGTTTGTCCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTgccg                                                                           >  2:2233708/1‑150 (MQ=255)
                                                                         agcTGTCGGATATGGGCTACAAGTTCCAGTTCATCACCCTGGCAGGTATCCACAGCATGTGGTTCAACATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTgccg                                                                           >  2:2938173/1‑150 (MQ=255)
                                                                           cTGTCGGATATGGGCTACAAGTTCCAGTTCGTCACCCTGGCAGGTATCCACAGCATGTGGTTCAACATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTgccgcc                                                                         <  2:21170/150‑1 (MQ=255)
                                                                           cTGTCGGATATGGGCTACAAGTTCCAGTTCATCACCCTGGCAGGTATCCACAGCATGTGGTTCAACATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTgccgcc                                                                         <  2:234082/150‑1 (MQ=255)
                                                                            tGTCGGATATGGGCTACAAGTTCCAGTTCATCACCCTGGCAGGTATCCACAGCATGTGGTTCAACATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGccgccg                                                                        >  2:3474650/1‑150 (MQ=255)
                                                                                ggATATGGGCTACAAGTTCCAGTTCATCACCCTGGCAGGTATCCACAGCATGTGGTTCAACATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGaa                                                                    >  2:4684165/1‑150 (MQ=255)
                                                                                  atatGGGCTACAAGTTCCAGTTCATCACCCTGGCAGGTATCCACAGCATGTGGTTCAACATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAg                                                                  <  2:850430/150‑1 (MQ=255)
                                                                                  atatGGGCTACAAGTTCCAGTTCATCACCCTGGCAGGTATCCACAGCATGTGGTTCAACATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAg                                                                  >  2:1561269/1‑150 (MQ=255)
                                                                                     tGGGCTACAAGTTCCAGTTCATCACCCTGGCAGGTATCCACAGCATGTGGTTCAACATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATg                                                               <  2:139498/150‑1 (MQ=255)
                                                                                      gggCTACAAGTTCCAGTTCATCACCCTGGCAGGTATCCACAGCATGTGGTTCAACATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATgg                                                              <  2:4493629/150‑1 (MQ=255)
                                                                                      gggCTACAAGTTCCAGTTCATCACCCTGGCAGGTATCCACAGCATGTGGTTCAACATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGGATGAAGCACTACGGTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATgg                                                              >  2:580773/1‑150 (MQ=255)
                                                                                         cTACAAGTTCCAGTTCATCACCCTGGCAGGTATCCACAGCATGTGGTTCAACATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCta                                                           >  1:715389/1‑150 (MQ=255)
                                                                                         cTACAAGTTCCAGTTCATCACCCTGGCAGGTATCCACAGCATGTGGTTCAACATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCta                                                           >  1:4821423/1‑150 (MQ=255)
                                                                                         cTACAAGTTCCAGTTCATCACCCTGGCAGGTATCCACAGCATGTGGTTCAACATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCta                                                           <  2:886538/150‑1 (MQ=255)
                                                                                          tACAAGTTCCAGTTCATCACCCTGGCAGGTATCCACAGCATGTGGTTCAACATGTTTGTCCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCtat                                                          <  1:1331979/150‑1 (MQ=255)
                                                                                          tACAAGTTCCAGTTCATCACCCTGGCAGGTATCCACAGCATGTGGTTCAACATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCtat                                                          >  1:4711799/1‑150 (MQ=255)
                                                                                                ttCCAGTTCATCACCCTGGCAGGTATCCACAGCATGTGGTTCAACATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTc                                                    >  1:3924973/1‑150 (MQ=255)
                                                                                                  ccAGTTCATCACCCTGGCAGGTATCCACAGCATGTGGTTCAACATGTTTGTCCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGt                                                  <  2:1127685/150‑1 (MQ=255)
                                                                                                   cAGTTCATCACCCTGGCAGGTATCCACAGCATGTGGTTCAACATGTTTGTCCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTa                                                 <  2:195690/150‑1 (MQ=255)
                                                                                                   cAGTTCATCACCCTGGCAGGTATCCACAGCATGTGGTTCAACATGTTTGTCCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTa                                                 <  1:4305011/150‑1 (MQ=255)
                                                                                                    aGTTCATCACCCTGGCAGGTATCCACAGCATGTGGTTCAACATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTAt                                                >  1:4543188/1‑150 (MQ=255)
                                                                                                    aGTTCATCACCCTGGCAGGTATCCACAGCATGTGGTTCAACATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTGGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTAt                                                >  1:1423782/1‑150 (MQ=255)
                                                                                                     gTTCATCACCCTGGCAGGTATCCACAGCATGTGGTTCAACATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTAtc                                               >  1:1409742/1‑150 (MQ=255)
                                                                                                         atcaCCCTGGCAGGTATCCACAGCATGTGGTTCAACATGTTTGTCCTGGCAAACGCCTATGCCCAGGGCGAGGGTCTGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGACAGATGGCTATCCCTGCGTCTCTCAc                                           >  2:4547783/1‑150 (MQ=255)
                                                                                                         atcaCCCTGGCAGGTATCCACAGCATGTGGTTCAACATGTTTGTCCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGACAGATGGCTATACCTTCGTATCTCAc                                           <  1:4138163/150‑1 (MQ=255)
                                                                                                         atcaCCCTGGCAGGTATCCACAGCATGTGGTTCAACATGTTTGTCCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCAc                                           <  1:1158657/150‑1 (MQ=255)
                                                                                                          tcaCCCTGGCAGGTATCCACAGCATGTGGTTCAACATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCAcc                                          <  1:1604075/150‑1 (MQ=255)
                                                                                                           caCCCTGGCAGGTATCCACAGCATGTGGTTCAACATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACca                                         <  2:907883/150‑1 (MQ=255)
                                                                                                               cTGGCAGGTATCCACAGCATGTGGTTCAACATGTTTGTCCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACcagcag                                     >  1:3172856/1‑150 (MQ=255)
                                                                                                                    aGGTATCCACAGCATGTGGTTCAACATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGt                                >  1:2186944/1‑150 (MQ=255)
                                                                                                                       tATCCACAGCATGTGGTTCAACATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTggg                             >  1:437386/1‑150 (MQ=255)
                                                                                                                       tATCCACAGCATGTGGTTCAACATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTggg                             >  1:3783166/1‑150 (MQ=255)
                                                                                                                        aTCCACAGCATGTGGTTCAACATGTTTGTCCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGt                            <  1:4707447/150‑1 (MQ=255)
                                                                                                                         tCCACAGCATGTGGTTCAACATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTa                           <  2:1452583/150‑1 (MQ=255)
                                                                                                                         tCCACAGCATGTGGTTCAACATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTa                           <  1:4013088/150‑1 (MQ=255)
                                                                                                                          ccACAGCATGTGGTTCAACATGTTTGTCCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTAc                          >  1:1516148/1‑150 (MQ=255)
                                                                                                                           cacaGCATGTGGTTCAACATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACa                         >  1:3246268/1‑150 (MQ=255)
                                                                                                                            acaGCATGTGGTTCAACATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAg                        <  2:2059297/150‑1 (MQ=255)
                                                                                                                             caGCATGTGGTTCAACATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAgg                       >  2:4520912/1‑150 (MQ=255)
                                                                                                                              aGCATGTGGTTCAACATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGt                      <  1:3577744/150‑1 (MQ=255)
                                                                                                                              aGCATGTGGTTCAACATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGt                      <  2:52516/150‑1 (MQ=255)
                                                                                                                                cATGTGGTTCAACATGTTTGTCCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTa                    >  1:2816497/1‑150 (MQ=255)
                                                                                                                                cATGTGGTTCAACATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGTAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTCTCTCACCAGCGGGTAGTGGGTACAGGTTa                    >  2:4012556/1‑150 (MQ=255)
                                                                                                                                cATGTGGTTCAACATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTa                    >  2:2623590/1‑150 (MQ=255)
                                                                                                                                 aTGTGGTTCAACATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTAc                   >  1:1455146/1‑150 (MQ=255)
                                                                                                                                  tgtgGTTCAACATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACt                  <  2:1071661/150‑1 (MQ=255)
                                                                                                                                  tgtgGTTCAACATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACt                  >  2:1218130/1‑150 (MQ=255)
                                                                                                                                  tgtgGTTCAACATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACt                  <  2:4331379/150‑1 (MQ=255)
                                                                                                                                    tgGTTCAACATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTc                >  2:3679413/1‑150 (MQ=255)
                                                                                                                                     ggTTCAACATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCg               <  2:1872047/150‑1 (MQ=255)
                                                                                                                                     ggTTCAACATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCg               >  1:1926956/1‑150 (MQ=255)
                                                                                                                                       ttCAACATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGAt             >  1:2857134/1‑150 (MQ=255)
                                                                                                                                          aaCATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATaaa          <  2:271646/150‑1 (MQ=255)
                                                                                                                                           aCATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAg         >  2:4435690/1‑150 (MQ=255)
                                                                                                                                             aTGTTTGTCCTGGCAAACGCCTATGCCAAGGGCGAGGGTATGAAGAACTACGTTGAGAAAGTGCAGAAGCCGGAATATGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTg       >  2:2261582/1‑150 (MQ=255)
                                                                                                                                              tGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTga      <  2:3573593/150‑1 (MQ=255)
                                                                                                                                              tGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTga      <  1:308131/150‑1 (MQ=255)
                                                                                                                                                tttGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTgagg    <  2:3077401/150‑3 (MQ=255)
                                                                                                                                                  tGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGACGTGGGTACAGGTTACTTCGATAAAGTgacgac  >  1:105452/1‑150 (MQ=255)
                                                                                                                                                  tGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTgacgac  <  2:1260968/150‑1 (MQ=255)
                                                                                                                                                  tGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTgacgac  >  2:3166571/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GCTTATAACTGCTCGCCGTCGTTCAACTGGCAGAAAAACCTCGACGACAAAACTATTGCCAGCTTCCAGCAGCAGCTGTCGGATATGGGCTACAAGTTCCAGTTCATCACCCTGGCAGGTATCCACAGCATGTGGTTCAACATGTTTGACCTGGCAAACGCCTATGCCCAGGGCGAGGGTATGAAGCACTACGTTGAGAAAGTGCAGCAGCCGGAATTTGCCGCCGCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGAC  >  CP009273/4207979‑4208274

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: ‑