Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 614,003 T→A 9.1% I98I (ATT→ATA fepE → regulator of length of O‑antigen component of lipopolysaccharide chains

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273614,0030TA 9.1% 921.1 / 57.4 307I98I (ATT→ATAfepEregulator of length of O‑antigen component of lipopolysaccharide chains
Reads supporting (aligned to +/- strand):  ref base T (119/160);  new base A (15/12);  total (135/172)
Fisher's exact test for biased strand distribution p-value = 2.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

Only 100 of 312 total aligned reads displayed.

GCGTTTGCCTGCGCAGGCTTGCTGATCTCTTTCATCCTGCCGCAAAAATGGACCAGCGCGGCGGTTGTCACGCCTCCAGAACCTGTTCAGTGGCAAGAGTTGGAGAAATCATTCACTAAGCTTCGTGTGCTGGATCTGGATATCAAAATTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAGGCGAAAATCGACGAACTGGATTTGCATCGCGCAATTGTCGCATT  >  CP009273/613854‑614152
                                                                                                                                                     |                                                                                                                                                     
   tttGCCTTCTCAGGCTTGCTGATCTCTTTCATCCTGCCGCAAAAATGGACCAGCGCGGCGGTTGTCACGCCTCCAGAACCTGTTCAGTGGCAAGAGTTGGAGAAATCATTCACTAAGCTTCGTGTGCTGGATATGGATATCAAAATTGAt                                                                                                                                                    <  2:4186016/150‑1 (MQ=255)
    ttGCCTGCGCAGGCTTGCTGATCTCTTTCATCCTGCCGCAAAAATGGACCAGCGCGGCGGTTGTCACGCCTCCAGAACCTGTTCAGTGGCAAGAGTTGGAGAAATCATTCACTAAGCTTCGTGTGCTGGATCTGGATATCAAAATTGATc                                                                                                                                                   <  1:3924034/150‑1 (MQ=255)
    ttGCCTGCGCAGGCTTGCTGATCTCTTTCATCCTGCCGCAAAAATGGACCAGCGCGGCGGTTGTCACGCCTCCAGAACCTGTTCAGTGGCAAGAGTTGGAGAAATCATTCACTAAGCTTCGTGTGCTGGATCTGGATATCAAAATTGATc                                                                                                                                                   <  1:3655843/150‑1 (MQ=255)
    ttGCCTGCGCAGGCTTGCTGATCTCTTTCATCCTGCCGCAAAAATGGACCAGCGCGGCGGTTGTCACGCCTCCAGAACCTGTTCAGTGGCAAGAGTTGGAGAAATCATTCACTAAGCTTCGTGTGCTGGATCTGGATATCAAAATTGATc                                                                                                                                                   <  2:685909/150‑1 (MQ=255)
      gCCTGCGCAGGCTTGCTGATCTCTTTCATCCTGCCGCAAAAATGGACCAGCGCGGCGGTTGTCACGCCTCCAGAACCTGTTCAGTGGCAAGAGTTGGAGAAATCATTCACTAAGCTTCGTGTGCTGGATCTGGATATCAAAATTGATCGt                                                                                                                                                 <  1:4524980/150‑1 (MQ=255)
       ccTGCGCAGGCTTGCTGATCTCTTTCATCCTGCCGCAAAAATGGACCAGCGCGGCGGTTGTCACGCCTCCAGAACCTGTTCAGTGGCAAGAGTTGGAGAAATCATTCACTAAGCTTCGTGTGCTGGATCTGGATATCAAAATTGATCGTa                                                                                                                                                <  1:5157946/150‑1 (MQ=255)
       ccTGCGCAGGCTTGCTGATCTCTTTCATCCTGCCGCAAAAATGGACCAGCGCGGCGGTTGTCACGCCTCCAGAACCTGTTCAGTGGCAAGAGTTGGAGAAATCATTCACTAAGCTTCGTGTGCTGGATCTGGATATCAAAATAGATCGTa                                                                                                                                                <  1:1311252/150‑1 (MQ=255)
        cTGCGCAGGCTTGCTGATCTCTTTCATCCTGCCGCAAAAATGGACCAGCGCGGCGGTTGTCACGCCTCCAGAACCTGTTCAGTGGCAAGAGTTGGAGAAATCATTCACTAAGCTTCGTGTGCTGGATCTGGATATCAAAATTGATCGTAc                                                                                                                                               <  2:2603644/150‑1 (MQ=255)
        cTGCGCAGGCTTGCTGATCTCTTTCATCCTGCCGCAAAAATGGACCAGCGCGGCGGTTGTCACGCCTCCAGAACCTGTTCAGTGGCAAGAGTTGGAGAAATCATTCACTAAGCTTCGTGTGCTGGATCTGGATATCAAAATTGATCGTAc                                                                                                                                               <  1:3134755/150‑1 (MQ=255)
          gcgcAGGCTTGCTGATCTCTTTCATCCTGCCGCAAAAATGGACCAGCGCGGCGGTTGTCACGCCTCCAGAACCTGTTCAGTGGCAAGAGTTGGAGAAATCATTCACTAAGCTTCGTGTGCTGGATCTGGATATCAAAATTGATCGTACAg                                                                                                                                             <  2:5130673/150‑1 (MQ=255)
           cgcAGGCTTGCTGATCTCTTTCATCCTGCCGCAAAAATGGACCAGCGCGGCGGTTGTCACGCCTCCAGAACCTGTTCAGTGGCAAGAGTTGGAGAAATCATTCACTAAGCTTCGTGTGCTGGATCTGGATATCAAAATTGATCGTACAGa                                                                                                                                            <  2:3330847/150‑1 (MQ=255)
              aGGCTTGCTGATCTCTTTCATCCTGCCGCAAAAATGGACCAGCGCGGCGGTTGTCACGCCTCCAGAACCTGTTCAGTGGCAAGAGTTGGAGAAATCATTCACTAAGCTTCGTGTGCTGGATCTGGATATCAAAATTGATCGTACAGAAGc                                                                                                                                         >  2:743906/1‑150 (MQ=255)
                 cTTGCTGATCTCTTTCATCCTGCCGCAAAAATGGACCAGCGCGGCGGTTGTCACGCCTCCAGAACCTGTTCAGTGGCAAGAGTTGGAGAAATCATTCACTAAGCTTCGTGTGCTGGATCTGGATATCAAAATTGATCGTACAGAAGCAtt                                                                                                                                      <  2:2817605/150‑1 (MQ=255)
                 cTTGCTGATCTCTTTCATCCTGCCGCAAAAATGGACCAGCGCGGCGGTTGTCACGCCTCCAGAACCTGTTCAGTGGCAAGAGTTGGAGAAATCATTCACTAAGCTTCGTGTGCTGGATCGGGATATCAAAATAGATCGTACAGAAGCAtt                                                                                                                                      >  2:1058528/1‑150 (MQ=255)
                  ttGCTGATCTCTTTCATCCTGCCGCAAAAATGGACCAGCGCGGCGGTTGTCACGCCTCCAGAACCTGTTCAGTGGCAAGAGTTGGAGAAATCATTCACTAAGCTTCGTGTGCTGGATCTGGATATCAAAATTGATCGTACAGAAGCAttt                                                                                                                                     <  2:439254/150‑1 (MQ=255)
                   tGCTGATCTCTTTCATCCTGCCGCAAAAATGGACCAGCGCGGCGGTTGTCACGCCTCCAGAACCTGTTCAGTGGCAAGAGTTGGAGAAATCATTCACTAAGCTTCGTGTGCTGGATCTGGATATCAAAATTGATCGTACAGAAGCATTTa                                                                                                                                    <  1:5095753/150‑1 (MQ=255)
                   tGCTGATCTCTTTCATCCTGCCGCAAAAATGGACCAGCGCGGCGGTTGTCACGCCTCCAGAACCTGTTCAGTGGCAAGAGTTGGAGAAATCATTCACTAAGCTTCGTGTGCTGGATCTGGATATCAAAATTGATCGTACAGAAGCATTTa                                                                                                                                    >  1:3323379/1‑150 (MQ=255)
                   tGCTGATCTCTTTCATCCTGCCGCAAAAATGGACCAGCGCGGCGGTTGTCACGCCTCCAGAACCTGTTCAGTGGCAAGAGTTGGAGAAATCATTCACTAAGCTTCGTGTGCTGGATCTGGATATCAAAATTGATCGTACAGAAGCATTTa                                                                                                                                    >  1:1707627/1‑150 (MQ=255)
                    gCTGATCTCTTTCATCCTGCCGCAAAAATGGACCAGCGCGGCGGTTGTCACGCCTCCAGAACCTGTTCAGTGGCAAGAGTTGGAGAAATCATTCACTAAGCTTCGTGTGCTGGATCTGGATATCAAAATTGATCGTACAGAAGCATTTaa                                                                                                                                   >  1:3418794/1‑150 (MQ=255)
                      tGATCTCTTTCATCCTGCCGCAAAAATGGACCAGCGCGGCGGTTGTCACGCCTCCAGAACCTGTTCAGTGGCAAGAGTTGGAGAAATCATTCACTAAGCTTCGTGTGCTGGATCTGGATATCAAAATTGATCGTACAGAAGCATTTAAcc                                                                                                                                 <  1:2225351/150‑1 (MQ=255)
                      tGATCTCTTTCATCCTGCCGCAAAAATGGACCAGCGCGGCGGTTGTCACGCCTCCAGAACCTGTTCAGTGGCAAGAGTTGGAGAAATCATTCACTAAGCTTCGTGTGCTGGATCTGGATATCAAAATAGATCGTACAGAAGCATTTAAcc                                                                                                                                 <  2:3803326/150‑1 (MQ=255)
                       gATCTCTTTCATCCTGCCGCAAAAATGGACCAGCGCGGCGGTTGTCACGCCTCCAGAACCTGTTCAGTGGCAAGAGTTGGAGAAATCATTCACTAAGCTTCGTGTGCTGGATCTGGATATCAAAATTGATCGTACAGAAGCATTTAACCt                                                                                                                                <  2:4613057/150‑1 (MQ=255)
                          ctctTTCATCCTGCCGCAAAAATGGACCAGCGCGGCGGTTGTCACGCCTCCAGAACCTGTTCAGTGGCAAGAGTTGGAGAAATCATTCACTAAGCTTCGTGTGCTGGATCTGGATATCAAAATTGATCGTACAGAAGCATTTAACCTGtt                                                                                                                             <  1:1606683/150‑1 (MQ=255)
                           tctTTCATCCTGCCGCAAAAATGGACCAGCGCGGCGGTTGTCACGCCTCCAGAACCTGTTCAGTGGCAAGAGTTGGAGAAATCATTCACTAAGCTTCGTGTGCTGGATCTGGATATCAAAATTGATCGTACAGAAGCATTTAACCTGttt                                                                                                                            <  2:5292781/150‑1 (MQ=255)
                           tctTTCATCCTGCCGCAAAAATGGACCAGCGCGGCGGTTGTCACGCCTCCAGAACCTGTTCAGTGGCAAGAGTTGGAGAAATCATTCACTAAGCTTCGTGTGCTGGATCTGGATATCAAAATAGATCGTACAGAAGCATTTAACCTGttt                                                                                                                            <  1:2407022/150‑1 (MQ=255)
                             tttCATCCTGCCGCAAAAATGGACCAGCGCGGCGGTTGTCACGCCTCCAGAACCTGTTCAGTGGCAAGAGTTGGAGAAATCATTCACTAAGCTTCGTGTGCTGGATCTGGATATCAAAATTGATCGTACAGAAGCATTTAACCTGTTTAt                                                                                                                          >  2:3278064/1‑150 (MQ=255)
                             tttCATCCTGCCGCAAAAATGGACCAGCGCGGCGGTTGTCACGCCTCCAGAACCTGTTCAGTGGCAAGAGTTGGAGAAATCATTCACTAAGCTTCGTGTGCTGGATCTGGATATCAAAATTGATCGTACAGAAGCATTTAACCTGTTTAt                                                                                                                          >  1:1890333/1‑150 (MQ=255)
                             tttCATCCTGCCGCAAAAATGGACCAGCGCGGCGGTTGTCACGCCTCCAGAACCTGTTCAGTGGCAAGAGTTGGAGAAATCATTCACTAAGCTTCGTGTGCTGGATCTGGATATCAAAATTGATCGTACAGAAGCATTTAACCTGTTTAt                                                                                                                          <  2:3889415/150‑1 (MQ=255)
                               tCATCCTGCCGCAAAAATGGACCAGCGCGGCGGTTGTCACGCCTCCAGAACCTGTTCAGTGGCAAGAGTTGGAGAAATCATTCACTAAGCTTCGTGTGCTGGATCTGGATATCAAAATTGATCGTACAGAAGCATTTAACCTGTTTATCa                                                                                                                        >  2:4447542/1‑150 (MQ=255)
                                       ccGCAAAAATGGACCAGCGCGGCGGTTGTCACGCCTCCAGAACCTGTTCAGTGGCAAGAGTTGGAGAAATCATTCACTAAGCTTCGTGTGCTGGATCTGGATATCAAAATTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGttt                                                                                                                <  1:3881026/150‑1 (MQ=255)
                                       ccGCAAAAATGGACCAGCGCGGCGGTTGTCACGCCTCCAGAACCTGTTCAGTGGCAAGAGTTGGAGAAATCATTCACTAAGCTTCGTGTGCTGGATCTGGATATCAAAATTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGttt                                                                                                                <  2:1330450/150‑1 (MQ=255)
                                            aaaaTGGACCAGCGCGGCGGTTGTCACGCCTCCAGAACCTGTTCAGTGGCAAGAGTTGGAGAAATCATTCACTAAGCTTCGTGTGCTGGATCTGGATATCAAAATTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTc                                                                                                           >  2:1951512/1‑150 (MQ=255)
                                            aaaaTGGACCAGCGCGGCGGTTGTCACGCCTCCAGAACCTGTTCAGTGGCAAGAGTTGGAGAAATCATTCACTAAGCTTCGTGTGCTGGATCTGGATATCAAAATTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTc                                                                                                           <  1:1012269/150‑1 (MQ=255)
                                             aaaTGGACCAGCGCGGCGGTTGTCACGCCTCCAGAACCTGTTCAGTGGCAAGAGTTGGAGAAATCATTCACTAAGCTTCGTGTGCTGGATCTGGATATCAAAATTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCg                                                                                                          >  2:1009723/1‑150 (MQ=255)
                                             aaaTGGACCAGCGCGGCGGTTGTCACGCCTCCAGAACCTGTTCAGTGGCAAGAGTTGGAGAAATCATTCACTAAGCTTCGTGTGCTGGATCTGGATATCAAAATTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCg                                                                                                          >  1:1931926/1‑150 (MQ=255)
                                             aaaTGGACCAGCGCGGCGGTTGTCACGCCACCAGAAACAGTTCAGTGGCAAGAGTTGTAGAAAACATTCAATAATCATAGTGTGCTGGATCTGGATATCAAAAATTAACGTACAGAAGCATATAAACTGTTTATAAAGAAGTTTCAGTCg                                                                                                          >  1:2270754/1‑150 (MQ=255)
                                                 ggACCAGCGCGGCGGTTGTCACGCCTCCAGAACCTGTTCAGTGGCAAGAGTTTGAGAAATCATTCACTAAGCTTCGTGTGCAGGATCAGGATATAAAAATTGATCGTACAGAATCATTAAACCTGTTTATCAAGAAGTTTCAGTCGGTTa                                                                                                      >  1:3675326/1‑150 (MQ=255)
                                                 ggACCAGCGCGGCGGTTGTCACGCCTCCAGAACCTGTTCAGTGGCAAGAGTTGGAGAAATCATTCACTAAGCTTCGTGTGCTGGATCTGGATATCAAAATAGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTa                                                                                                      <  1:845522/150‑1 (MQ=255)
                                                    ccAGCGCGGCGGTTGTCACGCCTCCAGAACCTGTTCAGTGGCAAGAGTTGGAGAAATCATTCACTAAGCTTCGTGTGCTGGATCTGGATATCAAAATTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCt                                                                                                   <  1:3173629/150‑1 (MQ=255)
                                                      aGCGCGGCGGTTGTCACGCCTCCAGAACCTGTTCAGTGGCAAGAGTTGGAGAAATCATTCACTAAGCTTCGTGTGCTGGATCTGGATATCAAAATTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTg                                                                                                 <  1:3820665/150‑1 (MQ=255)
                                                         gcggcgGTTTTCACGCATCCAGAAACACTTCAGTGGCAAGAGTTGGAGAAATCATTCACTAAGAAAAGTGTGCTGTATCATGATATCAAAATTGATCGTACATAACCATTTAACCTGATTATCAAGAAGTTTCAGTCGGTTAGAtttatg                                                                                              >  2:1631334/1‑146 (MQ=255)
                                                         gcggcgGTTGTCACGCCTCCAGAACCTGTTCAGTGGCAAGAGTTGGAGAAATCATTCACTAAGCTTCGTGTGCTGGATCTGGATATCAAAATAGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGTTTGCTg                                                                                              >  1:26971/1‑150 (MQ=255)
                                                             cggTTGTCACGCCTCCAGAACCTGTTCAGTGGCAAGAGTTGGAGAAATCATTCACTAAGCTTCGTGTGCTGGATCTGGATATCAAAATAGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAag                                                                                          >  1:4907769/1‑150 (MQ=255)
                                                              ggTTGTCACGCCTCCAGAACCTGTTCAGTGGCAAGAGTTGGAGAAATCATTCACTAAGCTTCGTGTGCTGGATCTGGATATCAAAATTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAaga                                                                                         <  2:3512607/150‑1 (MQ=255)
                                                                    cACGCCTCCAGAACCTGTTCAGTGGCAAGAGTTGGAGAAATCATTCACTAAGCTTCGTGTGCTGGATCTGGATATCAAAATTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCt                                                                                   <  1:5298524/150‑1 (MQ=255)
                                                                         cTCCAGAACCTGTTCAGTGGCAAGAGTTGGAGAAATCATTCACTAAGCTTCGTGTGCTGGATCTGGATATCAAAATTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGtt                                                                              <  2:429907/150‑1 (MQ=255)
                                                                          tCCAGAACCTGTTCAGTGGCAAGAGTTGGAGAAATCATTCACTAAGCTTCGTGTGCTGGATCTGGATATCAAAATTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTtc                                                                             <  2:5234667/150‑1 (MQ=255)
                                                                           ccAGAACCTGTTCAGTGGCAAGAGTTGGAGAAATCATTCACTAAGCTTCGTGTGCTGGATCTGGATATCAAAATTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTtca                                                                            >  1:4308938/1‑150 (MQ=255)
                                                                                  cTGTTCAGTGGCAAGAGTTGGAGAAATCATTCACTAAGCTTCGTGTGCTGGATCTGGATATCAAAATTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCtt                                                                     <  2:3093526/150‑1 (MQ=255)
                                                                                       cAGTGGCAAGAGTTGGAGAAATCATTCACTAAGCTTCGTGTGCTGGATCTGGATATCAAAATTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTAtgtg                                                                >  1:590210/1‑150 (MQ=255)
                                                                                        aGTGGCAAGAGTTGGAGAAATCATTCACTAAGCTTCGTGTGCTGGATCTGGATATCAAAATTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGa                                                               <  1:708271/150‑1 (MQ=255)
                                                                                         gTGGCAAGAGTTGGAGAAATCATTCACTAAGCTTCGTGTGCTGGATCTGGATATCAAAATTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGAt                                                              <  1:1596908/150‑1 (MQ=255)
                                                                                         gTGGCAAGAGTTGGAGAAATCATTCACTAAGCTTCGTGTGCTGGATCTGGATATCAAAATTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGAt                                                              >  1:133574/1‑150 (MQ=255)
                                                                                          tGGCAAGAGTTGGAGAAATCATTCACTAAGCTTCGTGTGCTTGATCTGGATATCAAAATTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATg                                                             <  1:1531688/150‑1 (MQ=255)
                                                                                          tGGCAAGAGTTGGAGAAATCATTCACTAAGCTTCGTGTGCTGGATCTGGATATCAAAATTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATg                                                             <  1:2119905/150‑1 (MQ=255)
                                                                                             cAAGAGTTGGAGAAATCATTCACTAAGCTTCGTGTGCTGGATCTGGATATCAAAATTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGAc                                                          <  2:1770914/150‑1 (MQ=255)
                                                                                              aaGAGTTGGAGAAATCATTCACTAAGCTTCGTGTGCTGGATCTGGATATCAAAATTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGAcc                                                         >  1:1789460/1‑150 (MQ=255)
                                                                                                 agTTGGAGAAATCATTCACTAAGCTTCGTGTGCTGGATCTGGATATCAAAATTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAAt                                                      <  2:874973/150‑1 (MQ=255)
                                                                                                 agTTGGAGAAATCATTCACTAAGCTTCGTGTGCTGGATCTGGATATCAAAATTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAAt                                                      <  2:3584332/150‑1 (MQ=255)
                                                                                                  gTTGGAGAAATCATTCACTAAGCTTCGTGTGCTGGATCTGGATATCAAAATTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAAtt                                                     <  2:2784921/150‑1 (MQ=255)
                                                                                                  gTTGGAGAAATCATTCACTAAGCTTCGTGTGCTGGATCTGGATATCAAAATTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAAtt                                                     <  1:4101452/150‑1 (MQ=255)
                                                                                                    tGGAGAAATCATTCACTAAGCTTCGTGTGCTGGATCTGGATATCAAAATTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTaa                                                   >  2:2854/1‑150 (MQ=255)
                                                                                                    tGGAGAAATCATTCACTAAGCTTCGTGTGCTGGATCTGGATATCAAAATTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTaa                                                   >  2:3152154/1‑150 (MQ=255)
                                                                                                      gagaAATCATTCACTAAGCTTCGTGTGCTGGATCTGGATATCAAAATTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGTGTTCATCACCTTATGTGATGGACCAATTaaaa                                                 >  2:2035175/1‑150 (MQ=255)
                                                                                                      gagaAATCATTCACTAAGCTTCGTGTGCTGGATCTGGATATCAAAATTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTaaaa                                                 >  1:2897133/1‑150 (MQ=255)
                                                                                                       agaAATCATTCACTAAGCTTCGTGTGCTGGATCTGGATATCAAAATTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGATGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAag                                                >  2:1160649/1‑150 (MQ=255)
                                                                                                        gaAATCATTCACTAAGCTTCGTGTGCTGGATCTGGATATCAAAATTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGTTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAaga                                               <  2:5184938/150‑1 (MQ=255)
                                                                                                         aaaTCATTCACTAAGCTTCGTGTGCTGGATCTGGATATCAAAATTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAagag                                              <  1:4351636/150‑1 (MQ=255)
                                                                                                         aaaTCATTCACTAAGCTTCGTGTGCTGGATCTGGATATCAAAATTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAagag                                              >  1:1770319/1‑150 (MQ=255)
                                                                                                         aaaTCATTCACTAAGCTTCGTGTGCTGGATCTGGATATCAAAATTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAagag                                              <  1:2962644/150‑1 (MQ=255)
                                                                                                         aaaTCATTCACTAAGCTTCGTGTGCTGGATCTGGATATCAAAATTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAagag                                              >  2:4597492/1‑150 (MQ=255)
                                                                                                         aaaTCATTCACTAAGCTTCGTGTGCTGGATCTGGATATCAAAATTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAagag                                              >  2:2123766/1‑150 (MQ=255)
                                                                                                          aaTCATTCACTAAGCTTCGTGTGCTGGATCTGGATATCAAAATTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAgg                                             >  2:2503883/1‑150 (MQ=255)
                                                                                                            tCATTCACTACGCTTCGTGTGCTGGATCTGGATATCAAAATTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAGGCg                                           <  2:3986065/150‑1 (MQ=255)
                                                                                                             cATTCACTAAGCTTCGTGTGCTGGATCTGGATATCAAAATTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCa                                                        <  1:4966011/136‑1 (MQ=255)
                                                                                                                tCACTAAGCTTCGTGTGCTGGATCTGGATATCAAAATTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAGGCGaaaa                                       >  2:4912814/1‑150 (MQ=255)
                                                                                                                 cACTAAGCTTCGTGTGCTGGATCTGGATATCAAAATTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAGGCGAAAAt                                      >  1:2507547/1‑150 (MQ=255)
                                                                                                                  aCTAAGCTTCGTGTGCTGGATCTTGATATCAAAATTGATAGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGTAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAGGCGAAAATc                                     <  1:4041328/150‑1 (MQ=255)
                                                                                                                      aGCTTCGTGTGCTGGATCTGGATATCAAAATTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAGGCGAAAATcgact                                 >  2:5077388/1‑149 (MQ=255)
                                                                                                                      aGCTTCGTGTGCTGGATCTGGATATCAAAATTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAGGCGAAAATcgacg                                 <  1:4840095/150‑1 (MQ=255)
                                                                                                                        cTTCGTGTGCTGGATCTGGATATCAAAATTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAGGCGAAAATCGACGaa                               >  1:777423/1‑150 (MQ=255)
                                                                                                                           cGTGTGCTGGATCTGGATATCAAAATTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAGGCGAAAATCGACGAACTg                            <  2:2217/150‑1 (MQ=255)
                                                                                                                              gtgCTGGATCTGGATATCAAAATTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAGGCGAAAATCGACGAACTGGAt                         >  2:1502146/1‑150 (MQ=255)
                                                                                                                              gtgCTGGATCTGGATATCAAAATTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAGGCGAAAATCGACGAACTGGAt                         <  1:4395984/150‑1 (MQ=255)
                                                                                                                                    gatctggatATCAAAATTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAGGCGAAAATCGACGAACTGGATTTGCAt                   <  2:4235474/150‑1 (MQ=255)
                                                                                                                                     atctggatATCAAAATTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAGGCGAAAATCGACGAACTGGATTTGCATc                  <  1:166591/150‑1 (MQ=255)
                                                                                                                                       ctggatATCAAAATTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAGGCGAAAATCGACGAACTGGATTTGCATcgc                >  1:1193924/1‑150 (MQ=255)
                                                                                                                                          gatATCAAAATTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAGGCGAAAATCGACGAACTGGATTTGCATCGCGCa             >  2:4524198/1‑150 (MQ=255)
                                                                                                                                           atatCAAAATTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAGGCGAAAATCGACGAACTGGATTTGCATCGCGCaa            >  2:2904325/1‑150 (MQ=255)
                                                                                                                                           atatCAAAATTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAGGCGAAAATCGACGAACTGGATTTGCATCGCGCaa            <  2:321485/150‑1 (MQ=255)
                                                                                                                                             atCAAAATTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAGGCGAAAATCGACGAACTGGATTTGCATCGCGCAAtt          <  1:3103511/150‑1 (MQ=255)
                                                                                                                                             atCAAAATTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAGGCGAAAATCGACGAACTGGATTTGCATCGCGCAAtt          <  2:469158/150‑1 (MQ=255)
                                                                                                                                              tCAAAATTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAGGCGAAAATCGACGAACTGGATTTGCATCGCGCAATTg         <  2:243392/150‑1 (MQ=255)
                                                                                                                                              tCAAAATTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAGGCGAAAATCGACGAACTGGATTTGCATCGCGCAATTg         >  1:4601361/1‑150 (MQ=255)
                                                                                                                                               cAAAATTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAGGCGAAAATCGACGAACTGGATTTGCATCGCGCAATTGt        <  2:3391077/150‑1 (MQ=255)
                                                                                                                                                aaaaTTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAGGCGAAAATCGACGAACTGGATTTGCATCGCGCAATTGTc       >  1:3463975/1‑150 (MQ=255)
                                                                                                                                                 aaaTTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAGGCGAAAATCGACGAACTGGATTTGCATCGCGCAATTGTCg      >  1:1450249/1‑150 (MQ=255)
                                                                                                                                                  aaTTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAGGCGAAAATCGACGAACTGGATTTGCATCGCGCAATTGTCGc     >  1:3606601/1‑150 (MQ=255)
                                                                                                                                                  aaTTGATCATACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAGGCGAAAATCGACGAACTGGATTTGCATCGCGCAATTGTCGc     >  1:4739442/1‑150 (MQ=255)
                                                                                                                                                     tGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAGGCGAAAATCGACGAACTGGATTTGCATCGCGCAATTGTCGCAtt  <  1:5272979/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCGTTTGCCTGCGCAGGCTTGCTGATCTCTTTCATCCTGCCGCAAAAATGGACCAGCGCGGCGGTTGTCACGCCTCCAGAACCTGTTCAGTGGCAAGAGTTGGAGAAATCATTCACTAAGCTTCGTGTGCTGGATCTGGATATCAAAATTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAGGCGAAAATCGACGAACTGGATTTGCATCGCGCAATTGTCGCATT  >  CP009273/613854‑614152

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: