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,8410AC98.2% 1016.9 / ‑1.0 254K3T (AAA→ACA) yjbEextracellular polysaccharide production threonine‑rich protein
Reads supporting (aligned to +/- strand):  ref base A (4/1);  new base C (140/109);  total (144/110)
Fisher's exact test for biased strand distribution p-value = 3.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-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 260 total aligned reads displayed.

TTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCGCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAAAAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAG  >  CP009273/4225695‑4225985
                                                                                                                                                  |                                                                                                                                                
tttGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAg                                                                                                                                                >  2:505141/1‑149 (MQ=255)
 ttGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGt                                                                                                                                               <  1:2740130/149‑1 (MQ=255)
 ttGATATTCATTCCAATTATCCTAAAACGCCATCGCTAAGTCCCCTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGt                                                                                                                                               <  2:3735975/149‑1 (MQ=255)
    atatTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCt                                                                                                                                            <  1:3123452/149‑1 (MQ=255)
      attcattcCAATTATCCTAAAACGCCATCGCTAATTCCCCTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGt                                                                                                                                          <  1:1561418/149‑1 (MQ=255)
          attcCAATTATCCTAAAACGCCATCGCTAATTCCCCTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATgg                                                                                                                                      >  2:2026367/1‑149 (MQ=255)
             ccAATTATCCTAAAACGCCATCGCTAATTCCCCTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCAt                                                                                                                                   <  2:4079195/149‑1 (MQ=255)
               aaTTATCCTAAAACGCCATCGCTAATTCCCCTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCAttt                                                                                                                                 >  2:3515350/1‑149 (MQ=255)
                aTTATCCTAAAACGCCATCGCTAATTCCCCTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCAtttt                                                                                                                                <  2:1707447/149‑1 (MQ=255)
                  tATCCTAAAACGCCATCGCTAATTCCCCTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTg                                                                                                                              <  1:1618891/149‑1 (MQ=255)
                   aTCCTAAAACGCCATCGCTAATTCCCCTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGc                                                                                                                             >  1:2981000/1‑149 (MQ=255)
                    tCCTAAAACGCCATCGCTAATTCCCCTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGcc                                                                                                                            >  2:162803/1‑149 (MQ=255)
                     ccTAAAACGCCATCGCTAATTCCCCTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCa                                                                                                                           <  2:1743487/149‑1 (MQ=255)
                      cTAAAACGCCATCGCTAATTCCCCTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCat                                                                                                                          >  1:1306243/1‑149 (MQ=255)
                         aaaCGCCATCGCTAATTCCCCTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATc                                                                                                                       <  1:3794222/149‑1 (MQ=255)
                         aaaCGCCATCGCTAATTCCCCTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATc                                                                                                                       <  1:2974529/149‑1 (MQ=255)
                          aaCGCCATCGCTAATTCCCCTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCt                                                                                                                      <  1:1639283/149‑1 (MQ=255)
                            cGCCATCGCTAATTCCCCTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTgc                                                                                                                    <  2:2806081/149‑1 (MQ=255)
                               cATCGCTAATTCCCCTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCt                                                                                                                 <  1:528571/149‑1 (MQ=255)
                               cATCGCTAATTCCCCTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCt                                                                                                                 <  1:792063/149‑1 (MQ=255)
                               cATCGCTAATTCCCCTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCt                                                                                                                 <  1:3555638/149‑1 (MQ=255)
                                aTCGCTAATTCCCCTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCtt                                                                                                                >  1:2590416/1‑149 (MQ=255)
                                aTCGCTAATTCCCCTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCtt                                                                                                                <  1:361783/149‑1 (MQ=255)
                                aTCGCTAATTCCCCTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCtt                                                                                                                <  2:1596719/149‑1 (MQ=255)
                                  cGCTAATTCCCCTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTgc                                                                                                              <  1:2633347/149‑1 (MQ=255)
                                   gCTAATTCCCCTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTgcg                                                                                                             >  1:1153110/1‑149 (MQ=255)
                                   gCTAATTCCCCTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTgcg                                                                                                             >  1:1774477/1‑149 (MQ=255)
                                    cTAATTCCCCTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTgcgg                                                                                                            >  2:2271952/1‑149 (MQ=255)
                                     tAATTCCCCTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTgcggc                                                                                                           >  2:3544880/1‑149 (MQ=255)
                                      aaTTCCCCTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTgcggcg                                                                                                          >  1:2046724/1‑149 (MQ=255)
                                      aaTTCCCCTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTgcggcg                                                                                                          >  1:1991197/1‑149 (MQ=255)
                                        ttCCCCTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAATTTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGAc                                                                                                        <  2:1471778/149‑1 (MQ=255)
                                        ttCCCCTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAATTTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGAc                                                                                                        <  2:2063190/149‑1 (MQ=255)
                                         tCCCCTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACt                                                                                                       <  1:3074701/149‑1 (MQ=255)
                                           cccTCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTc                                                                                                     <  2:1575236/149‑1 (MQ=255)
                                                   gTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGTCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTggg                                                                                             >  1:2808358/1‑149 (MQ=255)
                                                   gTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTggg                                                                                             >  1:2582824/1‑149 (MQ=255)
                                                            aTGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTg                                                                                    >  1:1856681/1‑149 (MQ=255)
                                                             tGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGt                                                                                   >  1:3311035/1‑149 (MQ=255)
                                                             tGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGt                                                                                   >  1:2444526/1‑149 (MQ=255)
                                                               cTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTAc                                                                                 >  2:3312729/1‑149 (MQ=255)
                                                               cTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTAc                                                                                 >  1:2074108/1‑149 (MQ=255)
                                                                 ttAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACGGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAg                                                                               <  2:447890/149‑1 (MQ=255)
                                                                     ttGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTgg                                                                           >  1:820727/1‑149 (MQ=255)
                                                                        tgtATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCg                                                                        >  2:3503717/1‑149 (MQ=255)
                                                                         gtATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGa                                                                       <  1:452133/149‑1 (MQ=255)
                                                                               tCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGc                                                                 <  2:3025138/149‑1 (MQ=255)
                                                                               tCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGc                                                                 <  1:2238928/149‑1 (MQ=255)
                                                                                  aTCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAgg                                                              >  2:3119395/1‑149 (MQ=255)
                                                                                  aTCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAgg                                                              >  2:2479806/1‑149 (MQ=255)
                                                                                  aTCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAgg                                                              >  2:1885767/1‑149 (MQ=255)
                                                                                   tCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAggg                                                             <  2:3417379/149‑1 (MQ=255)
                                                                                   tCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAggg                                                             <  2:4001172/149‑1 (MQ=255)
                                                                                    cccgcccgAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGt                                                            >  1:399690/1‑149 (MQ=255)
                                                                                     ccgcccgAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTc                                                           <  2:2422840/149‑1 (MQ=255)
                                                                                        cccgAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGgc                                                        <  2:1059601/149‑1 (MQ=255)
                                                                                        cccgAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGgc                                                        <  1:3063140/149‑1 (MQ=255)
                                                                                             aaTGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCaa                                                   <  2:954629/149‑1 (MQ=255)
                                                                                              aTGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAAc                                                  >  1:3426149/1‑149 (MQ=255)
                                                                                              aTGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAAc                                                  >  2:2679404/1‑149 (MQ=255)
                                                                                               tGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCACCg                                                 >  2:3179046/1‑149 (MQ=255)
                                                                                                 tttttGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTc                                               <  1:3912234/149‑1 (MQ=255)
                                                                                                  ttttGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCg                                              <  2:2456723/149‑1 (MQ=255)
                                                                                                   tttGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCgg                                             >  2:3067482/1‑149 (MQ=255)
                                                                                                    ttGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGt                                            <  1:3892282/149‑1 (MQ=255)
                                                                                                     tGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGtt                                           <  1:2379826/149‑1 (MQ=255)
                                                                                                      gggTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGttt                                          <  1:1138113/149‑1 (MQ=255)
                                                                                                         tAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTcgg                                       <  2:1371225/149‑1 (MQ=255)
                                                                                                           aaTCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTcggcg                                     >  2:3262532/1‑149 (MQ=255)
                                                                                                            aTCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTcggcgg                                    >  1:335846/1‑149 (MQ=255)
                                                                                                                 cattcattCAATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGt                               <  2:2444525/149‑1 (MQ=255)
                                                                                                                    tcattcaATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTtcc                            <  1:286039/149‑1 (MQ=255)
                                                                                                                    tcattcaATGAAGGGAAATTGTTATGAAAAAAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTtcc                            >  2:3979059/1‑149 (MQ=255)
                                                                                                                      attcaATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTtcctc                          <  1:2520808/149‑1 (MQ=255)
                                                                                                                         caATGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTTCCTCcgc                       >  1:2497078/1‑149 (MQ=255)
                                                                                                                            tGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGAc                    >  2:2616259/1‑149 (MQ=255)
                                                                                                                            tGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGAc                    >  2:2643612/1‑149 (MQ=255)
                                                                                                                            tGAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGAc                    >  2:583245/1‑149 (MQ=255)
                                                                                                                             gAAGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGAcc                   <  1:3243122/149‑1 (MQ=255)
                                                                                                                              aaGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCa                  >  1:1346423/1‑149 (MQ=255)
                                                                                                                              aaGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCa                  >  1:4039704/1‑149 (MQ=255)
                                                                                                                              aaGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCa                  >  1:1689823/1‑149 (MQ=255)
                                                                                                                              aaGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCa                  >  1:861534/1‑149 (MQ=255)
                                                                                                                               aGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCAg                 >  2:753812/1‑149 (MQ=255)
                                                                                                                               aGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCAg                 >  2:1271951/1‑149 (MQ=255)
                                                                                                                               aGGGAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCAg                 >  2:2843345/1‑149 (MQ=255)
                                                                                                                                  gAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGt              >  2:745744/1‑149 (MQ=255)
                                                                                                                                  gAAATTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGt              >  2:3804990/1‑149 (MQ=255)
                                                                                                                                  gAAATTGTTATGAAAAAAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGt              >  2:2652051/1‑149 (MQ=255)
                                                                                                                                    aaTTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCa            <  2:1306243/149‑1 (MQ=255)
                                                                                                                                    aaTTGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCa            >  2:2106956/1‑149 (MQ=255)
                                                                                                                                      ttGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGc          >  1:3605157/1‑149 (MQ=255)
                                                                                                                                      ttGTTATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGc          >  2:2979711/1‑149 (MQ=255)
                                                                                                                                         ttATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCAcc       >  2:2492002/1‑149 (MQ=255)
                                                                                                                                         ttATGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCAcc       >  1:957764/1‑149 (MQ=255)
                                                                                                                                            tGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTa    >  2:2510457/1‑149 (MQ=255)
                                                                                                                                            tGAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTa    <  1:446696/149‑1 (MQ=255)
                                                                                                                                             gAAAACAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTaa   >  1:193080/1‑149 (MQ=255)
                                                                                                                                              aaaacaGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAg  <  2:3604207/149‑1 (MQ=255)
                                                                                                                                              aaaacaGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAg  >  2:600747/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
TTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCGCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAAAAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAG  >  CP009273/4225695‑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: