Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,401,206 G→T 18.7% Q283K (CAG→AAG)  ydaM ← diguanylate cyclase, csgD regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,401,2060GT18.7% 635.2 / 104.4 219Q283K (CAG→AAG) ydaMdiguanylate cyclase, csgD regulator
Reads supporting (aligned to +/- strand):  ref base G (80/97);  new base T (21/20);  total (102/117)
Fisher's exact test for biased strand distribution p-value = 4.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 219 total aligned reads displayed.

CCAACGAAACACCAAATCGCCTTTGCGAGCGCAACTTTCGAGGGTGCGGGCGAGGGCGCATAACACCTCATCACCTTTAGAATGCCCATAGAGATCGTTAATGTGTTTAAAACGATCGGTGTCGATGAGCAACAAGCTGTAATCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAAT  >  CP009273/1401060‑1401350
                                                                                                                                                  |                                                                                                                                                
  aaCGAAACACCAAATCGCCTTTGCGAGCGCAACTTTCGAGGGTGCGGGCGAGGGCGCATAACACCTCATCACCTTTAGAATGCCCATAGAGATCGTTAATGTGTTTAAAACGATCGGTGTCGATGAGCAACAAGCTGTAATCCTGAGCGa                                                                                                                                             >  1:1297705/1‑150 (MQ=255)
    cGAAACACCAAATCGCCTTTGCGAGCGCAACTTTCGAGGGTGCGGGCGAGGGCGCATAACACCTCATCACCTTTAGAATGCCCATAGAGATCGTTAATGTGTTTAAAACGATCGGTGTCGATGAGCAACAAGCTGTAATCCTGAGCGATg                                                                                                                                           >  2:415478/1‑150 (MQ=255)
    cGAAACACCAAATCGCCTTTGCGAGCGCAACTTTCGAGGGTGCGGGCGAGGGCGCATAACACCTCATCACCTTTAGAATGCCCATAGAGATCGTTAATGTGTTTAAAACGATCGGTGTCGATGAGCAACAAGCTGTAATCCTGAGCGATg                                                                                                                                           >  1:774911/1‑150 (MQ=255)
     gAAACACCAAATCGCCTTTGCGAGCGCAACTTTCGAGGGTGCGGGCGAGGGCGCATAACACCTCATCACCTTTAGAATGCCCATAGAGATCGTTAATGTGTTTAAAACGATCGGTGTCGATGAGCAACAAGCTGTAATCCTGAGCGATgg                                                                                                                                          >  1:3861360/1‑150 (MQ=255)
      aaaCACCAAATCGCCTTTGCGAGCGCAACTTTCGAGGGTGCGGGCGAGGGCGCATAACACCTCATCACCTTTAGAATGCCCATAGAGATCGTTAATGTGTTTAAAACGATCGGTGTCGATGAGCAACAAGCTGTAATCCTGAGCGATGGc                                                                                                                                         >  2:2507632/1‑150 (MQ=255)
      aaaCACCAAATCGCCTTTGCGAGCGCAACTTTCGAGGGTGCGGGCGAGGGCGCATAACACCTCATCACCTTTAGAATGCCCATAGAGATCGTTAATGTGTTTAAAACGATCGGTGTCGATGAGCAACAAGCTGTAATCCTGAGCGATGGc                                                                                                                                         >  2:2121825/1‑150 (MQ=255)
       aaCACCAAATCGCCTTTGCGAGCGCAACTTTCGAGGGTGCGGGCGAGGGCGCATAACACCTCATCACCTTTAGAATGCCCATAGAGATCGTTAATGTGTTTAAAACGATCGGTGTCGATGAGCAACAAGCTGTAATCCTGAGCGATGGCg                                                                                                                                        <  1:2290064/150‑1 (MQ=255)
            cAAATCGCCTTTGCGAGCGCAACTTTCGAGGGTGCGGGCGAGGGCGCATAACACCTCATCACCTTTAGAATGCCCATAGAGATCGTTAATGTGTTTAAAACGATCGGTGTCGATGAGCAACAAGCTGTAATCCTTAGCGATGGCGAGAtg                                                                                                                                   >  2:2350156/1‑150 (MQ=255)
              aaTCGCCTTTGCGAGCGCAACTTTCGAGGGTGCGGGCGAGGGCGCATAACACCTCATCACCTTTAGAATGCCCATAGAGATCGTTAATGTGTTTAAAACGATCGGTGTCGATGAGCAACAAGCTGTAATCCTTAGCGATGGCGAGAtgct                                                                                                                                 <  2:3063137/150‑1 (MQ=255)
                   ccTTTGCGAGCGCAACTTTCGAGGGTGCGGGCGAGGGCGCATAACACCTCATCACCTTTAGAATGCCCATAGAGATCGTTAATGTGTTTAAAACGATCGGTGTCGATGAGCAACAAGCTGTAATCCTGAGCGATGGCGAGATGCTGCAtt                                                                                                                            <  1:4017885/150‑1 (MQ=255)
                   ccTTTGCGAGCGCAACTTTCGAGGGTGCGGGCGAGGGCGCATAACACCTCATCACCTTTAGAATGCCCATAGAGATCGTTAATGTGTTTAAAACGATCGGTGTCGATGAGCAACAAGCTGTAATCCTGAGCGATGGCGAGATGCTGCAtt                                                                                                                            <  1:1557170/150‑1 (MQ=255)
                    cTTTGCGAGCGCAACTTTCGAGGGTGCGGGCGAGGGCGCATAACACCTCATCACCTTTAGAATGCCCATAGAGATCGTTAATGTGTTTAAAACGATCGGTGTCGATGAGCAACAAGCTGTAATCCTTAGCGATGGCGAGATGCTGCAttt                                                                                                                           >  2:4472174/1‑150 (MQ=255)
                         cGAGCGCAACTTTCGAGGGTGCGGGCGAGGGCGCATAACACCTCATCACCTTTAGAATGCCCATAGAGATCGTTAATGTGTTTAAAACGATCGGTGTCGATGAGCAACAAGCTGTAATCCTTAGCGATGGCGAGATGCTGCATTTGGCCt                                                                                                                      <  2:3806377/150‑1 (MQ=255)
                         cGAGCGCAACTTTCGAGGGTGCGGGCGAGGGCGCATAACACCTCATCACCTTTAGAATGCCCATAGAGATCGTTAATGTGTTTAAAACGATCGGTGTCGATGAGCAACAAGCTGTAATCCTTAGCGATGGCGAGATGCTGCATTTGGCCt                                                                                                                      <  2:1351715/150‑1 (MQ=255)
                         cGAGCGCAACTTTCGAGGGTGCGGGCGAGGGCGCATAACACCTCATCACCTTTAGAATGCCCATAGAGATCGTTAATGTGTTTAAAACGATCGGTGTCGATGAGCAACAAGCTGTAATCCTGAGCGATGGCGAGATGCTGCATTTGGCCt                                                                                                                      <  2:4061638/150‑1 (MQ=255)
                          gAGCGCAACTTTCGAGGGTGCGGGCGAGGGCGCATAACACCTCATCACCTTTAGAATGCCCATAGAGATCGTTAATGTGTTTAAAACGATCGGTGTCGATGAGCAACAAGCTGTAATCCTTAGCGATGGCGAGATGCTGCATTTGGCCTg                                                                                                                     <  1:3511158/150‑1 (MQ=255)
                          gAGCGCAACTTTCGAGGGTGCGGGCGAGGGCGCATAACACCTCATCACCTTTAGAATGCCCATAGAGATCGTTAATGTGTTTAAAACGATCGGTGTCGATGAGCAACAAGCTGTAATCCTTAGCGATGGCGAGATGCTGCATTTGGCCTg                                                                                                                     <  1:3942289/150‑1 (MQ=255)
                          gAGCGCAACTTTCGAGGGTGCGGGCGAGGGCGCATAACACCTCATCACCTTTAGAATGCCCATAGAGATCGTTAATGTGTTTAAAACGATCGGTGTCGATGAGCAACAAGCTGTAATCCTGAGCGATGGCGAGATGCTGCATTTGGCCTg                                                                                                                     <  1:3039688/150‑1 (MQ=255)
                          gAGCGCAACTTTCGAGGGTGCGGGCGAGGGCGCATAACACCTCATCACCTTTAGAATGCCCATAGAGATCGTTAATGTGTTTAAAACGATCGGTGTCGATGAGCAACAAGCTGTAATCCTGAGCGATGGCGAGATGCTGCATTTGGCCTg                                                                                                                     >  2:4737504/1‑150 (MQ=255)
                           aGCGCAACTTTCGAGGGTGCGGGCGAGGGCGCATAACACCTCATCACCTTTAGAATGCCCATAGAGATCGTTAATGTGTTTAAAACGATCGGTGTCGATGAGCAACAAGCTGTAATCCTTAGCGATGGCGAGATGCTGCATTTGGCCTgg                                                                                                                    >  1:702281/1‑150 (MQ=255)
                             cgcAACTTTCGAGGGTGCGGGCGAGGGCGCATAACACCTCATCACCTTTAGAATGCCCATAGAGATCGTTAATGTGTTTAAAACGATCGGTGTCGATGAGCAACAAGCTGTAATCCTTAGCGATGGCGAGATGCTGCATTTGGCCTGGtt                                                                                                                  <  1:1414394/150‑1 (MQ=255)
                             cgcAACTTTCGAGGGTGCGGGCGAGGGCGCATAACACCTCATCACCTTTAGAATGCCCATAGAGATCGTTAATGTGTTTAAAACGATCGGTGTCGATGAGCAACAAGCTGTAATCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGtt                                                                                                                  >  1:3668682/1‑150 (MQ=255)
                                aaCTTTCGAGGGTGCGGGCGAGGGCGCATAACACCTCATCACCTTTAGAATGCCCATAGAGATCGTTAATGTGTTTAAAACGATCGGTGTCGATGAGCAACAAGCTGTAATCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCg                                                                                                               <  1:4107786/150‑1 (MQ=255)
                                    ttCGAGGGTGCGGGCGAGGGCGCATAACACCTCATCACCTTTAGAATGCCCATAGAGATCGTTAATGTGTTTAAAACGATCGGTGTCGATGAGCAACAAGCTGTAATCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAAt                                                                                                           <  2:4702800/150‑1 (MQ=255)
                                    ttCGAGGGTGCGGGCGAGGGCGCATAACACCTCATCACCTTTAGAATGCCCATAGAGATCGTTAATGTGTTTAAAACGATCGGTGTCGATGAGCAACAAGCTGTAATCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAAt                                                                                                           <  1:3496624/150‑1 (MQ=255)
                                      cGAGGGTGCGGGCGAGGGCGCATAACACCTCATCACCTTTAGAATGCCCATAGAGATCGTTAATGTGTTTAAAACGATCGGTGTCGATGAGCAACAAGCTGTAATCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAAtgt                                                                                                         >  1:2602022/1‑150 (MQ=255)
                                          ggTGCGGGCGAGGGCGCATAACACCTCATCACCTTTAGAATGCCCATAGAGATCGTTAATGTGTTTAAAACGATCGGTGTCGATGAGCAACAAGCTGTAATCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATa                                                                                                     >  2:3345231/1‑150 (MQ=255)
                                               gggCGAGGGCGCATAACACCTCATCACCTTTAGAATGCCCATAGAGATCGTTAATGTGTTTAAAACGATCGGTGTCGATGAGCAACAAGCTGTAATCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACt                                                                                                >  1:4679379/1‑150 (MQ=255)
                                                ggCGAGGGCGCATAACACCTCATCACCTTTAGAATGCCCATAGAGATCGTTAATGTGTTTAAAACGATCGGTGTCGATGAGCAACAAGCTGTAATCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTg                                                                                               >  1:2611253/1‑150 (MQ=255)
                                                  cGAGGGCGCATAACACCTCATCACCTTTAGAATGCCCATAGAGATCGTTAATGTGTTTAAAACGATCGGTGTCGATGAGCAACAAGCTGTAATCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTc                                                                                             <  1:511000/150‑1 (MQ=255)
                                                     gggCGCATAACACCTCATCACCTTTAGAATGCCCATAGAGATCGTTAATGTGTTTAAAACGATCGGTGTCGATGAGCAACAAGCTGTAATCCTTAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGcc                                                                                          >  1:2213061/1‑150 (MQ=255)
                                                     gggCGCATAACACCTCATCACCTTTAGAATGCCCATAGAGATCGTTAATGTGTTTAAAACGATCGGTGTCGATGAGCAACAAGCTGTAATCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGcc                                                                                          >  1:4711288/1‑150 (MQ=255)
                                                          cATAACACCTCATCACCTTTAGAATGCCCATAGAGATCGTTAATGTGTTTAAAACGATCGGTGTCGATGAGCAACAAGCTGTAATCCTTAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTc                                                                                     >  2:3673508/1‑150 (MQ=255)
                                                          cATAACACCTCATCACCTTTAGAATGCCCATAGAGATCGTTAATGTGTTTAAAACGATCGGTGTCGATGAGCAACAAGCTGTAATCCTTAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTc                                                                                     >  2:350573/1‑150 (MQ=255)
                                                           aTAACACCTCATCACCTTTAGAATGCCCATAGAGATCGTTAATGTGTTTAAAACGATCGGTGTCGATGAGCAACAAGCTGTAATCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCa                                                                                    <  2:3633279/150‑1 (MQ=255)
                                                           aTAACACCTCATCACCTTTAGAATGCCCATAGAGATCGTTAATGTGTTTAAAACGATCGGTGTCGATGAGCAACAAGCTGTAATCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCa                                                                                    <  2:1171220/150‑1 (MQ=255)
                                                             aaCACCTCATCACCTTTAGAATGCCCATAGAGATCGTTAATGTGTTTAAAACGATCGGTGTCGATGAGCAACAAGCTGTAATCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGt                                                                                  >  2:787808/1‑150 (MQ=255)
                                                              acacCTCATCACCTTTAGAATGCCCATAGAGATCGTTAATGTGTTTAAAACGATCGGTGTCGATGAGCAACAAGCTGTAATCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTa                                                                                 <  1:209359/150‑1 (MQ=255)
                                                              acacCTCATCACCTTTAGAATGCCCATAGAGATCGTTAATGTGTTTAAAACGATCGGTGTCGATGAGCAACAAGCTGTAATCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTa                                                                                 >  2:1093938/1‑150 (MQ=255)
                                                                 ccTCATCACCTTTAGAATGCCCATAGAGATCGTTAATGTGTTTAAAACGATCGGTGTCGATGAGCAACAAGCTGTAATCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATc                                                                              <  2:717456/150‑1 (MQ=255)
                                                                  cTCATCACCTTTAGAATGCCCATAGAGATCGTTAATGTGTTTAAAACGATCGGTGTCGATGAGCAACAAGCTGTAATCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATcc                                                                             >  1:3883162/1‑150 (MQ=255)
                                                                     atcaCCTTTAGAATGCCCATAGAGATCGTTAATGTGTTTAAAACGATCGGTGTCGATGAGCAACAAGCTGTAATCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGt                                                                          >  2:290545/1‑150 (MQ=255)
                                                                         ccTTTAGAATGCCCATAGAGATCGTTAATGTGTTTAAAACGATCGGTGTCGATGAGCAACAAGCTGTAATCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATc                                                                      >  1:3992910/1‑150 (MQ=255)
                                                                          cTTTAGAATGCCCATAGAGATCGTTAATGTGTTTAAAACGATCGGTGTCGATGAGCAACAAGCTGTAATCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATcg                                                                     >  1:4816946/1‑150 (MQ=255)
                                                                           tttAGAATGCCCATAGAGATCGTTAATGTGTTTAAAACGATCGGTGTCGATGAGCAACAAGCTGTAATCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATcgc                                                                    >  1:2825146/1‑150 (MQ=255)
                                                                           tttAGAATGCCCATAGAGATCGTTAATGTGTTTAAAACGATCGGTGTCGATGAGCAACAAGCTGTAATCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATcgc                                                                    >  1:3289844/1‑150 (MQ=255)
                                                                              aGAATGCCCATAGAGATCGTTAATGTGTTTAAAACGATCGGTGTCGATGAGCAACAAGCTGTAATCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTc                                                                 >  1:1018772/1‑150 (MQ=255)
                                                                              aGAATGCCCATAGAGATCGTTAATGTGTTTAAAACGATCGGTGTCGATGAGCAACAAGCTGTAATCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTc                                                                 <  2:2755314/150‑1 (MQ=255)
                                                                               gAATGCCCATAGAGATCGTTAATGTGTTTAAAACGATCGGTGTCGATGAGCAACAAGCTGTAATCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCa                                                                >  2:2820915/1‑150 (MQ=255)
                                                                                aaTGCCCATAGAGATCGTTAATGTGTTTAAAACGATCGGTGTCGATGAGCAACAAGCTGTAATCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCAt                                                               <  2:3271416/150‑1 (MQ=255)
                                                                                aaTGCCCATAGAGATCGTTAATGTGTTTAAAACGATCGGTGTCGATGAGCAACAAGCTGTAATCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCAt                                                               <  2:4169617/150‑1 (MQ=255)
                                                                                 aTGCCCATAGAGATCGTTAATGTGTTTAAAACGATCGGTGTCGATGAGCAACAAGCTGTAATCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATg                                                              <  1:2681832/150‑1 (MQ=255)
                                                                                 aTGCCCATAGAGATCGTTAATGTGTTTAAAACGATCGGTGTCGATGAGCAACAAGCTGTAATCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATg                                                              <  1:2681903/150‑1 (MQ=255)
                                                                                  tGCCCATAGAGATCGTTAATGTGTTTAAAACGATCGGTGTCGATGAGCAACAAGCTGTAATCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATgg                                                             >  1:2835270/1‑150 (MQ=255)
                                                                                     ccATAGAGCTCGTTAATGTGTTTAAAACGATCGGTGTCGATGAGCAACCAGCTGTAATCCTGCGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGa                                                          <  2:1675856/150‑1 (MQ=255)
                                                                                      cATAGAGATCGTTAATGTGTTTAAAACGATCGGTGTCGATGAGCAACAAGCTGTAATCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGag                                                         <  2:3704533/150‑1 (MQ=255)
                                                                                       aTAGAGATCGTTAATGTGTTTAAAACGATCGGTGTCGATGAGCAACAAGCTGTAATCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGagc                                                        >  2:1786614/1‑150 (MQ=255)
                                                                                        tAGAGATCGTTAATGTGTTTAAAACGATCGGTGTCGATGAGCAACAAGCTGTAATCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGagca                                                       <  2:1899486/150‑1 (MQ=255)
                                                                                         agagATCGTTAATGTGTTTAAAACGATCGGTGTCGATGAGCAACAAGCTGTAATCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGagcag                                                      <  2:4146549/150‑1 (MQ=255)
                                                                                          gagaTCGTTAATGTGTTTAAAACGATCGGTGTCGATGAGCAACAAGCTGTAATCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGagcagc                                                     >  1:4010486/1‑150 (MQ=255)
                                                                                            gaTCGTTAATGTGTTTAAAACGATCGGTGTCGATGAGCAACAAGCTGTAATCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCAt                                                   >  2:430181/1‑150 (MQ=255)
                                                                                            gaTCGTTAATGTGTTTAAAACGATCGGTGTCGATGAGCAACAAGCTGTAATCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCAt                                                   <  1:4284375/150‑1 (MQ=255)
                                                                                            gaTCGTTAATGTGTTTAAAACGATCGGTGTCGATGAGCAACAAGCTGTAAACCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCAt                                                   >  1:2719061/1‑150 (MQ=255)
                                                                                             aTCGTTAATGTGTTTAAAACGATCGGTGTCGATGAGCAACAAGCTGTAATCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATg                                                  <  2:1062877/150‑1 (MQ=255)
                                                                                              tCGTTAATGTGTTTAAAACGATCGGTGTCGATGAGCAACAAGCTGTAATCCTGCGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGt                                                 <  2:390677/150‑1 (MQ=255)
                                                                                                 ttAATGTGTTTAAAACGATCGGTGTCGATGAGCAACAAGCTGTAATCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTcc                                              >  1:1560904/1‑150 (MQ=255)
                                                                                                  tAATGTGTTTAAAACGATCGGTGTCGATGAGCAACAAGCTGTAATCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTACTGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCa                                             >  2:328178/1‑150 (MQ=255)
                                                                                                  tAATGTGTTTAAAACGATCGGTGTCGATGAGCAACAAGCTGTAATCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCa                                             >  1:2594345/1‑150 (MQ=255)
                                                                                                  tAATGTGTTTAAAACGATCGGTGTCGATGAGCAACAAGCTGTAATCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCa                                             <  1:1093938/150‑1 (MQ=255)
                                                                                                     tgtgTTTAAAACGATCGGTGTCGATGAGCAACAAGCTGTAATCCTTAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAgct                                          <  1:2114416/150‑1 (MQ=255)
                                                                                                     tgtgTTTAAAACGATCGGTGTCGATGAGCAACAAGCTGTAATCCTTAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAgct                                          <  1:928516/150‑1 (MQ=255)
                                                                                                      gtgtTTAAAACGATCGGTGTCGATGAGCAACAAGCTGTAATCCTTAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAgctg                                         >  2:607005/1‑150 (MQ=255)
                                                                                                       tgtTTAAAACGATCGGTGTCGATGAGCAACAAGCTGTAATCCTTAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAgctgc                                        >  1:4797326/1‑150 (MQ=255)
                                                                                                       tgtTTAAAACGATCGGTGTCGATGAGCAACAAGCTGTAATCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAgctgc                                        <  2:4416398/150‑1 (MQ=255)
                                                                                                            aaaaCGATCGGTGTCGATGAGCAACAAGCTGTAATCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGctcctc                                   >  2:1447684/1‑150 (MQ=255)
                                                                                                             aaaCGATCGGTGTCGATGAGCAACAAGCTGTAATCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCtcctcc                                  <  1:4737504/150‑1 (MQ=255)
                                                                                                              aaCGATCGGTGTCGATGAGCAACAAGCTGTAATCCTTAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCa                                 >  1:301205/1‑150 (MQ=255)
                                                                                                              aaCGATCGGTGTCGATGAGCAACAAGCTGTAATCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCa                                 <  2:1746830/150‑1 (MQ=255)
                                                                                                              aaCGATCGGTGTCGATGAGCAACAAGCTGTAATCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCa                                 <  2:4341958/150‑1 (MQ=255)
                                                                                                              aaCGATCGGTGTCGATGAGCAACAAGCTGTAATCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCa                                 <  2:1746931/150‑1 (MQ=255)
                                                                                                                 gATCGGTGTCGATGAGCAACAAGCTGTAATCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGcc                              >  1:1390672/1‑150 (MQ=255)
                                                                                                                  aTCGGTGTCGATGAGCAACAAGCTGTAATCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCg                             >  2:2528725/1‑150 (MQ=255)
                                                                                                                    cGGTGTCGATGAGCAACAAGCTGTAATCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGtt                           <  2:3469865/150‑1 (MQ=255)
                                                                                                                     ggTGTCGATGAGCAACAAGCTGTAATCCTTAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGttt                          <  2:1890836/150‑1 (MQ=255)
                                                                                                                      gtgtCGATGAGCAACAAGCTGTAATCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGtttt                         <  2:4046884/150‑1 (MQ=255)
                                                                                                                       tgtCGATGAGCAACAAGCTGTAATCCTTAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGttttt                        <  1:1971708/150‑1 (MQ=255)
                                                                                                                        gtCGATGAGCAACAAGCTGTAATCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTg                       >  1:2036778/1‑150 (MQ=255)
                                                                                                                          cGATGAGCAACAAGCTGTAATCCTGAGCGAGGGCGAGATGCTGCATTTGGCCTGGGTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCt                     <  2:2066064/150‑1 (MQ=255)
                                                                                                                               aGCAACAAGCTGTAATCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTa                >  1:4374996/1‑150 (MQ=255)
                                                                                                                                  aactaGCTGTAATCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAata             <  1:4330547/150‑1 (MQ=255)
                                                                                                                                  aacaaGCTGTAATCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAata             <  1:616321/150‑1 (MQ=255)
                                                                                                                                  aacaaGCTGTAATCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTGTTTGCGAAGTAata             >  1:3393772/1‑150 (MQ=255)
                                                                                                                                   acaaGCTGTAATCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAatat            <  2:3125934/150‑1 (MQ=255)
                                                                                                                                     aaGCTGTAATCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTGTTTTGTCAGTAATATCa          >  1:3706115/1‑150 (MQ=255)
                                                                                                                                      aGCTGTAATCCTTAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCAt         >  1:2968997/1‑150 (MQ=255)
                                                                                                                                       gCTGTAATCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTTCTCAGTAATATCATg        >  2:2933051/1‑150 (MQ=255)
                                                                                                                                       gCTGTAATCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATg        <  2:2973869/150‑1 (MQ=255)
                                                                                                                                       gCTGTAATCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATg        <  1:116676/150‑1 (MQ=255)
                                                                                                                                            aaTCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATTCACaa   >  1:196569/1‑150 (MQ=255)
                                                                                                                                            aaTCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTGTTTTCTCAGTAATATAATGCACaa   >  1:197699/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
CCAACGAAACACCAAATCGCCTTTGCGAGCGCAACTTTCGAGGGTGCGGGCGAGGGCGCATAACACCTCATCACCTTTAGAATGCCCATAGAGATCGTTAATGTGTTTAAAACGATCGGTGTCGATGAGCAACAAGCTGTAATCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAAT  >  CP009273/1401060‑1401350

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: