Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,510,854 G→T 94.4% 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,8540GT94.4% 887.1 / 38.3 304A340A (GCG→GCTpatDgamma‑aminobutyraldehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (10/7);  new base T (156/130);  total (166/138)
Fisher's exact test for biased strand distribution p-value = 8.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.93e-01

Only 100 of 304 total aligned reads displayed.

CATTTACGATACGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGCCCTT  >  CP009273/1510707‑1511000
                                                                                                                                                   |                                                                                                                                                  
cATTTACGATACGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTAc                                                                                                                                                  >  1:1344677/1‑150 (MQ=255)
 aTTTACGATACGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACa                                                                                                                                                 <  2:4488575/150‑1 (MQ=255)
  tttACGATACGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAg                                                                                                                                                <  2:753481/150‑1 (MQ=255)
  tttACGATACGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAg                                                                                                                                                >  2:2007037/1‑150 (MQ=255)
   ttACGATACGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAgg                                                                                                                                               <  1:454646/150‑1 (MQ=255)
   ttACGATACGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAgg                                                                                                                                               <  2:3450562/150‑1 (MQ=255)
   ttACGATACGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAgg                                                                                                                                               >  1:3857547/1‑150 (MQ=255)
   ttACGATACGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAgg                                                                                                                                               >  1:3521091/1‑150 (MQ=255)
    tACGATACGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAggg                                                                                                                                              <  1:2432600/150‑1 (MQ=255)
    tACGATACGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAggg                                                                                                                                              >  1:2467547/1‑150 (MQ=255)
      cGATACGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGca                                                                                                                                            <  1:221787/150‑1 (MQ=255)
       gATACGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGcac                                                                                                                                           >  2:1401621/1‑150 (MQ=255)
        aTACGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGcaca                                                                                                                                          >  2:4032651/1‑150 (MQ=255)
         tACGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACAt                                                                                                                                         >  2:1472085/1‑150 (MQ=255)
          aCGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATc                                                                                                                                        >  2:3629504/1‑150 (MQ=255)
             cTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCaaa                                                                                                                                     >  2:4815243/1‑150 (MQ=255)
              tggtggAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAg                                                                                                                                    >  1:101774/1‑150 (MQ=255)
                  ggAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGAt                                                                                                                                >  1:2648509/1‑150 (MQ=255)
                     aaaaCTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCAc                                                                                                                             >  1:732898/1‑150 (MQ=255)
                      aaaCTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGTCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACt                                                                                                                            >  1:3943746/1‑150 (MQ=255)
                      aaaCTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACt                                                                                                                            >  2:1278744/1‑150 (MQ=255)
                       aaCTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTg                                                                                                                           >  1:3969964/1‑150 (MQ=255)
                          tGGGTGCTGCGGTGGCACCGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCg                                                                                                                        >  1:2906517/1‑150 (MQ=255)
                          tGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCg                                                                                                                        >  2:1451549/1‑150 (MQ=255)
                               gctgcGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGaa                                                                                                                   <  1:2691026/150‑1 (MQ=255)
                                 tgcGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGaaaa                                                                                                                 >  2:4190696/1‑150 (MQ=255)
                                 tgcGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGaaaa                                                                                                                 >  2:4546506/1‑150 (MQ=255)
                                  gcGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCAACAGGGCATATCAAAGTGATCACTGGCGGTGAAAAg                                                                                                                <  1:652318/150‑1 (MQ=255)
                                   cGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAgc                                                                                                               <  2:3006394/150‑1 (MQ=255)
                                      tGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCa                                                                                                            <  1:1485887/150‑1 (MQ=255)
                                       ggCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACCTCCAAGTGATCACTGGCGGTGAAAAGCGCaa                                                                                                           >  2:67608/1‑150 (MQ=255)
                                       ggCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCaa                                                                                                           >  1:417306/1‑150 (MQ=255)
                                         cAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAgg                                                                                                         <  2:4120592/150‑1 (MQ=255)
                                             gTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTaa                                                                                                     >  1:3926496/1‑150 (MQ=255)
                                             gTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTaa                                                                                                     >  2:393580/1‑150 (MQ=255)
                                                aaaaTCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATgg                                                                                                  <  2:4950930/150‑1 (MQ=255)
                                                   aTCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTa                                                                                               <  2:4695757/150‑1 (MQ=255)
                                                         tGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTa                                                                                         <  2:783089/150‑1 (MQ=255)
                                                         tGCGCAAGATTACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGAAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCATAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTa                                                                                         <  2:3860126/150‑1 (MQ=255)
                                                          gcgcCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTAt                                                                                        >  1:3745400/1‑150 (MQ=255)
                                                            gcCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATgc                                                                                      <  2:107795/150‑1 (MQ=255)
                                                               aGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGcc                                                                                   >  1:2295981/1‑150 (MQ=255)
                                                                 aTGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGa                                                                                 >  1:4972164/1‑150 (MQ=255)
                                                                     tgAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACgct                                                                             >  2:2450041/2‑150 (MQ=255)
                                                                     cGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCNCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACgct                                                                             >  1:2673938/1‑150 (MQ=255)
                                                                      gAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACgctg                                                                            >  1:2023013/1‑150 (MQ=255)
                                                                          cTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTgctgg                                                                        <  1:4333644/150‑1 (MQ=255)
                                                                           tACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTgctggc                                                                       >  2:2381432/1‑150 (MQ=255)
                                                                              ggAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGctggctgg                                                                    >  1:1896187/1‑150 (MQ=255)
                                                                               gAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACTCTGCtggctggc                                                                   <  1:1586964/150‑1 (MQ=255)
                                                                                aGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGgcg                                                                  <  1:698767/150‑1 (MQ=255)
                                                                                  cTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCa                                                                <  2:1866824/150‑1 (MQ=255)
                                                                                     ggACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTa                                                             <  1:3730184/150‑1 (MQ=255)
                                                                                     ggACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTa                                                             <  1:2561749/150‑1 (MQ=255)
                                                                                      gACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTAc                                                            >  1:4613249/1‑150 (MQ=255)
                                                                                       aCCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACa                                                           >  1:2088768/1‑150 (MQ=255)
                                                                                       aCCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACa                                                           <  2:550472/150‑1 (MQ=255)
                                                                                        ccTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAg                                                          <  1:239272/150‑1 (MQ=255)
                                                                                         cTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAgg                                                         <  2:2276599/150‑1 (MQ=255)
                                                                                          tttAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGa                                                        >  1:2892664/1‑150 (MQ=255)
                                                                                            tAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACg                                                      >  1:1159464/1‑150 (MQ=255)
                                                                                             aaGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCCAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGa                                                     >  2:1506874/1‑150 (MQ=255)
                                                                                              aGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGAt                                                    <  1:1759274/150‑1 (MQ=255)
                                                                                                 tCGCTGGCGCATCTCGAACTCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGTCGGTGAAAAGCGCAAGGGTACTGGCTATTACAATGCGCCGACGCTGCTGGCTGGCGCATTCCCGGACGATGcc                                                 >  1:4341172/1‑150 (MQ=255)
                                                                                                  cGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCa                                                <  2:868629/150‑1 (MQ=255)
                                                                                                    cTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATc                                              >  1:3064398/1‑150 (MQ=255)
                                                                                                      ggCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGt                                            <  2:2821492/150‑1 (MQ=255)
                                                                                                       gcgcAGCCCGACCGCGTCGCCAAGGCAGCAGAAGAGGCGCAAGCTACAGGGCACATCCACGGTATCACGTGCGGTGAAAAGCGCAAGGGTAAGGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTg                                           <  2:4838499/150‑1 (MQ=255)
                                                                                                         gcATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTCATGGCTAGTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCAGCGTGcc                                         >  2:4756083/1‑149 (MQ=255)
                                                                                                         gcATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCa                                         >  2:616990/1‑150 (MQ=255)
                                                                                                         gcATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCa                                         >  2:5060656/1‑150 (MQ=255)
                                                                                                              tcGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAag                                    >  2:1684960/1‑150 (MQ=255)
                                                                                                               cGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAaga                                   >  1:1910853/1‑150 (MQ=255)
                                                                                                                gAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAagag                                  >  1:83964/1‑150 (MQ=255)
                                                                                                                gAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAagag                                  <  2:3745400/150‑1 (MQ=255)
                                                                                                                 aaCGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAgg                                 <  1:1664504/150‑1 (MQ=255)
                                                                                                                  aCGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGt                                <  1:99669/150‑1 (MQ=255)
                                                                                                                    gcgTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTAt                              <  2:5023710/150‑1 (MQ=255)
                                                                                                                     cgTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTAtt                             <  2:3501608/150‑1 (MQ=255)
                                                                                                                     cgTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTAtt                             <  2:3146608/150‑1 (MQ=255)
                                                                                                                      gTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTAttt                            >  2:2337650/1‑150 (MQ=255)
                                                                                                                      gTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTAttt                            <  2:2089902/150‑1 (MQ=255)
                                                                                                                      gTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTAttt                            >  1:4411438/1‑150 (MQ=255)
                                                                                                                        cGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTgg                          >  1:4724491/1‑150 (MQ=255)
                                                                                                                         ggCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGt                         >  2:3189285/1‑150 (MQ=255)
                                                                                                                          gCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCTCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTc                        <  1:1132066/150‑1 (MQ=255)
                                                                                                                          gCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTc                        <  2:3842103/150‑1 (MQ=255)
                                                                                                                                cAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCagtag                  >  2:4507209/1‑150 (MQ=255)
                                                                                                                                  gTAGAAGAGGCTAATGCTACTGGGCACCTCAAAGTGATCACTGGCGGTGAATAGCGCAAGTGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACATGCCGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTagtg                <  2:2606229/150‑1 (MQ=255)
                                                                                                                                  gTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTagtg                >  2:4639836/1‑150 (MQ=255)
                                                                                                                                    agaagaGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTagtgag              <  1:89510/150‑1 (MQ=255)
                                                                                                                                     gaagaGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTagtgagt             >  1:3373284/1‑150 (MQ=255)
                                                                                                                                      aagaGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTagtgagtg            >  1:3647507/1‑150 (MQ=255)
                                                                                                                                       agagGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAgtgt           <  1:4835535/150‑1 (MQ=255)
                                                                                                                                        gagGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGtt          <  2:2648509/150‑1 (MQ=255)
                                                                                                                                        gagGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGtt          <  1:1401621/150‑1 (MQ=255)
                                                                                                                                         agGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTa         >  1:302722/1‑150 (MQ=255)
                                                                                                                                          ggCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTAc        >  1:2354460/1‑150 (MQ=255)
                                                                                                                                             gAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGcc     >  1:5123686/1‑150 (MQ=255)
                                                                                                                                              aaaGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGccc    >  1:4670764/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CATTTACGATACGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGCCCTT  >  CP009273/1510707‑1511000

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: