Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,969,710 G→T 100% G25C (GGT→TGT)  wzxE → O‑antigen translocase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,969,7100GT99.6% 730.6 / ‑4.3 220G25C (GGT→TGT) wzxEO‑antigen translocase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (112/107);  total (112/108)
Fisher's exact test for biased strand distribution p-value = 4.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 220 total aligned reads displayed.

CGCTGTTCTACAACCTGTCGCCCGTCAATCAGCGTACGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATC  >  CP009273/3969566‑3969854
                                                                                                                                                |                                                                                                                                                
cGCTGTTCTACAACCTGTCGCCCGTCAATCAGCGTACGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTa                                                                                                                                               <  1:616453/148‑1 (MQ=255)
cGCTGTTCTACAACCTGTCGCCCGTCAATCAGCGTACGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTa                                                                                                                                               <  1:272952/148‑1 (MQ=255)
 gCTGTTCTACAACCTGTCGCCCGTCAATCAGCGTACGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTaa                                                                                                                                              >  2:1205202/1‑148 (MQ=255)
 gCTGTTCTACAACCTGTCGCCCGTCAATCAGCGTACGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTaa                                                                                                                                              <  2:1393678/148‑1 (MQ=255)
  cTGTTCTACAACCTGTCGCCCGTCAATCAGCGTACGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAATAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAg                                                                                                                                             >  1:958318/1‑148 (MQ=255)
  cTGTTCTACAACCTGTCGCCCGTCAATCAGCGTACGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAg                                                                                                                                             <  1:2647315/148‑1 (MQ=255)
     ttCTACAACCTGTCGCCCGTCAATCAGCGTACGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTg                                                                                                                                          <  2:2981553/148‑1 (MQ=255)
          cAACCTGTCGCCCGTCAATCAGCGTACGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGc                                                                                                                                     >  2:180542/1‑148 (MQ=255)
          cAACCTGTCGCCCGTCAATCAGCGTACGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGc                                                                                                                                     <  2:1948303/148‑1 (MQ=255)
           aaCCTGTCGCCCGTCAATCAGCGTACGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCg                                                                                                                                    >  1:1749651/1‑148 (MQ=255)
             ccTGTCGCCCGTCAATCAGCGTACGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGgt                                                                                                                                  >  1:1082807/1‑148 (MQ=255)
             ccTGTCGCCCGTCAATCAGCGTACGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGgt                                                                                                                                  >  2:2011573/1‑148 (MQ=255)
                  cGCCCGTCAATCAGCGTACGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCAt                                                                                                                             <  2:1496653/148‑1 (MQ=255)
                           aTCAGCGTACGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCcgg                                                                                                                    >  2:2237705/1‑148 (MQ=255)
                                 gtccGGTAATTGCGACTTTGTTGAACGACATTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGc                                                                                                              <  2:1781598/145‑1 (MQ=255)
                                 gTACGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGc                                                                                                              <  2:549825/148‑1 (MQ=255)
                                  tACGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCt                                                                                                             <  1:1230005/148‑1 (MQ=255)
                                   aCGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCtt                                                                                                            <  2:3240617/148‑1 (MQ=255)
                                   aCGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCtt                                                                                                            <  1:2087452/148‑1 (MQ=255)
                                      gTAATTGAGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAATAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTggg                                                                                                         <  1:1690485/148‑1 (MQ=255)
                                       tAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGc                                                                                                        >  2:1224746/1‑148 (MQ=255)
                                       tAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGc                                                                                                        <  1:2094065/148‑1 (MQ=255)
                                              gACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTggcggc                                                                                                 <  1:1869565/148‑1 (MQ=255)
                                               aCTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCa                                                                                                >  1:434062/1‑148 (MQ=255)
                                               aCTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCa                                                                                                >  1:340800/1‑148 (MQ=255)
                                               aCTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCa                                                                                                >  1:1529343/1‑148 (MQ=255)
                                                 tttGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTAGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCaaa                                                                                              >  1:1380788/1‑148 (MQ=255)
                                                    gttgAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAAttt                                                                                           >  1:1556806/1‑148 (MQ=255)
                                                         actactTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGcc                                                                                      <  1:2560442/148‑1 (MQ=255)
                                                         actactTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGcc                                                                                      >  1:1596892/1‑148 (MQ=255)
                                                          ctactTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCa                                                                                     >  2:2512031/1‑148 (MQ=255)
                                                          ctactTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCa                                                                                     >  1:2661167/1‑148 (MQ=255)
                                                          ctactTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCa                                                                                     >  1:1006005/1‑148 (MQ=255)
                                                           tactTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTTGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAg                                                                                    >  2:586812/1‑148 (MQ=255)
                                                           tactTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAg                                                                                    >  1:2439942/1‑148 (MQ=255)
                                                           tactTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAg                                                                                    <  1:486557/148‑1 (MQ=255)
                                                           tactTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAg                                                                                    >  2:1745939/1‑148 (MQ=255)
                                                            actTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGt                                                                                   >  1:2826539/1‑148 (MQ=255)
                                                              ttttCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTg                                                                                 <  1:1635320/148‑1 (MQ=255)
                                                                ttCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGAt                                                                               <  2:1749651/148‑1 (MQ=255)
                                                                ttCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGAt                                                                               <  2:731547/148‑1 (MQ=255)
                                                                ttCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGAt                                                                               <  2:3101149/148‑1 (MQ=255)
                                                                 tCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGAtt                                                                              <  2:1855984/148‑1 (MQ=255)
                                                                 tCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGAtt                                                                              <  2:1262552/148‑1 (MQ=255)
                                                                    tGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTAcc                                                                           >  1:1338584/1‑148 (MQ=255)
                                                                    tGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTAcc                                                                           >  1:1338454/1‑148 (MQ=255)
                                                                     gATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCg                                                                          >  1:1265724/1‑148 (MQ=255)
                                                                     gATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCg                                                                          <  2:3050436/148‑1 (MQ=255)
                                                                       tatGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTg                                                                        <  2:2496149/148‑1 (MQ=255)
                                                                       tatGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTg                                                                        <  2:3216547/148‑1 (MQ=255)
                                                                       tatGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTg                                                                        >  2:1120600/1‑148 (MQ=255)
                                                                        atGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGc                                                                       >  1:2437280/1‑148 (MQ=255)
                                                                          gTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTc                                                                     <  2:2613885/148‑1 (MQ=255)
                                                                          gTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTc                                                                     >  2:908869/1‑148 (MQ=255)
                                                                          gTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTc                                                                     >  1:778823/1‑148 (MQ=255)
                                                                            cGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCgg                                                                   <  2:1827750/148‑1 (MQ=255)
                                                                                ggCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTg                                                               <  2:881474/148‑1 (MQ=255)
                                                                                  cAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCt                                                             >  1:1255505/1‑148 (MQ=255)
                                                                                  cAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCt                                                             <  1:3225057/148‑1 (MQ=255)
                                                                                    aaaGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTg                                                           >  2:485214/1‑148 (MQ=255)
                                                                                    aaaGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTg                                                           >  2:2359308/1‑148 (MQ=255)
                                                                                       gCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCg                                                        >  2:1542083/1‑148 (MQ=255)
                                                                                           ccTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGc                                                    <  1:884143/148‑1 (MQ=255)
                                                                                           ccTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGc                                                    <  1:884787/148‑1 (MQ=255)
                                                                                           ccTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGc                                                    >  2:1933920/1‑148 (MQ=255)
                                                                                           ccTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTAACGAGCTCGGCGTGCTTGCCGGGGc                                                    >  1:301302/1‑148 (MQ=255)
                                                                                            cTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCt                                                   <  2:1571268/148‑1 (MQ=255)
                                                                                                 ggACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCAt                                              <  1:2143154/148‑1 (MQ=255)
                                                                                                 ggACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCAt                                              <  2:605291/148‑1 (MQ=255)
                                                                                                  gACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATc                                             <  1:1914299/148‑1 (MQ=255)
                                                                                                  gACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATc                                             <  1:3056053/148‑1 (MQ=255)
                                                                                                  gACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATc                                             <  2:430889/148‑1 (MQ=255)
                                                                                                  gACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATGACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATc                                             <  2:1706303/148‑1 (MQ=255)
                                                                                                       cggcCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTaa                                        <  2:1596575/148‑1 (MQ=255)
                                                                                                         gcCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACg                                      <  1:2545847/148‑1 (MQ=255)
                                                                                                          ccAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACgg                                     >  2:2090123/1‑148 (MQ=255)
                                                                                                           cAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGgt                                    <  2:848936/148‑1 (MQ=255)
                                                                                                            aGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGgtg                                   <  1:580117/148‑1 (MQ=255)
                                                                                                              tACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTa                                 <  2:2903937/148‑1 (MQ=255)
                                                                                                                cacTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAAc                               >  1:1767243/1‑148 (MQ=255)
                                                                                                                      tCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATa                         >  2:2593874/1‑148 (MQ=255)
                                                                                                                         aGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGt                      <  1:1574938/148‑1 (MQ=255)
                                                                                                                         aGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGt                      >  2:3208149/1‑148 (MQ=255)
                                                                                                                         aGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGt                      >  2:3286378/1‑148 (MQ=255)
                                                                                                                         aGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGt                      >  1:2750879/1‑148 (MQ=255)
                                                                                                                           aTTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTg                    <  1:2233524/148‑1 (MQ=255)
                                                                                                                           aTTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTg                    <  1:2499748/148‑1 (MQ=255)
                                                                                                                            ttGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGc                   >  1:2983123/1‑148 (MQ=255)
                                                                                                                               gTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCa                >  1:3163321/1‑148 (MQ=255)
                                                                                                                                tGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAg               >  1:980822/1‑148 (MQ=255)
                                                                                                                                 gCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGt              <  2:2204004/148‑1 (MQ=255)
                                                                                                                                 gCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGt              <  2:841797/148‑1 (MQ=255)
                                                                                                                                  ccGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTa             <  2:174022/148‑1 (MQ=255)
                                                                                                                                   cGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTAc            >  2:3025925/1‑148 (MQ=255)
                                                                                                                                     ggTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCa          >  2:1055941/1‑148 (MQ=255)
                                                                                                                                     ggTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCa          >  2:90197/1‑148 (MQ=255)
                                                                                                                                         aCTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGAt      >  1:2400213/1‑148 (MQ=255)
                                                                                                                                          cTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATa     >  1:1506675/1‑148 (MQ=255)
                                                                                                                                           tGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATaa    >  1:1996492/1‑148 (MQ=255)
                                                                                                                                            ggTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAAt   <  2:2618795/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                
CGCTGTTCTACAACCTGTCGCCCGTCAATCAGCGTACGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATC  >  CP009273/3969566‑3969854

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: