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,8410AC97.8% 972.5 / 1.4 242K3T (AAA→ACA) yjbEextracellular polysaccharide production threonine‑rich protein
Reads supporting (aligned to +/- strand):  ref base A (3/2);  new base C (135/101);  total (138/104)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 249 total aligned reads displayed.

TTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCGCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAAAAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGC  >  CP009273/4225696‑4225986
                                                                                                                                                 |                                                                                                                                                 
ttGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAg                                                                                                                                                 >  2:2824180/1‑148 (MQ=255)
ttGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAg                                                                                                                                                 <  2:1468365/148‑1 (MQ=255)
ttGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAg                                                                                                                                                 <  2:1107987/148‑1 (MQ=255)
ttGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAg                                                                                                                                                 <  2:876078/148‑1 (MQ=255)
 tGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGt                                                                                                                                                <  1:3017367/148‑1 (MQ=255)
 tGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGt                                                                                                                                                <  1:22967/148‑1 (MQ=255)
  gATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGtt                                                                                                                                               <  2:3182176/148‑1 (MQ=255)
   atatTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTc                                                                                                                                              <  2:3069100/148‑1 (MQ=255)
    tatTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCGCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAAAAGTTCt                                                                                                                                             <  2:2785250/148‑1 (MQ=255)
     attcattcCAATTATCCTAAAACTCAATCGCTAATTACCCTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCTACCGAAATGTTTTTGGTTAAATCTCCATTCATTCAATGAAGCGAAATTGTTATGAAAACAGTTCTg                                                                                                                                            >  2:3281129/1‑148 (MQ=255)
     attcattcCAATTATCCTAAAACGCCATCGCTAATTCCCCTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTCTGAAAACAGTTCTg                                                                                                                                            >  2:1620741/1‑148 (MQ=255)
         attcCAATTATCCTAAAACGCCATCGCTAATTCCCCTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATg                                                                                                                                        >  2:1298008/1‑148 (MQ=255)
          ttcCAATTATCCTAAAACGCCATCGCTAATTCCCCTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATgg                                                                                                                                       >  1:1104174/1‑148 (MQ=255)
           tcCAATTATCCTAAAACGCCATCGCTAATTCCCCTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGc                                                                                                                                      >  2:2402773/1‑148 (MQ=255)
             cAATTATCCTAAAACGCCATCGCTAATTCCCCTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCAt                                                                                                                                    >  2:561722/1‑148 (MQ=255)
                ttATCCTAAAACGCCATCGCTAATTCCCCTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCAtttt                                                                                                                                 <  1:1024447/148‑1 (MQ=255)
                ttATCCTAAAACGCCATCGCTAATTCCCCTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCAtttt                                                                                                                                 >  1:700314/1‑148 (MQ=255)
                    ccTAAAACGCCATCGCTAATTCCCCTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGcc                                                                                                                             <  1:1367082/148‑1 (MQ=255)
                    ccTAAAACGCCATCGCTAATTCCCCTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGcc                                                                                                                             <  2:2319425/148‑1 (MQ=255)
                       aaaaCGCCATCGCTAAGGCCCCTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCata                                                                                                                          <  2:94489/148‑1 (MQ=255)
                        aaaCGCCATCGCTAATTCCCCTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCatat                                                                                                                         <  1:2994855/148‑1 (MQ=255)
                        aaaCGCCATCGCTAATTCCCCTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCatat                                                                                                                         <  2:3491931/148‑1 (MQ=255)
                        aaaCGCCATCGCTAATTCCCCTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCatat                                                                                                                         <  1:2455359/148‑1 (MQ=255)
                            gCCATCGCTAATTCCCCTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTgc                                                                                                                     <  2:2657492/148‑1 (MQ=255)
                             ccATCGCTAATGCCCCTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTgcg                                                                                                                    <  2:2046933/148‑1 (MQ=255)
                              cATCGCTAATTCCCCTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTgcgc                                                                                                                   <  1:3050267/148‑1 (MQ=255)
                                tCGCTAATTCCCCTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCtt                                                                                                                 <  1:949211/148‑1 (MQ=255)
                                 cGCTAATTCCCCTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTg                                                                                                                >  1:1047300/1‑148 (MQ=255)
                                    gaaTTCCCCTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCAGTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTgcgg                                                                                                             <  2:2967321/147‑1 (MQ=255)
                                       ttCCCCTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGa                                                                                                          >  1:52396/1‑148 (MQ=255)
                                          cccTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACtt                                                                                                       >  1:1193066/1‑148 (MQ=255)
                                          cccTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACtt                                                                                                       <  1:2011598/148‑1 (MQ=255)
                                             tcgcCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCGTTTTTGCCATATCTGCGCTTGCGGCGACTTCTg                                                                                                    >  2:723489/2‑148 (MQ=255)
                                                 cGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTg                                                                                                >  1:275192/1‑148 (MQ=255)
                                                 cGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTg                                                                                                <  1:1525177/148‑1 (MQ=255)
                                                   tAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATATCCATTCATTCAATGAAGGGAAATTGTTATGAAAAAAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTggg                                                                                              >  1:176050/1‑148 (MQ=255)
                                                    aaTTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGc                                                                                             <  2:1030378/148‑1 (MQ=255)
                                                    aaTTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAAAAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGc                                                                                             <  2:2595313/148‑1 (MQ=255)
                                                      ttCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTg                                                                                           <  1:3476950/148‑1 (MQ=255)
                                                      ttCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTg                                                                                           <  2:1046387/148‑1 (MQ=255)
                                                          cATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCAcc                                                                                       >  1:494315/1‑148 (MQ=255)
                                                             gCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGt                                                                                    >  1:2850143/1‑148 (MQ=255)
                                                               tttAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTAc                                                                                  <  2:289670/148‑1 (MQ=255)
                                                                ttAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACa                                                                                 >  2:2522075/1‑148 (MQ=255)
                                                                 tAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGACATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAg                                                                                >  2:1553418/1‑148 (MQ=255)
                                                                  aGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAgg                                                                               <  1:2157664/148‑1 (MQ=255)
                                                                     tgtgTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTgg                                                                            >  1:2063890/1‑148 (MQ=255)
                                                                        gtATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCg                                                                         >  1:321449/1‑148 (MQ=255)
                                                                         tataCTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGa                                                                        >  1:1567243/1‑148 (MQ=255)
                                                                           taCTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAg                                                                      <  2:1980368/148‑1 (MQ=255)
                                                                            aCTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGc                                                                     >  1:80936/1‑148 (MQ=255)
                                                                             cTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCg                                                                    >  1:664774/1‑148 (MQ=255)
                                                                              tctaTCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCgg                                                                   >  2:358375/4‑148 (MQ=255)
                                                                               cGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGc                                                                  >  1:1242716/1‑148 (MQ=255)
                                                                               cGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGc                                                                  >  1:1440768/1‑148 (MQ=255)
                                                                                 aTCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAg                                                                >  2:388343/1‑148 (MQ=255)
                                                                                  tCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAgg                                                               >  1:1834241/1‑148 (MQ=255)
                                                                                  tCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAgg                                                               >  2:700218/1‑148 (MQ=255)
                                                                                   cccgcccgAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAggg                                                              <  1:1472139/148‑1 (MQ=255)
                                                                                   cccgcccgAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAggg                                                              <  2:3451530/148‑1 (MQ=255)
                                                                                          gAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGgcag                                                       >  1:60934/1‑148 (MQ=255)
                                                                                             aTGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCaa                                                    <  1:2352084/148‑1 (MQ=255)
                                                                                             aTGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCaa                                                    >  1:2400309/1‑148 (MQ=255)
                                                                                             aTGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCaa                                                    >  2:114116/1‑148 (MQ=255)
                                                                                              tGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAAc                                                   >  1:2037877/1‑148 (MQ=255)
                                                                                               gTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACg                                                  <  1:2934974/148‑1 (MQ=255)
                                                                                                tttttGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGt                                                 <  2:1783387/148‑1 (MQ=255)
                                                                                                 ttttGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCGGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTc                                                <  1:2424631/148‑1 (MQ=255)
                                                                                                  tttGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCg                                               <  1:2884513/148‑1 (MQ=255)
                                                                                                  tttGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCg                                               <  1:1985328/148‑1 (MQ=255)
                                                                                                   ttGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCgg                                              <  2:1590420/148‑1 (MQ=255)
                                                                                                        tAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTcg                                         >  1:2917692/1‑148 (MQ=255)
                                                                                                          aaTCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTcggc                                       <  2:3400146/148‑1 (MQ=255)
                                                                                                               ccATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGa                                  >  1:1933445/1‑148 (MQ=255)
                                                                                                               ccATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGa                                  <  2:500080/148‑1 (MQ=255)
                                                                                                                 attcattcAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGt                                <  1:1584837/148‑1 (MQ=255)
                                                                                                                 attcattcAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGt                                <  1:3321170/148‑1 (MQ=255)
                                                                                                                  ttcattcaATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTCCAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGtt                               >  2:3357976/1‑148 (MQ=255)
                                                                                                                   tcattcaATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTtc                              <  1:2824180/148‑1 (MQ=255)
                                                                                                                     attcaATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTtcct                            >  1:2168243/1‑148 (MQ=255)
                                                                                                                      ttcaATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTtcctc                           >  1:1932714/1‑148 (MQ=255)
                                                                                                                             aaGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGAcc                    <  1:2553084/148‑1 (MQ=255)
                                                                                                                             aaGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGAcc                    >  2:1323319/1‑148 (MQ=255)
                                                                                                                             aaGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGAcc                    <  1:358375/148‑1 (MQ=255)
                                                                                                                              aGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCa                   >  1:2490706/1‑148 (MQ=255)
                                                                                                                              aGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCa                   >  1:3153497/1‑148 (MQ=255)
                                                                                                                               gggAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCAg                  >  2:2847531/1‑148 (MQ=255)
                                                                                                                                ggAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCAGc                 <  1:821386/148‑1 (MQ=255)
                                                                                                                                 gAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCg                <  1:103233/148‑1 (MQ=255)
                                                                                                                                  aaaTTGTTATTAATATAGTTCTTTATGGCATTTTTGCAATATCTGAGCTTGCGGCGACTTATGCGTGGGATGCAAGTGTACTGGTGGGCGTATCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGt               <  1:3281129/148‑1 (MQ=255)
                                                                                                                                  aaaTTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGt               <  2:1242716/148‑1 (MQ=255)
                                                                                                                                   aaTTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTc              >  2:2529234/1‑148 (MQ=255)
                                                                                                                                    aTTGTTATGAAAACAGTTCTGTATGTCATTTTTGCCATATCTGAGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTTACTCCGCGACCAGCGTCa             >  1:921813/1‑148 (MQ=255)
                                                                                                                                     ttGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAg            >  2:1106011/1‑148 (MQ=255)
                                                                                                                                      tGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGc           >  1:1853187/1‑148 (MQ=255)
                                                                                                                                        ttATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCAc         >  1:56132/1‑148 (MQ=255)
                                                                                                                                          aTGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCACCg       >  1:429055/1‑148 (MQ=255)
                                                                                                                                            gAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTa     >  2:300317/1‑148 (MQ=255)
                                                                                                                                            gAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTa     >  2:257867/1‑148 (MQ=255)
                                                                                                                                              aaacaGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAg   <  2:2108431/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
TTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCGCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAAAAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGC  >  CP009273/4225696‑4225986

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: