Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,510,179 G→T 5.2% 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,1790GT 5.2% 1260.4 / 32.3 326A115A (GCG→GCTpatDgamma‑aminobutyraldehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (144/165);  new base T (10/7);  total (154/172)
Fisher's exact test for biased strand distribution p-value = 4.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01

Only 100 of 349 total aligned reads displayed.

TGTCTGCTGAAACTGGCTGATGTTATCGAAGAAAATGGTCAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGCTGATGATGGCCGCGTGGAAACTTGCTCCGGCGCTGGCGGC  >  CP009273/1510030‑1510328
                                                                                                                                                     |                                                                                                                                                     
tGTCTGCTGAAACTGGCTGATGTTATCGAAGAAAATGGTCAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTgcggcg                                                                                                                                                       <  2:3556905/150‑1 (MQ=255)
   ctgctgAAACTGGCTGATGTTATCGAAGAAAATGGTCAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGgcgcgc                                                                                                                                                    <  2:5663352/150‑1 (MQ=255)
    tgctgAAACTGGCTGATGTTATCGAAGAAAATGGTCAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCt                                                                                                                                                   <  2:1390356/150‑1 (MQ=255)
     gctgAAACTGGCTGATGTTATCGAAGAAAATGGTCAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTg                                                                                                                                                  <  1:4522747/150‑1 (MQ=255)
      ctgAAACTGGCTGATGTTATCGAAGAAAATGGTCAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGt                                                                                                                                                 <  1:3798903/150‑1 (MQ=255)
      ctgAAACTGGCTGATGTTATCGAAGAAAATGGTCAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGt                                                                                                                                                 >  1:604844/1‑150 (MQ=255)
       tgAAACTGGCTGATGTTATCGAAGAAAATGGTCAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTc                                                                                                                                                >  1:3383277/1‑150 (MQ=255)
         aaaCTGGCTGATGTTATCGAAGAAAATGGTCAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTg                                                                                                                                              >  2:2418148/1‑150 (MQ=255)
          aaCTGGCTGATGTTATCGAAGAAAATGGTCAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGa                                                                                                                                             <  2:2846173/150‑1 (MQ=255)
             tGGCTGATGTTATCGAAGAAAATGGTCAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATg                                                                                                                                          <  2:5059227/150‑1 (MQ=255)
               gCTGATGTTATCGAAGAAAATGGTCAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCTCGCTGTCTGAATGGt                                                                                                                                        <  2:4543892/150‑1 (MQ=255)
               gCTGATGTTATCGAAGAAAATGGTCAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGt                                                                                                                                        >  2:4078753/1‑150 (MQ=255)
               gCTGATGTTATCGAAGAAAATGGTCAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGt                                                                                                                                        >  1:1142988/1‑150 (MQ=255)
                cTGATGTTATCGAAGAAAATGGTCAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTc                                                                                                                                       <  2:540050/150‑1 (MQ=255)
                  gATGTTATCGAAGAAAATGGTCAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTg                                                                                                                                     >  2:4887203/1‑150 (MQ=255)
                  gATGTTATCGAAGAAAATGGTCAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTg                                                                                                                                     >  1:246042/1‑150 (MQ=255)
                    tGTTATCGAAGAAAATGGTCAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTggc                                                                                                                                   <  1:200305/150‑1 (MQ=255)
                     gTTATCGAAGAAAATGGTCAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTggcg                                                                                                                                  >  1:5602623/1‑150 (MQ=255)
                      ttATCGAAGAAAATGGTCAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCTCGCTGTCTGAATGGTCTggcgg                                                                                                                                 >  2:4185333/1‑150 (MQ=255)
                      ttATCGAAGAAAATGGTCAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTggcgg                                                                                                                                 <  1:4381925/150‑1 (MQ=255)
                        aTCGAAGAAAATGGTCAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCa                                                                                                                               <  2:835723/150‑1 (MQ=255)
                         tCGAAGAAAATGGTCAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAg                                                                                                                              <  2:726167/150‑1 (MQ=255)
                          cGAAGAAAATGGTCAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCTCGCTGTCTGAATGGTCTGGCGGCAgg                                                                                                                             >  1:1037207/1‑150 (MQ=255)
                           gaagaaAATGGTCAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGt                                                                                                                            <  1:902951/150‑1 (MQ=255)
                              gaaAATGGTCAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGaa                                                                                                                         <  1:1191436/150‑1 (MQ=255)
                              gaaAATGGTCAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGaa                                                                                                                         <  2:2592800/150‑1 (MQ=255)
                               aaaaTGGTCAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAat                                                                                                                        >  1:2322364/1‑150 (MQ=255)
                                aaaTGGTCAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAata                                                                                                                       >  1:2955262/1‑150 (MQ=255)
                                aaaTGGTCAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAata                                                                                                                       <  1:2215927/150‑1 (MQ=255)
                                   tGGTCAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCt                                                                                                                    >  1:1279059/1‑150 (MQ=255)
                                       cAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGCTTGTCGATGTTTTTCGCTTTTTCGAGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGAGAATATCTTgtc                                                                                                                >  2:2969111/1‑148 (MQ=255)
                                       cAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGaa                                                                                                                <  1:1618730/150‑1 (MQ=255)
                                       cAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGaa                                                                                                                <  2:4817455/150‑1 (MQ=255)
                                       cAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGaa                                                                                                                >  2:2315988/1‑150 (MQ=255)
                                          gTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGt                                                                                                             <  1:4168002/150‑1 (MQ=255)
                                                 ccGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACtt                                                                                                      >  2:4666494/1‑150 (MQ=255)
                                                  cGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTc                                                                                                     >  1:4073318/1‑150 (MQ=255)
                                                   gAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGGCATACTTCg                                                                                                    >  2:1872952/1‑150 (MQ=255)
                                                    aaCTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCga                                                                                                   >  2:478017/1‑150 (MQ=255)
                                                      cTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCgatg                                                                                                 <  1:4737566/150‑1 (MQ=255)
                                                      cTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCgatg                                                                                                 <  2:5073946/150‑1 (MQ=255)
                                                        ggAGTCCCGTAAGTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCTCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCgatgat                                                                                               <  2:5228210/150‑1 (MQ=255)
                                                          aGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATcc                                                                                             >  2:4923033/1‑150 (MQ=255)
                                                            tCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGt                                                                                           <  1:2688757/150‑1 (MQ=255)
                                                            tCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGAGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGt                                                                                           <  1:5398057/150‑1 (MQ=255)
                                                             cccGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTc                                                                                          >  2:2749676/1‑150 (MQ=255)
                                                              ccGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTcg                                                                                         <  1:4893232/150‑1 (MQ=255)
                                                              ccGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTcg                                                                                         >  1:1318575/1‑150 (MQ=255)
                                                               cGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCTCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTcgc                                                                                        >  2:3916960/1‑150 (MQ=255)
                                                               cGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTcgc                                                                                        <  2:445140/150‑1 (MQ=255)
                                                                gTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTcgcg                                                                                       >  1:3812761/1‑150 (MQ=255)
                                                                gTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTcgcg                                                                                       >  1:2677570/1‑150 (MQ=255)
                                                                    ttGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATcc                                                                                   <  1:3295697/150‑1 (MQ=255)
                                                                      gtgGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGt                                                                                 >  1:1961173/1‑150 (MQ=255)
                                                                         gCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTgg                                                                              <  1:1715872/150‑1 (MQ=255)
                                                                         gCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTgg                                                                              >  2:3383685/1‑150 (MQ=255)
                                                                           aaaCCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTgggg                                                                            <  1:5291977/150‑1 (MQ=255)
                                                                               cGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCg                                                                        >  1:391580/1‑150 (MQ=255)
                                                                                gCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGt                                                                       <  1:3438025/150‑1 (MQ=255)
                                                                                     aTAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCtt                                                                  >  2:382957/1‑150 (MQ=255)
                                                                                     aTAGTGAGTTCAATGATGAAATACCGGAGATTGTCGATGTTTTTCGCTTTTTAGCGGGTGCGGAGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTAGATTATCCGTCGCGATCAGTTGGGGGTCGAGGCtt                                                                  >  2:1133474/1‑150 (MQ=255)
                                                                                           cGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCg                                                            <  1:5044776/150‑1 (MQ=255)
                                                                                           cGTTCAATGATGAAATCCCGGCGAGTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCg                                                            <  1:3895128/150‑1 (MQ=255)
                                                                                            gTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGc                                                           >  1:2127965/1‑150 (MQ=255)
                                                                                             ttCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCa                                                          <  1:1445720/150‑1 (MQ=255)
                                                                                             ttCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCa                                                          <  1:4509441/150‑1 (MQ=255)
                                                                                              tCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCTCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCAc                                                         <  1:1868707/150‑1 (MQ=255)
                                                                                               cAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCGGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCAcc                                                        <  1:4291176/150‑1 (MQ=255)
                                                                                                     tgaAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGaa                                                  <  2:3759402/150‑1 (MQ=255)
                                                                                                     tgaAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGaa                                                  >  1:2415615/1‑150 (MQ=255)
                                                                                                      gaAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAAt                                                 >  1:3104509/1‑150 (MQ=255)
                                                                                                         aTCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTAt                                              <  2:5538619/150‑1 (MQ=255)
                                                                                                          tCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATc                                             <  2:1634915/150‑1 (MQ=255)
                                                                                                           cccGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATcc                                            >  1:570137/1‑150 (MQ=255)
                                                                                                             cGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGc                                          >  2:3662025/1‑150 (MQ=255)
                                                                                                             cGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGc                                          >  2:4002427/1‑150 (MQ=255)
                                                                                                                  aTTGTCGATGTTTTTCGCTTTTTCGAGGGAGATGCGCTCTGTCTCAAAGGACTGGCGGCAGGTGAATAAAATGAAGGTCATACAACGATGATCCGACGCGATCCGGTGGGGGTCGAGGCTTATATAGCTCCGTGGAATTATCCGCtgatg                                     >  2:4348904/1‑150 (MQ=255)
                                                                                                                   ttGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGCtgatga                                    <  2:2645510/150‑1 (MQ=255)
                                                                                                                     gTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGCTGatgatg                                  <  1:5369434/150‑1 (MQ=255)
                                                                                                                     gTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGCTGatgatg                                  >  1:1442722/1‑150 (MQ=255)
                                                                                                                      tCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGCTGATGATgg                                 >  2:3152375/1‑150 (MQ=255)
                                                                                                                      tCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGCTGATGATgg                                 <  1:3910346/150‑1 (MQ=255)
                                                                                                                        gATGTTTTTCGCTTTTTCGCGGGTGCGGCTCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGCTGATGATGGcc                               >  1:4844496/1‑150 (MQ=255)
                                                                                                                        gATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGCTGATGATGGcc                               >  2:4611599/1‑150 (MQ=255)
                                                                                                                        gATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGCTGATGATGGcc                               >  1:3969861/1‑150 (MQ=255)
                                                                                                                         aTGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGCTGATGATGGCcg                              <  2:1181318/150‑1 (MQ=255)
                                                                                                                                tcgctttttcgcGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGCTGATGATGGCCGCGTGGaa                       <  2:1569205/150‑1 (MQ=255)
                                                                                                                                tcgctttttcgcGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGCTGATGATGGCCGCGTGGaa                       >  2:4946478/1‑150 (MQ=255)
                                                                                                                                   ctttttcgcGGGTGCGGCTCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGCTGATGATGGCCGCGTGGAAACt                    >  1:1919107/1‑150 (MQ=255)
                                                                                                                                    tttttcgcGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGCTGATGATGGCCGCGTGGAAACtt                   <  1:5221873/150‑1 (MQ=255)
                                                                                                                                    tttttcgcGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGCTGATGATGGCCGCGTGGAAACtt                   >  2:615524/1‑150 (MQ=255)
                                                                                                                                           cgGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGCTGATGATGGCCGCGTGGAAACTTGCTCCgg            <  1:349601/150‑1 (MQ=255)
                                                                                                                                           cgGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGCTGATGATGGCCGCGTGGAAACTTGCTCCgg            >  2:2033531/1‑150 (MQ=255)
                                                                                                                                            gggTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGCTGATGATGGCCGCGTGGAAACTTGCTCCGgc           >  2:770795/1‑150 (MQ=255)
                                                                                                                                             ggTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGCTGATGATGGCCGCGTGGAAACTTGCTCCGgcg          >  1:3642581/1‑150 (MQ=255)
                                                                                                                                                  ggcgCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGCTGATGATGGCCGCGTGGAAACTTGCTCCGGCGCTggc     >  2:1038049/1‑150 (MQ=255)
                                                                                                                                                   gcgcgcTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGCTGATGATGGCCGCGTGGAAACTTGCTCCGGCGCTggcg    <  1:1564731/150‑1 (MQ=255)
                                                                                                                                                   gcgcgcTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGCTGATGATGGCCGCGTGGAAACTTGCTCCGGCGCTggcg    <  2:1279059/150‑1 (MQ=255)
                                                                                                                                                   gcgcgcTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGCTGATGATGGCCGCGTGGAAACTTGCTCCGGCGCTggcg    <  2:3383277/150‑1 (MQ=255)
                                                                                                                                                     gcgcTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGCTGATGATGGCCGCGTGGAAACTTGCTCCGGCGCTggcggc  >  1:4211194/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGTCTGCTGAAACTGGCTGATGTTATCGAAGAAAATGGTCAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGCTGATGATGGCCGCGTGGAAACTTGCTCCGGCGCTGGCGGC  >  CP009273/1510030‑1510328

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: