Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,077,302 C→A 100% D657Y (GAT→TAT)  speA ← biosynthetic arginine decarboxylase, PLP‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,077,3020CA99.1% 1082.1 / ‑1.3 337D657Y (GAT→TAT) speAbiosynthetic arginine decarboxylase, PLP‑binding
Reads supporting (aligned to +/- strand):  ref base C (4/0);  new base A (150/183);  total (154/183)
Fisher's exact test for biased strand distribution p-value = 4.27e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 337 total aligned reads displayed.

TGGATGGACTCCTATTACAAAGTCGATATAAAAAAACCCTTCCGCGTCGTTAGGCCCGACGAGGAAGGGTTGGATTTGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCACTTGATCGCGGAACTGGGTTAACAGCGTTTTCGGATCGAGCTGTACATATTGCAGCATGTC  >  CP009273/3077156‑3077449
                                                                                                                                                  |                                                                                                                                                   
 ggatggaCTCCTATTACAAAGTCGATATAAAAAAACCCTTCCGCGTCGTTAGGCCCGACGAGGAAGGGTTGGATTTGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCa                                                                                                                                                 >  2:1246358/1‑150 (MQ=255)
   atggaCTCCTATTACAAAGTCGATATAAAAAAACCCTTCCGCGTCGTTAGGCCCGACGAGGAAGGGTTGGATTTGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAg                                                                                                                                               >  1:4941370/1‑150 (MQ=255)
   atggaCTCCTATTACAAAGTCGAAATAAAAAAAACCTACCGCGTCGATACGCCCGACGAGGAAGGCTAGGATTTTTAACAAAAAATATTGGCGGAATATAACCGCTAAATTAAACCGAATCAAGTAACACACGACTTACTAATATTCAAg                                                                                                                                               >  1:1847212/1‑150 (MQ=255)
    tggaCTCCTATTACAAAGTCGATATAAAAAAACCCTTCCGCGTCGTTAGGCCCGACGAGGAAGGGTTGGATTTGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGa                                                                                                                                              >  1:2771869/1‑150 (MQ=255)
     ggaCTCCTATTACAAAGTCGATATAAAAAAACCCGTCCGCGTCGTTAGGCCCGACGAGGAAGGGTTGGATTTGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGAt                                                                                                                                             <  1:1504717/150‑1 (MQ=255)
        cTCCTATTACAAAGTCGATATAAAAAAACCCTTCCGCGTCGTTAGGCCCGACGAGGAAGGGTTGGATTTGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAg                                                                                                                                          >  2:200204/1‑150 (MQ=255)
         tCCTATTACAAAGTCGATATAAAAAAACCCTTCCGCGTCGTTAGGCCCGACGAGGAAGGGTTGGATTTGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGt                                                                                                                                         >  2:1356552/1‑150 (MQ=255)
           cTATTACAAAGTCGATATAAAAAAACCCTTCCGCGTCGTTAGGCCCGACGAGGAAGGGTTGGATTTGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGtat                                                                                                                                       >  2:4677854/1‑150 (MQ=255)
               tACAAAGTCGATATAAAAAAACCCTTCCGCGTCGTTAGGCCCGACGAGGAAGGGTTGGATTTGTCACTATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAAcc                                                                                                                                   <  1:3420362/150‑1 (MQ=255)
                aCAAAGTCGATATAAAAAAACCCTTCCGCGTCGTTAGGCCCGACGAGGAAGGGTTGGATTTGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCg                                                                                                                                  >  1:576347/1‑150 (MQ=255)
                aCAAAGTCGATATAAAAAAACCCTTCCGCGTCGTTAGGCCCGACGAGGAAGGGTTGGATTTGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCg                                                                                                                                  <  1:3289317/150‑1 (MQ=255)
                  aaaGTCGATATAAAAAAACCCTTCCGCGTCGTTAGGCCCGACGAGGAAGGGTTGGATTTGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTa                                                                                                                                <  1:4721335/150‑1 (MQ=255)
                   aaGTCGATATAAAAAAACCCTTCCGCGTCGTTAGGCCCGACGAGGAAGGGTTGGATTTGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTAc                                                                                                                               <  2:4755501/150‑1 (MQ=255)
                    aGTCGATATAAAAAAACCCTTCCGCGTCGTTAGGCCCGACGAGGAAGGGTTGGATTTGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACa                                                                                                                              >  2:2528254/1‑150 (MQ=255)
                       cGATATAAAAAAACCCTTCCGCGTCGTTAGGCCCGACGAGGAAGGGTTGGATTTGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAAc                                                                                                                           >  1:1271258/1‑150 (MQ=255)
                        gATATAAAAAAACCCTTCCGCGTCGTTAGGCCCGACGAGGAAGGGTTGGATTTGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAGAcc                                                                                                                          >  1:203377/1‑150 (MQ=255)
                         atatAAAAAAACCCTTCCGCGTCGTTAGGCCCGACGAGGAAGGGTTGGATTTGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCt                                                                                                                         <  1:964849/150‑1 (MQ=255)
                           ataAAAAAACCCTTCCGCGTCGTTAGGCCCGACGAGGAAGGGTTGGATTTGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGc                                                                                                                       >  2:1390190/1‑150 (MQ=255)
                             aaaaaaaCCCTTCCGCGTCGTTAGGCCCGACGAGGAAGGGTTGGATTTGTCACAATAAATTGTGGCGGATTATCACCGCGAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCt                                                                                                                     >  2:321149/1‑150 (MQ=255)
                              aaaaaaCCCTTCCGCGTCGTTAGGCCCGACGAGGAAGGGTTGGATTTGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTc                                                                                                                    >  1:4336441/1‑150 (MQ=255)
                              aaaaaaCCCTTCCGCGTCGTTAGGCCCGACGAGGAAGGGTTGGATTTGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTc                                                                                                                    <  2:2964991/150‑1 (MQ=255)
                              aaaaaaCCCTTCCGCGTCGTTAGGCCCGACGAGGAAGGGTTGGATTTGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTc                                                                                                                    <  2:280801/150‑1 (MQ=255)
                              aaaaaaCCCGTCCGCGTCGTTAGGCCCGACGAGGAAGGGTTGGATTTGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTc                                                                                                                    <  2:3366741/150‑1 (MQ=255)
                                aaaaCCCTTCCGCGTCGTTAGGCCCGACGAGGAAGGGTTGGATTTGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGa                                                                                                                  >  2:284045/1‑150 (MQ=255)
                                   aCCCTTCCGCGTCGTTAGGCCCGACGAGGAAGGGTTGGATTTGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAAct                                                                                                               >  2:3969759/1‑150 (MQ=255)
                                       ttCCGCGTCGTTAGGCCCGACGAGGAAGGGTTGGATTTGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTc                                                                                                           <  1:2795512/150‑1 (MQ=255)
                                        tCCGCGTCGTTAGGCCCGACGAGGAAGGGTTGGATTTGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCa                                                                                                          <  2:2481864/150‑1 (MQ=255)
                                            cgtcgtTAGGCCCGACGAGGAAGGGTTGGATTTGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCAGCCTCGAACTCTTCAAGGa                                                                                                      >  2:4137263/1‑150 (MQ=255)
                                              tcgtTAGGCCCGACGAGGAAGGGTTGGATTTGTCACAATAAATTGTGGCGGATTATCAACGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAAACTGCCTCGAACTCTTCAAGGAAc                                                                                                    <  2:643683/150‑1 (MQ=255)
                                                gtTAGGCCCGACGAGGAAGGGTTGGATTTGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACtg                                                                                                  >  2:3435145/1‑150 (MQ=255)
                                                  tAGGCCCGACGAGGAAGGGTTGGATTTGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACtgtt                                                                                                >  1:4508187/1‑150 (MQ=255)
                                                  tAGGCCCGACGAGGAAGGGTTGGATTTGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACtgtt                                                                                                <  2:5270623/150‑1 (MQ=255)
                                                   aGGCCCGACGAGGAAGGGTTGGATTTGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACtgttg                                                                                               >  1:3301161/1‑150 (MQ=255)
                                                   aGGCCCGACGAGGAAGGGTTGGATTTGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACtgttg                                                                                               >  1:1018626/1‑150 (MQ=255)
                                                        cgacgaGGAAGGGTTGGATTTGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCa                                                                                          <  2:1672760/150‑1 (MQ=255)
                                                            gaGGAAGGGTTGGATTTGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTc                                                                                      >  1:4990566/1‑150 (MQ=255)
                                                             aGGAAGGGTTGGATTTGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCa                                                                                     >  1:1783001/1‑150 (MQ=255)
                                                              ggAAGGGTTGGATTTGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAg                                                                                    >  1:1295484/1‑150 (MQ=255)
                                                                 aGGGTTGGATTTGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCAt                                                                                 >  2:1913751/1‑150 (MQ=255)
                                                                  gggTTGGATTAGTCACAATAAATTGTGGCGGATTATCACCGCGAAAGTAATCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATc                                                                                >  2:1120485/1‑150 (MQ=255)
                                                                       ggATTTGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGAt                                                                           >  2:4459826/1‑150 (MQ=255)
                                                                        gATTTGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATc                                                                          <  2:2350085/150‑1 (MQ=255)
                                                                         aTTTGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCg                                                                         <  2:226147/150‑1 (MQ=255)
                                                                         aTTTGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCg                                                                         >  2:747164/1‑150 (MQ=255)
                                                                         aTTTGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCg                                                                         <  2:1073108/150‑1 (MQ=255)
                                                                         aTTTGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCg                                                                         >  2:1608505/1‑150 (MQ=255)
                                                                                   aTAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCa                                                               >  1:2319215/1‑150 (MQ=255)
                                                                                     aaaTTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCACt                                                             <  2:2739871/150‑1 (MQ=255)
                                                                                      aaTTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCACtt                                                            <  1:840350/150‑1 (MQ=255)
                                                                                          gtgGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCACTTGATc                                                        <  1:69248/150‑1 (MQ=255)
                                                                                          gtgGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAAATCGGTTTTCTTCACTTGATc                                                        <  2:5071101/150‑1 (MQ=255)
                                                                                          gtgGCGGATTATCACCGCTAAAGTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCACTTGATc                                                        <  2:3453156/150‑1 (MQ=255)
                                                                                            ggCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCACTTGATcgc                                                      <  1:3741910/150‑1 (MQ=255)
                                                                                             gCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCACTTGATcgcg                                                     >  2:4214072/1‑150 (MQ=255)
                                                                                                gATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCACTTGATCGCGGaa                                                  <  2:2171052/150‑1 (MQ=255)
                                                                                                gATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCACTTGATCGCGGaa                                                  <  2:574482/150‑1 (MQ=255)
                                                                                                 aTTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCACTTGATCGCGGAAc                                                 >  1:4916303/1‑150 (MQ=255)
                                                                                                 aTTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCACTTGATCGCGGAAc                                                 >  2:5178432/1‑150 (MQ=255)
                                                                                                  ttATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCACTTGATCGCGGAACt                                                <  1:2612815/150‑1 (MQ=255)
                                                                                                   tATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCACTTGATCGCGGAACTg                                               <  1:747164/150‑1 (MQ=255)
                                                                                                   tATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCACTTGATCGCGGAACTg                                               <  1:5287961/150‑1 (MQ=255)
                                                                                                   tATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCACTTGATCGCGGAACTg                                               >  1:5135769/1‑150 (MQ=255)
                                                                                                     tCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCACTTGATCGCGGAACTggg                                             >  2:4078493/1‑150 (MQ=255)
                                                                                                        ccGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCACTTGATCGCGGAACTGGGTTa                                          <  2:1446366/150‑1 (MQ=255)
                                                                                                         cgcgAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGATGCAGTTCAGCATCAAGATCGGAATTCTTCACTTGATCGCGGAACTGGGTTaa                                         >  2:1550536/1‑150 (MQ=255)
                                                                                                         cGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCACTTGATCGCGGAACTGGGTTaa                                         <  2:1658558/150‑1 (MQ=255)
                                                                                                             aaaTTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCACTTGATCGCGGAACTGGGTTAACAGc                                     >  2:4334389/1‑150 (MQ=255)
                                                                                                              aaTTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCACTTGATCGCGGAACTGGGTTAACAGCg                                    <  1:986233/150‑1 (MQ=255)
                                                                                                              aaTTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAAACGTACAAACCTGCCTCGAACTCTTCAAGGAACAGGGGAGGCAGTTCAGCATCAAGATCGGTGTACTTCACTTGATAGCGTAACTGTGATAACACCg                                    >  1:2418307/1‑150 (MQ=255)
                                                                                                                 tAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCACTTGATCGCGGAACTGGGTTAACAGCGttt                                 <  1:2290693/150‑1 (MQ=255)
                                                                                                                  aaGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCACTTGATCGCGGACCTGGGTTAACAGCGtttt                                >  2:3545093/1‑150 (MQ=255)
                                                                                                                   aGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCACTTGATCGCGGAACTGGGTTAACAGCGTTTTc                               <  2:1755050/150‑1 (MQ=255)
                                                                                                                   aGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCACTTGATCGCGGAACTGGGTTAACAGCGTTTTc                               <  1:615896/150‑1 (MQ=255)
                                                                                                                   aGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCACTTGATCGCGGAACTGGGTTAACAGCGTTTTc                               <  1:398604/150‑1 (MQ=255)
                                                                                                                   aGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCACTTGATCGCGGAACTGGGTTAACAGCGTTTTc                               >  1:1056623/1‑150 (MQ=255)
                                                                                                                     cGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCACTTGATCGCGGAACTGGGTTAACAGCGTTTTCgg                             >  2:4301155/1‑150 (MQ=255)
                                                                                                                     cGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCACTTGATCGCGGAACTGGGTTAACAGCGTTTTCgg                             >  1:2145964/1‑150 (MQ=255)
                                                                                                                       gATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCACTTGATCGCGGAACTGGGTTAACAGCGTTTTCGGAt                           >  2:955095/1‑150 (MQ=255)
                                                                                                                       gATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCACTTGATCGCGGAACTGGGTTAACAGCGTTTTCGGAt                           <  2:1301859/150‑1 (MQ=255)
                                                                                                                        aTTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCACTTGATCGCGGAACTGGGTTAACAGCGTTTTCGGATc                          <  1:1623432/150‑1 (MQ=255)
                                                                                                                        aTTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCACTTGATCGCGGAACTGGGTTAACAGCGTTTTCGGATc                          <  2:1269655/150‑1 (MQ=255)
                                                                                                                         ttCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCACTTGATCGCGGAACTGGGTTAACAGCGTTTTCGGATCg                         <  2:3300701/150‑1 (MQ=255)
                                                                                                                          tCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCACTTGATCGCGGAACTGGGTTAACAGCGTTTTCGGATCGa                        <  1:4459826/150‑1 (MQ=255)
                                                                                                                           cAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCACTTGATCGCGGAACTGGGTTAACAGCGTTTTCGGATCGAg                       <  1:2857776/150‑1 (MQ=255)
                                                                                                                            aaGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCACTTGATCGCGGAACTGGGTTAACAGCGTTTTCGGATCGAGc                      <  1:1083597/150‑1 (MQ=255)
                                                                                                                               tAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCACTTGATCGCGGAACTGGGTTAACAGCGTTTTCGGATCGAGCTGt                   <  1:1245794/150‑1 (MQ=255)
                                                                                                                                aaCACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCACTTGATCGCGGAACTGGGTTAACAGCGTTTTCGGATCGAGCTGTa                  <  2:3728353/150‑1 (MQ=255)
                                                                                                                                  cacaGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCACTTGATCGCGGAACTGGGTTAACAGCGTTTTCGGATCGAGCTGTACa                <  2:4243484/150‑1 (MQ=255)
                                                                                                                                     aGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCACTTGATCGCGGAACTGGGTTAACAGCGTTTTCGGATCGAGCTGTACatat             <  1:4151948/150‑1 (MQ=255)
                                                                                                                                      ggACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCACTTGATCGCGGAACTGGGTTAACAGCGTTTTCGGATCGAGCTGTACATAtt            >  2:4903362/1‑150 (MQ=255)
                                                                                                                                         cTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCACTTGATCGCGGAACTGGGTTAACAGCGTTTTCGGATCGAGCTGTACATATTgca         <  1:3809126/150‑1 (MQ=255)
                                                                                                                                           tACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCACTTGATCGCGGAACTGGGTTAACAGCGTTTTCGGATCGAGCTGTACATATTgcagc       <  2:4582333/150‑1 (MQ=255)
                                                                                                                                           tACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCACTTGATCGCGGAACTGGGTTAACAGCGTTTTCGGATCGAGCTGTACATATTgcagc       <  1:3563964/150‑1 (MQ=255)
                                                                                                                                           tACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCACTTGATCGCGGAACTGGGTTAACAGCGTTTTCGGATCGAGCTGTACATATTgcagc       <  1:1445477/150‑1 (MQ=255)
                                                                                                                                            aCTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCACTTGATCGCGGAAATGGGTTAACAGCGTTTTCTGATCGAGCTGTACATATTgcagca      >  2:4982203/1‑150 (MQ=255)
                                                                                                                                             cTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCACTTGATCGCGGAACTGGGTTAACAGCGTTTTCGGATCGAGCTGTACATATTGCAGCAt     >  1:3084210/1‑150 (MQ=255)
                                                                                                                                             cTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCACTTGATCGCGGAACTGGGTTAACAGCGTTTTCGGATCGAGCTGTACATATTGCAGCAt     <  1:2161694/150‑1 (MQ=255)
                                                                                                                                              tCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCACTTGATCGCGGAACTGGGTTAACAGCGTTTTCGGATCGAGCTGTACATATTGCAGCATg    <  2:812716/150‑1 (MQ=255)
                                                                                                                                               cATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCACTTGATCGCGGAACTGGGTTAACAGCGTTTTCGGATCGAGCTGTACATATTGCAGCATGt   >  1:3958973/1‑150 (MQ=255)
                                                                                                                                                aTCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCACTTGATCGCGGAACTGGGTTAACAGCGTTTTCGGATCGAGCTGTACATATTGCAGCATGTc  >  2:208948/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
TGGATGGACTCCTATTACAAAGTCGATATAAAAAAACCCTTCCGCGTCGTTAGGCCCGACGAGGAAGGGTTGGATTTGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCACTTGATCGCGGAACTGGGTTAACAGCGTTTTCGGATCGAGCTGTACATATTGCAGCATGTC  >  CP009273/3077156‑3077449

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: