Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,513,935 A→T 100% Y73N (TAC→AAC)  gltX ← glutamyl‑tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,513,9350AT96.1% 961.1 / 18.4 296Y73N (TAC→AAC) gltXglutamyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base A (7/5);  new base T (148/136);  total (155/141)
Fisher's exact test for biased strand distribution p-value = 7.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 302 total aligned reads displayed.

CATAACGCGGCTTCTCACCTTTCGCCATTTGCTCTTCGCGCAGCGCTTCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTACGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCTCAAGATCGGTGTCTTCAATACGCAGCACGAACTCACCGCCGTGGTTACGTGCAAAAAGCCAGGAGTAAA  >  CP009273/2513788‑2514081
                                                                                                                                                   |                                                                                                                                                  
cATAACGCGGCTTCTCACCTTTCGCCATTTGCTCTTCGCGCAGCGCTTCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCg                                                                                                                                                  <  2:1439424/150‑1 (MQ=255)
cATAACGCGGCTTCTCACCTTTCGCCATTTGCTCTTCGCGCAGCGCTTCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCg                                                                                                                                                  <  2:770256/150‑1 (MQ=255)
      gcgGCTTCTCACCTTTCGCCATTTGCTCTTCGCGCAGCGCTTCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCtt                                                                                                                                            >  1:1323688/1‑150 (MQ=255)
       cgGCTTCTCACCTTTCGCCATTTGCTCTTCGCGCAGCGCTTCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTc                                                                                                                                           >  2:12052/1‑150 (MQ=255)
       cgGCTTCTCACCTTTCGCCATTTGCTCTTCGCGCAGCGCTTCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTc                                                                                                                                           >  1:1966993/1‑150 (MQ=255)
       cgGCTTCTCACCTTTCGCCATTTGCTCTTCGCGCAGCGCTTCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTc                                                                                                                                           >  2:3938619/1‑150 (MQ=255)
         gCTTCTCACCTTTCGCCATTTGCTCTTCGCGCAGCGCTTCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCAt                                                                                                                                         <  1:5011843/150‑1 (MQ=255)
              tcACCTTTCGCCATTTGCTCTTCGCGCAGCGCTTCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCAc                                                                                                                                    >  1:86538/1‑150 (MQ=255)
               cACCTTTCGCCATTTGCTCTTCGCGCAGCGCTTCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACt                                                                                                                                   >  1:4571329/1‑150 (MQ=255)
                  cTTTCGCCATTTGCTCTTCGCGCAGCGCTTCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCa                                                                                                                                >  1:3751998/1‑150 (MQ=255)
                      cGCCATTTGCTCTTCGCGCAGCGCTTCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCt                                                                                                                            >  1:3107854/1‑150 (MQ=255)
                      cGCCATTTGCTCTTCGCGCAGCGCTTCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCt                                                                                                                            >  1:2914595/1‑150 (MQ=255)
                             tGCTCTTCGCGCAGCGCTTCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAg                                                                                                                     >  2:4273709/1‑150 (MQ=255)
                             tGCTCTTCGCGCAGCGCTTCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAg                                                                                                                     >  2:1695805/1‑150 (MQ=255)
                             tGCTCTTCGCGCAGCGCTTCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTACCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAg                                                                                                                     <  2:3276526/150‑1 (MQ=255)
                                tctTCGCGCAGCGCTTCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTc                                                                                                                  <  1:2092593/150‑1 (MQ=255)
                                   tCGCGCAGCGCTTCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATg                                                                                                               <  2:4788123/150‑1 (MQ=255)
                                    cgcgCAGCGCTTCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGc                                                                                                              >  1:1461434/1‑150 (MQ=255)
                                     gcgcAGCGCTTCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGcc                                                                                                             >  1:3779868/1‑150 (MQ=255)
                                      cgcAGCGCTTCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTATCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGcca                                                                                                            >  1:1227153/1‑150 (MQ=255)
                                       gcAGCGCTTCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTACGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGccat                                                                                                           >  1:58050/1‑150 (MQ=255)
                                        cAGCGCTTCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGccatc                                                                                                          >  2:3656336/1‑150 (MQ=255)
                                           cgcTTCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGccatccat                                                                                                       <  1:3273147/150‑1 (MQ=255)
                                           cgcTTCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGccatccat                                                                                                       >  1:1684230/1‑150 (MQ=255)
                                           cgcTTCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGccatccat                                                                                                       <  2:377/150‑1 (MQ=255)
                                            gcTTCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATa                                                                                                      >  2:3972482/1‑150 (MQ=255)
                                             cTTCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATaa                                                                                                     <  2:276259/150‑1 (MQ=255)
                                             cTTCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATaa                                                                                                     >  2:1778960/1‑150 (MQ=255)
                                               tCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATg                                                                                                   <  2:2999860/150‑1 (MQ=255)
                                               tCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATg                                                                                                   <  2:4071807/150‑1 (MQ=255)
                                                 cAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGc                                                                                                 <  2:2240123/150‑1 (MQ=255)
                                                    gcgcTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTc                                                                                              <  2:1904070/150‑1 (MQ=255)
                                                    gcgcTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTc                                                                                              <  1:518568/150‑1 (MQ=255)
                                                     cgcTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCg                                                                                             <  2:1819445/150‑1 (MQ=255)
                                                        tctTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATa                                                                                          <  1:1741345/150‑1 (MQ=255)
                                                        tctTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATa                                                                                          <  2:4713755/150‑1 (MQ=255)
                                                         ctTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGCCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAg                                                                                         >  2:4936228/1‑150 (MQ=255)
                                                            tAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCtt                                                                                      >  2:2827389/1‑150 (MQ=255)
                                                             agagCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCGTCGATAGCTTc                                                                                     >  2:512155/1‑150 (MQ=255)
                                                               agCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCg                                                                                   >  1:4838939/1‑150 (MQ=255)
                                                               agCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCg                                                                                   <  2:4170939/150‑1 (MQ=255)
                                                               agCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCg                                                                                   >  1:661834/1‑150 (MQ=255)
                                                                  aaTAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCg                                                                                <  2:3867809/150‑1 (MQ=255)
                                                                   aTAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGt                                                                               <  2:2380468/150‑1 (MQ=255)
                                                                    tAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTg                                                                              <  1:4029053/150‑1 (MQ=255)
                                                                         tttATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACg                                                                         <  1:4709310/150‑1 (MQ=255)
                                                                         tttATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACg                                                                         <  2:2117566/150‑1 (MQ=255)
                                                                           tctaAGCAGGGCCCTCGTCCAGCATCGGATCGAGCACCGCGTTGTAGCGATCAAAACGTATGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCt                                                                       <  2:4416331/148‑1 (MQ=255)
                                                                           tataAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCt                                                                       <  1:733162/150‑1 (MQ=255)
                                                                            ataAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCTc                                                                      >  1:4808832/1‑150 (MQ=255)
                                                                            ataAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCTc                                                                      >  1:5001543/1‑150 (MQ=255)
                                                                             taAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCTCa                                                                     >  2:3423761/1‑150 (MQ=255)
                                                                             taAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCTCa                                                                     <  1:38335/150‑1 (MQ=255)
                                                                               aGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCTCAAg                                                                   >  2:1959639/1‑150 (MQ=255)
                                                                                  aGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCTCAAGATc                                                                >  1:4401807/1‑150 (MQ=255)
                                                                                    tGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCTCAAGATCgg                                                              >  2:1154485/1‑150 (MQ=255)
                                                                                     gCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTACGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCTCAAGATCGgt                                                             <  2:1295956/150‑1 (MQ=255)
                                                                                      cccTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCTCAAGATCGgtg                                                            <  1:4469123/150‑1 (MQ=255)
                                                                                       ccTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCTCAAGATCGgtgt                                                           >  1:975491/1‑150 (MQ=255)
                                                                                        ctctTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCTCAAGATCGGTGTc                                                          <  2:522894/150‑1 (MQ=255)
                                                                                          ctTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCTCAAGATCGGTGTCtt                                                        <  2:283737/150‑1 (MQ=255)
                                                                                           ttCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCTCAAGATCGGTGTCTTc                                                       <  1:3438195/150‑1 (MQ=255)
                                                                                            tCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCTCAAGATCGGTGTCTTCa                                                      >  2:1977982/1‑150 (MQ=255)
                                                                                             ccAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCTCAAGATCGGTGTCTTCaa                                                     >  1:334050/1‑150 (MQ=255)
                                                                                               aGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCTCAAGATCGGTGTCTTCAATa                                                   <  1:2843560/150‑1 (MQ=255)
                                                                                                   tCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCTCAAGATCGGTGTCTTCAATACgca                                               <  2:593276/150‑1 (MQ=255)
                                                                                                     tGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCTCAAGATCGGTGTCTTCAATACgcagc                                             <  2:4391151/150‑1 (MQ=255)
                                                                                                      gatcgatcACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCTCAAGATCGGTGTCTTCAATACgcagca                                            >  2:4198596/1‑150 (MQ=255)
                                                                                                         cgatcACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCTCAAGATCGGTGTCTTCAATACGCAGCACGa                                         <  2:2576832/150‑1 (MQ=255)
                                                                                                          gatcACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCTCAAGATCGGTGTCTTCAATACGCAGCACGaa                                        <  2:5014128/150‑1 (MQ=255)
                                                                                                          gatcACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCTCAAGATCGGTGTCTTCAATACGCAGCACGaa                                        >  2:1092286/1‑150 (MQ=255)
                                                                                                           atcACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCTCAAGATCGGTGTCTTCAATACGCAGCACGAAc                                       <  1:2261481/150‑1 (MQ=255)
                                                                                                           atcACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCTCAAGATCGGTGTCTTCAATACGCAGCACGAAc                                       <  1:4171556/150‑1 (MQ=255)
                                                                                                            tcACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCTCAAGATCGGTGTCTTCAATACGCAGCACGAACt                                      <  1:2522392/150‑1 (MQ=255)
                                                                                                              aCCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCTCAAGATCGGTGTCTTCAATACGCAGCACGAACTCa                                    >  1:212473/1‑150 (MQ=255)
                                                                                                               ccGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCTCAAGATCGGTGTCTTCAATACGCAGCACGAACTCAc                                   >  2:3195040/1‑150 (MQ=255)
                                                                                                                cgcgTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCTCAAGATCGGTGTCTTCAATACGCAGCACGAACTCAcc                                  >  2:2875939/1‑150 (MQ=255)
                                                                                                                 gcgTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCTCAAGATCGGTGTCTTCAATACGCAGCACGAACTCAccg                                 <  2:1509498/150‑1 (MQ=255)
                                                                                                                  cgTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCTCAAGATCGGTGTCTTCAATACGCAGCACGAACTCAccgc                                <  2:2655312/150‑1 (MQ=255)
                                                                                                                   gTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCTCAAGATCGGTGTCTTCAATACGCAGCACGAACTCAccgcc                               <  1:4238356/150‑1 (MQ=255)
                                                                                                                   gTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCTCAAGATCGGTGTCTTCAATACGCAGCACGAACTCAccgcc                               >  1:2441/1‑150 (MQ=255)
                                                                                                                     tGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCTCAAGATCGGTGTCTTCAATACGCAGCACGAACTCACCGCCGt                             >  2:856328/1‑150 (MQ=255)
                                                                                                                      gTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCTCAAGATCGGTGTCTTCAATACGCAGCACGAACTCACCGCCGTg                            <  2:2805251/150‑1 (MQ=255)
                                                                                                                      gTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCTCAAGATCGGTGTCTTCAATACGCAGCACGAACTCACCGCCGTg                            <  2:3107854/150‑1 (MQ=255)
                                                                                                                       tAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCTCAAGATCGGTGTCTTCAATACGCAGCACGAACTCACCGCCGTgg                           >  2:752546/1‑150 (MQ=255)
                                                                                                                        aGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCTCAAGATCGGTGTCTTCAATACGCAGCACGAACTCACCGCCGTGGt                          >  1:611847/1‑150 (MQ=255)
                                                                                                                         gCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCTCAAGATCGGTGTCTTCAATACGCAGCACGAACTCACCGCCGTGGtt                         <  1:12052/150‑1 (MQ=255)
                                                                                                                         gCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCTCAAGATCGGTGTCTTCAATACGCAGCACGAACTCACCGCCGTGGtt                         <  1:2358417/150‑1 (MQ=255)
                                                                                                                           gATCAAAACGTTTGGTCTGGTAGTACGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCTCAAGATCGGTGTCTTCAATACGCAGCACGAACTCACCGCCGTGGTTAc                       >  1:1126551/1‑150 (MQ=255)
                                                                                                                            aTCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCTCAAGATCGGTGTCTTCAATACGCAGCACGAACTCACCGCCGTGGTTACg                      <  2:2662623/150‑1 (MQ=255)
                                                                                                                              cAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCTCAAGATCGGTGTCTTCAATACGCAGCACGAACTCACCGCCGTGGTTACGTg                    >  2:3340277/1‑150 (MQ=255)
                                                                                                                              cAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCTCAAGATCGGTGTCTTCAATACCCAGCACGAACTCACCGCCGTGGTTACGTg                    <  2:87347/150‑1 (MQ=255)
                                                                                                                                 aaCGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCTCAAGATCGGTGTCTTCAATACGCAGCACGAACTCACCGCCGTGGTTACGTGCaa                 >  2:2468604/1‑150 (MQ=255)
                                                                                                                                   cGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCTCAAGATCGGTGTCTTCAATACGCAGCACGAACTCACCGCCGTGGTTACGTGCaaaa               >  2:4529014/1‑150 (MQ=255)
                                                                                                                                        ggTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCTCAAGATCGGTGTCTTCAATACGCAGCACGAACTCACCGCCGTGGTTACGTGCAAAAAGCCa          <  2:3694679/150‑1 (MQ=255)
                                                                                                                                           cTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCTCAAGATCGGTGTCTTCAATACGCAGCACGAACTCACCGCCGTGGTTACGTGCAAAAAGCCAGGa       >  1:3651908/1‑150 (MQ=255)
                                                                                                                                           cTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCTCAAGATCGGTGTCTTCAATACGCAGCACGAACTCACCGCCGTGGTTACGTGCAAAAAGCCAGGa       >  1:2733446/1‑150 (MQ=255)
                                                                                                                                            tGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCTCAAGATCGGTGTCTTCAATACGCAGCACGAACTCACCGCCGTGGTTACGTGCAAAAAGCCAGGAg      >  2:3495053/1‑150 (MQ=255)
                                                                                                                                            tGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCTCAAGATCGGTGTCTTCAATACGCAGCACGAACTCACCGCCGTGGTTACGTGCAAAAAGCCAGGAg      <  1:4106548/150‑1 (MQ=255)
                                                                                                                                                agttCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCTCAAGATCGGTGTCCTCAATACGCAGCACGAACTCACCGCCGTGGTTACGTGCAAAAAGCCAGGAGTaaa  >  1:184878/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CATAACGCGGCTTCTCACCTTTCGCCATTTGCTCTTCGCGCAGCGCTTCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTACGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCTCAAGATCGGTGTCTTCAATACGCAGCACGAACTCACCGCCGTGGTTACGTGCAAAAAGCCAGGAGTAAA  >  CP009273/2513788‑2514081

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: ‑