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% 737.2 / ‑2.4 229D657Y (GAT→TAT) speAbiosynthetic arginine decarboxylase, PLP‑binding
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base A (121/106);  total (122/107)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01
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 229 total aligned reads displayed.

ATGGATGGACTCCTATTACAAAGTCGATATAAAAAAACCCTTCCGCGTCGTTAGGCCCGACGAGGAAGGGTTGGATTTGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCACTTGATCGCGGAACTGGGTTAACAGCGTTTTCGGATCGAGCTGTACATATTGCAGCATG  >  CP009273/3077155‑3077447
                                                                                                                                                   |                                                                                                                                                 
atggatggACTCCTATTACAAAGTCGATATAAAAAAACCCTTCCGCGTCGTTAGGCCCGACGAGGAAGGGTTGGATTTGTCACAATAAATAGTGTCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATAtt                                                                                                                                                 >  2:4725508/1‑150 (MQ=255)
atggatggACTCCTATTACAAAGTCGATATAAAAAAACCCATCCGCGTCGTTAGGCCCGACGAGGAAGGGTTGGATTTGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATAtt                                                                                                                                                 >  2:3185049/1‑150 (MQ=255)
  ggatggaCTCCTATTACAAAGTCGATATAAAAAAACCCTTCCGCGTCGTTAGGCCCGACGAGGAAGGGTTGGATTTGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCa                                                                                                                                               >  2:20401/1‑150 (MQ=255)
       gaCTCCTATTACAAAGTCGATATAAAAAAACCCTTCCGCGTCGTTAGGCCCGACGAGGAAGGGTTGGATTTGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATa                                                                                                                                          >  1:3217133/1‑150 (MQ=255)
       gaCTCCTATTACAAAGTCGATATAAAAAAACCCTTCCGCGTCGTTAGGCCCGACGAGGAAGGGTTGGATTTGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATa                                                                                                                                          >  1:3472254/1‑150 (MQ=255)
        aCTCCTATTACAAAGTCGATATAAAAAAACCCTTCCGCGTCGTTAGGCCCGACGAGGAAGGGTTGGATTTGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATaa                                                                                                                                         >  1:877115/1‑150 (MQ=255)
           ccTATTACAAAGTCGATATAAAAAAACCCTTCCGCGTCGTTAGGCCCGACGAGGAAGGGTTGGATTTGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGta                                                                                                                                      <  2:5007818/150‑1 (MQ=255)
            cTATTACAAAGTCGATATAAAAAAACCCTTCCGCGTCGTTAGGCCCGACGAGGAAGGGTTGGATTTGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGtat                                                                                                                                     <  1:5065057/150‑1 (MQ=255)
             tATTACAAAGTCGATATAAAAAAACCCTTCCGCGTCGTTAGGCCCGACGAGGAAGGGTTGGATTTGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGtata                                                                                                                                    <  2:437320/150‑1 (MQ=255)
             tATTACAAAGTCGATATAAAAAAACCCTTCCGCGTCGTTAGGCCCGACGAGGAAGGGTTGGATTTGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGtata                                                                                                                                    <  1:675157/150‑1 (MQ=255)
                tACAAAGTCGATATAAAAAAACCCTTCCGCGTCGTTAGGCCCGACGAGGAAGGGTTGGATTTGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAAcc                                                                                                                                 >  1:3939609/1‑150 (MQ=255)
                 aCAAAGTCGATATAAAAAAACCCTTCCGCGTCGTTAGGCCCGACGAGGAAGGGTTGGATTTGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCg                                                                                                                                >  1:2664182/1‑150 (MQ=255)
                  cAAAGTCGATATAAAAAAACCCTTCCGCGTCGTTAGGCCCGACGAGGAAGGGTTGGATTTGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGt                                                                                                                               <  2:2930652/150‑1 (MQ=255)
                  cAAAGTCGATATAAAAAAACCCTTCCGCGTCGTTAGGCCCGACGAGGAAGGGTTGGATTTGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGt                                                                                                                               >  1:2619168/1‑150 (MQ=255)
                    aaGTCGATATAAAAAAACCCTTCCGCGTCGTTAGGCCCGACGAGGAAGGGTTGGATTTGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTAc                                                                                                                             >  1:2490612/1‑150 (MQ=255)
                     aGTCGATATAAAAAAACCCTTCCGCGTCGTTAGGCCCGACGAGGAAGGGTTGGATTTGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACa                                                                                                                            <  2:1557763/150‑1 (MQ=255)
                     aGTCGATATAAAAAAACCCTTCCGCGTCGTTAGGCCCGACGAGGAAGGGTTGGATTTGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACa                                                                                                                            >  2:1645583/1‑150 (MQ=255)
                      gTCGATATAAAAAAACCCTTCCGCGTCGTTAGGCCCGACGAGGAAGGGTTGGATTTGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACaa                                                                                                                           <  2:2488134/150‑1 (MQ=255)
                         gATATAAAAAAACCCTTCCGCGTCGTTAGGCCCGACGAGGAAGGGTTGGATTTGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAAcc                                                                                                                        <  2:3491834/150‑1 (MQ=255)
                           tataAAAAAACCCTTCCGCGTCGTTAGGCCCGACGAGGAAGGGTTGGATTTGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTg                                                                                                                      <  1:3225209/150‑1 (MQ=255)
                           tataAAAAAACCCTTCCGCGTCGTTAGGCCCGACGAGGAAGGGTTGGATTTGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTg                                                                                                                      <  1:1952633/150‑1 (MQ=255)
                           tataAAAAAACCCTTCCGCGTCGTTAAGCCCGACGAGGAAGGGTTGGATTTGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTg                                                                                                                      <  1:639279/150‑1 (MQ=255)
                            ataAAAAAACCCTTCCGCGTCGTTAGGCCCGACGAGGACGGGTTGGATTTGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGc                                                                                                                     >  1:4111028/1‑150 (MQ=255)
                              aaaaaaaCCCTTCCGCGTCGTTAGGCCCGACGAGGAAGGGTTGGATTTGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCt                                                                                                                   >  1:1251059/1‑150 (MQ=255)
                              aaaaaaaCCCTTCCGCGTCGTTAGGCCCGACGAGGAAGGGTTGGATTTGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCt                                                                                                                   >  2:1857206/1‑150 (MQ=255)
                              aaaaaaaCCCTTCCGCGTCGTTAGGCCCGACGAGGAAGGGTTGGATTTGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCt                                                                                                                   >  1:3830639/1‑150 (MQ=255)
                                aaaaaCCCTTCCGCGTCGTTAGGCCCGACGAGGAAGGGTTGGATTTGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCg                                                                                                                 >  2:2123525/1‑150 (MQ=255)
                                 aaaaCCCTTCCGCGTCGTTAGGCCCGACGAGGAAGGGTTGGATTTGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGa                                                                                                                >  2:169809/1‑150 (MQ=255)
                                    aCCCTTCCGCGTCGTTAGGCCCGACGAGGAAGGGTTGGATTTGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACCCAGGACTTACTCATATTCAAGATAAGTATAACAGTACAAACCTGCCTCGaaat                                                                                                             >  2:3386391/1‑148 (MQ=255)
                                     cccTTCCGCGTCGTTAGGCCCGACGAGGAAGGGTTGGATTTGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAActc                                                                                                            <  1:3585459/150‑1 (MQ=255)
                                     cccTTCCGCGTCGTTAGGCCCGACGAGGAAGGGTTGGATTTGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAActc                                                                                                            >  2:3647049/1‑150 (MQ=255)
                                      ccTTCCGCGTCGTTAGGCCCGACGAGGAAGGGTTGGATTTGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAActct                                                                                                           >  2:3059489/1‑150 (MQ=255)
                                         tCCGCGTCGTTAGGCCCGACGAGGAAGGGTTGGATTTGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCa                                                                                                        >  1:3737402/1‑150 (MQ=255)
                                         tCCGCGTCGTTAGGCCCGACGAGGAAGGGTTGGATTTGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCa                                                                                                        >  2:2341447/1‑150 (MQ=255)
                                           cgcgTCGTTAGGCCCGACGAGGAAGGGTTGGATTTGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAg                                                                                                      <  2:2823890/150‑1 (MQ=255)
                                           cgcgTCGTTAGGCCCGACGAGGAAGGGTTGGATTTGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAg                                                                                                      >  2:2089872/1‑150 (MQ=255)
                                           cgcgTCGTTAGGCCCGACGAGGAAGGGTTGGATTTGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAg                                                                                                      >  1:5067381/1‑150 (MQ=255)
                                                   tAGGCCCGACGAGGAAGGGTTGGATTTGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACtgtt                                                                                              <  1:1090826/150‑1 (MQ=255)
                                                   tAGGCCCGACGAGGAAGGGTTGGATTTGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACtgtt                                                                                              <  1:4536001/150‑1 (MQ=255)
                                                    aGGCCCGACGAGGAAGGGTTGGATTTGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACtgttg                                                                                             >  2:4183843/1‑150 (MQ=255)
                                                     ggCCCGACGAGGAAGGGTTGGATTTGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACtgttgt                                                                                            >  2:2583600/1‑150 (MQ=255)
                                                         cgacgaGGAAGGGTTGGATTTGTCACAATAAATAGTGGCGGATTATCACCTCTAAATTAAGCGGATTCAAGTAACACAGGACTTAATAATATACAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGttgttgaa                                                                                        >  1:4637204/1‑148 (MQ=255)
                                                          gacgaGGAGGGGTTGGATTTGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAg                                                                                       >  2:587584/1‑150 (MQ=255)
                                                             gaGGAAGGGTTGGATTTGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTc                                                                                    >  1:961117/1‑150 (MQ=255)
                                                             gaGGAAGGGTTGGATTTGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTc                                                                                    >  1:689844/1‑150 (MQ=255)
                                                              aGGAAGGGTTGGATTTGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCa                                                                                   >  2:616835/1‑150 (MQ=255)
                                                               ggAAGGGTTGGATTTGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAg                                                                                  >  2:3419998/1‑150 (MQ=255)
                                                                      ttGGATTTGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAg                                                                           >  2:4210651/1‑150 (MQ=255)
                                                                       tGGATTTGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGa                                                                          >  2:3219982/1‑150 (MQ=255)
                                                                        ggATTTGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGAt                                                                         <  2:1085844/150‑1 (MQ=255)
                                                                           tttGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCgg                                                                      <  1:2204501/150‑1 (MQ=255)
                                                                           tttGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTAATCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCgg                                                                      <  1:1825744/150‑1 (MQ=255)
                                                                                cacaATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTttc                                                                 <  1:5111651/150‑1 (MQ=255)
                                                                                   aaTAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTttcttc                                                              >  2:4425263/1‑150 (MQ=255)
                                                                                   aaTAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTttcttc                                                              >  1:4887644/1‑150 (MQ=255)
                                                                                    aTAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCa                                                             >  2:1575241/1‑150 (MQ=255)
                                                                                        aTTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCACTTg                                                         <  1:1382157/150‑1 (MQ=255)
                                                                                           gtgGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCACTTGATc                                                      <  1:3055376/150‑1 (MQ=255)
                                                                                             ggCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCACTTGATcgc                                                    >  2:3207412/1‑150 (MQ=255)
                                                                                              gCGGCTGATCACCGCTAAATTAAGCGGATTCCAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCACTTGATcgcg                                                   <  1:921721/150‑1 (MQ=255)
                                                                                              gCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCACTTGATcgcg                                                   <  1:1915124/150‑1 (MQ=255)
                                                                                                 gATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCACTTGATCGCGGaa                                                <  1:471550/150‑1 (MQ=255)
                                                                                                   ttATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCACTTGATCGCGGAACt                                              >  1:4340435/1‑150 (MQ=255)
                                                                                                     aTCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCACTTGATCGCGGAACTgg                                            >  1:4994752/1‑150 (MQ=255)
                                                                                                     aTCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCACTTGATCGCGGAACTgg                                            >  2:5128035/1‑150 (MQ=255)
                                                                                                            cTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCACTTGATCGCGGAACTGGGTTAACa                                     <  1:2089872/150‑1 (MQ=255)
                                                                                                               aaTTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCGCTTGATCGCGGAACTGGGTTAACAGCg                                  >  2:3811797/1‑150 (MQ=255)
                                                                                                               aaTTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCACTTGATCGCGGAACTGGGTTAACAGCg                                  >  1:3429028/1‑150 (MQ=255)
                                                                                                               aaTTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTAATCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTGTCTTCACTTGATCTCGGAACTGGGTTAACAGCg                                  >  2:5154562/1‑150 (MQ=255)
                                                                                                                aTTAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCACTTGATCGCGGAACTGGGTTAACAGCGt                                 >  1:440593/1‑150 (MQ=255)
                                                                                                                  tAAGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCACTTGATCGCGGAACTGGGTTAACAGCGttt                               >  1:100725/1‑150 (MQ=255)
                                                                                                                   aaGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCACTTGATCGCGGAACTGGGTTAACAGCGtttt                              <  2:2635725/150‑1 (MQ=255)
                                                                                                                   aaGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCACTTGATCGCGGAACTGGGTTAACAGCGtttt                              <  1:2881462/150‑1 (MQ=255)
                                                                                                                    aGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCACTTGATCGCGGAACTGGGTTAACAGCGTTTTc                             <  1:2257237/150‑1 (MQ=255)
                                                                                                                    aGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCACTTGATCGCGGAACTGGGTTAACAGCGTTTTc                             <  1:811555/150‑1 (MQ=255)
                                                                                                                    aGCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCACTTGATCGCGGAACTGGGTTAACAGCGTTTTc                             <  1:3185049/150‑1 (MQ=255)
                                                                                                                     gCGGATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCACTTGATCGCGGACCTGGGTTAACAGCGTTTTCg                            >  2:2522776/1‑150 (MQ=255)
                                                                                                                       ggATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCACTTGATCGCGGAACTGGGTTAACAGCGTTTTCGGa                          >  2:1173316/1‑150 (MQ=255)
                                                                                                                        gATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCACTTGATCGCGGAACTGGGTTAACAGCGTTTTCGGAt                         >  2:274521/1‑150 (MQ=255)
                                                                                                                        gATTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCACTTGATCGCGGAACTGGGTTAACAGCGTTTTCGGAt                         >  2:2860292/1‑150 (MQ=255)
                                                                                                                         aTTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCACTTGATCGCGGAACTGGGTTAACAGCGTTTTCGGATc                        <  2:4416575/150‑1 (MQ=255)
                                                                                                                         aTTCAAGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCACTTGATCGCGGAACTGGGTTAACAGCGTTTTCGGATc                        <  1:4007195/150‑1 (MQ=255)
                                                                                                                            cAAGGAACACAGGACTTACTCATATTCAAGCTAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCACTTGATCGCGGAACTGGGTTAACAGCGTTTTCGGATCGAg                     <  2:4972244/150‑1 (MQ=255)
                                                                                                                             aaGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCACTTGATCGCGGAACTGGGTTAACAGCGTTTTCGGATCGAGc                    <  2:1251059/150‑1 (MQ=255)
                                                                                                                              aGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCACTTGATCGCGGAACTGGGTTAACAGCGTTTTCGGATCGAGCt                   >  1:3537702/1‑150 (MQ=255)
                                                                                                                              aGTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCACTTGATCGCGGAACTGGGTTAACAGCGTTTTCGGATCGAGCt                   <  1:3688568/150‑1 (MQ=255)
                                                                                                                               gTAACACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCACTTGATCGCGGAACTGGGTTAACAGCGTTTTCGGATCGAGCTg                  <  2:186362/150‑1 (MQ=255)
                                                                                                                                 aaCACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCACTTGATCGCGGAACTGGGTTAACAGCGTTTTCGGATCGAGCTGTa                >  2:1052906/1‑150 (MQ=255)
                                                                                                                                 aaCACAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCAGCCTCGAACTCTACAAGGAACTGTTGATGCCGTTCAGCAACAAGATCGTTTTTCTTCACTTGATCCCGGAACGGGGTTAACAGCGTGTTCTGATCGAGCTGTa                >  2:864589/1‑150 (MQ=255)
                                                                                                                                  acacAGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCACTTGATCGCGGAACTGGGTTAACAGCGTTTTCGGATCGAGCTGTAc               >  2:1404024/1‑150 (MQ=255)
                                                                                                                                   cacaGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCACTTGATCGCGGAACTGGGTTAACAGCGTTTTCGGATCGAGCTGTACa              >  1:3255369/1‑150 (MQ=255)
                                                                                                                                     caGGACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCACTTGATCGCGGAACTGGGTTAACAGCGTTTTCGGATCGAGCTGTACata            >  2:4005797/1‑150 (MQ=255)
                                                                                                                                       ggACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCACTTGATCGCGGAACTGGGTTAACAGCGTTTTCGGATCGAGCTGTACATAtt          >  2:765585/1‑150 (MQ=255)
                                                                                                                                        gACTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCACTTGATCGCGGAACTGGGTTAACAGCGTTTTCGGATCGAGCTGTACATATTg         <  1:2379516/150‑1 (MQ=255)
                                                                                                                                         aCTTACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCACTTGATCGCGGAACTGGGTTAACAGCGTTTTCGGATCGAGCTGTACATATTgc        >  1:4191576/1‑150 (MQ=255)
                                                                                                                                            tACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCACTTGATCGCGGAACTGGGTTAACAGCGTTTTCGGATCGAGCTGTACATATTgcagc     <  2:4111028/150‑1 (MQ=255)
                                                                                                                                            tACTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCACTTGATCGCGGAACTGGGTTAACAGCGTTTTCGGATCGAGCTGTACATATTgcagc     <  2:1311004/150‑1 (MQ=255)
                                                                                                                                             aCTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCACTTGATCGCGGAACTGGGTTAACAGCGTTTTCGGATCGAGCTGTACATATTgcagca    <  2:4862162/150‑1 (MQ=255)
                                                                                                                                              cTCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCACTTGATCGCGGAACTGGGTTAACAGCGTTTTCGGATCGAGCTGTACATATTGCAGCAt   <  2:3859434/150‑1 (MQ=255)
                                                                                                                                               tCATATTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCACTTGATCGCGGAACTGGGTTAACAGCGTTTTCGGATCGAGCTGTACATATTGCAGCATg  <  1:2858965/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
ATGGATGGACTCCTATTACAAAGTCGATATAAAAAAACCCTTCCGCGTCGTTAGGCCCGACGAGGAAGGGTTGGATTTGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATCGGTTTTCTTCACTTGATCGCGGAACTGGGTTAACAGCGTTTTCGGATCGAGCTGTACATATTGCAGCATG  >  CP009273/3077155‑3077447

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: