Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,208,975 C→A 7.1% 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,9750CA 7.1% 1192.6 / 40.3 322S151Y (TCT→TAT) aceKisocitrate dehydrogenase kinase/phosphatase
Reads supporting (aligned to +/- strand):  ref base C (124/175);  new base A (12/11);  total (136/186)
Fisher's exact test for biased strand distribution p-value = 3.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

Only 100 of 332 total aligned reads displayed.

TTTTTTAACTCCGTGTACTGTCGGTTATTTGACCACCGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACAACCTCGCGGAAAGTCATTTACAGGTGGCGAACGAACTGT  >  CP009273/4208827‑4209124
                                                                                                                                                    |                                                                                                                                                     
ttttttAACTCCGTGTACTGTCGGTTATTTGACCACCGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGAGAAAGACTTTCACCCCGATCACGGCTGGGAAtct                                                                                                                                                      >  2:3149379/1‑150 (MQ=255)
 tttttAACTCCGTGTACTGTCGGTTATTTGACCACCGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAAtctc                                                                                                                                                     <  1:3717755/150‑1 (MQ=255)
 tttttAACTCCGTGTACTGTCGGTTATTTGACCACCGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAAtctc                                                                                                                                                     <  1:1878669/150‑1 (MQ=255)
   tttAACTCCGTGTACTGTCGGTTATTTGACCACCGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTa                                                                                                                                                   <  2:1466188/150‑1 (MQ=255)
        cTCCGTGTACTGTCGGTTATTTGACCACCGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGAGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGAt                                                                                                                                              <  1:1174840/150‑1 (MQ=255)
         tCCGTGTACTGTCGGTTATTTGACCACCGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGATAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATg                                                                                                                                             <  2:2811640/150‑1 (MQ=255)
           cGTGTACTGTCGGTTATTTGACCACCGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATgcg                                                                                                                                           >  1:2605869/1‑150 (MQ=255)
               tACTGTCGGTTATTTGACCACCGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGAGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGtt                                                                                                                                       >  2:2213638/1‑150 (MQ=255)
                aCTGTCGGTTATTTTACCACCGCTCGCTTACTCCCGAGCGGCTTTTTATATTTAGCACTCAGCCAGAGCGCCGCGTTCGTACCATTACCCCCCCCCAGGCGAAAAACATTAACCCCCAAAACCGCTGGGAAACTCAACAGATGCGCGTTa                                                                                                                                      >  1:467660/1‑150 (MQ=255)
                aCTGTCGGTTATTTGACCACCGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTa                                                                                                                                      <  1:4037773/150‑1 (MQ=255)
                 cTGTCGGTTATTTGACCACCGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTAt                                                                                                                                     >  1:3494336/1‑150 (MQ=255)
                 cTGTCGGTTATTTGACCACCGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTAt                                                                                                                                     >  2:680047/1‑150 (MQ=255)
                   gTCGGTTATTTGACCACCGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGAGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCa                                                                                                                                   <  2:3132857/150‑1 (MQ=255)
                      ggTTATTTGGCCACCGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGAGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCg                                                                                                                                >  1:579396/1‑150 (MQ=255)
                        ttATTTGACCACCGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGAc                                                                                                                              >  2:4038569/1‑150 (MQ=255)
                              gACCACCGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCg                                                                                                                        <  2:425807/150‑1 (MQ=255)
                              gACCACCGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGAGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCg                                                                                                                        <  1:1682864/150‑1 (MQ=255)
                              gACCACCGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCAGTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCg                                                                                                                        <  2:425867/150‑1 (MQ=255)
                                 caccGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGctg                                                                                                                     >  2:867089/1‑150 (MQ=255)
                                     gctcgctTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGAGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGctgcgcc                                                                                                                 >  1:592313/1‑150 (MQ=255)
                                      ctcgctTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGAGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGctgcgcct                                                                                                                >  2:1535428/1‑150 (MQ=255)
                                        cgctTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCAGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGctgcgcctgc                                                                                                              >  1:3306624/1‑150 (MQ=255)
                                          ctTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGAGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGctgcgcctgcgc                                                                                                            >  2:364863/1‑150 (MQ=255)
                                          ctTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGAGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGctgcgcctgcgc                                                                                                            >  2:4229237/1‑150 (MQ=255)
                                           ttACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGAGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCt                                                                                                           >  1:284503/1‑150 (MQ=255)
                                               tCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCa                                                                                                       <  1:4182107/150‑1 (MQ=255)
                                                  cGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGAGAAAGACTTTCACCCCGATCACGGCTGGGAATATCTACTGATGCGCGTTATCAGCGACCTCCCGCTGCGCCTGCGCTGGCAGaa                                                                                                    >  2:649050/1‑150 (MQ=255)
                                                    aGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATa                                                                                                  >  2:3741442/1‑150 (MQ=255)
                                                      cGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATaaa                                                                                                <  1:657448/150‑1 (MQ=255)
                                                      cGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGAGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATaaa                                                                                                <  2:87214/150‑1 (MQ=255)
                                                        gCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAg                                                                                              >  1:3730029/1‑150 (MQ=255)
                                                         cattttATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGATAGACTTTCACCCCGATCTCGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGc                                                                                             <  1:1114810/148‑1 (MQ=255)
                                                         cattttATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGAGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGc                                                                                             <  1:1900965/148‑1 (MQ=255)
                                                         cTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGAGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGc                                                                                             <  2:4341999/150‑1 (MQ=255)
                                                           ttttATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGAGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCg                                                                                           >  1:72417/1‑150 (MQ=255)
                                                           ttttATCTTGAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCg                                                                                           <  2:1880268/150‑1 (MQ=255)
                                                             ttATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTg                                                                                         <  2:2353847/150‑1 (MQ=255)
                                                             ttATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTg                                                                                         <  1:1319723/150‑1 (MQ=255)
                                                             ttATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGAGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTg                                                                                         >  2:2964807/1‑150 (MQ=255)
                                                                  tttAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGAGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATc                                                                                    <  2:2123931/150‑1 (MQ=255)
                                                                     aGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGAGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCcat                                                                                 <  2:3335087/150‑1 (MQ=255)
                                                                      gCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGAGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCcatt                                                                                <  2:3423383/150‑1 (MQ=255)
                                                                       ctctCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGAGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCcatta                                                                               <  1:1604792/150‑1 (MQ=255)
                                                                        tctcAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGAGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCcattac                                                                              <  2:540722/150‑1 (MQ=255)
                                                                            aGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGACTAAAAGCCGTGACATCcattacatta                                                                          >  1:2528666/1‑150 (MQ=255)
                                                                             gCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACattat                                                                         <  1:2006010/150‑1 (MQ=255)
                                                                              ccAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACattatt                                                                        <  1:2269426/150‑1 (MQ=255)
                                                                               cAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGAGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTc                                                                       <  2:4005744/150‑1 (MQ=255)
                                                                                agagCGCCGCTTTCGTACCATTCCCCGCCCGCTGGAGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCg                                                                      >  1:3397974/1‑150 (MQ=255)
                                                                                agagCGCCGCTTTCGTACCATTCCCCGCCCGCTGGAGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCg                                                                      >  1:4042116/1‑150 (MQ=255)
                                                                                agagCGCCGCTTTCGTACCATTCCCCGCCCGCTGGAGAAAGACTTTCACCCCGATCACGGCTGGGAATATCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCg                                                                      <  1:67133/150‑1 (MQ=255)
                                                                                 gagCGCCGCTTTCGTACCATTCCCCGCCCGCTGGAGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGc                                                                     <  1:3536040/150‑1 (MQ=255)
                                                                                        gcTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGa                                                              <  2:3522151/150‑1 (MQ=255)
                                                                                        gcTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGa                                                              <  2:3522190/150‑1 (MQ=255)
                                                                                        gcTTTCGTACCATTCCCCGCCCGCTGGAGAAAGACTTTCACCCCGATCACGGCTGGGAATATCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGa                                                              <  1:3212133/150‑1 (MQ=255)
                                                                                         cTTTCGTACCATTCCCCTCCCGCTGGAGAAAGACTTTCACCCCGATCACGTCTGGGAATCTCTACTGATGCGCGTTATCAGAGACCTACCGCTGCGCCTGAGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGAc                                                             <  2:4339044/150‑1 (MQ=255)
                                                                                         cTTTCGTACCATTCCCCGCCCGCTGGAGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGAc                                                             <  2:2093449/150‑1 (MQ=255)
                                                                                         cTTTCGTACCATTCCCCGCCCGCTGGAGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGAc                                                             <  1:4348645/150‑1 (MQ=255)
                                                                                              gTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGaaa                                                        >  2:2956346/1‑150 (MQ=255)
                                                                                              gTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGaaa                                                        >  1:1006278/1‑150 (MQ=255)
                                                                                               tACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAAc                                                       <  2:2852162/150‑1 (MQ=255)
                                                                                               tACCATTCCCCGCCCGCTGGAGAAAGACTTTCACCCCGATCACGGCTGGGAATATCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAAc                                                       <  2:1708221/150‑1 (MQ=255)
                                                                                                aCCATTCCCCGCCCGCTGGAGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACg                                                      <  1:2976852/150‑1 (MQ=255)
                                                                                                  cATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCt                                                    >  2:41679/1‑150 (MQ=255)
                                                                                                     tCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTggg                                                 <  2:1866295/150‑1 (MQ=255)
                                                                                                     tCCCCGCCCGCTGGAGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGGTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGACTAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTggg                                                 >  2:3272513/1‑150 (MQ=255)
                                                                                                      ccccGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTgggg                                                >  1:2383081/1‑150 (MQ=255)
                                                                                                      ccccGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTgggg                                                >  1:2383099/1‑150 (MQ=255)
                                                                                                       cccgcccgCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGga                                               >  2:1269601/1‑150 (MQ=255)
                                                                                                           cccgcTGGAGAAAGACTTTCACCCCGATCACGGCTGGGAATATCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGgacaga                                           >  2:1801830/1‑150 (MQ=255)
                                                                                                             cgcTGGAGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGgacagaca                                         <  2:2504537/150‑1 (MQ=255)
                                                                                                             cgcTGGAGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGgacagaca                                         <  2:4212520/150‑1 (MQ=255)
                                                                                                             cgcTGGAGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGgacagaca                                         >  2:3119411/1‑150 (MQ=255)
                                                                                                             cgcTGGAGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGgacagaca                                         <  2:205406/150‑1 (MQ=255)
                                                                                                                 ggCCAAACACATCCACCCCCAACACCCCTCGGAATCACCCCTGCTGCGCGAACCCACCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACAACCt                                     <  2:3761242/150‑1 (MQ=255)
                                                                                                                   cGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACAACCTcg                                   <  2:4461152/150‑1 (MQ=255)
                                                                                                                   agAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACAACCTcg                                   <  2:3098327/149‑1 (MQ=255)
                                                                                                                    gAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACAACCTcgc                                  <  1:1505305/150‑1 (MQ=255)
                                                                                                                     aaaGACTTTCACCCCGATCACTGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACAACCTcgcg                                 >  1:1233109/1‑150 (MQ=255)
                                                                                                                     aaaGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACAACCTcgcg                                 >  2:4390242/1‑150 (MQ=255)
                                                                                                                            ttCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACAACCTCGCGGAAAGTc                          <  1:205671/150‑1 (MQ=255)
                                                                                                                             tCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACAACCTCGCGGAAAGTCa                         <  1:1548715/150‑1 (MQ=255)
                                                                                                                              cACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACAACCTCGCGGAAAGTCAt                        <  2:2281919/150‑1 (MQ=255)
                                                                                                                              cACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACAACCTCGCGGAAAGTCAt                        <  2:467660/150‑1 (MQ=255)
                                                                                                                               aCCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACAACCTCGCGGAAAGTCAtt                       <  1:1806902/150‑1 (MQ=255)
                                                                                                                               aCCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACAACCTCGCGGAAAGTCAtt                       <  1:3816319/150‑1 (MQ=255)
                                                                                                                                   cGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACAACCTCGCGGAAAGTCATTTACa                   >  1:2669025/1‑150 (MQ=255)
                                                                                                                                    gATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACAACCTCGCGGAAAGTCATTTACAg                  >  2:409177/1‑150 (MQ=255)
                                                                                                                                    gATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACAACCTCGCGGAAAGTCATTTACAg                  <  1:4085517/150‑1 (MQ=255)
                                                                                                                                    gATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACAACCTCGCGGAAAGTCATTTACAg                  >  1:3608344/1‑150 (MQ=255)
                                                                                                                                     aTCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACAACCTCGCGGAAAGTCATTTACAgg                 <  2:2399630/150‑1 (MQ=255)
                                                                                                                                         cGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACAACCTCGCGGAAAGTCATTTACAGGTGGc             >  2:1005073/1‑150 (MQ=255)
                                                                                                                                          ggCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACAACCTCGCGGAAAGTCATTTACAGGTGGcg            <  2:3494281/150‑1 (MQ=255)
                                                                                                                                          ggCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACAACCTCGCGGAAAGTCATTTACAGGTGGcg            >  2:2162376/1‑150 (MQ=255)
                                                                                                                                           gCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACAACCTCGCGGAAAGTCATTTACAGGTGGcga           <  1:810009/150‑1 (MQ=255)
                                                                                                                                           gCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACAACCTCGCGGAAAGTCATTTACAGGTGGcga           <  1:1271929/150‑1 (MQ=255)
                                                                                                                                            cTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACAACCTCGCGGAAAGTCATTTACAGGTGGcgaa          <  2:2701194/150‑1 (MQ=255)
                                                                                                                                                gAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACAACCTCGCGGAAAGTCATTTACAGGTGGcgaacgaa      >  2:3972524/1‑150 (MQ=255)
                                                                                                                                                  aTCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACAACCTCGCGGAAAGTCATTTACAGGTGGCGAACGAACt    <  1:3916926/150‑1 (MQ=255)
                                                                                                                                                  aTCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACAACCTCGCGGAAAGTCATTTACAGGTGGCGAACGAACt    >  2:392183/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TTTTTTAACTCCGTGTACTGTCGGTTATTTGACCACCGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACAACCTCGCGGAAAGTCATTTACAGGTGGCGAACGAACTGT  >  CP009273/4208827‑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: