Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,860,159 C→A 100% P33P (CCG→CCT ibpB ← heat shock chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,860,1590CA95.5% 781.9 / 26.6 271P33P (CCG→CCTibpBheat shock chaperone
Reads supporting (aligned to +/- strand):  ref base C (4/8);  new base A (110/147);  total (116/155)
Fisher's exact test for biased strand distribution p-value = 5.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.62e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 275 total aligned reads displayed.

TTTTTTCTCTTCTTTTGGCTGCTCCGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGCGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTAAGAAGCGAGTAAGTACCTGCAAATCCGAAGA  >  CP009273/3860012‑3860305
                                                                                                                                                   |                                                                                                                                                  
ttttttCTCTTCTTTTGGCTGCTCCGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTAcggagg                                                                                                                                                  <  1:88582/150‑1 (MQ=255)
     tctcTTCTTTTGGCTGCTCCGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAAATTTGTACGGAGGGAAGc                                                                                                                                             >  1:2905947/1‑150 (MQ=255)
        cttcttTTGGCTGCTCCGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGcgct                                                                                                                                          >  1:1151677/1‑150 (MQ=255)
          tcttTTGGCTGCTCCGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTgg                                                                                                                                        <  2:4285545/150‑1 (MQ=255)
               tGGCTGCTCCGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTc                                                                                                                                   >  1:2670341/1‑150 (MQ=255)
               tGGCTGCTCCGGCGTGCATTTTACGCTCTGGCGCGTACCTTACAGTTGAATCTATAAATATTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTTGTTATCGTCGCTTTTCTCAATTTTGTACGGAGGGAAGCTCTGGCTTTc                                                                                                                                   <  1:2978702/150‑1 (MQ=255)
                ggCTGCTCCGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCa                                                                                                                                  >  1:1475237/1‑150 (MQ=255)
                  ctgctCCGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGTATCTCTAAATCTTCCTGACGGAAACATGCAAGCGCAAGGTTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAATCTCTGTCTTTCaac                                                                                                                                <  2:2097099/150‑3 (MQ=255)
                  ctgctCCGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCCGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCAcc                                                                                                                                >  2:4155166/1‑150 (MQ=255)
                   tgctCCGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCg                                                                                                                               <  1:666569/150‑1 (MQ=255)
                      tCCGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCg                                                                                                                            >  2:1740647/1‑150 (MQ=255)
                          gCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGtttt                                                                                                                        <  2:4599650/150‑1 (MQ=255)
                           cGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTg                                                                                                                       <  1:1374839/150‑1 (MQ=255)
                              gCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAg                                                                                                                    <  1:147435/150‑1 (MQ=255)
                                cTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTg                                                                                                                  >  2:1356111/1‑150 (MQ=255)
                                 ttttACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGc                                                                                                                 <  1:2297418/150‑1 (MQ=255)
                                    tACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGtt                                                                                                              >  2:577229/1‑150 (MQ=255)
                                    tACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGtt                                                                                                              <  1:178193/150‑1 (MQ=255)
                                     aCGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTg                                                                                                             <  1:3969322/150‑1 (MQ=255)
                                          cAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAg                                                                                                        <  1:4047006/150‑1 (MQ=255)
                                             gcgcGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGttt                                                                                                     >  2:3078789/1‑150 (MQ=255)
                                             gcgcGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGttt                                                                                                     >  2:4535665/1‑150 (MQ=255)
                                             gcgcGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGttt                                                                                                     <  1:1318092/150‑1 (MQ=255)
                                             gcgcGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGttt                                                                                                     <  2:4050494/150‑1 (MQ=255)
                                              cgcgTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTg                                                                                                    >  2:2279725/1‑150 (MQ=255)
                                              cgcgTACCTTCCAGATGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCATTTCTCAATGTTGTACGTAGGGAAGCACTGGCATTCACCGGCGTTTTGCAGTGCGTTTGCCAGtttt                                                                                                    >  2:2280006/1‑149 (MQ=255)
                                               gcgTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGt                                                                                                   <  1:2718929/150‑1 (MQ=255)
                                                cgTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTc                                                                                                  <  2:4586917/150‑1 (MQ=255)
                                                cgTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTc                                                                                                  <  2:2971473/150‑1 (MQ=255)
                                                 gTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGATGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCa                                                                                                 >  1:2281385/1‑150 (MQ=255)
                                                 gTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTATCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGGTTTGCAGTGCGTTGGCCAGTTTTTCa                                                                                                 >  1:4056602/1‑150 (MQ=255)
                                                  tACCTTCCAGTTGAATCTCTAACTCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGCGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCaa                                                                                                <  2:1666979/150‑1 (MQ=255)
                                                  tACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGCGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCaa                                                                                                <  2:3995910/150‑1 (MQ=255)
                                                     cTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTGAATGTTGTAAGGAGGGAAGCTCTGGCGTTCACCGGCGGTTTGCAGTGCGGTGGCCAGTTTGTCAAAAc                                                                                             >  1:2405921/1‑150 (MQ=255)
                                                     cTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAAc                                                                                             <  1:1243045/150‑1 (MQ=255)
                                                       tCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCAACGGCGTTTTGCAGTGCGTTGGCCAGGTTGTAAAAACCg                                                                                           >  2:3819218/1‑150 (MQ=255)
                                                        ccAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGa                                                                                          <  2:508186/150‑1 (MQ=255)
                                                         cAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGAt                                                                                         <  2:2789593/150‑1 (MQ=255)
                                                          aGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATc                                                                                        <  1:1965204/150‑1 (MQ=255)
                                                           gTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATcc                                                                                       >  2:1797456/1‑150 (MQ=255)
                                                             tGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCAt                                                                                     <  1:3894449/150‑1 (MQ=255)
                                                             tGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCAt                                                                                     <  1:2575299/150‑1 (MQ=255)
                                                             tGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCAt                                                                                     <  2:1959522/150‑1 (MQ=255)
                                                              gAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCAtt                                                                                    <  2:1475271/150‑1 (MQ=255)
                                                              gAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCAtt                                                                                    <  2:1475232/150‑1 (MQ=255)
                                                               aaTCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTAAATGTTGTAAGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGGGTTGGCCAGTTTGTCAAAAACGATCCATTg                                                                                   >  1:1674346/1‑150 (MQ=255)
                                                                 tctcTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGAc                                                                                 <  2:3129702/150‑1 (MQ=255)
                                                                  ctctAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACg                                                                                <  2:4454467/150‑1 (MQ=255)
                                                                     tAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCAt                                                                             >  1:2657065/1‑150 (MQ=255)
                                                                     tAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTAACCGGCGTTTTGCAGTGCGTTGGCCCGTTTGTCAAAACCGATCCCTTGACGCAt                                                                             >  2:3555898/1‑150 (MQ=255)
                                                                      aaaTCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATc                                                                            <  2:2837158/150‑1 (MQ=255)
                                                                      aaaTCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATc                                                                            <  2:2059996/150‑1 (MQ=255)
                                                                       aaTCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTATCAATGTTGTACGGAGGGAAGCGCTGGCTTTCACCGGCGGTTTGCAGTGCGTTGGCCCGTTTGTCAAAACCGATACATTGTCGCATCa                                                                           >  2:3561098/1‑150 (MQ=255)
                                                                              cTGACGGAAACCTGCCAGCGCACGGGTAATGCGGTAGTGGTTATCGTCGCGTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCATTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGa                                                                    >  2:3302101/1‑150 (MQ=255)
                                                                               tGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGAt                                                                   <  2:2550132/150‑1 (MQ=255)
                                                                               tGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGAt                                                                   >  2:1262575/1‑150 (MQ=255)
                                                                               tGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGAt                                                                   <  1:976936/150‑1 (MQ=255)
                                                                                   ggAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAAt                                                               <  2:19236/150‑1 (MQ=255)
                                                                                   ggAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAAt                                                               <  1:3328391/150‑1 (MQ=255)
                                                                                   ggAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAAt                                                               <  1:542894/150‑1 (MQ=255)
                                                                                   ggAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTATCAATGTTTTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATaaaa                                                               >  2:4091464/1‑149 (MQ=255)
                                                                                      aaCCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGGAGCTCTGGCGGTCACCGGCGGTTTGCAGTGCGTTGGCCAGGTTGTCAAAACCGATCCCTTGACGCATCAGTGGGGAGAAAACGa                                                            >  2:1673128/1‑150 (MQ=255)
                                                                                        ccTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAg                                                          <  1:3803000/150‑1 (MQ=255)
                                                                                         cTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGt                                                         <  2:1656582/150‑1 (MQ=255)
                                                                                          tGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTGTCAATGTTGTACGGGGGGTAGCTCTGGCTTTCACCGGCGTTTTTCAGTGCGTTGGGCACTTTGTCAAAAACGATCCAGTGACGCATGAGTGGGGATAAATAGAAGtt                                                        >  1:789258/1‑150 (MQ=255)
                                                                                          tGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTTTACGGAGGGAAGCTCTGGCGTTCACCGGCGTTTTGCAGTTCGTTTGCCAGTTTGTCAAAAACGATCCATTGAAGCATAAGGGGGGATAAATAGAAGtt                                                        >  1:156578/1‑150 (MQ=255)
                                                                                          tGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGtt                                                        <  1:3483304/150‑1 (MQ=255)
                                                                                           gCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTTCGTTTGCCAGTTTTTCAAAACCGATCCAGTGACGCATCAGTGGGGATAAATAGAAGTTa                                                       >  1:3739884/1‑150 (MQ=255)
                                                                                             cAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCCTTGACGCATCAGTGGGGATAAATCGAAGTTACg                                                     >  2:239479/1‑150 (MQ=255)
                                                                                             cAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACg                                                     <  1:2267428/150‑1 (MQ=255)
                                                                                             cAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACg                                                     <  1:1504408/150‑1 (MQ=255)
                                                                                                cgcAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGCGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCAt                                                  <  1:3618050/150‑1 (MQ=255)
                                                                                                   aaGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTAAAAGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAAACCATACAATGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGt                                               >  2:3709045/1‑150 (MQ=255)
                                                                                                     gggTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCa                                             >  2:693008/1‑150 (MQ=255)
                                                                                                        tAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTTCAGTGCGTTGGCCAGTTTGTCAAAAACGATCCAATGACGCAACAGTGGGGATAAATCGAAGTTACGCATAGTCAttt                                          >  2:3544030/1‑150 (MQ=255)
                                                                                                             cGGTAGTGGTTATCGTCGCTTTTCTCAATTTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCt                                     >  2:1999467/1‑150 (MQ=255)
                                                                                                             cGGTAGTGGTTATCGTCGCTTTTATCAATGTTGTACGGAGGGAAGCACTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTTTCAAAAACGATCCAATGACGCATCAGTGGGGATAAATAGAAGTTACGCATAGTCAtttataca                                     >  1:3703594/1‑145 (MQ=255)
                                                                                                              ggTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCtt                                    <  2:3146868/150‑1 (MQ=255)
                                                                                                              ggTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCtt                                    >  2:3355911/1‑150 (MQ=255)
                                                                                                              ggTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCtt                                    <  2:2311792/150‑1 (MQ=255)
                                                                                                              ggTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAAACGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCtt                                    >  2:4532749/1‑150 (MQ=255)
                                                                                                                   tGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTaag                               <  1:4392028/150‑1 (MQ=255)
                                                                                                                     gTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTaagaa                             <  1:4136220/150‑1 (MQ=255)
                                                                                                                     gTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTaagaa                             <  2:4189186/150‑1 (MQ=255)
                                                                                                                      ttATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTaagaag                            <  1:2355661/150‑1 (MQ=255)
                                                                                                                         tcgtcgCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTAAGAAGCGa                         <  1:2206060/150‑1 (MQ=255)
                                                                                                                         tcgtcgCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTAAGAAGCGa                         <  2:3944177/150‑1 (MQ=255)
                                                                                                                              gCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTAAGAAGCGagtaag                    <  1:3371157/150‑1 (MQ=255)
                                                                                                                               cTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTAAGAAGCGagtaagt                   <  2:4426373/150‑1 (MQ=255)
                                                                                                                                   tctcAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTAAGAAGCGAGTAAGTACCt               <  1:1460383/150‑1 (MQ=255)
                                                                                                                                    ctcAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTAAGAAGCGAGTAAGTACCTg              <  1:1408396/150‑1 (MQ=255)
                                                                                                                                      cAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTAAGAAGCGAGTAAGTACCTGCa            >  1:2488718/1‑150 (MQ=255)
                                                                                                                                       aaTGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTAAGAAGCGAGTAAGTACCTGCaa           >  1:1025898/1‑150 (MQ=255)
                                                                                                                                        aTGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTAAGAAGCGAGTAAGTACCTGCaaa          >  1:3116300/1‑150 (MQ=255)
                                                                                                                                           ttgtACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTAAGAAGCGAGTAAGTACCTGCAAATcc       >  1:2518949/1‑150 (MQ=255)
                                                                                                                                            tgtACGGCGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTAAGAAGCGAGTAAGTACCTGCAAATCCg      <  1:842804/150‑1 (MQ=255)
                                                                                                                                              tACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTAAGAAGCGAGTAAGTACCTGCAAATCCGaa    >  1:732768/1‑150 (MQ=255)
                                                                                                                                              tACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTAAGAAGCGAGTAAGTACCTGCAAATCCGaa    >  2:2970451/1‑150 (MQ=255)
                                                                                                                                               aCGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTAAGAAGCGAGTAAGTACCTGCAAATCCGAAg   <  1:2055757/150‑1 (MQ=255)
                                                                                                                                                cggaggGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTAAGAAGCGAGTAAGTACCTGCAAATCCGAAGa  >  2:2520251/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TTTTTTCTCTTCTTTTGGCTGCTCCGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGCGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTAAGAAGCGAGTAAGTACCTGCAAATCCGAAGA  >  CP009273/3860012‑3860305

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: