Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009273 2,579,392 T→A H101Q (CAT→CAA)  narQ → sensory histidine kinase in two‑component regulatory system with NarP (NarL)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,579,3920TA98.9% 292.9 / NA 93H101Q (CAT→CAA) narQsensory histidine kinase in two‑component regulatory system with NarP (NarL)
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base A (45/47);  minor base C (1/0);  total (46/47)
Fisher's exact test for biased strand distribution p-value = 4.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGTTACCGCCTGGGCTACGACTTGCAAAGTGGCAGTCCACAACTCAATGCACATCGCCAGTTATTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGTATGTGCCAGAAGCAGTAAAAACTCGCTAT‑GCGCATCTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTGCCGTGGTATCAGGCCAATATTAATAATTATGTTAATCAGATAGACCTGTTCGTACTGGCTTTACAGCACTACGCTGAACGCAAAATGCTGCTGG  >  CP009273/2579246‑2579537
                                                                                                                                                   |                                                                                                                                                 
agTTACCGCCTGGGCTACGACTTGCAAAGTGGCAGTCCACAACTCAATGCACATCGCCAGTTATTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCCAGGTCCGGCATGTGCCAGAAGCAGTAAAAACTCGCTAT‑GCGCAACt                                                                                                                                                 <  1:1061280/149‑1 (MQ=255)
    aCCGCCTGGGCTACGACTTGCAAAGTGGCAGTCCACAACTCAATGCACATCGCCAGTTATTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGTATGTGCCAGAAGCAGTAAAAACTCGCTAT‑GCGCAACTGAAt                                                                                                                                             <  1:649807/149‑1 (MQ=255)
      cGCCTGGGCTACGACTTGCAAAGTGGCAGTCCACAACTCAATGCACATCGCCAGTTATTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGTATGTGCCAGAAGCAGTAAAAACTCGCTAT‑GCGCAACTGAATGc                                                                                                                                           >  1:512568/1‑149 (MQ=255)
        ccagggCTACGACTTGCAAAGGGGCAGTGCACAACTCAATGCACATCGCCAGTTATTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGTATGTGCCAGAAGCAGTAAAAACTCGCTAT‑GCGCAACTGAATGCCa                                                                                                                                         <  1:27325/146‑1 (MQ=255)
         cTGGGCTACGACTTGCAAAGTGGCAGTCCACAACTCAATGCACATCGCCAGTTATTTCAGCAGGCACTGCATTCCCCGGTATTAACCAACCTCAACGTCTGGTATGTGCCAGAAGCAGTAAAAACTCGCTAT‑GCGCAACTGAATGCCa                                                                                                                                         >  2:1007279/1‑148 (MQ=255)
         cTGGGCTACGACTTGCAAAGTGGCAGTCCACAACTCAATGCACATCGCCAGTTATTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGTATGTGCCAGAAGCAGTAAAAACTCGCTAT‑GCGCAACTGAATGCCaa                                                                                                                                        <  1:248477/149‑1 (MQ=255)
              cTACGACTTGCAAAGTGGCAGTCCACAACTCAATGCACATCGCCAGTTATTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGTATGTGCCAGAAGCAGTAAAAACTCGCTAT‑GCGCAACTGAATGCCAActggc                                                                                                                                   <  1:972958/149‑1 (MQ=255)
                 cGACTTGCAAAGTGGCAGTCCACAACTCAATGCACATCGCCAGTTATTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGTATGTGCCAGAAGCAGTAAAAACTCGCTAT‑GCGCAACTGAATGCCAActggctg                                                                                                                                 >  2:550129/1‑148 (MQ=255)
                 cGACTTGCAAAGTGGCAGTCCACAACTCAATGCACATCGCCAGTTATTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGTATGTGCCAGAAGCAGTAAAAACTCGCTAT‑GCGCAACTGAATGCCAActggctg                                                                                                                                 <  2:853356/148‑1 (MQ=255)
                 cGACTTGCAAAGCGGCAGGTCACAACTCAATGCACATCGCCAGTTATTTCAGCAGGCACTGCATTCACCGGTATTCACCAACCTCAACGTCTGGTATGTGCCAGAAGCAGTAAAAACTCGCTAT‑GCGCAACTGAATGCCAActggctg                                                                                                                                 <  2:281302/148‑1 (MQ=255)
                  gACTTGCAAAGTGGCAGTGCCCAACTCAATGCACATCGCCAGTTATTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGTATGTGCCAGAAGCAGTAAAAACTCGCTAT‑GCGCAACTGAATGCCAACTGGCTGga                                                                                                                               <  1:876416/149‑1 (MQ=255)
                  gACTTGCAAAGTGGCAGTCCACAACTCAATGCACATCGCCAGTTATTTAAGCACGCACTGCAAGCACGGGCATTAACCAACCCCAACGTCTGGCATGTGCCAGAAGCAGTAAAAACTCGCTAT‑GCGCAACTGAATGCCAAcaggccga                                                                                                                                >  2:268000/1‑145 (MQ=255)
                     ttGCAAAGTGGCAGTCCACAACTCAATGCACATCGCCAGTTATTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGTATGTGCCAGAAGCAGTAAAAACTCGCTAT‑GCGCAACTGAATGCCAACTGGCTGGAGAt                                                                                                                            >  1:720025/1‑149 (MQ=255)
                     ttGCAAAGTGGCAGTCCACAACTCAATGCACATCGCCAGTTATTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGTATGTGCCAGAAGCAGTAAAAACTCGCTAT‑GCGCAACTGAATGCCAACTGGCTGGAGAt                                                                                                                            >  1:1338416/1‑149 (MQ=255)
                     ttGCAAAGTGGCAGTCCACAACTCAATGCACATCGCCAGTTATTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGTATGTGCCAGAAGCAGTAAAAACTCGCTAT‑GCGCAACTGAATGCCAACTGGCTGGAGAt                                                                                                                            >  1:290427/1‑149 (MQ=255)
                       gCAAAGTGGCAGTCCACAACTCAATGCACATCGCCAGTTATTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGTATGTGCCAGAAGCAGTAAAAACTCGCTAT‑GCGCAACTGAATGCCAACTGGCTGGAGATg                                                                                                                           >  2:200750/1‑148 (MQ=255)
                        cAAAGCGGCAGTCCACAACGCAATGCACATCGCCAGTTATTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGTATGTGCCAGAAGCAGTAAAAACTCGCTAT‑GCGCAACCGAATGCCAACTGGCTGGAGATGa                                                                                                                          <  2:848516/148‑1 (MQ=255)
                                   tCCACAACTCAATGCACATCGCCAGTTATTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGTATGTGCCAGAAGCAGTAAAAACTCGCTAT‑GCGCAACTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCa                                                                                                              >  1:979362/1‑149 (MQ=255)
                                   tCCACAACCCAATGCACATCGCCAGTTATTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGTATGTGCCAGAAGCAGTAAAAACTCGCTAT‑GCGCAACTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCa                                                                                                              <  1:642489/149‑1 (MQ=255)
                                     cacaACTCAATGCACATCGCCAGTTATTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGTATGTGCCAGAAGCAGTAAAAACTCGCTAT‑GCGCAACTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGc                                                                                                            >  1:385333/1‑149 (MQ=255)
                                     cacaACTCAATGCACATCGCCAGTTATTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGTATGTGCCAGAAGCAGTAAAAACTCGCTAT‑GCGCAACTGAATGCCAACTGGCTGGAGATGAAAAACCGGCTCAg                                                                                                             >  2:579535/1‑148 (MQ=255)
                                     cacaACTCAATGCACATCGCCAGTTATTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGTATGTGCCAGAAGCAGTAAAAACTCGCTAT‑GCGCAACGGAATGCGAACTGGCTGGAGAGGAATAATCGGCTCAg                                                                                                             >  2:364803/1‑148 (MQ=255)
                                      acacCTCAATGCACAGCGCCAGTTATTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGTATGTGCCAGAAGCAGTAAAAACTCGCTAT‑GCGCAACTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGc                                                                                                            <  2:195639/148‑1 (MQ=255)
                                       caACTCAATGCACATCGCCAGTTATTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGTATGTGCCAGAAGCAGTAAAAACTCGCTAT‑GCGCAACTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCa                                                                                                           >  2:369307/1‑148 (MQ=255)
                                        aaCTCAATGCACATCGCCAGTTATTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGTATGTGCCAGAAGCAGT‑AAACTTCTGTGTCGCGCAACTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAg                                                                                                         >  1:935949/1‑149 (MQ=255)
                                          cgcAACGCACCTCGCCAGTCATTTACGTAGGCACTCCATTACCCGGTATGAACCAACCTCAACGTCTGGTATGTGCCAGAAGCAGTAAAAACTCGCTAT‑GCGCAACTGGATGCCAACTGACTTGAGATGAATAATCGTCTCAGCAAgg                                                                                                        <  2:274624/146‑1 (MQ=255)
                                          cTCAATGCACATCGCCAGTTATTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGTATGTGCCAGAAGCAGTAAAAGCTCGCTAT‑GCGCAACTGAATGCCAACTGGCTGGAGATGAATAACCGGCTCAGCAAgg                                                                                                        >  2:133840/1‑148 (MQ=255)
                                           tgaaTGCACATCGCCAGGTATTTCAGCAGGCACTGCATTCACCGGTATTATCCAACCTCACCGTCTGGTATGGGCCAGAAGCAGTAAATACGCGCTAT‑GCGCAACTGAATGCCACCTGCCTGGAGATGAATAATCGGCTCTGCAAggg                                                                                                       <  2:23582/146‑1 (MQ=255)
                                           tCAATGCACATCGCCAGTTATTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGTATGTGCCAGAAGCAGTAAAAACTCGCTAT‑GCGCAACTGAATGCCAACTGGCTGGCGATGAATAATCGGCTCAGCAAGGGc                                                                                                      >  1:244200/1‑149 (MQ=255)
                                           tCAATGCACATCGCCAGTTATTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGTATGTGCCAGAAGCAGTAAAAACTCGCTAT‑GCGCAACTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAggg                                                                                                       >  2:69494/1‑148 (MQ=255)
                                                  acaTCGCCAGTTATTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGTATGTGCCAGAAGCAGTAAAAACTCGCTAT‑GCGCAACTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTg                                                                                                <  2:532004/148‑1 (MQ=255)
                                                     tggCCAGTTATTTCACCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGTATGTGCCAGAAGCAGGAAAAACTCGCTAC‑GCACAACGTAATGCCAACGGGCTGTAGATGAATAATCGGCTCAGCAAGGGCGATTTGCCg                                                                                             <  2:554009/146‑1 (MQ=255)
                                                     tCGCCAGTTATTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGTATGTGCCAGAAGCAGTAAAAACTCGCTAT‑GCGCAACTGAATGCCAACTGGCTGGAGGTGAATAATCGGCTCAGCCAGGGCGCTTTGcgg                                                                                             >  2:606881/1‑146 (MQ=255)
                                                       gCCAGTTATTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGTATGTGCCAGAAGCAGTAAAAACTCGCTAT‑GCGCAACTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTGCCGTg                                                                                           <  2:128554/148‑1 (MQ=255)
                                                        ccAGTTATTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGTATGTGCCAGAAGCAGTAAAAACTCGCTAT‑GCGCAACTGAATGCCGACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTGCCGTgg                                                                                          >  2:754518/1‑148 (MQ=255)
                                                         cAGTTATTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGTATGTGCCAGAAGCAGTAAAAACTCGCTAT‑GCGCAACTGAATGCCAACTGGCTGGAGATGAATAATCGGGTCAGAAAGGGCGATTTGCCGGGGTa                                                                                        >  1:1255931/1‑149 (MQ=255)
                                                           gTTATTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGTATGTGCCAGAAGCAGTAAAAACTCGCTAT‑GCGCAACTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGAGATTTGCCGTGGTAt                                                                                       >  2:252670/1‑148 (MQ=255)
                                                             tATTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGTATGTGCCAGAAGCAGTAAAAACTCGCTAT‑GCGCAACTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTGCCGTGGTATCAg                                                                                    <  1:55598/149‑1 (MQ=255)
                                                             tATTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGTATGTGCCAGAAGCAGTAAAAACTCGCTAT‑GCGCAACTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTGCCGTGGTATCAg                                                                                    <  1:407338/149‑1 (MQ=255)
                                                              aTTTCCGCAGGCACTGCATTCACCGGTATTCACCAACCTCCACGTCGTGTATGTGCCAGAAGGAGTAAAAACCCGCTAT‑GCGCAACTGAATGCCAACGGGCTGGTGATGAGTAATGGGCTCAGCAATGGCGATTTGCCGTGGTATCAg                                                                                    <  2:517014/148‑1 (MQ=255)
                                                               tttCAGCAGGCACTTCATTCACCCGTATTAACCAACCTCAACTTCTGCTATTTGCCAGAAGCAGTAAAAACTCGCTAT‑GCGCAACTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTGCCGTGGTATCAgg                                                                                   <  2:733844/148‑1 (MQ=255)
                                                                    gcagGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGTATGTGCCAGAAGCAGTAAAAACTCGCTAT‑GCGCAACTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTGCCGTGTTATCAGGCCAata                                                                             <  1:937705/149‑1 (MQ=255)
                                                                       ggcagtgaatgtaCCGGCATTAGCGAACCGCAACGTCTGGGAGGTAAGCTAAGCAGTAAAAGCGCTTTAT‑GCGCAACTGAATGCCAACTTGCGGGAGATGAATAATCGGCTCAGCAAGGGCGATTTGCCGTGGTATCAGGCCAATATta                                                                          <  1:1199410/137‑1 (MQ=255)
                                                                       ggCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGTATGTGCCAGAAGCAGTAAAAACTCGCTAT‑GCGCAACTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTGCCGTGGTATCAGGCCAATATta                                                                          <  1:258233/149‑1 (MQ=255)
                                                                        gCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGTGGGCTGCAGAAGCAGTAAAAACTCGCTAT‑GCGCAACTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTGCCGTGGTATCAGGCCAATATtaa                                                                         <  1:379201/149‑1 (MQ=255)
                                                                         cACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGTATGTGCCAGAAGCAGTAAAAACTCGCTAT‑GCGCAACTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTGCCGTGGTATCAGGCCAATATtaat                                                                        >  1:1197688/1‑149 (MQ=255)
                                                                          agctgaTTCACCGGTATTAACCAACCTCAACGTCTGGTATGTGCCAGAAGCAGTAAAAACTCGCTCT‑GCGCACCTGAATGCCAACTGGCTGGAAATGAAAAATCGGCTCACCAAGGGCGAGTTGCCGTGGTAACCGGCCAATAGtaat                                                                        >  2:145165/6‑148 (MQ=255)
                                                                          aCTGCATTCACCGGTATTAACCAACCTCAACGTCTGGTATGTGCCAGAAGCAGTAAAAACTCGCTAT‑GCGCAACTGAATGCCAACTGGCTGGAGCTGAATAATCGGCTCAGCAAGGGCGATTTGCCGTGGTATCAGGCCAATATtaat                                                                        <  2:663576/148‑1 (MQ=255)
                                                                               aTTCACCGGTATTAACCAACCTCAACGTCGGGGATGGTCCAGAGTCAGTAAAAACTCGCTAT‑GCGCAACTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTGCCGTGGTATCAGGCCAATATtaatacgta                                                                   <  2:262967/148‑5 (MQ=255)
                                                                                ttCACCGGTATTAACCAACCTCATCGTCTGGTATGTGCCAGAAGCAGTAAAAACTCGCTAT‑GCGCAACTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTGCCGTGGTATCAGGCCAATATTAATAATTAt                                                                  <  2:851039/148‑1 (MQ=255)
                                                                                ttCACCGGTATTAACCAACCTCAATGTCTGGTATGTGCCAGAAGCAGTAAAAACTCGCTAT‑GCGCAACTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTGCCGTGGTATCAGGCCAATATTAATAATTAt                                                                  <  2:350159/148‑1 (MQ=255)
                                                                                   aCCGGTATTAACCAACCTCAACGTCTGGTATGTGCCAGAAGCAGTAAAAACTCGCTAT‑GCGCAACTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTGCCGTGGTATCAGGCCAGTATTAATAATCATGtt                                                               >  2:1042259/1‑148 (MQ=255)
                                                                                    ccGGTATTAACCAACCTCAACGTCTGGTATGGGCCAGAAGCAGTAAAAACTCGCTAT‑GCGCAACTGAATGCCAACTGGATGGTGATGAATAATCGGCTCAGCAAGGGCGATTTGCCGTGGTATCAGGCCAATATTAATAATTATGTTa                                                              <  2:65853/148‑1 (MQ=255)
                                                                                     cGGTATTAAACCACCTCAACGTCTGGTATGTGCCAGAAGCAGTAAAAACTCGCTAT‑GCGCAACTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTGCCGTGGTATCAGGCCAATATTAATAATTATGTTAAt                                                            <  1:177095/149‑1 (MQ=255)
                                                                                       gTATTAACCAACCTCAACGTCTGGTATGTGCCAGAAGCAGTAAAAACTCGCTAT‑GCGCAACTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTGCCGTGGTATCAGGCCAATATTAATAATTATGTTAATCa                                                          >  1:495882/1‑149 (MQ=255)
                                                                                            aaccaaccTCAACGTATGGTATGTGCCAGAAGCAGTAAAAACTCGCTAT‑GCGCAACTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTGCCGTGGTATCAGGCCAATATTAATAATTATGTTAATCAGATa                                                      >  2:686415/1‑148 (MQ=255)
                                                                                              ccaaccTCAACGTCTGGTATGTGCCAGAAGCAGTAAAAACTCGCTAT‑GCGCAACTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTGCCGTGGTATCAGGCCAATATTAATAATTATGTTAATCAGATAGa                                                    <  2:1250346/148‑1 (MQ=255)
                                                                                              ccaaccTCAACGTCTGGTATGTGCCAGAAGCAGTAAAAACTCGCTAT‑GCGCAACTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTGCCGTGGTATCAGGCCAATATTAATAATTATGTTAATCAGATAGAc                                                   >  1:200559/1‑149 (MQ=255)
                                                                                               caaccTCAACGTCTGGTATGTGCCAGAAGCAGTAAAAACTCGCTAT‑GCGCAACTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTGCCGTGGTATCAGGCCAATATTAATAATTATGTTAATTAGATAGAcc                                                  >  1:299858/1‑149 (MQ=255)
                                                                                               caaccTCAACGTCTGGTATGTGCCAGAAGCAGTAAAAACTCGCTAT‑GCGCAACTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTGCCGTGGTATCAGGCCAATATTAATAATTATGTTAATCAGATAGAcc                                                  >  1:298755/1‑149 (MQ=255)
                                                                                                aaccTCAACGTCTGGTATGTGCCAGAAGCAGTAAAAACTCGCTAT‑GCGCAACTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTGCCGTGGTATCAGGCCAATAGTGATAAATATGTTAATCAGATAGACCt                                                 >  1:1226287/1‑149 (MQ=255)
                                                                                                aaccTCAACGTCTGGTATGTGCCAGAAGCAGTAAAAACTCGCTAT‑GAGCAACTGAATGCCAACTGGCTGGAGATGAATAATCGCCTCAGCAGGGGCGATTTGCCGTGGTATCAGGCCAATATTAATACTTCTGTTAAtgtgatgaacc                                                  >  2:748209/1‑138 (MQ=255)
                                                                                                   cTCAACGTCTGGTATGTGCCAGAAGCAGTAAAAACTCGCTAT‑GCGCAACTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTGCCGTGGTATCAGGCCAATATTAATAATTATGTTAATCAGATAGACCTGtt                                              <  1:684817/149‑1 (MQ=255)
                                                                                                   cTCAACGTCTGGTATGTGCCAGAAGCAGTAAAAACTCGCTAT‑GCGCAACTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTGCCGTGGTATCAGGCCAATATTAATAATTATGTTAATCAGATAGACCTGt                                               <  2:647188/148‑1 (MQ=255)
                                                                                                   cTCAACGTCTGGTATGTGCCAGAAGCAGTAAAAACTCGCTAT‑GCGCAACTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTGCCGTGGTATCAGGCCAATATTAATAATTATGTTAATCAGATAGACCTGt                                               <  2:624617/148‑1 (MQ=255)
                                                                                                   cTCAACGTCTGGTATGTGCCAGAAGCAGTAAAAACTCGCTAT‑GCGCAACTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTGACCAGGTAGTATGCCAATATTAATAATTCTCTTGATAAGGTAGAACTGt                                               >  2:1064724/1‑148 (MQ=255)
                                                                                                    tCAACGTCTGGTATGTGCCAGAAGCAGTAAAAACTCGCTAT‑GCGCAACTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTGCCGTGGTATCAGGCCAATATTAATAATTATGTTAATCAGATAGACCTGtt                                              <  2:1111783/148‑1 (MQ=255)
                                                                                                     caccGTCTGGTATGTGGCAGAAGCAGTAAAAAGTCGCTGT‑GCGCAACTGAAGGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTGCCGTGGTATCAGGCCAATATTAATAATTATGTTAATCAGATAGACCTGTTCg                                            <  1:930625/146‑1 (MQ=255)
                                                                                                     cAACGTCTGGTAGGTGCCAGAAGCAGTAAAAACTCGCTAT‑GCGCAACTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTGCCGTGGTATCAGGCCAATATTAATAATTATGTTAATCAGATAGACTTGTTc                                             <  2:872892/148‑1 (MQ=255)
                                                                                                            tGGTATGTGCCAGAAGCAGTAAAAACTCGCTAT‑GCGCAACTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTGCCGTGGTATCAGGCCAATATTAATAATTATGTTAATCAGATAGACCTGTTCGTACTGGc                                     >  1:761468/1‑149 (MQ=255)
                                                                                                              gTATGTGCCAGAAGCAGTAAAAACTCGCTAT‑GCGCAACTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTGCCGTGGTATCAGGCCGAGATTAATAATTATGTTAATCAGATAGACCTGTTCGTACTGGCtt                                   >  1:117877/1‑149 (MQ=255)
                                                                                                                  tggcaaGAAGCAGTAAAACCTCGCTAT‑GCGCAACTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTGCCGTGGTATCAGGCCAATATTAATAATTATGTTAATCAGATAGACCTGTTCGTACTGGCTTTAc                                <  2:492128/143‑1 (MQ=255)
                                                                                                                   tgCCAGAAGCAGTAAAAACTCGCTAT‑GCGCGACTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCCAGGGCGATTTGCGGTGGTAGCCGGCCACGATTAATAATTATGTTAATCAGATAGACCCATTCGTACTGGCtatcca                               >  2:670900/1‑143 (MQ=255)
                                                                                                                    gCCAGAAGCAGTAAAAACTCGCTAT‑GCGCAACTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTGCCGTGGTATCAGGCCAATATTAATCATTATGTTAATCAGATAGACCTGTTCGTACTGGCATTACAg                              >  2:1227306/1‑148 (MQ=255)
                                                                                                                     ccggAAGCAGTAAAAACTCGCTAT‑GCGCAACTGAGTGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTGCCGTGGTATCAGGCCAATATTAATAATTATGTTAATCAGATAGACCTGTTCGTACTGGCttattatg                             <  2:627173/145‑7 (MQ=255)
                                                                                                                     ccAGAAGCAGTAAAAACTCGCTAT‑GCGCAACTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTGCCGTGGTATCAGGCCAATATTAATAATTATGTTAATCAGATAGACCTGTTCGTACTGGCTTTACAGc                             >  2:861103/1‑148 (MQ=255)
                                                                                                                     caaGAAACAGTAAAACCTCGCTAC‑GCGCAACTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGACTTGCCGTGGTATCAGGCCAATATTAATAATTATGTTAATCAGATAGACCTGTTCGTACTGGCTTTACAGc                             <  2:59872/146‑1 (MQ=255)
                                                                                                                       aGAAGCAGTAAAAACTCGCTAT‑GCGCAACTGAATGCCAACTGGCAGGACATGAATAATCGGCTACGCAAGGGCGATTTGCCGCGGGAACAGGCCAATGGTTCAAATTATGTTTATCATATAGCCCTGTTCGTCCTGGCCTTACAGCAc                           >  2:110111/1‑148 (MQ=255)
                                                                                                                           gCAGTAAACACTCGCTAT‑GCGCAACTGAATGCCAACTGGCTGGAGATGACTAATCGGCTCAGCAAGGGCGATTTGCCGTGGTATCAGGCCAATATTAATAATTATCTTAATCAGATAGACCTGTTCGTACTGGCTTTACAGCACTACg                       <  2:952244/148‑1 (MQ=255)
                                                                                                                               tAAAAACTCGCTAT‑GCGCAACTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTGCCGTGGTATCAGGCCAATATTAATAATTATGTTAATCAGATAGACCTGTTCGTACTGGCTTTACAGCACTACGCTGa                   >  2:193193/1‑148 (MQ=255)
                                                                                                                               tAAAAACTCGCTAT‑GCGCAACTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTGCCGTGGTATCAGGCCAATATTAATAATTATGTTAAACAGATAGACCTGTTCGTACTGGCTTTACAGCACTACAGTGaa                  >  1:246767/1‑149 (MQ=255)
                                                                                                                                aaaaaCTCGCTAT‑GCGCAACTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTGCCGTGGTATCAGGCCAATATTAATAATTATGTTAATCAGATAGACCTGTTCGTACTGGCTTTACAGCACTACGCTGAAc                 <  1:77200/149‑1 (MQ=255)
                                                                                                                                   aaCTCGCTAT‑GCGCAACTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTGCCGTGGTATCAGGCTAATATTAATAATTATGTTAATCAGATAGACCTGTTCGTCCTGGCTTTACAGCACTACGCTGAACGCa              <  1:825593/149‑1 (MQ=255)
                                                                                                                                    aCTCGCTAT‑GCGCAACTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTGCCGTGGTATCAGGCCAATATTAATAATTATGTTAATCAGATAGACCTGTTCGTACTGGCTTTACAGCACTACGCTAAACGCaa             >  1:584267/1‑149 (MQ=255)
                                                                                                                                      tCGCTAT‑GCGCAACTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTGCCGTGGTATCAGGCCCATATTAATAATTATGTTAATAATATAGAACTGTTCTCAGTGGCTTGTCAGCACTACGCTGAAC‑Caaaa           >  2:439695/1‑148 (MQ=255)
                                                                                                                                      tCGCTAT‑GCGCAACTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTGCCGTGGTATCAGGCCAATATTAATAATTATGTTAATCAGATAGACCTGTTCGTACTGGCTTTACAGCACTACGCTGAACGCaaaa           >  1:711687/1‑149 (MQ=255)
                                                                                                                                      tCGCTAT‑GCGCAACTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTGCCGTGGTATCAGGCCAATATTAATAATTATGTTAATCAGATAGACCTGTTCGTACTGGCTTTACAGCACTACGCTGAACGCaaaa           <  1:4786/149‑1 (MQ=255)
                                                                                                                                      tCGCTAT‑GCGCAACAGAATGCCAACTGGCTGGACATGAATAATCGGCTCAGCAAGGGCGATTTGCCGTGGTATCAGGCCAAAATTAATAATTATGATAATCAGATAGACCTGTTCGTACTGGCTTTACAGCACTACGCTGAACGCaaaa           <  1:834310/149‑1 (MQ=255)
                                                                                                                                       cGCTAT‑GCGCAACTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTGCCGTGGTATCAGGCCAATATTAATAATTATGTTAATCGGATAGACCTGTTCGTACTGGCTTTACAGCACTACGCTGAACGCAAAAt          <  1:1188492/149‑1 (MQ=255)
                                                                                                                                          tAT‑GCGCAACTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAGGGGCGATTTGCCGTGGTATCAGGCCAATATTAATAATTATGTTAATCAGATAGACCTGTTCGTACTGGCTTTACAGCACTCCGCTGAACGCACAAtgc        >  2:241496/1‑148 (MQ=255)
                                                                                                                                          tAT‑GCGCAACTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTGCCGTGGTATCAGGCCAATATTAATAATTATGTTAATCAGATAGACCTGTTCGTACTGGCTTTACAGCACTACGCTGAACGCAAAAtgct       >  1:1257240/1‑149 (MQ=255)
                                                                                                                                               cgcAACTGAATGCCAACTGGCGGGAGATGAATAATCGGCTCAGCAAGGGCGATTTGCCGTGGTATCAGGCCAATATTAATACTTATGTTAATCAGATACACCTGTTCGTGCTGGCTTTACAGCACTACGCTGAACGCAAAATgctgct    <  2:357200/148‑1 (MQ=255)
                                                                                                                                                gcAACTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCAATTTGCCGTGGTATCAGGCCCATATTAATAATTATGTTAATCGGGTAGACCTCTTCGTCCTGGCTTTCCAGAACGACGCTGAACGCAAAATGCTGCggg  >  1:961105/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
AGTTACCGCCTGGGCTACGACTTGCAAAGTGGCAGTCCACAACTCAATGCACATCGCCAGTTATTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGTATGTGCCAGAAGCAGTAAAAACTCGCTAT‑GCGCATCTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTGCCGTGGTATCAGGCCAATATTAATAATTATGTTAATCAGATAGACCTGTTCGTACTGGCTTTACAGCACTACGCTGAACGCAAAATGCTGCTGG  >  CP009273/2579246‑2579537

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑