Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,510,854 G→T 100% A340A (GCG→GCT)  patD → gamma‑aminobutyraldehyde dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,510,8540GT96.3% 574.5 / 12.0 190A340A (GCG→GCT) patDgamma‑aminobutyraldehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (4/3);  new base T (95/88);  total (99/91)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 190 total aligned reads displayed.

TTTACGATACGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGCCCT  >  CP009273/1510709‑1510999
                                                                                                                                                 |                                                                                                                                                 
tttACGATACGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTAc                                                                                                                                                 >  2:681815/1‑148 (MQ=255)
  tACGATACGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAg                                                                                                                                               <  2:3162311/148‑1 (MQ=255)
       tACGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGcac                                                                                                                                          <  1:2520721/148‑1 (MQ=255)
        aCGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGcaca                                                                                                                                         >  2:244941/1‑148 (MQ=255)
        aCGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGcaca                                                                                                                                         >  2:623993/1‑148 (MQ=255)
             ggtggAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCaaa                                                                                                                                    >  1:215631/1‑148 (MQ=255)
               tggAAAAACTGGGTGCTGCGGTGTCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGt                                                                                                                                  <  1:3419921/148‑1 (MQ=255)
                ggAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTg                                                                                                                                 <  1:2270301/148‑1 (MQ=255)
                 gAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGa                                                                                                                                >  2:776203/1‑148 (MQ=255)
                   aaaaCTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATc                                                                                                                              <  1:1444958/148‑1 (MQ=255)
                      aCTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACt                                                                                                                           >  2:2257680/1‑148 (MQ=255)
                      aCTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACt                                                                                                                           >  2:2121688/1‑148 (MQ=255)
                      aCTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGCAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACt                                                                                                                           >  2:2784012/1‑148 (MQ=255)
                       cTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTg                                                                                                                          <  2:2748003/148‑1 (MQ=255)
                         gggTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGc                                                                                                                        >  2:2266482/1‑148 (MQ=255)
                          ggTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCg                                                                                                                       >  1:785377/1‑148 (MQ=255)
                          ggTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCg                                                                                                                       >  1:148497/1‑148 (MQ=255)
                            tgctgcGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGt                                                                                                                     >  1:289186/1‑148 (MQ=255)
                            tgctgcGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGt                                                                                                                     >  1:2838639/1‑148 (MQ=255)
                               tgcGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGaa                                                                                                                  <  2:3477363/148‑1 (MQ=255)
                                gcGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGaaa                                                                                                                 <  2:30677/148‑1 (MQ=255)
                                gcGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGaaa                                                                                                                 <  2:917079/148‑1 (MQ=255)
                                 cGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGaaaa                                                                                                                <  1:884193/148‑1 (MQ=255)
                                   gTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAgc                                                                                                              >  1:2527008/1‑148 (MQ=255)
                                    tGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAgcg                                                                                                             <  2:237324/148‑1 (MQ=255)
                                       cAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCaa                                                                                                          <  2:2854680/148‑1 (MQ=255)
                                         aCGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAgg                                                                                                        <  2:1790413/148‑1 (MQ=255)
                                          cGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAggg                                                                                                       <  2:1684392/148‑1 (MQ=255)
                                            ttAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTa                                                                                                     <  2:2545773/148‑1 (MQ=255)
                                            ttAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTa                                                                                                     <  1:295147/148‑1 (MQ=255)
                                                aaTCTGGTGCGCCATATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCTGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATgg                                                                                                 <  1:1646174/148‑1 (MQ=255)
                                                  tCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCt                                                                                               <  2:2481393/148‑1 (MQ=255)
                                                  tCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCt                                                                                               >  1:3404995/1‑148 (MQ=255)
                                                   cTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTa                                                                                              >  2:2374062/1‑148 (MQ=255)
                                                     ggTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTAtt                                                                                            >  2:671762/1‑148 (MQ=255)
                                                        gcgcCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACt                                                                                         <  1:2960380/148‑1 (MQ=255)
                                                        gcgcCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACt                                                                                         <  2:752182/148‑1 (MQ=255)
                                                         cgcCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTa                                                                                        <  2:2762114/148‑1 (MQ=255)
                                                              gATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATgcgc                                                                                   <  1:2980832/148‑1 (MQ=255)
                                                                tGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCg                                                                                 >  1:909796/1‑148 (MQ=255)
                                                                 gACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGa                                                                                >  2:714528/1‑148 (MQ=255)
                                                                   cGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACg                                                                              >  1:2368200/1‑148 (MQ=255)
                                                                       tCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACgctgc                                                                          >  2:2302680/1‑148 (MQ=255)
                                                                        cTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACAACAAAGTGATCACTGGCGGTGAAAAGCGCAAGCGTAATGGCTATTACTATGCGCCGACgctgct                                                                         <  1:1109352/148‑1 (MQ=255)
                                                                          aCGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTgctgg                                                                       >  2:2087285/1‑148 (MQ=255)
                                                                              aGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGctggctgg                                                                   >  2:2842773/1‑148 (MQ=255)
                                                                                cTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGgcg                                                                 >  1:734057/1‑148 (MQ=255)
                                                                                  tGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCa                                                               <  1:3284131/148‑1 (MQ=255)
                                                                                      ccTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTAc                                                           <  1:1251742/148‑1 (MQ=255)
                                                                                       cTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTCGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACa                                                          >  2:2059362/1‑148 (MQ=255)
                                                                                       cTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACa                                                          >  2:1334774/1‑148 (MQ=255)
                                                                                       cTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACa                                                          >  2:1334781/1‑148 (MQ=255)
                                                                                       cTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACa                                                          >  1:3160563/1‑148 (MQ=255)
                                                                                          tAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGa                                                       <  1:1398933/148‑1 (MQ=255)
                                                                                             gCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGa                                                    <  2:2236698/148‑1 (MQ=255)
                                                                                             gCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGa                                                    <  1:2514471/148‑1 (MQ=255)
                                                                                             gCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGa                                                    >  1:2697438/1‑148 (MQ=255)
                                                                                                cGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGc                                                 >  2:1604057/1‑148 (MQ=255)
                                                                                                  cTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCa                                               <  1:2412777/148‑1 (MQ=255)
                                                                                                        aaTCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGc                                         >  2:520935/2‑148 (MQ=255)
                                                                                                               aaCGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAaga                                  >  1:1990310/1‑148 (MQ=255)
                                                                                                                aCGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAagag                                 >  1:1968839/1‑148 (MQ=255)
                                                                                                                aCGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAagag                                 <  2:2627886/148‑1 (MQ=255)
                                                                                                                aCGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAagag                                 <  1:1298936/148‑1 (MQ=255)
                                                                                                                aCGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAagag                                 >  1:1479244/1‑148 (MQ=255)
                                                                                                                aCGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAagag                                 <  1:634778/148‑1 (MQ=255)
                                                                                                                 cgcgTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAgg                                >  2:374755/1‑148 (MQ=255)
                                                                                                                  gcgTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGt                               <  2:2591424/148‑1 (MQ=255)
                                                                                                                    gTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTAt                             <  2:1111332/148‑1 (MQ=255)
                                                                                                                    gTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTAt                             <  2:2977055/148‑1 (MQ=255)
                                                                                                                    gTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACGGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTAt                             <  2:2052864/148‑1 (MQ=255)
                                                                                                                     tCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTAtt                            >  1:1906656/1‑148 (MQ=255)
                                                                                                                     tCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTAtt                            <  2:2887524/148‑1 (MQ=255)
                                                                                                                     tCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTAtt                            <  2:481843/148‑1 (MQ=255)
                                                                                                                     tCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTAtt                            <  2:481797/148‑1 (MQ=255)
                                                                                                                      cGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCCTTACAGGACGATGCCATCGTGCAAAAAGAGGTAttt                           >  2:1962364/1‑148 (MQ=255)
                                                                                                                        gCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTgg                         >  1:3032511/1‑148 (MQ=255)
                                                                                                                         cAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGt                        >  1:649131/1‑148 (MQ=255)
                                                                                                                         cAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGt                        >  1:649350/1‑148 (MQ=255)
                                                                                                                           aGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTcc                      >  1:3051711/1‑148 (MQ=255)
                                                                                                                           aGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTcc                      >  2:2018565/1‑148 (MQ=255)
                                                                                                                            ggCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCa                     >  1:588383/1‑148 (MQ=255)
                                                                                                                            ggCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCa                     >  1:1683704/1‑148 (MQ=255)
                                                                                                                              cAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCagt                   <  1:2607183/148‑1 (MQ=255)
                                                                                                                                 tAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCagtagt                >  1:1611603/1‑148 (MQ=255)
                                                                                                                                    aagaGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTagtgag             >  1:876780/1‑148 (MQ=255)
                                                                                                                                    aagaGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTagtgag             <  1:112134/148‑1 (MQ=255)
                                                                                                                                     agagGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCAAGTagtgagt            >  1:2721616/1‑148 (MQ=255)
                                                                                                                                      gagGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTagtgagtg           <  1:242933/148‑1 (MQ=255)
                                                                                                                                       agGCGACAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAgtgt          <  2:2158991/148‑1 (MQ=255)
                                                                                                                                       agGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAgtgt          >  2:2566466/1‑148 (MQ=255)
                                                                                                                                       agGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAgtgt          <  1:2787958/148‑1 (MQ=255)
                                                                                                                                       agGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAgtgt          >  2:2119714/1‑148 (MQ=255)
                                                                                                                                         gCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTa        >  1:270183/1‑148 (MQ=255)
                                                                                                                                            aaaGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGGATTTGGTCCAGTAGTGAGTGTTACGc     >  2:125020/1‑148 (MQ=255)
                                                                                                                                             aaGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGcc    >  1:712674/1‑148 (MQ=255)
                                                                                                                                              aGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGccc   <  2:1783848/148‑1 (MQ=255)
                                                                                                                                              aGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGccc   <  2:1683128/148‑1 (MQ=255)
                                                                                                                                              aGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGccc   <  2:145465/148‑1 (MQ=255)
                                                                                                                                               gCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGCCCt  <  2:2527008/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
TTTACGATACGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGCCCT  >  CP009273/1510709‑1510999

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: ‑