Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,820,443 T→A 9.9% T246T (ACA→ACT astE ← succinylglutamate desuccinylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,820,4430TA 9.9% 548.0 / 42.6 181T246T (ACA→ACTastEsuccinylglutamate desuccinylase
Reads supporting (aligned to +/- strand):  ref base T (73/90);  new base A (10/8);  total (83/98)
Fisher's exact test for biased strand distribution p-value = 4.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 194 total aligned reads displayed.

GTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGTCTGGCGCTG  >  CP009273/1820297‑1820590
                                                                                                                                                  |                                                                                                                                                   
gTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGgtg                                                                                                                                                    <  2:1888281/148‑1 (MQ=255)
 ttACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGgtgt                                                                                                                                                   <  2:1442202/148‑1 (MQ=255)
 ttACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGgtgt                                                                                                                                                   >  1:912343/1‑148 (MQ=255)
  tACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGtt                                                                                                                                                  <  2:1681386/148‑1 (MQ=255)
     ggTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTtctc                                                                                                                                               <  2:2705231/148‑1 (MQ=255)
     ggTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTtctc                                                                                                                                               <  1:1937787/148‑1 (MQ=255)
      gTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCa                                                                                                                                              <  1:2083092/148‑1 (MQ=255)
      gTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCa                                                                                                                                              >  1:1739724/1‑148 (MQ=255)
        aaaaCGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACg                                                                                                                                            <  1:362935/148‑1 (MQ=255)
         aaaCGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGa                                                                                                                                           >  2:945884/1‑148 (MQ=255)
         aaaCGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGa                                                                                                                                           <  1:1010504/148‑1 (MQ=255)
         aaaCGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGa                                                                                                                                           <  1:1638375/148‑1 (MQ=255)
          aaCGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGAt                                                                                                                                          >  1:2597254/1‑148 (MQ=255)
          aaCGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGAt                                                                                                                                          >  1:539923/1‑148 (MQ=255)
              ttCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCCGCGTGTCACTTGCCAGATGCAGTTCGAAGGACGGCGAGTGGCGAGTACTTTGCGACACCACCCGCGAACGGAGCGGCGGTGTGCTCACGATCCCg                                                                                                                                      >  2:1703303/1‑148 (MQ=255)
              ttCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCg                                                                                                                                      >  2:1496720/1‑148 (MQ=255)
              ttCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCg                                                                                                                                      >  1:2721549/1‑148 (MQ=255)
               tCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGa                                                                                                                                     >  2:2934819/1‑148 (MQ=255)
               tCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGa                                                                                                                                     >  1:1206703/1‑148 (MQ=255)
                 ctctCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCTAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGaca                                                                                                                                   <  1:423626/148‑1 (MQ=255)
                 ctctCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGaca                                                                                                                                   <  2:1379960/148‑1 (MQ=255)
                    tcCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACActc                                                                                                                                <  1:3204854/148‑1 (MQ=255)
                           tGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAg                                                                                                                         <  1:2429240/148‑1 (MQ=255)
                             cgcCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGAt                                                                                                                       >  1:2729879/1‑148 (MQ=255)
                             cgcCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGAt                                                                                                                       <  1:2490548/148‑1 (MQ=255)
                               ccAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATag                                                                                                                     <  1:1040429/148‑1 (MQ=255)
                                 aGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATagca                                                                                                                   >  1:2957133/1‑148 (MQ=255)
                                   cAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATagcagc                                                                                                                 >  2:1695522/1‑148 (MQ=255)
                                         tCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCgcagca                                                                                                           <  2:2866712/148‑1 (MQ=255)
                                          ccTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCaa                                                                                                          <  1:838145/148‑1 (MQ=255)
                                          ccTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCaa                                                                                                          <  1:196585/148‑1 (MQ=255)
                                          ccTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGAGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCaa                                                                                                          >  1:280600/1‑148 (MQ=255)
                                          ccTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGAGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCaa                                                                                                          >  1:280783/1‑148 (MQ=255)
                                             ttCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGAGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTg                                                                                                       <  2:809617/148‑1 (MQ=255)
                                               ctcAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTgcg                                                                                                     >  2:1400731/1‑148 (MQ=255)
                                               ctcAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTgcg                                                                                                     >  2:820157/1‑148 (MQ=255)
                                               ctcAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTgcg                                                                                                     >  2:3169729/1‑148 (MQ=255)
                                                              aTTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCa                                                                                      >  2:825654/1‑148 (MQ=255)
                                                               ttCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCaa                                                                                     <  2:3252342/148‑1 (MQ=255)
                                                               ttCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCaa                                                                                     >  2:1597105/1‑148 (MQ=255)
                                                                 cAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAAc                                                                                   <  2:2721635/148‑1 (MQ=255)
                                                                  aGCGTGTCACTTGCCATATGCATTTCGAGGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACt                                                                                  >  1:2181039/1‑148 (MQ=255)
                                                                  aGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACt                                                                                  >  1:519076/1‑148 (MQ=255)
                                                                  aGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACt                                                                                  <  1:1395200/148‑1 (MQ=255)
                                                                       gtCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGa                                                                             <  2:2774585/148‑1 (MQ=255)
                                                                        tCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGaa                                                                            >  2:2149592/1‑148 (MQ=255)
                                                                          acgtGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGa                                                                          <  1:442524/145‑1 (MQ=255)
                                                                          aCTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGa                                                                          <  1:638705/148‑1 (MQ=255)
                                                                            ttGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGAGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATc                                                                        <  1:1086936/148‑1 (MQ=255)
                                                                             tGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCg                                                                       >  2:630296/1‑148 (MQ=255)
                                                                             tGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCg                                                                       <  2:2229577/148‑1 (MQ=255)
                                                                             tGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCg                                                                       <  2:2624717/148‑1 (MQ=255)
                                                                              gCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGAGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGt                                                                      <  2:2534144/148‑1 (MQ=255)
                                                                                cATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGttt                                                                    <  1:2947539/148‑1 (MQ=255)
                                                                                 atatGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGtttt                                                                   >  1:2422999/1‑148 (MQ=255)
                                                                                   atGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGc                                                                 >  1:809437/1‑148 (MQ=255)
                                                                                      cATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGAGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCa                                                              >  2:1138127/1‑148 (MQ=255)
                                                                                        tttCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTTCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCaaa                                                            >  1:2971771/1‑148 (MQ=255)
                                                                                        tttCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCaaa                                                            <  1:717775/148‑1 (MQ=255)
                                                                                         ttCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAAATGGCGAAGATCGTTTTGCCCAAAg                                                           <  1:1258925/148‑1 (MQ=255)
                                                                                         ttCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAAATGGCGAAGATCGTTTTGCCCAAAg                                                           <  1:2907307/148‑1 (MQ=255)
                                                                                             aaGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCa                                                       <  2:3025204/148‑1 (MQ=255)
                                                                                              aGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCaa                                                      <  1:1292558/148‑1 (MQ=255)
                                                                                              aGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGAGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCaa                                                      <  2:1877529/148‑1 (MQ=255)
                                                                                                gACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACg                                                    >  1:1925867/1‑148 (MQ=255)
                                                                                                 aCGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGc                                                   <  2:2597254/148‑1 (MQ=255)
                                                                                                 aCGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGc                                                   <  2:539923/148‑1 (MQ=255)
                                                                                                 aCGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGc                                                   <  1:945884/148‑1 (MQ=255)
                                                                                                  cGGCGAGTGGCGAGTAATGTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCt                                                  <  2:2425961/148‑1 (MQ=255)
                                                                                                   ggcgagtggcgagtAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCtt                                                 <  1:3082468/148‑1 (MQ=255)
                                                                                                   ggcgagtggcgagtAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCtt                                                 <  2:694329/148‑1 (MQ=255)
                                                                                                    gcgagtggcgagtAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCANATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCttt                                                <  1:1680471/148‑1 (MQ=255)
                                                                                                     cgagtggcgagtAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTg                                               <  1:70671/148‑1 (MQ=255)
                                                                                                         tggcgagtAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCaa                                           >  1:2612138/1‑148 (MQ=255)
                                                                                                          ggcgagtAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCAAg                                          >  1:3277181/1‑148 (MQ=255)
                                                                                                          ggcgagtAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCAAg                                          <  2:1720206/148‑1 (MQ=255)
                                                                                                          ggcgagtAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCAAg                                          <  1:2539144/148‑1 (MQ=255)
                                                                                                                  aTTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCAAGTTCCAGCg                                  >  2:1912269/1‑148 (MQ=255)
                                                                                                                    ttGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCAAGTTCCAGCGTa                                >  1:1671056/1‑148 (MQ=255)
                                                                                                                        gAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCAAGTTCCAGCGTACAgg                            <  2:1206703/148‑1 (MQ=255)
                                                                                                                         aaaCCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCAAGTTCCAGCGTACAGGc                           <  2:2827384/148‑1 (MQ=255)
                                                                                                                         aaaCCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCAAGTTCCAGCGTACAGGc                           <  1:2344500/148‑1 (MQ=255)
                                                                                                                          aaCCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCAAGTTCCAGCGTACAGGcc                          >  2:2270114/1‑148 (MQ=255)
                                                                                                                                    tAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCaa                >  1:3285136/1‑148 (MQ=255)
                                                                                                                                    tAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCaa                >  1:3284997/1‑148 (MQ=255)
                                                                                                                                     aaCGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCaaa               >  1:2623727/1‑148 (MQ=255)
                                                                                                                                     aaCGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCaaa               <  2:2234768/148‑1 (MQ=255)
                                                                                                                                      aCGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCaaaa              >  2:2482534/1‑148 (MQ=255)
                                                                                                                                       cGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAAt             >  1:1989962/1‑148 (MQ=255)
                                                                                                                                         gAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGt           <  1:3169729/148‑1 (MQ=255)
                                                                                                                                          aGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGTc          >  2:2061111/1‑148 (MQ=255)
                                                                                                                                          aGCGGCGGAGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGTc          <  2:2415110/148‑1 (MQ=255)
                                                                                                                                           gcggcgGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGTCt         <  1:2131514/148‑1 (MQ=255)
                                                                                                                                            cggcggTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGTCTg        <  2:2366072/148‑1 (MQ=255)
                                                                                                                                             ggcggTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGTCTgg       <  2:465645/148‑1 (MQ=255)
                                                                                                                                                ggTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGTCTGgcgc    <  2:1308478/148‑1 (MQ=255)
                                                                                                                                                 gtgtTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGTCTGGCGCt   >  2:996135/1‑148 (MQ=255)
                                                                                                                                                  tgtTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGTCTGGCGCTg  >  1:3313153/1‑148 (MQ=255)
                                                                                                                                                  tgtTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGTCTGGCGCTg  >  1:1513457/1‑148 (MQ=255)
                                                                                                                                                  tgtTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGTCTGGCGCTg  >  2:2703019/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGTCTGGCGCTG  >  CP009273/1820297‑1820590

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: