Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,510,854 G→T 71.5% 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,8540GT71.5% 398.4 / 314.3 340A340A (GCG→GCTpatDgamma‑aminobutyraldehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (48/49);  new base T (107/136);  total (155/185)
Fisher's exact test for biased strand distribution p-value = 3.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

Only 100 of 342 total aligned reads displayed.

GGCATTTACGATACGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGCCCTTCGA  >  CP009273/1510705‑1511003
                                                                                                                                                     |                                                                                                                                                     
  cATTTACGATACGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTAc                                                                                                                                                     <  1:2780954/150‑1 (MQ=255)
      tACGATACGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAggg                                                                                                                                                 >  1:1492114/1‑150 (MQ=255)
          aTACGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGcaca                                                                                                                                             >  1:3574454/1‑150 (MQ=255)
           tACGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACAt                                                                                                                                            <  1:5208916/150‑1 (MQ=255)
            aCGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATc                                                                                                                                           >  2:1524048/1‑150 (MQ=255)
             cGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCa                                                                                                                                          >  2:1287932/1‑150 (MQ=255)
                tggtggAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAg                                                                                                                                       >  1:453457/1‑150 (MQ=255)
                 ggtggAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGt                                                                                                                                      <  1:895828/150‑1 (MQ=255)
                  gtggAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTg                                                                                                                                     >  2:4024582/1‑150 (MQ=255)
                  gtggAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTg                                                                                                                                     <  1:1798411/150‑1 (MQ=255)
                     gAATAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGCGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATc                                                                                                                                  <  2:5322747/150‑1 (MQ=255)
                      aaaaaCTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCa                                                                                                                                 >  1:501665/1‑150 (MQ=255)
                      aaaaaCTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCa                                                                                                                                 >  2:5128483/1‑150 (MQ=255)
                        aaaCTGGGTGCTGAGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACTAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACt                                                                                                                               >  1:3674980/1‑150 (MQ=255)
                          aCTGGGTGCTGCGGTGGCAACGTTAAATTCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTgg                                                                                                                             <  1:428213/150‑1 (MQ=255)
                           cTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGc                                                                                                                            <  1:3100017/150‑1 (MQ=255)
                            tGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCg                                                                                                                           <  2:1015969/150‑1 (MQ=255)
                                tgctgcGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGa                                                                                                                       <  1:218204/150‑1 (MQ=255)
                                     cGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAgc                                                                                                                  <  1:2838441/150‑1 (MQ=255)
                                     cGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAgc                                                                                                                  >  1:4110648/1‑150 (MQ=255)
                                      ggTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAgcg                                                                                                                 >  2:4066199/1‑150 (MQ=255)
                                       gTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTTAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGTGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAgcgc                                                                                                                <  1:3019801/150‑1 (MQ=255)
                                        tGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCa                                                                                                               >  1:789539/1‑150 (MQ=255)
                                        tGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCa                                                                                                               <  2:2256293/150‑1 (MQ=255)
                                        tGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCa                                                                                                               >  2:2781736/1‑150 (MQ=255)
                                          gCAACTTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAg                                                                                                             <  1:2042864/150‑1 (MQ=255)
                                            atcGTTAAAATCTGGTGCGCCAGTTGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTTGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAggg                                                                                                           <  1:2487777/148‑1 (MQ=255)
                                            aaCGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCACGGCAGTAGACGAGGCGAAAGCGACCGGGCACATCCAAGTGATCACTGGCGGTGAAAAGCGCACggg                                                                                                           >  2:4092315/1‑150 (MQ=255)
                                            aaCGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAggg                                                                                                           >  2:2003590/1‑150 (MQ=255)
                                            aaCGTAAAAATCTGGAGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAATCTCGATGGCTCATCTCGAACGCGTCGGCAATGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCATggg                                                                                                           >  2:4348596/1‑150 (MQ=255)
                                               gTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGACAAGCGCAAGGGTaa                                                                                                        >  2:3701997/1‑150 (MQ=255)
                                               gTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTaa                                                                                                        >  2:2235926/1‑150 (MQ=255)
                                               gTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTaa                                                                                                        >  2:191237/1‑150 (MQ=255)
                                                ttAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAAt                                                                                                       <  2:1813633/150‑1 (MQ=255)
                                                ttAAAATCTGGGGCGCCAGATGACGAGTCTCCGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCATTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAAt                                                                                                       <  1:554781/150‑1 (MQ=255)
                                                   aaaTCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGc                                                                                                    <  2:5040553/150‑1 (MQ=255)
                                                    aaTCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCt                                                                                                   >  2:5195189/1‑150 (MQ=255)
                                                       cTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACGGGCGGTGAAAAGCGCAAGGGTAATGGCTAtt                                                                                                <  1:321446/150‑1 (MQ=255)
                                                        tGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTa                                                                                               <  2:743677/150‑1 (MQ=255)
                                                          gTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACt                                                                                             >  1:3455329/1‑150 (MQ=255)
                                                            gcgcCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTAt                                                                                           <  2:3942075/150‑1 (MQ=255)
                                                             cgcCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATg                                                                                          <  1:2241392/150‑1 (MQ=255)
                                                             cgcCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATg                                                                                          <  1:1693937/150‑1 (MQ=255)
                                                             cgcCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATg                                                                                          <  2:4884527/150‑1 (MQ=255)
                                                               ccAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATgcg                                                                                        <  1:3847990/150‑1 (MQ=255)
                                                                cAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATgcgc                                                                                       >  1:2317257/1‑150 (MQ=255)
                                                                cAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATgcgc                                                                                       >  1:645879/1‑150 (MQ=255)
                                                                  gATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCg                                                                                     <  1:1215373/150‑1 (MQ=255)
                                                                        gAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACgctg                                                                               >  1:153864/1‑150 (MQ=255)
                                                                        gAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACgctg                                                                               <  2:346983/150‑1 (MQ=255)
                                                                            cTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTgctgg                                                                           >  2:1450854/1‑150 (MQ=255)
                                                                             tACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTgctggc                                                                          <  2:3574454/150‑1 (MQ=255)
                                                                               cGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTgctggctg                                                                        <  1:4456718/150‑1 (MQ=255)
                                                                                 gAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACTCTGCtggctggc                                                                      >  2:2639132/1‑150 (MQ=255)
                                                                                       ggACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACTCTGCTGGCTGGCGCATTa                                                                >  1:1054733/1‑150 (MQ=255)
                                                                                        taCCTTTAAGCACGCTGGCGAATCTCGAACGCGTCGGCAAGGATGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACTCTGCTGGCTGGCGCATTAc                                                               <  1:1884123/149‑1 (MQ=255)
                                                                                        gACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTAc                                                               >  2:1667886/1‑150 (MQ=255)
                                                                                             ttAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGAc                                                          <  1:2003590/150‑1 (MQ=255)
                                                                                              tAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACg                                                         >  1:5052327/1‑150 (MQ=255)
                                                                                               aaGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGa                                                        >  2:5166511/1‑150 (MQ=255)
                                                                                                tgCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACTCTGCTGGCTGGCGCATTACAGGACGAt                                                       <  2:4494891/149‑1 (MQ=255)
                                                                                                aGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGAt                                                       <  1:1252528/150‑1 (MQ=255)
                                                                                                     gCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACTCTGCTGGCTGGCGCATTACAGGACGATGCCAt                                                  <  1:3587735/150‑1 (MQ=255)
                                                                                                      cTGGCGCATCTCTAACGCATCGGCATGGCAGTAGATGAGGCTAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATc                                                 <  1:2354005/150‑1 (MQ=255)
                                                                                                       tGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCg                                                <  2:3319033/150‑1 (MQ=255)
                                                                                                        ggCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGt                                               >  2:856442/1‑150 (MQ=255)
                                                                                                           gcATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCa                                            <  2:4655634/150‑1 (MQ=255)
                                                                                                             ctctcGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACTCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCaaa                                          >  2:228380/2‑150 (MQ=255)
                                                                                                                 cGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCACAAaga                                      >  2:1763940/1‑150 (MQ=255)
                                                                                                                 cGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAaga                                      <  2:236370/150‑1 (MQ=255)
                                                                                                                  gAACGAGTCGGAAAGGCAGTATAAGATGCGAAAGCTACAGGGAACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGTCTATTACTATGCGCCGACGCTGCTGGCTGTCGCATTACATGACGATGCCATCGTGCAAAAagag                                     >  1:4815798/1‑150 (MQ=255)
                                                                                                                   aaCGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAgg                                    >  2:2155671/1‑150 (MQ=255)
                                                                                                                    aCGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGTCGATGCCATCGTGCAAAAAGAGGt                                   <  2:2944957/150‑1 (MQ=255)
                                                                                                                    aCGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGt                                   <  2:3521694/150‑1 (MQ=255)
                                                                                                                      gcgTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTAt                                 <  1:3949153/150‑1 (MQ=255)
                                                                                                                       cgTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACTCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTAtt                                <  2:763331/150‑1 (MQ=255)
                                                                                                                        gTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTAttt                               <  1:4591275/150‑1 (MQ=255)
                                                                                                                         tCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTg                              <  2:1098158/150‑1 (MQ=255)
                                                                                                                           ggCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGt                            >  1:4336919/1‑150 (MQ=255)
                                                                                                                           ggCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGt                            >  1:1844929/1‑150 (MQ=255)
                                                                                                                            gCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTc                           <  2:649937/150‑1 (MQ=255)
                                                                                                                             cAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCGATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTcc                          >  2:4234354/1‑150 (MQ=255)
                                                                                                                              aaGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCa                         >  1:2111999/1‑150 (MQ=255)
                                                                                                                              aaGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCa                         >  2:175036/1‑150 (MQ=255)
                                                                                                                                  cAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACTCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCagtag                     <  2:4369912/150‑1 (MQ=255)
                                                                                                                                   aGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCagtagt                    <  1:1618910/150‑1 (MQ=255)
                                                                                                                                   aGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACTCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCagtagt                    <  2:2508485/150‑1 (MQ=255)
                                                                                                                                     tAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTagtga                  <  1:606052/150‑1 (MQ=255)
                                                                                                                                        aagaGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACTCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTagtgagtg               <  2:4273520/150‑1 (MQ=255)
                                                                                                                                        aagaGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACTCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTagtgagtg               >  2:1379973/1‑150 (MQ=255)
                                                                                                                                         atagGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAgtgt              >  1:3259127/3‑150 (MQ=255)
                                                                                                                                             gCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACg          <  1:4782882/150‑1 (MQ=255)
                                                                                                                                             gCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACTCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACg          <  1:5128483/150‑1 (MQ=255)
                                                                                                                                              cGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGc         >  2:479086/1‑150 (MQ=255)
                                                                                                                                               gAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGcc        >  1:3366002/1‑150 (MQ=255)
                                                                                                                                               gAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCAAGTAGTGAGTGTTACGcc        >  2:1291835/1‑150 (MQ=255)
                                                                                                                                                 aaGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGCCCt      <  2:3667659/150‑1 (MQ=255)
                                                                                                                                                    cGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACTCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGCCCTTCg   <  2:940598/150‑1 (MQ=255)
                                                                                                                                                    cGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACTCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGCCCTTCg   <  2:4993187/150‑1 (MQ=255)
                                                                                                                                                     gACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACTCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGCCCTTCGa  >  2:3111155/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGCATTTACGATACGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGCCCTTCGA  >  CP009273/1510705‑1511003

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: