Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,733,753 C→T 100% E253K (GAG→AAG)  aroF ← 3‑deoxy‑D‑arabino‑heptulosonate‑7‑phosphate synthase, tyrosine‑repressible

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,733,7530CT93.0% 859.2 / 39.5 257E253K (GAG→AAG) aroF3‑deoxy‑D‑arabino‑heptulosonate‑7‑phosphate synthase, tyrosine‑repressible
Reads supporting (aligned to +/- strand):  ref base C (13/5);  new base T (135/104);  total (148/109)
Fisher's exact test for biased strand distribution p-value = 2.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01
Rejected as polymorphism: Polymorphic base substitution creates a homopolymer stretch.

Only 100 of 269 total aligned reads displayed.

TCATCAGACCAATAATTGAGCGATTGCCATCTTTGATTTGAGCAACCACGGATTCTGCCACCGCAGGCTGACGGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTCTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAGTTTGTAGCAACGCAACCTGCCCTGCCTGGTTAATGCCAACAAAACGGTGCGGCTGGGCG  >  CP009273/2733606‑2733898
                                                                                                                                                   |                                                                                                                                                 
tcatcaGACCAATAATTGAGCGATTGCCATCTTTGATTTGAGCAACCACGGATTCTGCCACCGCAGGCTGACGGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCtttt                                                                                                                                                 >  1:3252049/1‑150 (MQ=255)
 catcaGACCAATAATTGAGCGATTGCCATCTTTGATTTGAGCAACCACGGATTCTGCCACCGCAGGCTGACGGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTtttt                                                                                                                                                >  2:3065376/1‑150 (MQ=255)
    caGACCAATAATTGAGCGATTGCCATCTTTGATTTGAGCAACCACGGATTCTGCCACCGCAGGCTGACGGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTc                                                                                                                                             <  2:4225942/150‑1 (MQ=255)
     aGACCAATAATTGAGCGATTGCCATCTTTGATTTGAGCAACCACGGATTCTGCCACCGCAGGCTGACGGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTca                                                                                                                                            >  2:3969726/1‑150 (MQ=255)
     aGACCAATAATTGAGCGATTGCCATCTTTGATTTGAGCAACCACGGATTCTGCCACCGCAGGCTGACGGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTca                                                                                                                                            >  2:2838822/1‑150 (MQ=255)
       aCCAATAATTGAGCGATTGCCATCTTTGATTTGAGCAACCACGGATTCTGCCACCGCAGGCTGACGGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTcaca                                                                                                                                          <  2:144722/150‑1 (MQ=255)
         cAATAATTGAGCGATTGCCATCTTTGATTTGAGCAACCACGGATTCTGCCACCGCAGGCTGACGGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACAtt                                                                                                                                        <  1:3070642/150‑1 (MQ=255)
         cAATAATTGAGCGATTGCCATCTTTGATTTGAGCAACCACGGATTCTGCCACCGCAGGCTGACGGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACAtt                                                                                                                                        <  1:3511036/150‑1 (MQ=255)
         cAATAATTGAGCGATTGCCATCTTTGATTTGAGCAACCACGGATTCTGCCACCGCAGGCTGACGGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACAtt                                                                                                                                        <  1:3926934/150‑1 (MQ=255)
         cAATAATTGAGCGATTGCCATCTTTGATTTGAGCAACCACGGATTCTGCCACCGCAGGCTGACGGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACAtt                                                                                                                                        >  1:4649720/1‑150 (MQ=255)
           ataatTGAGCGATTGCCATCTTTGATTTGAGCAACCACGGATTCTGCCACCGCAGGCTGACGGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTgc                                                                                                                                      >  2:1030466/1‑150 (MQ=255)
             aatTGAGCGATTGCCATCTTTGATTTGAGCAACCACGGATTCTGCCACCGCAGGCTGACGGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTgcgc                                                                                                                                    >  1:3169663/1‑150 (MQ=255)
                 gAGCGATTGCCATCTTTGATTTGAGCAACCACGGATTCTGCCACCGCAGGCTGACGGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACa                                                                                                                                <  1:2101011/150‑1 (MQ=255)
                 gAGCGATTGCCATCTTTGATTTGAGCAACCACGGATTCTGCCACCGCAGGCTGACGGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACa                                                                                                                                <  1:1626278/150‑1 (MQ=255)
                    cGATTGCCATCTTTGATTTGAGCAACCACGGATTCTGCCACCGCAGGCTGACGGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATcc                                                                                                                             >  1:4736582/1‑150 (MQ=255)
                    cGATTGCCATCTTTGATTTGAGCAACCACGGATTCTGCCACCGCAGGCTGACGGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTCTTTTTCACATTGCGCAACATcc                                                                                                                             >  2:3953302/1‑150 (MQ=255)
                      aTTGCCATCTTTGATTTGAGCAACCACGGATTCTGCCACCGCAGGCTGACGGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGc                                                                                                                           >  2:1019790/1‑150 (MQ=255)
                       ttGCCATCTTTGATTTGAGCAACCACGGATTCTGCCACCGCAGGCTGACGGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCa                                                                                                                          >  1:2633882/1‑150 (MQ=255)
                       ttGCCATCTTTGATTTGAGCAACCACGGATTCTGCCACCGCAGGCTGACGGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCa                                                                                                                          >  1:1270882/1‑150 (MQ=255)
                        tGCCATCTTTGATTTGAGCAACCACGGATTCTGCCACCGCAGGCTGACGGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAg                                                                                                                         <  1:3736737/150‑1 (MQ=255)
                         gCCATCTTTGATTTGAGCAACCACGGATTCTGCCACCGCAGGCTGACGGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAgg                                                                                                                        <  2:211413/150‑1 (MQ=255)
                          ccATCTTTGATTTGAGCAACCACGGATTCTGCCACCGCAGGCTGACGGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAggg                                                                                                                       <  1:101352/150‑1 (MQ=255)
                           cATCTTTGATTTGAGCAACCACGGATTCTGCCACCGCAGGCTGACGGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGc                                                                                                                      >  2:1670032/1‑150 (MQ=255)
                            aTCTTTGATTTGAGCAACCACGGATTCTGCCACCGCAGGCTGACGGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCt                                                                                                                     <  2:1099603/150‑1 (MQ=255)
                            aTCTTTGATTTGAGCAACCACGGATTCTGCCACCGCAGGCTGACGGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCt                                                                                                                     <  2:911921/150‑1 (MQ=255)
                             tCTTTGATTTGAGCAACCACGGATTCTGCCACCGCAGGCTGACGGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCta                                                                                                                    <  2:3176686/150‑1 (MQ=255)
                                 tgatttgaGCAACCACGGATTCTGCCACCGCAGGCTGACGGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGt                                                                                                                >  2:2502882/1‑150 (MQ=255)
                                        aGCAACCACGGATTCTGCCACCGCAGGCTGACGGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGgcgc                                                                                                         >  2:4622722/1‑150 (MQ=255)
                                                cGGATTCTGCCACCGCAGGCTGACGGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTaccac                                                                                                 >  1:2952468/1‑150 (MQ=255)
                                                cGGATTCTGCCACCGCAGGCTGACGGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTaccac                                                                                                 >  1:60877/1‑150 (MQ=255)
                                                 ggATTCTGCCACCGCAGGCTGACGGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTaccacc                                                                                                <  1:360515/150‑1 (MQ=255)
                                                 ggATTCTGCCACCGCAGGCTGACGGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTaccacc                                                                                                <  1:638804/150‑1 (MQ=255)
                                                         ccACCGCAGGCTGACGGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCATATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGa                                                                                        <  1:379606/150‑1 (MQ=255)
                                                            ccGCAGGCTGACGGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATca                                                                                     <  1:2307562/150‑1 (MQ=255)
                                                                aGGCTGACGGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATg                                                                                 >  2:767002/1‑150 (MQ=255)
                                                                aGGCTGACGGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATg                                                                                 >  2:1763673/1‑150 (MQ=255)
                                                                  gCTGACGGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGc                                                                               >  1:727100/1‑150 (MQ=255)
                                                                   cTGACGGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGcc                                                                              <  1:3470723/150‑1 (MQ=255)
                                                                   cTGACGGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGcc                                                                              <  2:2714630/150‑1 (MQ=255)
                                                                   cTGACGGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGcc                                                                              <  2:4174967/150‑1 (MQ=255)
                                                                    tGACAGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCg                                                                             <  1:3317583/150‑1 (MQ=255)
                                                                      aCGGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTc                                                                           >  1:1092183/1‑150 (MQ=255)
                                                                       cGGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTcc                                                                          >  2:4048083/1‑150 (MQ=255)
                                                                       cGGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTcc                                                                          >  2:1128531/1‑150 (MQ=255)
                                                                       cGGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTcc                                                                          >  1:3866714/1‑150 (MQ=255)
                                                                       cGGCGATAATCTTTATTGGAATTAACGTGGCTGCAATCTACAATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTcc                                                                          >  2:2232317/1‑150 (MQ=255)
                                                                        ggCGAGAATCTTTATTGGAATTACCGTGGCTGCACTCTACCATCCGAGACGGGCGCAGTCCCGCCTGTTCCCTCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATCGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTcct                                                                         <  2:1805061/150‑2 (MQ=255)
                                                                           gATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGAt                                                                      >  2:246274/1‑150 (MQ=255)
                                                                            aTAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGAtt                                                                     >  1:2683327/1‑150 (MQ=255)
                                                                               aTCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTccc                                                                  >  1:1370946/1‑150 (MQ=255)
                                                                               aTCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTccc                                                                  >  2:3464424/1‑150 (MQ=255)
                                                                                tCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTcccc                                                                 >  1:1520880/1‑150 (MQ=255)
                                                                                    tATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAg                                                             >  2:2673530/1‑150 (MQ=255)
                                                                                     aTTGGAATTACCGTGGCTGCACTCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAGt                                                            >  1:3101161/1‑150 (MQ=255)
                                                                                     aTTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAGt                                                            <  2:1524362/150‑1 (MQ=255)
                                                                                     aTTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAGt                                                            >  2:720062/1‑150 (MQ=255)
                                                                                     aTTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAGt                                                            >  1:4049429/1‑150 (MQ=255)
                                                                                        ggAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAGTTTg                                                         >  1:4044397/1‑150 (MQ=255)
                                                                                          aaTTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAGTTTGTa                                                       >  2:4401916/1‑150 (MQ=255)
                                                                                          aaTTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAGTTTGTa                                                       >  2:362866/1‑150 (MQ=255)
                                                                                           aTTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAGTTTGTAg                                                      <  2:3370109/150‑1 (MQ=255)
                                                                                             tACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAGTTTGTAgca                                                    <  1:940559/150‑1 (MQ=255)
                                                                                              aCCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAGTTTGTAgcaa                                                   >  1:1691303/1‑150 (MQ=255)
                                                                                                cGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTCTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAGTTTGTAgcaacg                                                 >  2:1887314/1‑150 (MQ=255)
                                                                                                 gTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTCTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAGTTTGTAgcaacgc                                                >  2:4665522/1‑150 (MQ=255)
                                                                                                 gTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTCTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAGTTTGTAgcaacgc                                                >  2:2774228/1‑150 (MQ=255)
                                                                                                   ggCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAGTTTGTAgcaacgcaa                                              >  2:487023/1‑150 (MQ=255)
                                                                                                     cTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAGTTTGTAGCAACGCAAcc                                            >  1:3574849/1‑150 (MQ=255)
                                                                                                     cTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAGTTTGTAGCAACGCAAcc                                            <  2:542639/150‑1 (MQ=255)
                                                                                                      tGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAGTTTGTAGCAACGCAAcct                                           <  2:2917375/150‑1 (MQ=255)
                                                                                                       gCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAGTTTGTAGCAACGCAAcctg                                          <  2:2019354/150‑1 (MQ=255)
                                                                                                         aaTCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAGTTTGTAGCAACGCAAcctgcc                                        >  1:4365517/1‑150 (MQ=255)
                                                                                                          aTCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTCTAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAGTTTGTAGCAACGCAAcctgccc                                       >  1:1611295/1‑150 (MQ=255)
                                                                                                          aTCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAGTTTGTAGCAACGCAAcctgccc                                       >  1:23725/1‑150 (MQ=255)
                                                                                                           tCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAGTTTGTAGCAACGCAAcctgccct                                      >  1:4498496/1‑150 (MQ=255)
                                                                                                           tCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAGTTTGTAGCAACGCAAcctgccct                                      <  2:3252049/150‑1 (MQ=255)
                                                                                                            cTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAGTTTGTAGCAACGCAAcctgccctg                                     <  1:1739043/150‑1 (MQ=255)
                                                                                                               ccATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAGTTTGTAGCAACGCAACCTGCcctgcct                                  <  1:4530795/150‑1 (MQ=255)
                                                                                                                 aTCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAGTTTGTAGCAACGCAACCTGCCCTGCCTgg                                >  1:4349941/1‑150 (MQ=255)
                                                                                                                  tCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAGTTTGTAGCAACGCAACCTGCCCTGCCTGGt                               >  2:4214895/1‑150 (MQ=255)
                                                                                                                   cAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAGTTTGTAGCAACGCAACCTGCCCTGCCTGGtt                              >  2:1211784/1‑150 (MQ=255)
                                                                                                                   cAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAGTTTGTAGCAACGCAACCTGCCCTGCCTGGtt                              <  2:2197538/150‑1 (MQ=255)
                                                                                                                   cAGAGACGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAGTTTGTAGCAACGCAACCTGCCCTGCCTGGtt                              <  1:1165474/150‑1 (MQ=255)
                                                                                                                     gagaCGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAGTTTGTAGCAACGCAACCTGCCCTGCCTGGTTaa                            >  2:3862347/1‑150 (MQ=255)
                                                                                                                       gaCGGGCGCAGTCCCGCCTGTTCCATCTCTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAGTTTGTAGCAACGCAACCTGCCCTGCCTGGTCAATg                          <  2:2947962/150‑1 (MQ=255)
                                                                                                                        aCGGGCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAGTTTGTAGCAACGCAACCTGCCCTGCCTGGTTAATGc                         <  2:2592689/150‑1 (MQ=255)
                                                                                                                          gggCGCAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAGTTTGTAGCAACGCAACCTGCCCTGCCTGGTTAATGCca                       >  2:2240937/1‑150 (MQ=255)
                                                                                                                             cgcAGTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAGTTTGTAGCAACGCAACCTGCCCTGCCTGGTTAATGCcaaca                    <  2:3343399/150‑1 (MQ=255)
                                                                                                                               cAGTCCCGCCTGTTCCATCTCTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAGTTTGTAGCAACGCAACCTGCCCTGCCTGGTTAATGCcaacaac                  >  2:3692178/1‑149 (MQ=255)
                                                                                                                                 gTCCCGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAGTTTGTAGCAACGCAACCTGCCCTGCCTGGTTAATGCCAACAAAAc                >  2:2083203/1‑150 (MQ=255)
                                                                                                                                    ccGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAGTTTGTAGCAACGCAACCTGCCCTGCCTGGTTAATGCCAACAAAACGGt             <  1:2172768/150‑1 (MQ=255)
                                                                                                                                     cGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAGTTTGTAGCAACGCAACCTGCCCTGCCTGGTTAATGCCAACAAAACGGTg            >  1:2790976/1‑150 (MQ=255)
                                                                                                                                     cGCCTGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAGTTTGTAGCAACGCAACCTGCCCTGCCTGGTTAATGCCAACAAAACGGTg            >  2:2965715/1‑150 (MQ=255)
                                                                                                                                         tGTTCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAGTTTGTAGCAACGCAACCTGCCCTGCCTGGTTAATGCCAACAAAACGGTGCGGc        >  1:1979062/1‑150 (MQ=255)
                                                                                                                                           ttCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAGTTTGTAGCAACGCAACCTGCCCTGCCTGGTTAATGCCAACAAAACGGTGCGGCTg      >  1:106355/1‑150 (MQ=255)
                                                                                                                                            tCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAGTTTGTAGCAACGCACCCTGCCCTGCCTGGTTAATGCCAACAAAACGGTGCGGCTgg     <  1:2826876/150‑1 (MQ=255)
                                                                                                                                            tCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAGTTTGTAGCAACGCAACCTGCCCTGCCTGGTTAATGCCAACAAAACGGTGCGGCTgg     >  2:333625/1‑150 (MQ=255)
                                                                                                                                            tCCATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAGTTTGTAGCAACGCAACCTGCCCTGCCTGGTTAATGCCAACAAAACGGTGCGGCTgg     >  1:3373875/1‑150 (MQ=255)
                                                                                                                                              cATCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAGTTTGTAGCAACGCAACCTGCCCTGCCTGGTTAATGCCAACAAAACGGTGCGGCTGggc   >  2:1305671/1‑150 (MQ=255)
                                                                                                                                               aTCTTTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAGTTTGTAGCAACGCAACCTGCCCTGCCTGGTTAATGCCAACAAAACGGTGCGGCTGggcg  >  1:2142000/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
TCATCAGACCAATAATTGAGCGATTGCCATCTTTGATTTGAGCAACCACGGATTCTGCCACCGCAGGCTGACGGCGATAATCTTTATTGGAATTACCGTGGCTGCAATCTACCATCAGAGACGGGCGCAGTCCCGCCTGTTCCATCTCTTTTTCACATTGCGCAACATCCGCAGGGCTATAGTTCGGCGCTTTACCACCGCGCAGGATCACATGGCCGTCCGGATTCCCCTGAGTTTGTAGCAACGCAACCTGCCCTGCCTGGTTAATGCCAACAAAACGGTGCGGCTGGGCG  >  CP009273/2733606‑2733898

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