Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,484,072 G→A 49.7% pseudogene (150/252 nt) gapC ← pseudogene, GAP dehydrogenase; glyceraldehyde‑3‑phosphate dehydrogenase (second fragment)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,484,0720GA49.7% 15.4 / 243.1 145pseudogene (150/252 nt)gapCpseudogene, GAP dehydrogenase; glyceraldehyde‑3‑phosphate dehydrogenase (second fragment)
Reads supporting (aligned to +/- strand):  ref base G (30/43);  new base A (37/35);  total (67/78)
Fisher's exact test for biased strand distribution p-value = 2.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

Only 100 of 149 total aligned reads displayed.

TTAATGAGTGAAAGAGGCGGAGGTTTTTTCCTCCGCCTGTGCGCGTCAGAGTTTAGCGAATTTTTCGAGGGTGCGAATAAGCTGTGTGACGAAGCCATATTCGTTATCGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACGGCGGTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTC  >  CP009273/1483925‑1484219
                                                                                                                                                   |                                                                                                                                                   
ttAATGAGTGAAAGAGGCGGAGGTTTTTTCCTCCGCCTGTGCGCGTCAGAGTTTAGCGAATTTTTCGAGGGTGCGAATAAGCTGTGTGACGAAGCCATATTCGTTATCGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCAcggc                                                                                                                                                   >  2:348920/1‑150 (MQ=255)
 tAATGAGTGAAAGAGGCGGAGGTTTTTTCCTCCGCCTGTGCGCGTCAGAGTTTAGCGAATTTTTCGAGGGTGCGAATAAGCTGTGTGACGAAGCCATATTCGTTATCGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCAcggcg                                                                                                                                                  >  1:44897/1‑150 (MQ=255)
  aaaagagcgaccgaagtcgctctttttagatcacagtgtcaCGTCAGAGTTTAGCGAATTTTTCGAGGGTGCGAATAAGCTGCGTGACGAAGCCATATTCGTTATCGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCAcagcgg                                                                                                                                                 >  1:979254‑M1/42‑150 (MQ=255)
  aaTGAGTGAAAGAGGCGGAGGTTTTTTCCTCCGCCTGTGCGCGTCAGAGTTTAGCGAATTTTTCGAGGGTGCGAATAAGCTGTGTGACGAAGCCATATTCGTTATAGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCAcggcgg                                                                                                                                                 >  1:705627/1‑150 (MQ=255)
    tGAGTGAAAGAGGCGGAGGTTTTTTCCTCCGCCTGTGCGCGTCAGAGTTTAGCGAATTTTTCGAGGGTGCGAATAAGCTGTGTGACGAAGCCATATTCGTTATCGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACGGCGGTa                                                                                                                                               <  1:692036/150‑1 (MQ=255)
      gagcgaccgaagtcgctctttttagatcacagtgtcaCGTCAGAGTTTAGCGAATTTTTCGAGGGTGCGAATAAGCTGCGTGACGAAGCCATATTCGTTATCGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACAGCGGTAAt                                                                                                                                             >  1:764726‑M1/38‑150 (MQ=255)
         gAAAGAGGCGGAGGTTTTTTCCTCCGCCTGTGCGCGTCAGAGTTTAGCGAATTTTTCGAGGGTGCGAATAAGCTGTGTGACGAAGCCATATTCGTTATCGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACGGCGGTAAtttt                                                                                                                                          >  1:228095/1‑149 (MQ=255)
         gAAAGAGGCGGAGGTTTTTTCCTCCGCCTGTGCGCGTCAGAGTTTAGCGAATTTTTCGAGGGTGCGAATAAGCTGTGTGACGAAGCCATATTCGTTATCGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACGGCGGTAATTTc                                                                                                                                          >  2:660968/1‑150 (MQ=255)
           accgaagtcgctctttttagatcacagtgtcaCGTCAGAGTTTAGCGAATTTTTCGAGGGTGCGAATAAGCTGCGTGACGAAGCCATATTCGTTATCGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACAGCGGTAATTTCCg                                                                                                                                        >  1:326643‑M1/33‑150 (MQ=255)
             cgaagtcgctctttttagatcacagtgtcaCGTCAGAGTTTAGCGAATTTTTCGAGGGTGCGAATAAGCTGCGTGACGAAGCCATATTCGTTATCGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACAGCGGTAATTTCCGtt                                                                                                                                      <  1:338269‑M1/120‑1 (MQ=255)
             cgaagtcgctctttttagatcacagtgtcaCGTCAGAGTTTAGCGAATTTTTCGAGGGTGCGAATAAGCTGCGTGACGAAGCCATATTCGTTATCGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACAGCGGTAATTTCCGtt                                                                                                                                      <  1:85204‑M1/120‑1 (MQ=255)
              gaagtcgctctttttagatcacggtgtcaCGTCAGAGTTTAGCGAATTTTTCGAGGGTGCGAATAAGCTGCGTGACGAAGCCATATTCGTTATCGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACAGCGGTAATTTCCGttt                                                                                                                                     >  2:289264‑M1/30‑150 (MQ=255)
              gaagtcgctctttttagatcacagtgtcaCGTCAGAGTTTAGCGAATTTTTCGAGGGTGCGAATAAGCTGCGTGACGAAGCCATATTCGTTATCGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACAGCGGTAATTTCCGttt                                                                                                                                     <  1:790661‑M1/121‑1 (MQ=255)
              gaagtcgctctttttagatcacagtgtcaCGTCAGAGTTTAGCGAATTTTTCGAGGGTGCGAATAAGCTGCGTGACGAAGCCATATTCGTTATCGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACAGCGGTAATTTCCGttt                                                                                                                                     <  2:939658‑M1/121‑1 (MQ=255)
               aagtcgctctttttagatcacagtgtcaCGTCAGAGTTTAGCGAATTTTTCGAGGGTGCGAATAAGCTGCGTGACGAAGCCATATTCGTTATCGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACAGCGGTAATTTCCGTTTg                                                                                                                                    <  1:15371‑M1/122‑1 (MQ=255)
                 cGGAGGTTTTTTCCTCCGCCTGTGCGCGTCAGAGTTTAGCGAATTTTTCGAGGGTGCGAATAAGCTGTGTGACGAAGCCATATTCGTTATCGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACGGCGGTAATTTCCGTTTGCg                                                                                                                                  <  2:44897/150‑1 (MQ=255)
                  ggAGGTTTTTTCCTCCGCCTGTGCGCGTCAGAGTTTAGCGAATTTTTCGAGGGTGCGAATAAGCTGTGTGACGAAGCCATATTCGTTATCGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACGGCGGTAATTTCCGTTTGCGt                                                                                                                                 <  2:933987/150‑1 (MQ=255)
                    gctctttttagatcacagtgtcaCGTCAGAGTTTAGCGAATTTTTCGAGGGTGCGAATAAGCTGCGTGACGAAGCCATATTCGTTATCGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACAGCGGTAATTTCCGTTTGCGTg                                                                                                                                >  2:790661‑M1/24‑149 (MQ=255)
                     ctctttttagatcacagtgtcaCGTCAGAGTTTAGCGAATTTTTCGAGGGTGCGAATAAGCTGCGTGACGAAGCCATATTCGTTATCGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACAGCGGTAATTTCCGTTTGCGTgg                                                                                                                               <  1:531608‑M1/127‑1 (MQ=255)
                         ttttCCTCCGCCTGTGCGCGTCAGAGTTTAGCGAATTTTTCGAGGGTGCGAATAAGCTGTGTGACGAAGCCATATTCGTTATCGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACGGCGGTAATTTCCGTTTGCGTGGCATc                                                                                                                           <  2:54885/149‑1 (MQ=255)
                          tttagatcacagtgtcaCGTCAGAGTTTAGCGAATTTTTCGAGGGTGCGAATAAGCTGCGTGACGAAGCCATATTCGTTATCGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACAGCGGTAATTTCCGTTTGCGTGGCATCa                                                                                                                          <  1:567373‑M1/132‑1 (MQ=255)
                           ttCCTCCGCCTGTGCGCGTCAGAGTTTAGCGAATTTTTCGAGGGTGCGAATAAGCTGTGTGACGAAGCCATATTCGTTATCGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACGGCGGTAATTTCCGTTTGCGTGGCATCaaa                                                                                                                        >  1:24488/1‑150 (MQ=255)
                            tagatcacagtgtcaCGTCAGAGTTTAGCGAATTTTTCGAGGGTGCGAATAAGCTGCGTGACGAAGCCATATTCGTTATCGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACAGCGGTAATTTc                                                                                                                                          <  2:632503‑M1/116‑1 (MQ=255)
                             cctccGCCTGTGCGCGTCAGAGTTTAGCGAATTTTTCGAGGGTGCGAATAAGCTGTGTGACGAAGCCATATTCGTTATCGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACGGCGGTAATTTCCGTTTGCGTGGCATCAAaca                                                                                                                      >  2:586328/1‑150 (MQ=255)
                             agatcacagtgtcaCGTCAGAGTTTAGCGAATTTTTCGAGGGTGCGAATAAGCTGCGTGACGAAGCCATATTCGTTATCGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACAGCGGTAATTTcc                                                                                                                                         >  1:632503‑M1/15‑131 (MQ=255)
                              gatcacagtgtcaCGTCAGAGTTTAGCGAATTTTTCGAGGGTGCGAATAAGCTGCGTGACGAAGCCATATTCGTTATCGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACAGCGGTAATTTCCGTTTGCGTGGCATCAAacac                                                                                                                     <  2:600116‑M1/137‑1 (MQ=255)
                               tccGCCTGTGCGCGTCAGAGTTTAGCGAATTTTTCGAGGGTGCGAATAAGCTGTGTGACGAAGCCATATTCGTTATCGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACGGCGGTAATTTCCGTTTGCGTGGCATCAAACacc                                                                                                                    >  1:27045/1‑150 (MQ=255)
                                 cacagtgtcaCGTCAGAGTTTAGCGAATTTTTCGAGGGTGCGAATAAGCTGCGTGACGAAGCCATATTCGTTATCGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACAGCGGTAATTTCCGTTTGCGTGGCATCAAACaccg                                                                                                                   <  1:329114‑M1/139‑1 (MQ=255)
                                    agtgtcaCGTCAGAGTTTAGCGAATTTTTCGAGGGTGCGAATAAGCTGCGTGACGAAGCCATATTCGTTATCGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACAGCGGTAATTTCCGTTTGCGTGGCATCAAACaccgaacc                                                                                                               >  1:751381‑M1/8‑150 (MQ=255)
                                     gtgtcaCGTCAGAGTTTAGCGAATTTTTCGAGGGTGCGAATAAGCTGTGTGACGAAGCCATATTCGTTATCGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACGGCGGTAATTTCCGTTTGCGTGGCATCAAAAaccgaacc                                                                                                               <  2:715134‑M1/143‑1 (MQ=255)
                                         cgcgTCAGAGTTTAGCGAATTTTTCGAGGGTGCGAATAAGCTGTGTGACGAAGCCATATTCGTTATCGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACGGCGGTAATTTTCGTTTTCTTGGATCCAAACACCGAATCCaaaa                                                                                                          >  2:31844/1‑149 (MQ=255)
                                         cgcgTCAGAGTTTAGCGAATTTTTCGAGGGTGCGAATAAGCTGTGTGACGAAGCCATATTCGTTATCGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACGGCGGTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAAt                                                                                                          <  1:338006/150‑1 (MQ=255)
                                         cgcgTCAGAGTTTAGCGAATTTTTCGAGGGTGCGAATAAGCTGTGTGACGAAGCCATATTCGTTATCGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACGGCGGTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAAt                                                                                                          >  2:411274/1‑150 (MQ=255)
                                         cgcgTCAGAGTTTAGCGAATTTTTCGAGGGTGCGAATAAGCTGTGTGACGAAGCCATATTCGTTATCGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACGGCGGTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAAt                                                                                                          >  2:397614/1‑150 (MQ=255)
                                            gTCAGAGTTTAGCGAATTTTTCGAGGGTGCGAATAAGCTGCGTGACGAAGCCATATTCGTTATCGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACAGCGGTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGc                                                                                                       <  1:342808/150‑1 (MQ=255)
                                               agagTTTAGCGAATTTTTCGAGGGTGCGAATAAGCTGTGTGACGAAGCCATATTCGTTATCGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACGGCGGTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTg                                                                                                     <  2:443842/149‑1 (MQ=255)
                                                  gTTTAGCGAATTTTTCGAGGGTGCGAATAAGCTGTGTGACGAAGCCATATTCGTTATCGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACGGCGGTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCaa                                                                                                 >  1:129559/1‑150 (MQ=255)
                                                       gCGAATTTTTCGAGGGTGCGAATAAGCTGTGTGACGAAGCCATATTCGTTATCGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACGGCGGTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGata                                                                                            <  1:354903/150‑1 (MQ=255)
                                                         gAATTTTTCGAGGGTGCGAATAAGCTGTGTGACGAAGCCATATTCGTTATCGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACGGCGGTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATc                                                                                          <  1:130996/150‑1 (MQ=255)
                                                          aaTTTTTCGAGGGTGCGAATAAGCTGCGTGACGAAGCCATATTCGTTATCGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACAGCGGTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCg                                                                                         <  1:424107/150‑1 (MQ=255)
                                                           aTTTTTCGAGGGTGCGAATAAGCTGTGTGACGAAGCCATATTCGTTATCGCACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACGGCGGTAATTTCCGTTTGCGCGGCATCAAACACCGAACCGAAATGGCTGCCACTGATATCg                                                                                         >  1:133843/1‑149 (MQ=255)
                                                            tttttCGAGGGTGCGAATAAGCTGCGTGACGAAGCCATATTCGTTATCGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACAGCGGTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGa                                                                                       <  1:277090/150‑1 (MQ=255)
                                                               ttCGAGGGTGCGAATAAGCTGCGTGACGAAGCCATATTCGTTATCGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACAGCGGTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAAt                                                                                                >  1:929500/1‑138 (MQ=255)
                                                                 cGAGGGTGCGAATAAGCTGCGTGACGAAGCCATATTCGTTATCGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACAGCGGTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAagag                                                                                   <  1:383254/149‑1 (MQ=255)
                                                                 cGAGGGTGCGAATAAGCTGCGTGACGAAGCCATATTCGTTATCGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACAGCGGTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAgaga                                                                                  >  1:80860/1‑150 (MQ=255)
                                                                    gggTGCGAATAAGCTGCGTGACGAAGCCATATTCGTTATCGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACAGCGGTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTa                                                                               >  2:749360/1‑150 (MQ=255)
                                                                     ggTGCGAATAAGCTGTGTGACGAAGCCATATTCGTTATCGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACGGCGGTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTAt                                                                              <  2:350955/150‑1 (MQ=255)
                                                                         cGAATAAGCTGTGTGACGAAGCCATATTCGTTATCGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACGGCGGTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATttct                                                                          >  2:501287/1‑150 (MQ=255)
                                                                         cGAATAAGCTGTGTGACGAAGCCATATTCGTTATCGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACGGCGGTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATttc                                                                           <  1:715134/149‑1 (MQ=255)
                                                                            aTAAGCTGCGTGACGAAGCCATATTCGTTATCGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACAGCGGTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCa                                                                       >  2:330594/1‑150 (MQ=255)
                                                                              aaGCTGTGTGACGAAGCCATATTCGTTATCGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACGGCGGTAATTTCCGTTTGCGGGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATc                                                                     >  1:377088/1‑150 (MQ=255)
                                                                               aGCTGTGTGACGAAGCCATATTCGTTATCGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACGGCGGTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCg                                                                    <  2:981565/150‑1 (MQ=255)
                                                                                gCTGTGTGACGAAGCCATATTCGTTATCGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACGGCGGTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCgg                                                                   >  1:632979/1‑150 (MQ=255)
                                                                                gCTGCGTGACGAAGCCATATTCGTTATCGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACAGCGGTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCgg                                                                   >  1:785307/1‑150 (MQ=255)
                                                                                    cgtgACGAAGCCATATTCGTTATCGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACAGCGGTAATTTCCGTTTGCGTGGCATCAAACACCGATCCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGta                                                                 <  2:172578/147‑1 (MQ=255)
                                                                                    cgtgACGAAGCCATATTCGTTATCGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACAGCGGTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGtat                                                                <  1:617111/148‑1 (MQ=255)
                                                                                     gtgACGAAGCCATATTCGTTATCGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACAGCGGTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATaa                                                              >  1:161623/1‑150 (MQ=255)
                                                                                       aaCGAAGCCATATTCGTTATCGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACGGCGGTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAAc                                                             >  2:508512/2‑149 (MQ=255)
                                                                                        aCGAAGCCATATTCGTTATCGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACGGCGGTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAAcc                                                            >  2:302189/1‑149 (MQ=255)
                                                                                          gAAGCCATATTCGTTATCGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACAGCGGTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCaaa                                                         >  1:460532/1‑150 (MQ=255)
                                                                                           aaGCCATATTCGTTATCGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACGGCGGTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAAt                                                        <  2:752636/150‑1 (MQ=255)
                                                                                             gCCATATTCGTTATCGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACAGCGGTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGa                                                      >  2:89593/1‑150 (MQ=255)
                                                                                               cATATTCGTTATCGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACGGCGGTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACt                                                    <  2:370195/150‑1 (MQ=255)
                                                                                               cATATTCGTTATCGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACGGCGGTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACt                                                    >  1:24375/1‑150 (MQ=255)
                                                                                               cATATTCGTTATCGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACGGCGGTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACt                                                    <  2:911425/150‑1 (MQ=255)
                                                                                               cATATTCGTTATCGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACAGCGGTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACt                                                    >  2:270349/1‑150 (MQ=255)
                                                                                                atatTCGTTATCGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACGGCGGTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGGTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTc                                                   >  1:796014/1‑150 (MQ=255)
                                                                                                atatTCGTTATCGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACGGCGGTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTc                                                   <  1:716695/150‑1 (MQ=255)
                                                                                                atatTCGTTATCGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACAGCGGTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTc                                                   >  1:358307/1‑150 (MQ=255)
                                                                                                atatTCGTTATCGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACAGCGGTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTc                                                   <  2:15371/150‑1 (MQ=255)
                                                                                                atatTCGTTATCGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACAGCGGTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTc                                                   >  1:894143/1‑150 (MQ=255)
                                                                                                      gTTATCGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACGGCGGTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAGTGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGGCACGTTAt                                              <  2:566285/149‑1 (MQ=255)
                                                                                                       ttATCGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACAGCGGTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATtg                                            >  1:179956/1‑150 (MQ=255)
                                                                                                          tCGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACAGCGGTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATtggt                                          <  2:80860/149‑1 (MQ=255)
                                                                                                            gTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACGGCGGTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATtggtgg                                        <  1:397614/149‑1 (MQ=255)
                                                                                                            gTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACGGCGGTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTggtggt                                       <  2:27045/150‑1 (MQ=255)
                                                                                                            gTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACGGCGGTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTggtggt                                       <  2:24488/150‑1 (MQ=255)
                                                                                                                cAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACGGCGGTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCt                                   <  2:577838/150‑1 (MQ=255)
                                                                                                                cAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACAGCGGTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCt                                   >  2:836595/1‑150 (MQ=255)
                                                                                                                 aGGCGACCGTTTTCACCAGTTGTAAATCGCCCACAGCGGTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCtt                                  >  2:621267/1‑150 (MQ=255)
                                                                                                                  ggCGACCGTTTTCACCAGTTGTAAATCGCCCACAGCGGTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTg                                 <  1:221226/150‑1 (MQ=255)
                                                                                                                   gCGACCGTTTTCACCAGTTGTAAATCGCCCACGGCGGTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGt                                <  1:508512/150‑1 (MQ=255)
                                                                                                                       ccGTTTTCACCAGTTGTAAATCGCCCACAGCGGTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGttt                              >  1:314862/1‑148 (MQ=255)
                                                                                                                        cGTTTTCACCAGTTGTAAATCGCCCACGGCGGTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGGGGTTGCTTGTTTAAg                           >  2:149522/1‑150 (MQ=255)
                                                                                                                        cGTTTTCACCAGTTGTAAATCGCCCACAGCGGTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAg                           <  1:423889/150‑1 (MQ=255)
                                                                                                                         gTTTTCACCAGTTGTAAATCGCCCACGGCGGTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCCATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTTGTGGTTGCTTGTTTaatt                          >  2:601519/1‑148 (MQ=255)
                                                                                                                         gTTTTCACCAGTTGTAAATCGCCCACAGCGGTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGt                          <  2:969066/150‑1 (MQ=255)
                                                                                                                          ttttCACCAGTTGTAAATCGCCCACGGCGGTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTg                         >  1:888670/1‑150 (MQ=255)
                                                                                                                           tttCACCAGTTGTAAATCGCCCACGGCGGTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGc                        >  1:617053/1‑150 (MQ=255)
                                                                                                                           tttCACCAGTTGTAAATCGCCCACAGCGGTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGc                        >  1:298261/1‑150 (MQ=255)
                                                                                                                                 cAGTTGTAAATCGCCCACGGCGGTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTAtt                  <  1:586328/150‑1 (MQ=255)
                                                                                                                                        aaaTCGCCCACAGCGGTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACctct           >  1:307903/1‑150 (MQ=255)
                                                                                                                                         aaTCGCCCACAGCGGTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACctct           <  2:567373/149‑1 (MQ=255)
                                                                                                                                          aTCGCCCACAGCGGTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTc         >  1:289245/1‑150 (MQ=255)
                                                                                                                                             gCCCACGGCGGTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGc      <  2:77145/150‑1 (MQ=255)
                                                                                                                                             gCCCACGGCGGTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGc      <  1:186527/150‑1 (MQ=255)
                                                                                                                                             gCCCACAGCGGTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCgg       >  2:53848/1‑149 (MQ=255)
                                                                                                                                              cccACGGCGGTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCa     <  2:426234/150‑1 (MQ=255)
                                                                                                                                                cACGGCGGTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGt   <  1:933872/150‑1 (MQ=255)
                                                                                                                                                 aCGGCGGTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTc  <  2:129559/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TTAATGAGTGAAAGAGGCGGAGGTTTTTTCCTCCGCCTGTGCGCGTCAGAGTTTAGCGAATTTTTCGAGGGTGCGAATAAGCTGTGTGACGAAGCCATATTCGTTATCGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACGGCGGTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTC  >  CP009273/1483925‑1484219

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