Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,288,739 C→A 7.3% intergenic (‑361/‑244) hns ← / → tdk global DNA‑binding transcriptional dual regulator H‑NS/thymidine kinase/deoxyuridine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,288,7390CA 7.3% 1195.7 / 37.1 327intergenic (‑361/‑244)hns/tdkglobal DNA‑binding transcriptional dual regulator H‑NS/thymidine kinase/deoxyuridine kinase
Reads supporting (aligned to +/- strand):  ref base C (148/154);  new base A (10/15);  total (158/169)
Fisher's exact test for biased strand distribution p-value = 4.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 341 total aligned reads displayed.

CTGATTCATGTCACATTTTATGGGGAGATTATCGTAGGCTGACGACCTTTCAGTCTTCTGTATTAGTTGTGTTTACGAGAATTCCCTATTAAGCGAATGATGAAAAGTAGAACAGTCGCAATAAGAGCATGGACTTAGTATTGCACTATCTCCTGGAGGTCAACAGAGGGCTATTACTTGCGCAACAGGTTAAAGATTGTGAATAGTTACCAGCAGTCATTTACCCGCTTATAACAAGCGAGGCAGTTGTAATGATAGCTCAGAAGGATTATGCAAGGCTTCGTAAGGGAGAACGCATA  >  CP009273/1288590‑1288888
                                                                                                                                                     |                                                                                                                                                     
   aTTCATGTCACATTTTATGGGGAGATTATCGTAGGCTGACGACCTTTCAGTCTTCTGTATTAGTTGTGTTTACGAGTATTCCCTATTAAGCGAATGATGAAAAGTAGAACAGTCGCAATAAGAGCATGGACTTAGTATTGCACTATCTcc                                                                                                                                                    <  2:4753524/150‑1 (MQ=255)
   aTTCATGTCACATTTTATGGGGAGATTATCGTAGGCTGACGACCTTTCAGTCTTCTGTATTAGTTGTGTTTACGAGTATTCCCTATTAAGCGAATGATGAAAAGTAGAACAGTCGCAATAAGAGCATGGACTTAGTATTGCACTATCTcc                                                                                                                                                    >  2:3051201/1‑150 (MQ=255)
    ttCATGTCACATTTTATGGGGAGATTATCGTAGGCTGACGAACTTTCAGTCTTCTGTATTAGTTGTGTTTACGAGTATTCCCTATTAAGCGAATGATGAAAAGTAGAACAGTCGCAATAAGAGCATGGACTTAGTATTGCACTATATcca                                                                                                                                                   <  2:3973527/150‑2 (MQ=255)
      cATGTCACATTTTATGGGGAGATTATCGTAGGCTGACGACCTTTCAGTCTTCTGTATTAGTTGTGTTTACGAGTATTCCCTATTAAGCGAATGATGAAAAGTAGAACAGTCGCAATAAGAGCATGGACTTAGTATTGCACTATCTCCTgg                                                                                                                                                 >  1:4667559/1‑150 (MQ=255)
      cATGTCACATTTTATGGGGAGATTATCGTAGGCTGACGACCTTTCAGTCTTCTGTATTAGTTGTGTTTACGAGTATTCCCTATTAAGCGAATGATGAAAAGTAGAACAGTCGCAATAAGAGCATGGACTTAGTATTGCACTATCTCCTgg                                                                                                                                                 >  1:4413018/1‑150 (MQ=255)
       aTGTCACATTTTATGGGGAGATTATCGTAGGCTGACGACCTTTCAGTCTTCTGTATTAGTTGTGTTTACGAGTATTCCCTATTAAGCGAATGATGAAAAGTAGAACAGTCGCAATAAGAGCATGGACTTAGTATTGCACTATCTCCTGGa                                                                                                                                                >  2:5321409/1‑150 (MQ=255)
        tGTCACATTTTATGGGGAGATTATCGTAGGCTGACGACCTTTCAGTCTTCTGTATTAGTTGTGTTTACGAGTATTCCCTATTAAGCGAATGATGAAAAGTAGAACAGTCGCAATAAGAGCATGGACTTAGTATTGCACTATCTCCTGGAg                                                                                                                                               <  1:3291916/150‑1 (MQ=255)
          tCACATTTTATGGGGAGATTATCGTAGGCTGACGACCTTTCAGTCTTCTGTATTAGTTGTGTTTACGAGTATTCCCTATTAAGCGAATGATGAAAAGTAGAACAGTCGCAATAAGAGCATGGACTTAGTATTGCACTATCTCCTGGAGGt                                                                                                                                             >  2:2032510/1‑150 (MQ=255)
           cacaTTTTATGGGGAGATTATCGTAGGCTGACGACCTTTCAGTCTTCTGTATTAGTTGTGTTTACGAGTATTCCCTATTAAGCGAATGATGAAAAGTAGAACAGTCGCAATAAGAGCATGGACTTAGTATTGCACTATCTCCTGGAGGTc                                                                                                                                            >  2:1133152/1‑150 (MQ=255)
           cacaTTTTATGGGGAGATTATCGTAGGCTGACGACCTTTCAGTCTTCTGTATTAGTTGTGTTTACGAGTATTCCCTATTAAGCGAATGATGAAAAGTAGAACAGTCGCAATAAGAGCATGGACTTAGTATTGCACTATCTCCTGGAGGTc                                                                                                                                            <  1:3503878/150‑1 (MQ=255)
           cacaTTTTATGGGGAGATTATCGTAGGCTGACGACCTTTCAGTCTTCTGTATTAGTTGTGTTTACGAGTATTCCCTATTAAGCGAATGATGAAAAGTAGAACAGTCGCAATAAGAGCATGGACTTAGTATTGCACTATCTCCTGGAGGTc                                                                                                                                            <  1:2434732/150‑1 (MQ=255)
            acaTTTTATGGGGAGATTATCGTAGGCTGACGACCTTTCAGTCTTCTGTATTAGTTGTGTTTACGAGTATTCCCTATTAAGCGAATGATGAAAAGTAGAACTGTCGCAATAAGAGCATGGACTTAGTATTGCACTATCTCCTGGAGGTCa                                                                                                                                           <  2:889608/150‑1 (MQ=255)
              aTTTTATGGGGAGATTATCGTAGGCTGACGACCTTTCAGTCTTCTGTATTAGTTGTGTTTACGAGTATTCCCTATTAAGCGAATGATGAAAAGTAGAACAGTCGCAATAAGAGCATGGACTTAGTATTGCACTATCTCCTGGAGGTCAAc                                                                                                                                         <  2:4219000/150‑1 (MQ=255)
               ttttATGGGGAGATTATCGTAGGCTGACGACCTTTCAGTCTTCTGTATTAGTTGTGTTTACGAGTATTCCCTATTAAGCGAATGATGAAAAGTAGAACAGTCGCAATAAGAGCATGGACTTAGTATTGCACTATCTCCTGGAGGTCAACa                                                                                                                                        >  2:2295376/1‑150 (MQ=255)
                   aTGGGGAGATTATCGTAGGCTGACGACCTTTCAGTCTTCTGTATTAGTTGTGTTTACGAGTATTCCCTATTAAGCGAATGATGAAAAGTAGAACAGTCGCAATAAGAGCATGGACTTAGTATTGCACTATCTCCTGGAGGTCAACAGAgg                                                                                                                                    <  1:1170923/150‑1 (MQ=255)
                    tGGGGAGATTATCGTAGGCTGACGACCTTTCAGTCTTCTGTATTAGTTGTGTTTACGAGTATTCCCTATTAAGCGAATGATGAAAAGTAGAACAGTCGCAATAAGAGCATGGACTTAGTATTGCACTATCTCCTGGAGGTCAACAGAggg                                                                                                                                   <  1:1056373/150‑1 (MQ=255)
                     ggggAGATTATCGTAGGCTGACGACCTTTCAGTCTTCTGTATTAGTTGTGTTTACGAGTATTCCCTATTAAGCGAATGATGAAAAGTAGAACAGTCGCAATAAGAGCATGGACTTAGTATTGCACTATCTCCTGGAGGTCAACAGAGGGc                                                                                                                                  >  2:2804667/1‑150 (MQ=255)
                      gggAGATTATCGTAGGCTGACGACCTTTCAGTCTTCTGTATTAGTTGTGTTTACGAGTATTCCCTATTAAGCGAATGATGAAAAGTAGAACAGTCGCAATAAGAGCATGGACTTAGTATTGCACTATCTCCTGGAGGTCAACAGAGGGCt                                                                                                                                 <  1:2310302/150‑1 (MQ=255)
                      gggAGATTATCGTAGGCTGACGACCTTTCAGTCTTCTGTATTAGTTGTGTTTACGAGTATTCCCTATTAAGCGAATGATGAAAAGTAGAACAGTCGCAATAAGAGCATGGACTTAGTATTGCACTATCTCCTGGAGGTCAACAGAGGGCt                                                                                                                                 <  1:1391135/150‑1 (MQ=255)
                       ggAGATTATCGTAGGCTGACGACCTTTCAGTCTTCTGTATTAGTTGTGTTTACGAGTATTCCCTATTAAGCGAATGATGAAAAGTAGAACAGTCGCAATAAGAGCATGGACTTAGTATTGCACTATCTCCTGGAGGTCAACAGAGGGCTa                                                                                                                                <  2:1944800/150‑1 (MQ=255)
                       ggAGATTATCGTAGGCTGACGACCTTTCAGTCTTCTGTATTAGTTGTGTTTACGAGTATTCCCTATTAAGCGAATGATGAAAAGTAGAACAGTCGCAATAAGAGCATGGACTTAGTATTGCACTATCTCCTGGAGGTCAACAGAGGGCTa                                                                                                                                <  2:2699010/150‑1 (MQ=255)
                        gagaTTATCGTAGGCTGACGACCTTTCAGTCTTCTGTATTAGTTGTGTTTACGAGTATTCCCTATTAAGCGAATGATGAAAAGTAGAACAGTCGCAATAAGAGCATGGACTTAGTATTGCACTATCTCCTGGAGGTCAACAGAGGGCTAt                                                                                                                               <  1:1738287/150‑1 (MQ=255)
                        gagaTTATCGTAGGCTGACGACCTTTCAGTCTTCTGTATTAGTTGTGTTTACGAGTATTCCCTATTAAGCGAATGATGAAAAGTAGAACAGTCGCAATAAGAGCATGGACTTAGTATTGCACTATCTCCTGGAGGTCAACAGAGGGCTAt                                                                                                                               <  2:1285891/150‑1 (MQ=255)
                        gagaTTATCGTAGGCTGACGACCTTTCAGTCTTCTGTATTAGTTGTGTTTACGAGTATTCCCTATTAAGCGAATGATGAAAAGTAGAACAGTCGCAATAAGAGCATGGACTTAGTATTGCACTATCTCCTGGAGGTCAACAGAGGGCTAt                                                                                                                               >  2:2063052/1‑150 (MQ=255)
                          gaTTATCGTAGGCTGACGACCTTTCAGTCTTCTGTATTAGTTGTGTTTACGAGTATTCCCTATTAAGCGAATGATGAAAAGTAGAACAGTCGCAATAAGAGCATGGACTTAGTATTGCACTATCTCCTGGAGGTCAACAGAGGGCTATTa                                                                                                                             <  1:3704484/150‑1 (MQ=255)
                              aTCGTAGGCTGACGACCTTTCAGTCTTCTGTATTAGTTGTGTTTACGAGTATTCCCTATTAAGCGAATGATGAAAAGTAGAACAGTCGCAATAAGAGCATGGACTTAGTATTGCACTATCTCCTGGAGGTCAACAGAGGGCTATTACTTg                                                                                                                         <  2:4868159/150‑1 (MQ=255)
                               tCGTAGGCTGACGACCTTTCAGTATTCTGTATTAGTTGTGTTTACGAGTATTCCCTATTAAGCGAATGATGAAAAGTAGAACAGTCGCAATAAGAGCATGGACTTAGTATTGCACTATCTCCTGGAGGTCAACAGAGGGCTATTACTTgc                                                                                                                        <  2:3552342/150‑1 (MQ=255)
                                cGTAGGCTGACGACCTGTCAGTCTTCTGTATTAGTTGTGTTTACGAGTATTCCCTATTAAGCGAATGATGAAAAGTAGAACAGTCGCAATAAGAGCATGGACTTAGTATTGCACTATCTCCTGGAGGTCAACAGAGGGCTATTACTTgcg                                                                                                                       <  2:2537310/150‑1 (MQ=255)
                                     gCTGACGACCTTTCAGTCTTCTGTATTAGTTGTGTTTACGAGTATTCCCTATTAAGCGAATGATGAAAAGTAGAACAGTCGCAATAAGAGCATGGACTTAGTATTGCACTATCTCCTGGAGGTCAACAGAGGGCTATTACTTGCGCAACa                                                                                                                  <  2:798330/150‑1 (MQ=255)
                                       tGACGACCTTTCAGTCTTCTGTATTAGTTGTGTTTACGAGTATTCCCTATTAAGCGAATGATGAAAAGTAGAACAGTCGCAATAAGAGCATGGACTTAGTATTGCACTATCTCCTGGAGGTCAACAGAGGGCTATTACTTGCGCAACAgg                                                                                                                <  2:4411654/150‑1 (MQ=255)
                                       tGACGACCTTTCAGTCTTCTGTATTAGTTGTGTTTACGAGTATTCCCTATTAAGCGAATGATGAAAAGTAGAACAGTCGCAATAAGAGCATGGACTTAGTATTGCACTATCTCCTGGAGGTCAACAGAGGGCTATTACTTGCGCAACAgg                                                                                                                >  2:1937552/1‑150 (MQ=255)
                                        gacgacCTTTCAGTCTTCTGTATTAGTTGTGTTTACGAGTATTCCCTATTAAGCGAATGATGAAAAGTAGAACAGTCGCAATAAGAGCATGGACTTAGTATTGCACTATCTCCTGGAGGTCAACAGAGGGCTATTACTTGCGCAACAGGt                                                                                                               >  1:2002117/1‑150 (MQ=255)
                                        gacgacCTTTCAGTCTTCTGTATTAGTTGTGTTTACGAGTATTCCCTATTAAGCGAATGATGAAAAGTAGAACAGTCGCAATAAGAGCATGGACTTAGTATTGCACTATCTCCTGGAGGTCAACAGAGGGCTATTACTTGCGCAACAGGt                                                                                                               <  1:4250238/150‑1 (MQ=255)
                                         acgacCTTTCAGTCTTCTGTATTAGTTGTGTTTACGAGTATTCCCTATTAAGCGAATGATGAAAAGTAGAACAGTCGCAATAAGAGCATGGACTTAGTATTGCACTATCTCCTGGAGGTCAACAGAGGGCTATTACTTGCGCAACAGGtt                                                                                                              <  2:2225999/150‑1 (MQ=255)
                                         acgacCTTTCAGTCTTCTGTATTAGTTGTGTTTACGAGTATTCCCTATTAAGCGAATGATGAAAAGTAGAACAGTCGCAATAAGAGCATGGACTTAGTATTGCACTATCTCCTGGAGGTCAACAGAGGGCTATTACTTGCGCAACAGGtt                                                                                                              >  1:4560673/1‑150 (MQ=255)
                                           gacCTTTCAGTCTTCTGTATTAGTTGTGTTTACGAGTATTCCCTATTAAGCGAATGATGAAAAGTAGAACAGTCGCAATAAGAGCATGGACTTAGTATTGCACTATCTCCTGGAGGTCAACAGAGGGCTATTACTTGCGCAACAGGTTaa                                                                                                            >  2:188247/1‑150 (MQ=255)
                                            acCTTTCAGTCTTCTGTATTAGTTGTGTTTACGAGTATTCCCTATTAAGCGAATGATGAAAAGTAGAACAGTCGCAATAAGAGCATGGACTTAGTATTGCACTATCTCCTGGAGGTCAACAGAGGGCTATTACTTGCGCAACAGGTTaaa                                                                                                           >  2:1223787/1‑150 (MQ=255)
                                            acCTTTCAGTCTTCTGTATTAGTTGTGTTTACGAGTATTCCCTATTAAGCGAATGATGAAAAGTAGAACAGTCGCAATAAGAGCATGGACTTAGTATTGCACTATCTCCTGGAGGTCAACAGAGGGCTATTACTTGCGCAACAGGTTaaa                                                                                                           <  2:4070429/150‑1 (MQ=255)
                                             ccTTTCAGTCTTCTGTATTAGTTGTGTTTACGAGTATTCCCTATTAAGCGAATGATGAAAAGTAGAACAGTCGCAATAAGAGCATGGACTTAGTATTGCACTATCTCCTGGAGGTCAACAGAGGGCTATTACTTGCGCAACAGGTTAAAg                                                                                                          <  2:400074/150‑1 (MQ=255)
                                                ttCAGTCTTCTGTATTAGTTGTGTTTACGAGTATTCCCTATTAAGCGAATGATGAAAAGTAGAACAGTCGCAATAAGAGCATGGACTTAGTATTGCACTATCTCCTGGAGGTCAACAGAGGGCTATTACTTGCGCAACAGGTTAAAGAtt                                                                                                       <  2:1725279/150‑1 (MQ=255)
                                                ttCAGTCTTCTGTATTAGTTGTGTTTACGAGTATTCCCTATTAAGCGAATGATGAAAAGTAGAACAGTCGCAATAAGAGCATGGACTTAGTATTGCACTATCTCCTGGAGGTCAACAGAGGGCTATTACTTGCGCAACAGGTTAAAGAtt                                                                                                       >  1:74213/1‑150 (MQ=255)
                                                 tCAGTCTTCTGTATTAGTTGTGTTTACGAGTATTCCCTATTAAGCGAATGATGAAAAGTAGAACAGTCGCAATAAGAGCATGGACTTAGTATTGCACTATCTCCTGGAGGTCAACAGAGGGCTATTACTTGCGCAACAGGTTAAAGATtg                                                                                                      <  1:3802971/150‑1 (MQ=255)
                                                 tCAGTCTTCTGTATTAGTTGTGTTTACGAGTATTCCCTATTAAGCGAATGATGAAAAGTAGAACAGTCGCAATAAGAGCATGGACTTAGTATTGCACTATATCCTGGAGGTCAACAGAGGGCTATTACTTGCGCAACAGGTTAAAGATtg                                                                                                      <  1:1364436/150‑1 (MQ=255)
                                                   aGTCTTCTGTATTAGTTGTGTTTACGAGTATTCCCTATTAAGCGAATGATGAAAAGTAGAACAGTCGCAATAAGAGCATGGACTTAGTATTGCACTATCTCCTGGAGGTCAACAGAGGGCTATTACTTGCGCAACAGGTTAAAGATtgtg                                                                                                    >  1:796364/1‑150 (MQ=255)
                                                      cttctGTATTAGTTGTGTTTACGAGTATTCCCTATTAAGCGAATGATGAAAAGTAGAACAGTCGCAATAAGAGCATGGACTTAGTATTGCACTATCTCCTGGAGGTCAACAGAGGGCTATTACTTGCGCAACAGGTTAAAGATTGTGAAt                                                                                                 <  2:3722378/150‑1 (MQ=255)
                                                      cttctGTATTAGTTGTGTTTACGAGTATTCCCTATTAAGCGAATGATGAAAAGTAGAACAGTCGCAATAAGAGCATGGACTTAGTATTGCACTATCTCCTGGAGGTCAACAGAGGGCTATTACTTGCGCAACAGGTTAAAGATTGTGAAt                                                                                                 <  1:3085849/150‑1 (MQ=255)
                                                        tctGTATTAGTTGTGTTTACGAGTATTCCCTATTAAGCGAATGATGAAAAGTAGAACAGTCGCAATAAGAGCATGGACTTAGTATTGCACTATCTCCTGGAGGTCAACAGAGGGCTATTACTTGCGCAACAGGTTAAAGATTGTGAATAg                                                                                               <  1:1072385/150‑1 (MQ=255)
                                                         ctGTATTAGTTGTGTTTACGAGTATTCCCTATTAAGCGAATGATGAAAAGTAGAACAGTCGCAATAAGAGCATGGACTTAGTATTGCACTATCTCCTGGAGGTCAACAGAGGGCTATTACTTGCGCAACAGGTTAAAGATTGTGAATAGt                                                                                              <  1:484/150‑1 (MQ=255)
                                                          tGTATTAGTTGTGTTTACGAGTATTCCCTATTAAGCGAATGATGAAAAGTAGAACAGTCGCAATAAGAGCATGGACTTAGTATTGCACTATATCCTGGAGGTCAACAGAGGGCTATTACTTGCGCAACAGGTTAAAGATTGTGAATAGtt                                                                                             >  1:3154934/1‑150 (MQ=255)
                                                           gTATTAGTTGTGTTTACGAGTATTCCCTATTAAGCGAATGATGAAAAGTAGAACAGTCGCAATAAGAGCATGGACTTAGTATTGCACTATCTCCTGGAGGTCAACAGAGGGCTATTACTTGCGCAACAGGTTAAAGATTGTGAATAGTTa                                                                                            <  2:4606525/150‑1 (MQ=255)
                                                               tAGTTGTGTTTACGAGTATTCCCTATTAAGCGAATGATGAAAAGTAGAACAGTCGCAATAAGAGCATGGACTTAGTATTGCACTATCTCCTGGAGGTCAACAGAGGGCTATTACTTGCGCAACAGGTTAAAGATTGTGAATAGTTACcag                                                                                        >  2:1404546/1‑150 (MQ=255)
                                                               tAGTTGTGTTTACGAGTATTCCCTATTAAGCGAATGATGAAAAGTAGAACAGTCGCAATAAGAGCATGGACTTAGTATTGCACTATCTCCTGGAGGTCAACAGAGGGCTATTACTTGCGCAACAGGTTAAAGATTGTGAATAGTTACcag                                                                                        <  2:1595206/150‑1 (MQ=255)
                                                                 gTTGTGTTTACGAGTATTCCCTATTAAGCGAATGATGAAAAGTAGAACAGTCGCAATAAGAGCATGGACTTAGTATTGCACTATCTCCTGGAGGTCAACAGAGGGCTATTACTTGCGCAACAGGTTAAAGATTGTGAATAGTTACcagca                                                                                      >  1:4790674/1‑150 (MQ=255)
                                                                 gTTGTGTTTACGAGTATTCCCTATTAAGCGAATGATGAAAAGTAGAACAGTCGCAATAAGAGCATGGACTTAGTATTGCACTATCTCCTGGAGGTCAACAGAGGGCTATTACTTGCGCAACAGGTTAAAGATTGTGAATAGTTACcagca                                                                                      >  2:4251795/1‑150 (MQ=255)
                                                                    gtgtTTACGAGTATTCCCTATTAAGCGAATGATGAAAAGTAGAACAGTCGCAATAAGAGCATGGACTTAGTATTGCACTATCTCCTGGAGGTCAACAGAGGGCTATTACTTGCGCAACAGGTTAAAGATTGTGAATAGTTACCAGCAGTc                                                                                   <  1:1495060/150‑1 (MQ=255)
                                                                    gtgtTTACGAGTATTCCCTATTAAGCGAATGATGAAAAGTAGAACAGTCGCAATAAGAGCATGGACTTAGTATTGCACTATCTCCTGGAGGTCAACAGAGGGCTATTACTTGCGCAACAGGTTAAAGATTGTGAATAGTTACCAGCAGTc                                                                                   <  1:1290851/150‑1 (MQ=255)
                                                                     tgtTTACGAGTATTCCCTATTAAGCGAATGATGAAAAGTAGAACAGTCGCAATAAGAGCATGGACTTAGTATTGCACTATCTCCTGGAGGTCAACAGAGGGCTATTACTTGCGCAACAGGTTAAAGATTGTGAATAGTTACCAGCAGTCa                                                                                  <  1:4449144/150‑1 (MQ=255)
                                                                          aCGAGTATTCCCTATTAAGCGAATGATGAAAAGTAGAACAGTCGCAATAAGAGCATGGACTTAGTATTGCACTATCTCCTGGAGGTCAACAGAGGGCTATTACTTGCGCAACAGGTTAAAGATTGTGAATAGTTACCAGCAGTCATTTAc                                                                             <  1:841556/150‑1 (MQ=255)
                                                                             agtaTTCCCTATTAAGCGAATGATGAAAAGTAGAACAGTCGCAATAAGAGCATGGACTTAGTATTGCACTATCTCCTGGAGGTCAACAGAGGGCTATTACTTGCGCAACAGGTTAAAGATTGTGAATAGTTACCAGCAGTCATTTACCCg                                                                          <  2:1072177/147‑1 (MQ=255)
                                                                             agtaTTCCCTATTAAGCGAATGATGAAAAGTAGAACAGTCGCAATAAGAGCATGGACTTAGTATTGCACTATCTCCTGGAGGTCAACAGAGGGCTATTACTTGCGCAACAGGTTAAAGATTGTGAATAGTTACCAGCAGTCATTTACCCg                                                                          <  2:3235899/147‑1 (MQ=255)
                                                                             agaATTCCCTATTAAGCGAATGATGAAAAGTAGAACAGTCGCAATAAGAGCATGGACTTAGTATTGCACTATCTCCTGGAGGTCAACAGAGGGCTATTACTTGCGCAACAGGTTAAAGATTGTGAATAGTTACCAGCAGTCATTTACCCg                                                                          >  2:1019249/1‑150 (MQ=255)
                                                                                  tCCCTATTAAGCGAATGATGAAAAGTAGAACAGTCGCAATAAGAGCATGGACTTAGTATTGCACTATCTCCTGGAGGTCAACAGAGGGCTATTACTTGCGCAACAGGTTAAAGATTGTGAATAGTTACCAGCAGTCATTTACCCGCTtat                                                                     <  2:3646335/150‑1 (MQ=255)
                                                                                  tCCCTATTAAGCGAATGATGAAAAGTAGAACAGTCGCAATAAGAGCATGGACTTAGTATTGCACTATCTCCTGGAGGTCAACAGAGGGCTATTACTTGCGCAACAGGTTAAAGATTGTGAATAGTTACCAGCAGTCATTTACCCGCTtat                                                                     <  2:5144360/150‑1 (MQ=255)
                                                                                         tAAGCGTATGATGAAAAGTAGAACAGTCGCAATAAGAGCATGGACTTAGTATTGCACTATCTCCTGGAGGTCAACAGAGGGCTTTTACTTGCGCAACAGGTTAAAGATTGTGAATAGTTACCAGCAGTCATTTACCCGCTTATAACAAGc                                                              <  1:3760673/150‑1 (MQ=255)
                                                                                         tAAGCGAATGATGAAAAGTAGAACAGTCGCAATAAGAGCATGGACTTAGTATTGCACTATCTCCTGGAGGTCAACAGAGGGCTATTACTTGCGCAACAGGTTAAAGATTGTGAATAGTTACCAGCAGTCATTTACCCGCTTATAACAAGc                                                              <  1:2162865/150‑1 (MQ=255)
                                                                                         tAAGCGAATGATGAAAAGTAGAACAGTCGCAATAAGAGCATGGACTTAGTATTGCACTATCTCCTGGAGGTCAACAGAGGGCTATTACTTGCGCAACAGGTTAAAGATTGTGAATAGTTACCAGCAGTCATTTACCCGCTTATAACAAGc                                                              >  1:2674467/1‑150 (MQ=255)
                                                                                             cGAATGATGAAAAGTAGAACAGTCGCAATAAGAGCATGGACTTAGTATTGCACTATCTCCTGGAGGTCAACAGAGGGCTATTACTTGCGCAACAGGTTAAAGATTGTGAATAGTTACCAGCAGTCATTTACCCGCTTATAACAAGCGAgg                                                          >  2:5058595/1‑150 (MQ=255)
                                                                                             cGAATGATGAAAAGTAGAACAGTCGCAATAAGAGCATGGACTTAGTATTGCACTATCTCCTGGAGGTCAACAGAGGGCTATTACTTGCGCAACAGGTTAAAGATTGTGAATAGTTACCAGCAGTCATATACCCGCTTATAACAAGCGAgg                                                          >  2:2739165/1‑150 (MQ=255)
                                                                                                atgatgAAAAGTAGAACAGTCGCAATAAGAGCATGGACTTAGTATTGCACTATCTCCTGGAGGTCAACAGAGGGCTATTACTTGCGCAACAGGTTAAAGATTGTGAATAGTTACCAGCAGTCATTTACCCGCTTATAACAAGCGAGGCAg                                                       >  1:883947/1‑150 (MQ=255)
                                                                                                 tgatgaAAAGTAGAACAGTCGCAATAAGAGCATGGACTTAGTATTGCACTATCTCCTGGAGGTCAACAGAGGGCTATTACTTGCGCAACAGGTTAAAGATTGTGAATAGTTACCAGCAGTCATTTACCCGCTTATAACAAGCGAGGCAGt                                                      >  2:45166/1‑150 (MQ=255)
                                                                                                 tgatgaAAAGTAGAACAGTCGCAATAAGAGCATGGACTTAGTATTGCACTATCTCCTGGAGGTCAACAGAGGGCTATTACTTGCGCAACAGGTTAAAGATTGTGAATAGTTACCAGCAGTCATTTACCCGCTTATAACAAGCGAGGCAGt                                                      >  1:554390/1‑150 (MQ=255)
                                                                                                  gatgaAAAGTAGAACAGTCGCAATAAGAGCATGGACTTAGTATTGCACTATCTCCTGGAGGTCAACAGAGGGCTATTACTTGCGCAACAGGTTAAAGATTGTGAATAGTTACCAGCAGTCATTTACCCGCTTATAACAAGCGAGGCAGtt                                                     >  2:3880077/1‑150 (MQ=255)
                                                                                                   atgaAAAGTAGAACAGTCGCAATAAGAGCATGGACTTAGTATTGCACTATCTCCTGGAGGTCAACAGAGGGCTATTACTTGCGCAACAGGTTAAAGATTGTGAATAGTTACCAGCAGTCATTTACCCGCTTATAACAAGCGAGGCAGTTg                                                    >  2:3398670/1‑150 (MQ=255)
                                                                                                        aaGTAGAACAGTCGCAATAAGAGCATGGACTTAGTATTGCACTATCTCCTGGAGGTCAACAGAGGGCTATTACTTGCGCAACAGGTTAAAGATTGTGAATAGTTACCAGCAGTCATTTACCCGCTTATAACAAGCGAGGCAGTTGTAATg                                               <  1:2955034/150‑1 (MQ=255)
                                                                                                        aaGTAGAACAGTAGCAATAAGAGCATGGACTTAGTATTGCACTATCTCCTGGAGGTCAACAGAGGGCTATTACTTGCGCAACAGGTTAAAGATTGTGAATAGTTACCAGCAGTCATTTACCCGCTTATAACAAGCGAGGCAGTTGTAATg                                               <  1:3922531/150‑1 (MQ=255)
                                                                                                         aGTAGAACAGTCGCAATAAGAGCATGGACTTAGTATTGCACTATCTCCTGGAGGTCAACAGAGGGCTATTACTTGCGCAACAGGTTAAAGATTGTGAATAGTTACCAGCAGTCATTTACCCGCTTATAACAAGCGAGGCAGTTGTAATGa                                              >  1:2780475/1‑150 (MQ=255)
                                                                                                          gTAGAACAGTCGCAATAAGAGCATGGACTTAGTATTGCACTATCTCCTGGAGGTCAACAGAGGGCTATTACTTGCGCAACAGGTTAAAGATTGTGAATAGTTACCAGCAGTCATTTACCCGCTTATAACAAGCGAGGCAGTTGTAATGAt                                             <  2:3304108/150‑1 (MQ=255)
                                                                                                          gTAGAACAGTCGCAATAAGAGCATGGACTTAGTATTGCACTATCTCCTGGAGGTCAACAGAGGGCTATTACTTGCGCAACAGGTTAAAGATTGTGAATAGTTACCAGCAGTCATTTACCCGCTTATAACAAGCGAGGCAGTTGTAATGAt                                             <  1:2078605/150‑1 (MQ=255)
                                                                                                            aGAACAGTCGCAATAAGAGCATGGACTTAGTATTGCACTATCTCCTGGAGGTCAACAGAGGGCTATTACTTGCGCAACAGGTTAAAGATTGTGAATAGTTACCAGCAGTCATTTACCCGCTTATAACAAGCGAGGCAGTTGTAATGATAg                                           >  2:710063/1‑150 (MQ=255)
                                                                                                                 aGTCGCAATAAGAGCATGGACTTAGTATTGCACTATCTCCTGGAGGTCAACAGAGGGCTATTACTTGCGCAACAGGTTAAAGATTGTGAATAGTTACCAGCAGTCATTTACCCGCTTATAACAAGCGAGGCAGTTGTAATGATAGCTCAg                                      <  2:4616205/150‑1 (MQ=255)
                                                                                                                  gTCGCAATAAGAGCATGGACTTAGTATTGCACTATCTCCTGGAGGTCAACAGAGGGCTATTACTTGCGCAACAGGTTAAAGATTGTGAATAGTTACCAGCAGTCATTTACCCGCTTATAACAAGCGAGGCAGTTGTAATGATAGCTCAGa                                     >  1:2203262/1‑150 (MQ=255)
                                                                                                                     gCAATAAGAGCATGGACTTAGTATTGCACTATCTCCTGGAGGTCAACAGAGGGCTATTACTTGCGCAACAGGTTAAAGATTGTGAATAGTTACCAGCAGTCATTTACCCGCTTATAACAAGCGAGGCAGTTGTAATGATAGCTCAGAAgg                                  <  1:4064110/150‑1 (MQ=255)
                                                                                                                      cAATAAGAGCATGGACTTAGTATTGCACTATCTCCTGGAGGTCAACAGAGGGCTATTACTTGCGCAACAGGTTAAAGATTGTGAATAGTTACCAGCAGTCATTTACCCGCTTATAACAAGCGAGGCAGTTGTAATGATAGCTCAGAAGGa                                 >  1:238662/1‑150 (MQ=255)
                                                                                                                       aaTAAGAGCATGGACTTAGTATTGCACTATCTCCTGGAGGTCAACAGAGGGCTATTACTTGCGCAACAGGTTAAAGATTGTGAATAGTTACCAGCAGTCATTTACCCGCTTATAACAAGCGAGGCAGTTGTAATGATAGCTCAGAAGGAt                                <  2:482285/150‑1 (MQ=255)
                                                                                                                        aTAAGAGCATGGACTTAGTATTGCACTATCTCCTGGAGGTCAACAGAGGGCTATTACTTGCGCAACAGGTTAAAGATTGTGAATAGTTACCAGCAGTCATTTACCCGCTTATAACAAGCGAGGCAGTTGTAATGATAGCTCAGAAGGAtt                               <  2:74213/150‑1 (MQ=255)
                                                                                                                         tAAGAGCATGGACTTAGTATTGCACTATCTCCTGGAGGTCAACAGAGGGCTATTACTTGCGCAACAGGTTAAAGATTGTGAATAGTTACCAGCAGTCATTTACCCGCTTATAACAAGCGAGGCAGTTGTAATGATAGCTCAGAAGGATTa                              >  2:4375215/1‑150 (MQ=255)
                                                                                                                             agCATGGACTTAGTATTGCACTATCTCCTGGAGGTCAACAGAGGGCTATTACTTGCGCAACAGGTTAAAGATTGTGAATAGTTACCAGCAGTCATTTACCCGCTTATAACAAGCGAGGCAGTTGTAATGATAGCTCAGAAGGATTATGCa                          >  1:3327960/1‑150 (MQ=255)
                                                                                                                              gCATGGACTTAGTATTGCACTATCTCCTGGAGGTCAACAGAGGGCTATTACTTGCGCAACAGGTTAAAGATTGTGAATAGTTACCAGCAGTCATTTACCCGCTTATAACAAGCGAGGCAGTTGTAATGATAGCTCAGAAGGATTATGCaa                         <  2:1924221/150‑1 (MQ=255)
                                                                                                                                    aCTTAGTATTGCACTATCTCCTGGAGGTCAACAGAGGGCTATTACTTGCGCAACAGGTTAAAGATTGTGAATAGTTACCAGCAGTCATTTACCCGCTTATAACAAGCGAGGCAGTTGTAATGATAGCTCAGAAGGATTATGCAAGGCTTc                   <  1:952998/150‑1 (MQ=255)
                                                                                                                                    aCTTAGTATTGCACTATCTCCTGGAGGTCAACAGAGGGCTATTACTTGCGCAACAGGTTAAAGATTGTGAATAGTTACCAGCAGTCATTTACCCGCTTATAACAAGCGAGGCAGTTGTAATGATAGCTCAGAAGGATTATGCAAGGCTTc                   <  1:1927706/150‑1 (MQ=255)
                                                                                                                                    aCTGAGTATTGCACTATCTCCTGGAGGTCAACAGAGGGCTATTACTTGCGCAACAGGTTAAAGATTGTGAATAGTTACCAGCAGTCATTTACCCGCTTATAACAAGCGAGGCAGTTGTAATGATAGCTCAGAAGGATTATGCAAGGCTTc                   <  2:3915099/150‑1 (MQ=255)
                                                                                                                                         gTATTGCACTATCTCCTGGAGGTCAACAGAGGGATATTACTTGCGCAACAGGTTAAAGATTGTGAATAGTTACCAGCAGTCATTTACCCGCTTATAACAAGCGAGGCAGTTGTAATGATAGCTCAGAAGGATTATGCAAGGCTTCGTAAg              <  2:1553701/150‑1 (MQ=255)
                                                                                                                                              gCACTATCTCCTGGAGGTCAACAGAGGGCTATTACTTGCGCAACAGGTTAAAGATTGTGAATAGTTACCAGCAGTCATTTACCCGCTTATAACAAGCGAGGCAGTTGTAATGATAGCTCAGAAGGATTATGCAAGGCTTCGTAAGGgaga         >  1:149858/1‑150 (MQ=255)
                                                                                                                                               cACTATCTCCTGGAGGTCAACAGAGGGCTATTACTTGCGCAACAGGTTAAAGATTGTGAATAGTTACCAGCAGTCATTTACCCGCTTATAACAAGCGAGGCAGTTGTAATGATAGCTCAGAAGGATTATGCAAGGCTTCGTAAGGGAGaa        >  2:651379/1‑150 (MQ=255)
                                                                                                                                                aCTATCTCCTGGAGGTCAACAGAGGGCTATTACTTGCGCAACAGGTTAAAGATTGTGAATAGTTACCAGCAGTCATTTACCCGCTTATAACAAGCGAGGCAGTTGTAATGATAGCTCAGAAGGATTATGCAAGGCTTCGTAAGGGAGAAc       <  1:2656556/150‑1 (MQ=255)
                                                                                                                                                 cTATCTCCTGGAGGTCAACAGAGGGCTATTACTTGCGCAACAGGTTAAAGATTGTGAATAGTTACCAGCAGTCATTTACCCGCTTATAACAAGCGAGGCAGTTGTAATGATAGCTCAGAAGGATTATGCAAGGCTTCGTAAGGGAGAACg      >  1:1665265/1‑150 (MQ=255)
                                                                                                                                                  tATCTCCTGGAGGTCAACAGAGGGCTATTACTTGCGCAACAGGTTAAAGATTGTGAATAGTTACCAGCAGTCATTTACCCGCTTATAACAAGCGAGGCAGTTGTAATGATAGCTCAGAAGGATTATGCAAGGCTTCGTAAGGGAGAACGc     >  1:5197026/1‑150 (MQ=255)
                                                                                                                                                   aTCTCCTGGAGGTCAACAGAGGGCTATTACTTGCGCAACAGGTTAAAGATTGTGAATAGTTACCAGCAGTCATTTACCCGCTTATAACAAGCGAGGCAGTTGTAATGATAGCTCAGAAGGATTATGCAAGGCTTCGTAAGGGAGAACGCa    <  2:5968/150‑1 (MQ=255)
                                                                                                                                                    tctcCTGGAGGTCAACAGAGGGCTATTACTTGCGCAACAGGTTAAAGATTGTGAATAGTTACCAGCAGTCATTTACCCGCTTATAACAAGCGAGGCAGTTGTAATGATAGCTCAGAAGGATTATGCAAGGCTTCGTAAGGGAGAACGCat   >  1:5051851/1‑150 (MQ=255)
                                                                                                                                                    tctcCTGGAGGTCAACAGAGGGCTATTACTTGCGCAACAGGTTAAAGATTGTGAATAGTTACCAGCAGTCATTTACCCGCTTATAACAAGCGAGGCAGTTGTAATGATAGCTCAGAAGGATTATGCAAGGCTTCGTAAGGGAGAACGCat   >  1:4164928/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CTGATTCATGTCACATTTTATGGGGAGATTATCGTAGGCTGACGACCTTTCAGTCTTCTGTATTAGTTGTGTTTACGAGAATTCCCTATTAAGCGAATGATGAAAAGTAGAACAGTCGCAATAAGAGCATGGACTTAGTATTGCACTATCTCCTGGAGGTCAACAGAGGGCTATTACTTGCGCAACAGGTTAAAGATTGTGAATAGTTACCAGCAGTCATTTACCCGCTTATAACAAGCGAGGCAGTTGTAATGATAGCTCAGAAGGATTATGCAAGGCTTCGTAAGGGAGAACGCATA  >  CP009273/1288590‑1288888

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: