Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,208,942 C→A 18.2% A140E (GCG→GAG)  aceK → isocitrate dehydrogenase kinase/phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,208,9420CA18.2% 857.9 / 117.4 290A140E (GCG→GAG) aceKisocitrate dehydrogenase kinase/phosphatase
Reads supporting (aligned to +/- strand):  ref base C (113/123);  new base A (30/24);  total (143/147)
Fisher's exact test for biased strand distribution p-value = 3.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 291 total aligned reads displayed.

GGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGTTATTTGACCACCGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACAACCTCG  >  CP009273/4208796‑4209091
                                                                                                                                                  |                                                                                                                                                     
          cgcTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGTTATTTGACCACCGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCAc                                                                                                                                          >  2:929365/1‑150 (MQ=255)
          cgcTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGTTATTTGACCACCGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGAGAAAGACTTTCAc                                                                                                                                          <  1:3200305/150‑1 (MQ=255)
                 agaTTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGTTATTTGACCACCGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATc                                                                                                                                   >  2:2269234/1‑150 (MQ=255)
                 agaTTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGGTATTTGACCACCGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATc                                                                                                                                   >  2:513328/1‑150 (MQ=255)
                    ttGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGTTATTTGACCACCGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACg                                                                                                                                <  2:3027168/150‑1 (MQ=255)
                    ttGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGTTATTTGACCACCGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACg                                                                                                                                <  2:3729757/150‑1 (MQ=255)
                       cGGAGAGCTTTTTTAACTCCGTGTACTGTCGGTTATTTGACCACCGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGAGAAAGACTTTCACCCCGATCACGGCt                                                                                                                             <  1:4163423/150‑1 (MQ=255)
                        ggAGAGCTTTTTTAACTCCGTGTACTGTCGGTTATTTGACCACCGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTg                                                                                                                            <  1:2390305/150‑1 (MQ=255)
                            agCTTTTTTAACTCCGTGTACTGTCGGTTATTTGACCACCGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGaa                                                                                                                        <  2:4431629/150‑1 (MQ=255)
                             gCTTTTTTAACTCCGTGTACTGTCGGTTATTTGTCCACCGCTCGCTTACTCCCGAGCGGCTTTTTCTCTTTAGCTCTCAGCCTGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAAt                                                                                                                       <  2:279831/150‑1 (MQ=255)
                              cTTTTTTAACTCCGTGTACTGTCGGTTATTTGACCACCGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAAtc                                                                                                                      <  2:4200175/150‑1 (MQ=255)
                                 ttttAACTCCGTGTACTGTCGGTTATTTGACCACCGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATctct                                                                                                                   <  2:2382019/150‑1 (MQ=255)
                                   ttAACTCCGTGTACTGTCGGTTATTTGAACACCGCACGATTACTCCCGAGCGGCTTTTAATCTTTAGCTCTCAGCCAGAGAGCCGAATTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTAACACCGATCAAGGCTGGGAATCTCTAc                                                                                                                 >  2:3451191/1‑150 (MQ=255)
                                    tAACTCCGTGTACTGTCGGTTATTTGACCACCGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACt                                                                                                                <  2:2685347/150‑1 (MQ=255)
                                     aaCTCCGTGTACTGTCGGTTATTTGACCACCGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTg                                                                                                               <  1:3237828/150‑1 (MQ=255)
                                     aaCTCCGTGTACTGTCGGTTATTTGACCACCGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCCATCACGGCTGGGAATCTCTACTg                                                                                                               >  2:3258007/1‑150 (MQ=255)
                                      aCTCCGTGTACTGTCGGTTATTTGACCACCGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGa                                                                                                              >  1:2239618/1‑150 (MQ=255)
                                       cTCCGTGTACTGTCGGTTATTTGACCACCGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGAt                                                                                                             <  1:4628427/150‑1 (MQ=255)
                                        tCCGTGTACTGTCGGTTATTTGACCACCGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATg                                                                                                            <  1:1328208/150‑1 (MQ=255)
                                           gtgtACTGTCGGTTATTTGACCACCGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCCGAGCGCCCCATTCGTACCATTCCCCGCCCGCAGGAGAAAGACTTTAACCCCCATAACCGCGGGGAATATCTACAGATgcgc                                                                                                         >  1:248448/1‑150 (MQ=255)
                                              tACTGTCGGTTATTTGACCACCGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCTCTTTCGTACCATTCCCCGCCCGCTGGCGACAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGtt                                                                                                      <  2:61967/150‑1 (MQ=255)
                                              tACTGTCGGTTATTTGACCACCGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGtt                                                                                                      <  1:40087/150‑1 (MQ=255)
                                              tACTGTCGGTTATTTGACCACCGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATATCTACTGATGCGCGtt                                                                                                      >  2:2161932/1‑150 (MQ=255)
                                                cTGTCGGTTATTTGACCACCGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTAt                                                                                                    <  1:2356342/150‑1 (MQ=255)
                                                cTGTCGGTTATTTGACCACCGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTAt                                                                                                    >  2:2182257/1‑150 (MQ=255)
                                                  gTCGGTTATTTGACCACCGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCa                                                                                                  >  1:4587445/1‑150 (MQ=255)
                                                     ggTTATTTGACCACCGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGTCTGGGAATCTCTACTGATGCGCGTTATCAGCg                                                                                               >  1:1919265/1‑150 (MQ=255)
                                                      gTTATTTGACCACCGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGa                                                                                              <  2:3078731/150‑1 (MQ=255)
                                                       ttATTTGACCACCGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGAc                                                                                             <  1:1371916/150‑1 (MQ=255)
                                                        tATTTGACCACCGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGAcc                                                                                            <  2:198056/150‑1 (MQ=255)
                                                         aTTTGACCACCGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCt                                                                                           >  2:2886985/1‑150 (MQ=255)
                                                             gACCACCGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCTCCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCg                                                                                       >  2:181178/1‑150 (MQ=255)
                                                               ccaccGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGct                                                                                     >  1:4524391/1‑150 (MQ=255)
                                                               ccaccGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGAGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGct                                                                                     <  2:3576020/150‑1 (MQ=255)
                                                                caccGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGctg                                                                                    <  2:866937/150‑1 (MQ=255)
                                                                 accGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGctgc                                                                                   >  1:1036467/1‑150 (MQ=255)
                                                                  ccGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGctgcg                                                                                  >  1:4213526/1‑150 (MQ=255)
                                                                   cgctcgctTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGctgcgc                                                                                 <  1:388032/150‑1 (MQ=255)
                                                                     ctcgctTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGctgcgcct                                                                               <  1:1510320/150‑1 (MQ=255)
                                                                     ctcgctTACTCCCGAGCGGCTTTTTATATTTAGCACTCAGCAAGAGCGCAGCTTTCGTAACAATACCCGCCCGCTGGCGAAAGACTTTCAACCCGATCACGCCTGGGAATCTCTAAAGTAGATCGTTATAAGCGACCTAAAGctgcgcct                                                                               >  2:1935977/1‑150 (MQ=255)
                                                                      tcgctTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGctgcgcctg                                                                              >  2:4402507/1‑150 (MQ=255)
                                                                      tcgctTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGctgcgcctg                                                                              <  2:2406634/150‑1 (MQ=255)
                                                                      tcgctTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGctgcgcctg                                                                              <  1:2207720/150‑1 (MQ=255)
                                                                       cgctTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGAGAAAGACTTTCACCCCGATCACGGCTGGGAATATCTACTGATGCGCGTTATCAGCGACCTACCGctgcgcctgc                                                                             <  2:2563764/150‑1 (MQ=255)
                                                                         ctTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGAGAAAGACTTTCACCCCGATCACGGCTGGGAAAATATACTGATGCGCGTTATCAGCGACCAACCGctgcgccagcgc                                                                           >  1:1589395/1‑150 (MQ=255)
                                                                         ctTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTAGTAACATTACCCGCCCGCTGGAGAAAGAATATCACCCCGATCACGGCGGGGAATATATACTGATGCGCGTTATCAGCGACCTACCGctgcgcctgcgc                                                                           >  1:2525918/1‑150 (MQ=255)
                                                                          ttACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGAGAAAGACTTTCACCCCGATCACGGCTGGGAATATCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCt                                                                          <  2:49497/150‑1 (MQ=255)
                                                                           tACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTAACATTACCCGCCCGCTGGCGAAAGACTTTAACCCCGATCACGGCTGGGAATATCTACTGATGCGCGTTATCAGCGACCAACCGCAGCGCCAGCGCTg                                                                         >  1:2724254/1‑150 (MQ=255)
                                                                            aCTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTgg                                                                        >  2:1303260/1‑150 (MQ=255)
                                                                            aCTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTgg                                                                        <  1:4392630/150‑1 (MQ=255)
                                                                              tCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGAGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCa                                                                      <  2:4047820/150‑1 (MQ=255)
                                                                               cccGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAg                                                                     <  1:3136291/150‑1 (MQ=255)
                                                                                ccGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGAGAAAGACTTTCACCCCGATCACGGCTGGGAATATCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGa                                                                    <  2:1641155/150‑1 (MQ=255)
                                                                                  gAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAAt                                                                  <  1:1810358/150‑1 (MQ=255)
                                                                                   aGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATa                                                                 >  2:670059/1‑150 (MQ=255)
                                                                                      ggCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATaaaa                                                              >  2:3135151/1‑150 (MQ=255)
                                                                                      ggCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATaaaa                                                              >  1:2273403/1‑150 (MQ=255)
                                                                                      ggCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATaaaa                                                              >  2:553620/1‑150 (MQ=255)
                                                                                       gCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAg                                                             >  1:1712531/1‑150 (MQ=255)
                                                                                       gCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAg                                                             >  1:3654667/1‑150 (MQ=255)
                                                                                        cTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGc                                                            >  2:1449022/1‑150 (MQ=255)
                                                                                        cTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGc                                                            >  2:364425/1‑150 (MQ=255)
                                                                                         tttttATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGAGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGcc                                                           <  2:3291952/150‑1 (MQ=255)
                                                                                             tATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCAGTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGa                                                       <  1:838613/150‑1 (MQ=255)
                                                                                               tCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACa                                                     <  2:488851/150‑1 (MQ=255)
                                                                                                cTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACAt                                                    <  1:680519/150‑1 (MQ=255)
                                                                                                  ttAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATcc                                                  >  2:1172956/1‑150 (MQ=255)
                                                                                                    aGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCcat                                                <  1:899245/150‑1 (MQ=255)
                                                                                                          cAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCcattacatt                                          <  1:2716267/150‑1 (MQ=255)
                                                                                                          cAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCcattacatt                                          <  1:3377700/150‑1 (MQ=255)
                                                                                                          cAGCCAGAGCGCCGCTTTCGTACCAGTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCcattacatt                                          <  2:3064859/150‑1 (MQ=255)
                                                                                                           aGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCcattacatta                                         <  2:1486117/150‑1 (MQ=255)
                                                                                                           aGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGAGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCcattacatta                                         >  2:3708127/1‑150 (MQ=255)
                                                                                                            gCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACattat                                        >  2:2631535/1‑150 (MQ=255)
                                                                                                               agagCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGTCGTTCACCCCGATCACGGCTGGGACTCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCTGAATAAAAGCCGTGACATCCATTACATTATTCg                                     >  2:2755811/1‑150 (MQ=255)
                                                                                                               agagCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCg                                     >  1:1209086/1‑150 (MQ=255)
                                                                                                               agagCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCg                                     >  2:3490832/1‑150 (MQ=255)
                                                                                                                 agCGCCGCTTTCGTACCATTCCCCGCCCGCTGGAGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGcc                                   >  2:3859529/1‑150 (MQ=255)
                                                                                                                   cgccgcTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGGGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCAt                                 <  2:1876037/150‑1 (MQ=255)
                                                                                                                   cgccgcTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCAt                                 >  2:4139039/1‑150 (MQ=255)
                                                                                                                    gccgcTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATc                                >  1:2620043/1‑150 (MQ=255)
                                                                                                                    gccgcTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATc                                >  2:816133/1‑150 (MQ=255)
                                                                                                                     ccgcTTTCGTACCAGTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCt                               <  2:2300083/150‑1 (MQ=255)
                                                                                                                      cgcTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTg                              >  1:30200/1‑150 (MQ=255)
                                                                                                                           tCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCGACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGa                         <  1:3967139/150‑1 (MQ=255)
                                                                                                                            cGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGaa                        >  1:3776773/1‑150 (MQ=255)
                                                                                                                            cGTACCAGTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGaa                        <  2:820645/150‑1 (MQ=255)
                                                                                                                               aCCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACg                     <  2:3392762/150‑1 (MQ=255)
                                                                                                                                ccATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGc                    >  2:2005976/1‑150 (MQ=255)
                                                                                                                                 cATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCt                   >  1:2066250/1‑150 (MQ=255)
                                                                                                                                 cATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCt                   >  1:558697/1‑150 (MQ=255)
                                                                                                                                      cccgcccgCTGGAGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGga              >  1:4399113/1‑150 (MQ=255)
                                                                                                                                       ccgcccgcTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGgac             >  1:966539/1‑150 (MQ=255)
                                                                                                                                        cgcccgcTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGgaca            <  1:1157837/150‑1 (MQ=255)
                                                                                                                                          cccgcTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGgacaga          <  2:3766699/150‑1 (MQ=255)
                                                                                                                                             gcTGGCGAAAGACTTTCACCCCGATTACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACaa       >  2:4281112/1‑150 (MQ=255)
                                                                                                                                               tGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACAAcc     <  2:1290206/150‑1 (MQ=255)
                                                                                                                                                ggCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACAACCt    >  1:3229404/1‑150 (MQ=255)
                                                                                                                                                 gCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACAACCTc   <  2:4304162/150‑1 (MQ=255)
                                                                                                                                                  cGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACAACCTcg  <  2:2657012/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGTTATTTGACCACCGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACAACCTCG  >  CP009273/4208796‑4209091

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