Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,208,975 C→A 12.3% S151Y (TCT→TAT)  aceK → isocitrate dehydrogenase kinase/phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,208,9750CA12.3% 1085.6 / 80.4 316S151Y (TCT→TAT) aceKisocitrate dehydrogenase kinase/phosphatase
Reads supporting (aligned to +/- strand):  ref base C (138/138);  new base A (27/13);  total (165/151)
Fisher's exact test for biased strand distribution p-value = 4.28e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 325 total aligned reads displayed.

CTTTTTTAACTCCGTGTACTGTCGGTTATTTGACCACCGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACAACCTCGCGGAAAGTCATTTACAGGTGGCGAACGAACTGT  >  CP009273/4208826‑4209124
                                                                                                                                                     |                                                                                                                                                     
cTTTTTTAACTCCGTGTACTGTCGGTTATTTGACCACCGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAAtc                                                                                                                                                       <  2:4200175/150‑1 (MQ=255)
   ttttAACTCCGTGTACTGTCGGTTATTTGACCACCGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATctct                                                                                                                                                    <  2:2382019/150‑1 (MQ=255)
        aCTCCGTGTACTGTCGGTTATTTGACCACCGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGa                                                                                                                                               >  1:2239618/1‑150 (MQ=255)
        aCTCCGTGTACTGTCGGTTATTTGACCACCGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGa                                                                                                                                               <  2:3071630/150‑1 (MQ=255)
             gtgtACTGTCGGTTATTTGACCACCGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCCGAGCGCCCCATTCGTACCATTCCCCGCCCGCAGGAGAAAGACTTTAACCCCCATAACCGCGGGGAATATCTACAGATgcgc                                                                                                                                          >  1:248448/1‑150 (MQ=255)
                tACTGTCGGTTATTTGACCACCGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGtt                                                                                                                                       <  1:40087/150‑1 (MQ=255)
                tACTGTCGGTTATTTGACCACCGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGtt                                                                                                                                       <  2:61957/150‑1 (MQ=255)
                     tCGGTTATTTGACCACCGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAg                                                                                                                                  >  1:3923307/1‑150 (MQ=255)
                          tATTTGACCACCGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGAcc                                                                                                                             <  2:198056/150‑1 (MQ=255)
                          tATTTGACCACCGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGAcc                                                                                                                             >  1:122908/1‑150 (MQ=255)
                           aTTTGACCACCGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCt                                                                                                                            >  2:2886833/1‑150 (MQ=255)
                                 ccaccGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCAGTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGct                                                                                                                      <  2:2942170/150‑1 (MQ=255)
                                  caccGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCCGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCGCTACTGATGCGCGTTATCAGCGACCTACCGctg                                                                                                                     >  1:3690982/1‑150 (MQ=255)
                                  caccGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGctg                                                                                                                     <  2:866937/150‑1 (MQ=255)
                                  caccGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGctg                                                                                                                     <  1:986559/150‑1 (MQ=255)
                                     cgctcgctTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGctgcgc                                                                                                                  <  1:388032/150‑1 (MQ=255)
                                     cgctcgctTACTCCCGAGCGGCTTTTTATATTTAGCTCTCAGCCAGAGCGCCGCTTTCGTAACATTCCCCGCCCGCTGGCGAAAGACGTTCACCCCGATACCGGCTGGGAATCTATAATGATTCGCGTTAACAGCGACCTACCGctgcgc                                                                                                                  >  1:751559/1‑150 (MQ=255)
                                       ctcgctTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGctgcgcct                                                                                                                <  2:1465283/150‑1 (MQ=255)
                                         cgctTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGctgcgcctgc                                                                                                              >  2:175924/1‑150 (MQ=255)
                                         cgctTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGctgcgcctgc                                                                                                              <  2:2718105/150‑1 (MQ=255)
                                           ctTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGAGAAAGACTTTCACCCCCATCACGGCTGGGAATATATACTGATGCGCGTTATCAGCGACCAACCGctgcgccggcgc                                                                                                            >  1:1589408/1‑150 (MQ=255)
                                            ttACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGAGAAAGACTTTCACCCCGATCACGGCTGGGAATATCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCt                                                                                                           <  2:49434/150‑1 (MQ=255)
                                             tACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTg                                                                                                          >  1:2724100/1‑150 (MQ=255)
                                             tACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCCCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCGACCGCTGCGCCAGCGCTg                                                                                                          >  1:2724298/1‑150 (MQ=255)
                                              aCTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTgg                                                                                                         <  1:4392630/150‑1 (MQ=255)
                                               cTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGc                                                                                                        <  2:3824246/150‑1 (MQ=255)
                                               cTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCAGTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGc                                                                                                        <  2:1503654/150‑1 (MQ=255)
                                                      gCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATaa                                                                                                 <  2:850694/150‑1 (MQ=255)
                                                        ggCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATaaaa                                                                                               <  1:3109995/150‑1 (MQ=255)
                                                        ggCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATaaaa                                                                                               >  2:553327/1‑150 (MQ=255)
                                                        ggCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATaaaa                                                                                               >  1:4571689/1‑150 (MQ=255)
                                                         gCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAg                                                                                              >  1:3230757/1‑150 (MQ=255)
                                                           tttttATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGcc                                                                                            <  2:3463267/150‑1 (MQ=255)
                                                              ttATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTg                                                                                         >  1:151347/1‑150 (MQ=255)
                                                                    ttAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATcc                                                                                   >  2:1172956/1‑150 (MQ=255)
                                                                    ttAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATcc                                                                                   >  2:2661920/1‑150 (MQ=255)
                                                                    ttAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGAGAAAGACTTTCACCCCGATCACGGCTGGGAATATCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATcc                                                                                   <  1:2636090/150‑1 (MQ=255)
                                                                      aGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCcat                                                                                 <  1:899245/150‑1 (MQ=255)
                                                                            cAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCcattacatt                                                                           <  2:3064855/150‑1 (MQ=255)
                                                                             aGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCcattacatta                                                                          <  2:3720835/150‑1 (MQ=255)
                                                                               ccAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACattatt                                                                        <  2:2190622/150‑1 (MQ=255)
                                                                                 agagCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCg                                                                      >  2:3490859/1‑150 (MQ=255)
                                                                                 agagCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCg                                                                      >  1:2013679/1‑150 (MQ=255)
                                                                                  gagCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGc                                                                     <  2:798759/150‑1 (MQ=255)
                                                                                  gagCGCCGCTTTCGTACCATTCCCCGCCCGCTGGAGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGc                                                                     <  1:1824296/150‑1 (MQ=255)
                                                                                  gagCGCCGCTTTCGTACCATTCCCCGCCCGCTGGAGAAAGACTTTCACCCCGATCACGGCTGGGAATATCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGc                                                                     >  2:2277234/1‑150 (MQ=255)
                                                                                     cgccgcTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGGGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCAt                                                                  <  2:1876037/150‑1 (MQ=255)
                                                                                      gccgcTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATc                                                                 >  2:816191/1‑150 (MQ=255)
                                                                                      gccgcTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATc                                                                 >  1:4630408/1‑150 (MQ=255)
                                                                                           tttCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACg                                                            >  2:1699555/1‑150 (MQ=255)
                                                                                                 aCCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACg                                                      <  2:3392762/150‑1 (MQ=255)
                                                                                                  ccATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGc                                                     >  2:2005976/1‑150 (MQ=255)
                                                                                                   cATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCt                                                    >  1:558705/1‑150 (MQ=255)
                                                                                                   cATTCCCCGCCCGCTGGAGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCt                                                    >  1:1800963/1‑150 (MQ=255)
                                                                                                   cATTCCCCGCCCGCTGGAGAAAGACTTTCACCCCGATCACGGCTGGGAATATCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCt                                                    <  1:2294653/150‑1 (MQ=255)
                                                                                                        cccgcccgCTGGAGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGga                                               >  1:4399113/1‑150 (MQ=255)
                                                                                                          cgcccgcTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGgaca                                             <  1:1157837/150‑1 (MQ=255)
                                                                                                           gcccgcTGGAGAAAGACTTTCACCCCGATCACGGCTGGGAATATCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGgacag                                            <  1:3221423/150‑1 (MQ=255)
                                                                                                            cccgcTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGgacaga                                           <  2:3766699/150‑1 (MQ=255)
                                                                                                             ccgcTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGgacagac                                          >  2:3158965/1‑150 (MQ=255)
                                                                                                               gcTGGAGAAAGACTTTCACCCCGATCACGGCTGGGAATATCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACaa                                        >  2:3286960/1‑150 (MQ=255)
                                                                                                                cTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACAAc                                       <  1:1388604/150‑1 (MQ=255)
                                                                                                                  ggCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACAACCt                                     <  2:2651728/150‑1 (MQ=255)
                                                                                                                  ggCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACAACCt                                     <  2:4425018/150‑1 (MQ=255)
                                                                                                                   gagAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACAACCTc                                    <  1:842846/148‑1 (MQ=255)
                                                                                                                   gagAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACAACCTc                                    <  1:3053087/148‑1 (MQ=255)
                                                                                                                   gCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACAACCTc                                    <  2:4304179/150‑1 (MQ=255)
                                                                                                                   gCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACAACCTc                                    <  1:2126649/150‑1 (MQ=255)
                                                                                                                    agAAAGACTTTCACCCCGATCACGGCTGGGAATATCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACAACCTcg                                   >  2:4012044/2‑150 (MQ=255)
                                                                                                                     gAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACAACCTcgc                                  >  2:1382131/1‑150 (MQ=255)
                                                                                                                        aGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACAACCTCGCGGa                               >  1:337677/1‑150 (MQ=255)
                                                                                                                           cTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACAACCTCGCGGAAAg                            <  1:1332649/150‑1 (MQ=255)
                                                                                                                            tttCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACAACCTCGCGGAAAGt                           <  2:3648877/150‑1 (MQ=255)
                                                                                                                              tCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACAACCTCGCGGAAAGTCa                         <  2:1322079/150‑1 (MQ=255)
                                                                                                                              tCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACAACCTCGCGGAAAGTCa                         >  2:2135396/1‑150 (MQ=255)
                                                                                                                                 ccccGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACAACCTCGCGGAAAGTCAttt                      <  2:2424677/150‑1 (MQ=255)
                                                                                                                                  cccGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACAACCTCGCGGAAAGTCATTTa                     >  1:2362483/1‑150 (MQ=255)
                                                                                                                                     gATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACAACCTCGCGGAAAGTCATTTACAg                  <  2:4066522/150‑1 (MQ=255)
                                                                                                                                      aTCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACAACCTCGCGGAAAGTCATTTACAgg                 >  2:1279430/1‑150 (MQ=255)
                                                                                                                                       tCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACAACCTCGCGGAAAGTCATTTACAGGt                <  2:3319592/150‑1 (MQ=255)
                                                                                                                                       tCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACAACCTCGCGGAAAGTCATTTACAGGt                <  2:669605/150‑1 (MQ=255)
                                                                                                                                       tCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACAACCTCGCGGAAAGTCATTTACAGGt                <  2:670134/150‑1 (MQ=255)
                                                                                                                                       tCACGGCTGGGAATATCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACAACCTCGCGGAAAGTCATTTACAGGt                >  2:4276160/1‑150 (MQ=255)
                                                                                                                                        cACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACAACCTCGCGGAAAGTCATTTACAGGTg               >  2:2116688/1‑150 (MQ=255)
                                                                                                                                        cACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACAACCTCGCGGAAAGACATTTACAGTTg               >  2:1956880/1‑150 (MQ=255)
                                                                                                                                         aCGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACAACCTCGCGGAAAGTCATTTACAGGTgg              >  1:3535992/1‑150 (MQ=255)
                                                                                                                                          cGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACAACCTCGCGGAAAGTCATTTACAGGTGGc             <  2:3620609/150‑1 (MQ=255)
                                                                                                                                          cGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACAACCTCGCGGAAAGTCATTTACAGGTGGc             >  2:4288327/1‑150 (MQ=255)
                                                                                                                                          cGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACAACCTCGCGGAAAGTCATTTACAGGTGGc             <  2:2619582/150‑1 (MQ=255)
                                                                                                                                           ggCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACAACCTCGCGGAAAGTCATTTACAGGTGGcg            <  1:2532163/150‑1 (MQ=255)
                                                                                                                                             cTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACAACCTCGCGGAAAGTCATTTACAGGTGGcgaa          <  1:1034148/150‑1 (MQ=255)
                                                                                                                                              tGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACAACCTCGCGGAAAGTCATTTACAGGTGGcgaac         >  1:2410084/1‑150 (MQ=255)
                                                                                                                                              tGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACAACCTCGCGGAAAGTCATTTACAGGTGGcgaac         <  1:4148979/150‑1 (MQ=255)
                                                                                                                                              tGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACAACCTCGCGGAAAGTCATTTACAGGTGGcgaac         <  1:4212253/150‑1 (MQ=255)
                                                                                                                                                ggAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACAACCTCGCGGAAAGTCATTTACAGGTGGcgaacga       <  1:929971/150‑1 (MQ=255)
                                                                                                                                                 gAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACAACCTCGCGGAAAGTCATTTACAGGTGGcgaacgaa      <  2:452105/150‑1 (MQ=255)
                                                                                                                                                 gAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACAACCTCGCGGAAAGTCATTTACAGGTGGcgaacgaa      <  2:2541718/150‑1 (MQ=255)
                                                                                                                                                  aaTCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACAACCTCGCGGAAAGTCATTTACAGGTGGCgaacgaac     >  1:3391680/1‑150 (MQ=255)
                                                                                                                                                  aaTCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACAACCTCGCGGAAAGTCATTTACAGGTGGCgaacgaac     >  1:601761/1‑150 (MQ=255)
                                                                                                                                                     ctctACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACAACCTCGCGGAAAGTCATTTACAGGTGGCGAACGAACTGt  <  2:2167044/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CTTTTTTAACTCCGTGTACTGTCGGTTATTTGACCACCGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACAACCTCGCGGAAAGTCATTTACAGGTGGCGAACGAACTGT  >  CP009273/4208826‑4209124

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: