Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,510,179 G→T 26.0% A115A (GCG→GCT patD → gamma‑aminobutyraldehyde dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,510,1790GT26.0% 558.1 / 182.1 246A115A (GCG→GCTpatDgamma‑aminobutyraldehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (105/77);  new base T (25/39);  total (130/116)
Fisher's exact test for biased strand distribution p-value = 1.31e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 257 total aligned reads displayed.

GTCTGCTGAAACTGGCTGATGTTATCGAAGAAAATGGTCAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGCTGATGATGGCCGCGTGGAAACTTGCTCCGGCGCTGGC  >  CP009273/1510031‑1510325
                                                                                                                                                    |                                                                                                                                                  
gTCTGCTGAAACTGGCTGATGTTATCGAAGAAAATGGTCAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGgcgc                                                                                                                                                   <  1:695504/150‑1 (MQ=255)
gTCTGCTGAAACTGGCTGATGTTATCGAAGAAAATGGTCAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGgcgc                                                                                                                                                   <  1:4544356/150‑1 (MQ=255)
  ctgctgAAACTGGCTGATGTTATCGAAGAAAATGGTCAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGCTTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGgctcgc                                                                                                                                                 >  2:249259/1‑150 (MQ=255)
   tgctgAAACTGGCTGATGTTATCGAAGAAAATGGTCAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCTCGCt                                                                                                                                                >  2:4807433/1‑150 (MQ=255)
    gctgAAACTGGCTGATGTTATCGAAGAAAATGGTCAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTg                                                                                                                                               >  1:739825/1‑150 (MQ=255)
      tgAAACTGGCTGATGTTATCGAAGAAAATGGTCAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCTCGCTGTc                                                                                                                                             <  2:849959/150‑1 (MQ=255)
       gAAACTGGCTGATGTTATCGAAGAAAATGGTCAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCt                                                                                                                                            <  1:1805228/150‑1 (MQ=255)
       gAAACTGGCTGATGTTATCGAAGAAAATGGTCAGGTTTTTGCCGAACTGGAGTCCCGTAAGTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCTCGCTGTCt                                                                                                                                            <  1:4384365/150‑1 (MQ=255)
         aaCTGGCTGATGTTATCGAAGAAAATGGTCAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCTCGCTGTCTGa                                                                                                                                          >  1:1402896/1‑150 (MQ=255)
          aCTGGCTGATGTTATCGAAGAAAATGGTCAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGaa                                                                                                                                         <  2:3574042/150‑1 (MQ=255)
            tGGCTGATGTTATCGAAGAAAATGGTCAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATg                                                                                                                                       >  1:1301231/1‑150 (MQ=255)
            tGGCTGATGTTATCGAAGAAAATGGTCAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATg                                                                                                                                       >  1:2424432/1‑150 (MQ=255)
                   tGTTATCGAAGAAAATGGTCAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTggc                                                                                                                                >  1:3925539/1‑150 (MQ=255)
                    gTTATCGAAGAAAATGGTCAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTggcg                                                                                                                               <  2:4378132/150‑1 (MQ=255)
                        tCGAAGAAAATGGTCAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAg                                                                                                                           <  1:2221440/150‑1 (MQ=255)
                        tCGAAGAAAATGGTCAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAg                                                                                                                           <  1:3296352/150‑1 (MQ=255)
                            agaaAATGGTCAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGa                                                                                                                       <  1:3611618/150‑1 (MQ=255)
                             gaaAATGGTCAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGaa                                                                                                                      >  2:2895231/1‑150 (MQ=255)
                             gaaAATGGTCAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGaa                                                                                                                      >  2:1539718/1‑150 (MQ=255)
                              aaaaTGGTCAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAat                                                                                                                     >  2:1290580/1‑150 (MQ=255)
                               aaaTGGTCAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAata                                                                                                                    <  2:1983760/150‑1 (MQ=255)
                                aaTGGTCAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCTCGCTGTCTGAATGGTCTGGCGGCAGGTGAatat                                                                                                                   <  1:2114374/150‑1 (MQ=255)
                                aaTGGTCAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAatat                                                                                                                   >  2:3522616/1‑150 (MQ=255)
                                     tCAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGa                                                                                                              <  1:2732642/150‑1 (MQ=255)
                                     tCAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGa                                                                                                              <  2:1473314/150‑1 (MQ=255)
                                          tttttGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCTCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTc                                                                                                         >  1:3034740/1‑150 (MQ=255)
                                           ttttGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCTCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCa                                                                                                        <  2:450397/150‑1 (MQ=255)
                                              tGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATAc                                                                                                     <  1:791637/150‑1 (MQ=255)
                                              tGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATAc                                                                                                     >  2:1881054/1‑150 (MQ=255)
                                              tGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATAc                                                                                                     >  2:177053/1‑150 (MQ=255)
                                               gCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACt                                                                                                    <  1:2015020/150‑1 (MQ=255)
                                               gCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACt                                                                                                    <  2:1693513/150‑1 (MQ=255)
                                               gCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACt                                                                                                    <  1:1539718/150‑1 (MQ=255)
                                                ccGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCTCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACtt                                                                                                   >  1:2800415/1‑150 (MQ=255)
                                                ccGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACtt                                                                                                   >  1:2621271/1‑150 (MQ=255)
                                                ccGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACtt                                                                                                   >  2:2990893/1‑150 (MQ=255)
                                                 cGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTc                                                                                                  >  2:1013235/1‑150 (MQ=255)
                                                  gAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCg                                                                                                 <  2:167803/150‑1 (MQ=255)
                                                  gAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCg                                                                                                 <  2:2597242/150‑1 (MQ=255)
                                                  gAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCg                                                                                                 <  2:2597473/150‑1 (MQ=255)
                                                  gAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCg                                                                                                 >  1:1361190/1‑150 (MQ=255)
                                                    aCTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCTCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCgat                                                                                               <  1:3459722/150‑1 (MQ=255)
                                                    aCTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCgat                                                                                               <  2:3561797/150‑1 (MQ=255)
                                                      tGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCgatga                                                                                             >  1:1084720/1‑150 (MQ=255)
                                                      tGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCgatga                                                                                             <  2:3282071/150‑1 (MQ=255)
                                                          gTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCTCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCg                                                                                         <  2:4217429/150‑1 (MQ=255)
                                                             ccGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTcg                                                                                      >  2:4696986/1‑150 (MQ=255)
                                                               gTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCTCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTcgcg                                                                                    <  1:1761730/150‑1 (MQ=255)
                                                                tAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCTCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGa                                                                                   <  1:40548/150‑1 (MQ=255)
                                                                tAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCTCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGa                                                                                   <  1:3726315/150‑1 (MQ=255)
                                                                tAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGa                                                                                   <  1:883829/150‑1 (MQ=255)
                                                                  aTTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATc                                                                                 <  1:3863560/150‑1 (MQ=255)
                                                                   ttGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCTCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATcc                                                                                <  2:626097/150‑1 (MQ=255)
                                                                     gtgGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGt                                                                              >  2:4854641/1‑150 (MQ=255)
                                                                      tgGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATATTGAAGGTCATACTTCGATGATCCGTCGCGATACGtt                                                                             >  1:1269719/1‑150 (MQ=255)
                                                                       ggCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCTCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTg                                                                            <  2:2721827/150‑1 (MQ=255)
                                                                       ggCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTg                                                                            <  2:2663378/150‑1 (MQ=255)
                                                                         cAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCTCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTggg                                                                          <  1:525418/150‑1 (MQ=255)
                                                                          aaaCCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTgggg                                                                         <  2:3211707/150‑1 (MQ=255)
                                                                           aaCCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCTCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTggggg                                                                        >  1:293574/1‑150 (MQ=255)
                                                                            aCCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCTCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGt                                                                       >  2:4016767/1‑150 (MQ=255)
                                                                              cGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCTCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCg                                                                     <  1:1910603/150‑1 (MQ=255)
                                                                               gCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCTCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGt                                                                    >  1:4351426/1‑150 (MQ=255)
                                                                                  gCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGc                                                                 <  2:688363/150‑1 (MQ=255)
                                                                                     tAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACGTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTc                                                              >  1:4581981/1‑150 (MQ=255)
                                                                                       gTGCGTTCAATGATGACATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCTCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTa                                                            <  2:1508537/150‑1 (MQ=255)
                                                                                        tGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCTCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTAt                                                           <  2:753704/150‑1 (MQ=255)
                                                                                             tCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTGGGGGGTCGTGGCTTCTATCGCAc                                                      >  2:1789372/1‑150 (MQ=255)
                                                                                              cAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCTCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCAca                                                     <  2:2439667/150‑2 (MQ=255)
                                                                                              cAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCAcc                                                     >  2:4583365/1‑150 (MQ=255)
                                                                                              cAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCAcc                                                     >  2:2300086/1‑150 (MQ=255)
                                                                                                atgatgAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGt                                                   <  2:571863/150‑1 (MQ=255)
                                                                                                 tgatgaAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTg                                                  <  2:3135997/150‑1 (MQ=255)
                                                                                                     gaAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAAt                                              >  2:1844639/1‑150 (MQ=255)
                                                                                                        aTCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTAt                                           >  2:3365125/1‑150 (MQ=255)
                                                                                                        aTCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTAt                                           >  1:4280170/1‑150 (MQ=255)
                                                                                                          cccGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATcc                                         >  2:2410193/1‑150 (MQ=255)
                                                                                                            cGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCTCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGc                                       <  1:382613/150‑1 (MQ=255)
                                                                                                             ggCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGCt                                      <  2:2526517/150‑1 (MQ=255)
                                                                                                               cGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGCtga                                    >  2:1102852/1‑150 (MQ=255)
                                                                                                                  ttGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGCtgatga                                 >  1:1400003/1‑150 (MQ=255)
                                                                                                                   tGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGCTgatgat                                >  2:2175574/1‑150 (MQ=255)
                                                                                                                    gTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGCTGatgatg                               >  2:969967/1‑150 (MQ=255)
                                                                                                                    gTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGCTGatgatg                               >  2:367175/1‑150 (MQ=255)
                                                                                                                     tCGATGTTTTTCGCTTTTTCGCGGGTGCGGCTCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGACTTATCCGCTGATGATgg                              >  2:760121/1‑150 (MQ=255)
                                                                                                                     tCGATGTTATTCGCTTTTACGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATACGTCGCGATCCGTTGTGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGCTGATGATgg                              >  2:11817/1‑150 (MQ=255)
                                                                                                                       gATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGCTGATGATGGcc                            >  2:4374313/1‑150 (MQ=255)
                                                                                                                        aTGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGCTGATGATGGCcg                           <  2:3923791/150‑1 (MQ=255)
                                                                                                                               tcgctttttcgcGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGCTGATGATGGCCGCGTGGaa                    <  1:3141675/150‑1 (MQ=255)
                                                                                                                                  ctttttcgcGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGCTGATGATGGCCGCGTGGAAACt                 <  2:2089257/150‑1 (MQ=255)
                                                                                                                                   tttttcgcGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGCTGATGATGGCCGCGTGGAAACtt                <  1:4726850/150‑1 (MQ=255)
                                                                                                                                    ttttcgcGGGTGCGGCTCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGCTGATGATGGCCGCGTGGAAACTTg               >  2:3437482/1‑150 (MQ=255)
                                                                                                                                    ttttcgcGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGCTGATGATGGCCGCGTGGAAACTTg               >  1:3177872/1‑150 (MQ=255)
                                                                                                                                        cgcgGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGCTGATGATGGCCGCGTGGAAACTTGCTcc           >  2:2337809/1‑150 (MQ=255)
                                                                                                                                          cgGGTGCGGCTCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGCTGATGATGGCCGCGTGGAAACTTGCTCCgg         >  1:450326/1‑150 (MQ=255)
                                                                                                                                          cgGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGCTGATGATGGCCGCGTGGAAACTTGCTCCgg         <  1:4143256/150‑1 (MQ=255)
                                                                                                                                          cgGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGCTGATGATGGCCGCGTGGAAACTTGCTCCgg         <  1:4861225/150‑1 (MQ=255)
                                                                                                                                          cgGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGCTGATGATGGCCGCGTGGAAACTTGCTCCgg         >  2:4243787/1‑150 (MQ=255)
                                                                                                                                               gcggcgCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGCTGATGATGGCCGCGTGGAAACTTGCTCCGGCGCTg    >  1:3103609/1‑150 (MQ=255)
                                                                                                                                                 ggcgCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGCTGATGATGGCCGCGTGGAAACTTGCTCCGGCGCTggc  >  2:2367687/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GTCTGCTGAAACTGGCTGATGTTATCGAAGAAAATGGTCAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGCTGATGATGGCCGCGTGGAAACTTGCTCCGGCGCTGGC  >  CP009273/1510031‑1510325

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: