Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,510,854 G→T 64.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,8540GT64.4% 139.4 / 291.7 236A340A (GCG→GCTpatDgamma‑aminobutyraldehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (37/47);  new base T (74/78);  total (111/125)
Fisher's exact test for biased strand distribution p-value = 5.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

Only 100 of 237 total aligned reads displayed.

CATTTACGATACGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGCC  >  CP009273/1510707‑1510997
                                                                                                                                                   |                                                                                                                                               
  tttACGATACGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGCACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGAc                                                                                                                                               >  1:1567375/1‑148 (MQ=255)
   ttACGATACGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACa                                                                                                                                              <  1:2001671/148‑1 (MQ=255)
   ttACGATACGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACa                                                                                                                                              >  2:684809/1‑148 (MQ=255)
    tACGATACGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAg                                                                                                                                             <  2:309791/148‑1 (MQ=255)
      cGATACGCTGGTGGAATAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAggg                                                                                                                                           <  1:349364/148‑1 (MQ=255)
         tACGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGcac                                                                                                                                        >  1:469704/1‑148 (MQ=255)
         tACGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGcac                                                                                                                                        >  1:1673858/1‑148 (MQ=255)
          aCGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGcaca                                                                                                                                       >  2:2155901/1‑148 (MQ=255)
                 tggAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGt                                                                                                                                >  1:985425/1‑148 (MQ=255)
                 tggAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGt                                                                                                                                >  1:2972397/1‑148 (MQ=255)
                   gAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGa                                                                                                                              <  1:3206867/148‑1 (MQ=255)
                     aaaaCTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATc                                                                                                                            <  2:1640127/148‑1 (MQ=255)
                      aaaCTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCa                                                                                                                           >  1:2472822/1‑148 (MQ=255)
                      aaaCTGGGTGCTGCGGTGGAAACGTTAAAATCTGGGGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCa                                                                                                                           <  1:841622/148‑1 (MQ=255)
                       aaCTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCAc                                                                                                                          >  2:1759023/1‑148 (MQ=255)
                        aCTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTATAGGGCACATCAAAGTGATCACt                                                                                                                         >  2:3565399/1‑148 (MQ=255)
                        aCTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACt                                                                                                                         >  2:527574/1‑148 (MQ=255)
                         cTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGTTGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTg                                                                                                                        <  1:564156/148‑1 (MQ=255)
                             gTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCgg                                                                                                                    <  1:2017644/148‑1 (MQ=255)
                              tgctgcGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGt                                                                                                                   <  2:719438/148‑1 (MQ=255)
                              tgctgcGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGt                                                                                                                   <  2:137504/148‑1 (MQ=255)
                                   cGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGaaaa                                                                                                              <  1:1510799/148‑1 (MQ=255)
                                   cGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGaaaa                                                                                                              <  1:1510784/148‑1 (MQ=255)
                                       ggCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAgcgc                                                                                                          <  2:2510122/148‑1 (MQ=255)
                                        gCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCa                                                                                                         <  1:25810/148‑1 (MQ=255)
                                        gCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCa                                                                                                         <  1:2865799/148‑1 (MQ=255)
                                            cGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAggg                                                                                                     <  2:3231397/148‑1 (MQ=255)
                                              ttAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTa                                                                                                   <  2:2326371/148‑1 (MQ=255)
                                              ttAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTa                                                                                                   >  1:2477925/1‑148 (MQ=255)
                                              ttAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTa                                                                                                   <  1:877045/148‑1 (MQ=255)
                                              ttAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTa                                                                                                   >  1:968475/1‑148 (MQ=255)
                                               tAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTaa                                                                                                  >  2:633506/1‑148 (MQ=255)
                                                  aaTCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATgg                                                                                               <  2:861381/148‑1 (MQ=255)
                                                   aTCGGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGc                                                                                              <  1:394689/148‑1 (MQ=255)
                                                     cTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTa                                                                                            <  2:755092/148‑1 (MQ=255)
                                                       ggTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTAtt                                                                                          <  2:2856125/148‑1 (MQ=255)
                                                        gTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTa                                                                                         >  2:1137542/1‑148 (MQ=255)
                                                             ccAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATg                                                                                    >  1:1861365/1‑148 (MQ=255)
                                                                gATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATgcgc                                                                                 <  2:442769/148‑1 (MQ=255)
                                                                   gACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGa                                                                              >  2:3398590/1‑148 (MQ=255)
                                                                    aCGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGAc                                                                             >  1:1037401/1‑148 (MQ=255)
                                                                    aCGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGAc                                                                             >  1:2135795/1‑148 (MQ=255)
                                                                        ttttACGGAGATTGGACCTTTAAGCGCGGTGGCGCTTCTCGAACGCGTCGGCAAGGCAGTTGAAGAGGGGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACgctg                                                                         <  2:1077336/145‑1 (MQ=255)
                                                                        gTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACgctg                                                                         <  2:2782022/148‑1 (MQ=255)
                                                                        gTCTACGGAGCTTGGACATTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACgctg                                                                         >  2:242219/1‑148 (MQ=255)
                                                                        gTCTACGGAGCTTGGACATTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACgctg                                                                         >  2:1071647/1‑148 (MQ=255)
                                                                         tCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACgctgc                                                                        >  1:702779/1‑148 (MQ=255)
                                                                          cTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACgctgct                                                                       <  1:820246/148‑1 (MQ=255)
                                                                          cTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACgctgct                                                                       <  2:397171/148‑1 (MQ=255)
                                                                             cGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTgctggc                                                                    <  2:2477925/148‑1 (MQ=255)
                                                                                  cTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGgcg                                                               >  1:2221344/1‑148 (MQ=255)
                                                                                     ggACCTTTACGCGCGCTGGCTCATCTCGAACGCGCCGGCAAGGCAGTAGAAGAGGCTAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCAt                                                            <  2:1205106/148‑1 (MQ=255)
                                                                                       aCCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTa                                                          >  1:1153942/1‑148 (MQ=255)
                                                                                        ccTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTAc                                                         <  2:2713888/148‑1 (MQ=255)
                                                                                        ccTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTAc                                                         <  2:1645060/148‑1 (MQ=255)
                                                                                        catttAAGCACGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGTGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTAc                                                         <  1:712983/146‑1 (MQ=255)
                                                                                           ttAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAgg                                                      >  2:3349373/1‑148 (MQ=255)
                                                                                             aaGATCGCTGGCGCATCTCGAACGCGTGGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGGGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGAc                                                    <  2:1272862/148‑1 (MQ=255)
                                                                                               gCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGa                                                  >  2:3058132/1‑148 (MQ=255)
                                                                                               gCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGa                                                  >  1:2652233/1‑148 (MQ=255)
                                                                                                cTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGAt                                                 <  1:230512/148‑1 (MQ=255)
                                                                                                 tCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATg                                                <  2:1733625/148‑1 (MQ=255)
                                                                                                   gCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGcc                                              >  1:3087516/1‑148 (MQ=255)
                                                                                                    cTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCa                                             <  2:179549/148‑1 (MQ=255)
                                                                                                       gcgcATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCg                                          <  2:3357522/148‑1 (MQ=255)
                                                                                                       gcgcATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCg                                          >  2:962991/1‑148 (MQ=255)
                                                                                                       gagcATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAATCTACATGGCACATCTAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCg                                          <  2:3104224/146‑1 (MQ=255)
                                                                                                         gcATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTg                                        <  2:540395/148‑1 (MQ=255)
                                                                                                            tctcGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCaa                                     <  2:3512870/148‑1 (MQ=255)
                                                                                                             ctcGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCaaa                                    >  1:1201420/1‑148 (MQ=255)
                                                                                                              tcGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCaaaa                                   <  2:2908970/148‑1 (MQ=255)
                                                                                                              tcGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCaaaa                                   <  2:311547/148‑1 (MQ=255)
                                                                                                               cGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCaaaaa                                  >  1:1108320/1‑148 (MQ=255)
                                                                                                                gAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAag                                 >  2:2181642/1‑148 (MQ=255)
                                                                                                                gAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAag                                 <  1:1722955/148‑1 (MQ=255)
                                                                                                                  aCGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAagag                               <  1:2824386/148‑1 (MQ=255)
                                                                                                                  aCGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAagag                               <  1:636714/148‑1 (MQ=255)
                                                                                                                    gcgTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGt                             >  1:2802661/1‑148 (MQ=255)
                                                                                                                    gcgTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGt                             <  2:1567375/148‑1 (MQ=255)
                                                                                                                     cgTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTa                            <  1:458312/148‑1 (MQ=255)
                                                                                                                      gTCGGCACGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTAt                           <  1:1039666/148‑1 (MQ=255)
                                                                                                                       tCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTAtt                          >  2:3360042/1‑148 (MQ=255)
                                                                                                                       tCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTAtt                          <  2:1861365/148‑1 (MQ=255)
                                                                                                                           cAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGt                      <  1:3343201/148‑1 (MQ=255)
                                                                                                                             aGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTcc                    >  2:2360850/1‑148 (MQ=255)
                                                                                                                                cAGTAGAAGAGTCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCagt                 <  1:2160540/148‑1 (MQ=255)
                                                                                                                                cAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCagt                 >  2:1047620/1‑148 (MQ=255)
                                                                                                                                cAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCagt                 >  1:2110514/1‑148 (MQ=255)
                                                                                                                                  gTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCagtag               >  1:3367438/1‑148 (MQ=255)
                                                                                                                                  gTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCagtag               >  1:1720568/1‑148 (MQ=255)
                                                                                                                                    agaagaGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTagtg             <  2:1947042/148‑1 (MQ=255)
                                                                                                                                    agaagaGGCGAAAGCTACAGGGCACATAAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTagtg             <  2:2668120/148‑1 (MQ=255)
                                                                                                                                     gaagaGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTagtga            >  1:2422683/1‑148 (MQ=255)
                                                                                                                                       agagGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTagtgagt          >  1:137718/1‑148 (MQ=255)
                                                                                                                                         agGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAgtgt        <  2:1444425/148‑1 (MQ=255)
                                                                                                                                          ggCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGtt       >  2:2241159/1‑148 (MQ=255)
                                                                                                                                           gCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTa      <  1:2938016/148‑1 (MQ=255)
                                                                                                                                              aaaGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGc   >  1:1932743/1‑148 (MQ=255)
                                                                                                                                              aaaGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGc   >  1:2328642/1‑148 (MQ=255)
                                                                                                                                               aaGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGcc  >  2:579941/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
CATTTACGATACGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGCC  >  CP009273/1510707‑1510997

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: