Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,225,841 A→C 100% K3T (AAA→ACA)  yjbE → extracellular polysaccharide production threonine‑rich protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,225,8410AC96.7% 942.2 / 6.6 239K3T (AAA→ACA) yjbEextracellular polysaccharide production threonine‑rich protein
Reads supporting (aligned to +/- strand):  ref base A (3/5);  new base C (139/92);  total (142/97)
Fisher's exact test for biased strand distribution p-value = 2.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.37e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 244 total aligned reads displayed.

TTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCGCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAAAAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAG  >  CP009273/4225696‑4225985
                                                                                                                                                 |                                                                                                                                                
  gATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGtt                                                                                                                                              <  2:252470/148‑1 (MQ=255)
   atatTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTc                                                                                                                                             >  2:1031958/1‑148 (MQ=255)
   atatTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTc                                                                                                                                             <  2:1858208/148‑1 (MQ=255)
   atatTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTc                                                                                                                                             >  2:2812028/1‑148 (MQ=255)
   atatTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTc                                                                                                                                             <  2:3284137/148‑1 (MQ=255)
   atatTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTc                                                                                                                                             >  1:1904957/1‑148 (MQ=255)
    tatTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCTCGCCGTAATTCGCATGCCTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCt                                                                                                                                            >  1:1815589/1‑148 (MQ=255)
       tcattcCAATTATCCTAAAACGCCATCGCTAATTCCCCTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTa                                                                                                                                         >  2:2184330/1‑148 (MQ=255)
        cattcCAATTATCCTAAAACGCCATCGCTAATTCCCCTCGCCGTAATTCGCATGATTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTAt                                                                                                                                        >  2:2164054/1‑148 (MQ=255)
            ccAATTATCCTAAAACGCCATCGCTAATTCCCCTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCa                                                                                                                                    <  1:642997/148‑1 (MQ=255)
            ccAATTATCCTAAAACGCCATCGCTAATTCCCCTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCa                                                                                                                                    <  2:2278107/148‑1 (MQ=255)
              aaTTATCCTAAAACGCCATCGCTAATTCCCCTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCAtt                                                                                                                                  <  2:1260543/148‑1 (MQ=255)
                  aTCCTAAAACGCCATCGCTAATTCCCCTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTg                                                                                                                              >  2:799186/1‑148 (MQ=255)
                   tCCTAAAACGCCATCGCTAATTCCCCTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGc                                                                                                                             <  2:429239/148‑1 (MQ=255)
                    ccTAAAACGCCATCGCTAATTCCCCTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGcc                                                                                                                            >  2:3308508/1‑148 (MQ=255)
                    ccTAAAACGCCATCGCTAATTCCCCTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGcc                                                                                                                            <  2:395959/148‑1 (MQ=255)
                        aaaCGCCATCGCTAATTCCCCTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCatat                                                                                                                        <  2:3330886/148‑1 (MQ=255)
                        aaaCGCCATCGCTAATTCCCCTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCatat                                                                                                                        <  2:2555659/148‑1 (MQ=255)
                         aaCGCCATCGCTAATTCCCCTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATc                                                                                                                       >  2:1760212/1‑148 (MQ=255)
                             ccATCGCTAATTCCCCTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTgcg                                                                                                                   <  1:2214128/148‑1 (MQ=255)
                              cATCGCTAATTCCCCTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTgcgc                                                                                                                  <  1:3421282/148‑1 (MQ=255)
                              cATCGCTAATTCCCCTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTgcgc                                                                                                                  <  1:618257/148‑1 (MQ=255)
                              cATCGCTAATTCCCCTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTgcgc                                                                                                                  >  1:2518136/1‑148 (MQ=255)
                              cATCGCTAATTCCCCTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTgcgc                                                                                                                  >  1:240696/1‑148 (MQ=255)
                              cATCGCTAATTCCCCTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTgcgc                                                                                                                  <  1:618345/148‑1 (MQ=255)
                                tCGCTAATTCCCCTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCtt                                                                                                                >  2:860653/1‑148 (MQ=255)
                                 cGCTAATTCCCCTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTg                                                                                                               >  2:1470286/1‑148 (MQ=255)
                                 cGCTAATTCCCCTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTg                                                                                                               >  1:3122036/1‑148 (MQ=255)
                                  gCTAATTCCCCTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTgc                                                                                                              <  2:28915/148‑1 (MQ=255)
                                       ttCCCCTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGa                                                                                                         >  1:1737361/1‑148 (MQ=255)
                                       ttCCCCTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGa                                                                                                         >  2:1070715/1‑148 (MQ=255)
                                       ttCCCCTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGa                                                                                                         >  1:3460103/1‑148 (MQ=255)
                                       ttCCCCGCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAAAAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGa                                                                                                         <  1:374694/148‑1 (MQ=255)
                                        tCCCCTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGAc                                                                                                        >  1:1596495/1‑148 (MQ=255)
                                          cccGCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAAAAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACtt                                                                                                      <  1:3339821/148‑1 (MQ=255)
                                           cctcgcCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTc                                                                                                     <  1:1313220/145‑1 (MQ=255)
                                             tcgcCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTg                                                                                                   <  2:96258/147‑1 (MQ=255)
                                               gcCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCg                                                                                                 <  1:72210/148‑1 (MQ=255)
                                                   tAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTggg                                                                                             <  1:1006083/148‑1 (MQ=255)
                                                       tCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGc                                                                                         >  1:1700271/1‑148 (MQ=255)
                                                       tCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAAAAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGc                                                                                         >  1:2013058/1‑148 (MQ=255)
                                                         gCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCAc                                                                                       <  1:1518535/148‑1 (MQ=255)
                                                          cATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCAcc                                                                                      >  1:1068971/1‑148 (MQ=255)
                                                           aTGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCt                                                                                     <  2:3038242/148‑1 (MQ=255)
                                                             gCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGt                                                                                   >  1:3505105/1‑148 (MQ=255)
                                                             gCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGt                                                                                   >  2:2456483/1‑148 (MQ=255)
                                                               tttAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTAc                                                                                 <  1:2946177/148‑1 (MQ=255)
                                                               tttAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTAc                                                                                 >  2:3444806/1‑148 (MQ=255)
                                                               tttAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTAc                                                                                 <  1:1732710/148‑1 (MQ=255)
                                                                ttAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACa                                                                                >  2:1440687/1‑148 (MQ=255)
                                                                 tAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAg                                                                               >  1:3338679/1‑148 (MQ=255)
                                                                  aGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAAAAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAgg                                                                              <  1:1068728/148‑1 (MQ=255)
                                                                   gTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGt                                                                             <  2:1530042/148‑1 (MQ=255)
                                                                    ttGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTg                                                                            <  2:1973039/148‑1 (MQ=255)
                                                                      gtgtATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTggg                                                                          >  2:3475647/1‑148 (MQ=255)
                                                                      gtgtATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTggg                                                                          <  2:1140394/148‑1 (MQ=255)
                                                                       tgtATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGc                                                                         >  2:53876/1‑148 (MQ=255)
                                                                       tgtATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGc                                                                         >  2:54917/1‑148 (MQ=255)
                                                                          ataCTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGaa                                                                      >  2:461963/1‑148 (MQ=255)
                                                                           taCTCGATCCCGCCCGATATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAg                                                                     <  1:49238/148‑1 (MQ=255)
                                                                            aCTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGc                                                                    >  2:657547/1‑148 (MQ=255)
                                                                              tCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCgg                                                                  <  2:3316607/148‑1 (MQ=255)
                                                                                gATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCa                                                                >  1:2426534/1‑148 (MQ=255)
                                                                                 aTCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAg                                                               >  2:3155009/1‑148 (MQ=255)
                                                                                  tCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAgg                                                              <  2:3494493/148‑1 (MQ=255)
                                                                                             aTGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCaa                                                   <  1:3246312/148‑1 (MQ=255)
                                                                                              tGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAAc                                                  >  2:2720200/1‑148 (MQ=255)
                                                                                               gTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACg                                                 >  1:1291868/1‑148 (MQ=255)
                                                                                                 ttttGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTc                                               >  2:2397630/1‑148 (MQ=255)
                                                                                                  tttGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAAAAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCg                                              <  2:2013058/148‑1 (MQ=255)
                                                                                                    tGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGt                                            <  2:2654429/148‑1 (MQ=255)
                                                                                                     gggTAAATCTCCATTCATTCAATGAAGGGAAATTGTTTTGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGtt                                           <  1:3180190/148‑1 (MQ=255)
                                                                                                      ggTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGttt                                          <  2:3400438/148‑1 (MQ=255)
                                                                                                      ggTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGttt                                          <  1:799186/148‑1 (MQ=255)
                                                                                                       gTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTc                                         >  1:1511032/1‑148 (MQ=255)
                                                                                                        tAAATCTCCATTCATTCAATGAAGGGAACTTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTcg                                        >  1:1558182/1‑148 (MQ=255)
                                                                                                           aTCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTcggcg                                     <  2:1192715/148‑1 (MQ=255)
                                                                                                           aTCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTcggcg                                     >  1:1707253/1‑148 (MQ=255)
                                                                                                              tcCATTCATTCAATGAAGGGAAATTGTTATGAAAACTGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCCCCTGTACTGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCgggg                                  <  1:2716777/148‑1 (MQ=255)
                                                                                                                cattcattCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAg                                >  1:2151100/1‑148 (MQ=255)
                                                                                                                  ttcattcaATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGtt                              >  1:948647/1‑148 (MQ=255)
                                                                                                                    cattcaATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTtcc                            <  2:2315785/148‑1 (MQ=255)
                                                                                                                    cattcaATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTtcc                            >  1:1234115/1‑148 (MQ=255)
                                                                                                                     attcaATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTtcct                           <  2:2584727/148‑1 (MQ=255)
                                                                                                                     attcaATGAAGGGAAATTGTTATCAAAACAGTTCTGTATGGCATTTTAGCCATATCTGCGCTTGCGGCGACTACTGCGTGGGCTGCACCTGTACAGGTGGGAGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTtcct                           >  1:293360/1‑148 (MQ=255)
                                                                                                                      ttcaATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTtcctc                          <  2:3563313/148‑1 (MQ=255)
                                                                                                                        caATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTTCCTCcg                        <  1:1138780/148‑1 (MQ=255)
                                                                                                                           tGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGa                     >  2:2087237/1‑148 (MQ=255)
                                                                                                                             aaGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGAcc                   >  2:2958611/1‑148 (MQ=255)
                                                                                                                                 gAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCg               >  2:1391578/1‑148 (MQ=255)
                                                                                                                                  aaaTTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGt              >  2:2314804/1‑148 (MQ=255)
                                                                                                                                  aaaTTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGt              >  2:1318984/1‑148 (MQ=255)
                                                                                                                                    aTTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCa            >  1:867910/1‑148 (MQ=255)
                                                                                                                                     ttGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAg           >  1:2625166/1‑148 (MQ=255)
                                                                                                                                         tATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCAcc       >  1:1657250/1‑148 (MQ=255)
                                                                                                                                         tATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCAcc       >  1:1852764/1‑148 (MQ=255)
                                                                                                                                         tATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCAcc       >  1:2308296/1‑148 (MQ=255)
                                                                                                                                         tATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCAcc       >  1:52310/1‑148 (MQ=255)
                                                                                                                                              aaacaGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAg  <  2:2469623/148‑1 (MQ=255)
                                                                                                                                              aaacaGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAg  <  2:1794370/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
TTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCGCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAAAAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAG  >  CP009273/4225696‑4225985

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: ‑