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,1590CA99.1% 521.5 / ‑6.2 187P33P (CCG→CCTibpBheat shock chaperone
Reads supporting (aligned to +/- strand):  ref base C (3/0);  major base A (58/120);  minor base T (4/0);  total (67/120)
Fisher's exact test for biased strand distribution p-value = 1.26e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
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 188 total aligned reads displayed.

TTTTCTCTTCTTTTGGCTGCTCCGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGCGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTAAGAAGCGAGTAAGTACCTGCAAATCC  >  CP009273/3860014‑3860300
                                                                                                                                                 |                                                                                                                                               
ttttCTCTTCTTTTGGCTGCTCCGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTAcggagg                                                                                                                                               <  2:22855/148‑1 (MQ=255)
     tcttctTTTGGCTGCTCCGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTGTCAATGTTGTACGGAGGGCAGc                                                                                                                                          >  1:791325/1‑148 (MQ=255)
          ttttGGCTGCTCCGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTgg                                                                                                                                     <  1:1920771/148‑1 (MQ=255)
           tttGGCTGCTCCGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGGACGGAGGGAAGCTCTGGc                                                                                                                                    >  2:2531581/1‑148 (MQ=255)
            ttGGCTGCTCCGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCt                                                                                                                                   <  1:117382/148‑1 (MQ=255)
            ttGGCTGCTCCGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCt                                                                                                                                   <  1:1268957/148‑1 (MQ=255)
             tGGCTGCTCCGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCtt                                                                                                                                  <  1:594876/148‑1 (MQ=255)
               gctgctCCGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTc                                                                                                                                >  2:1756180/1‑148 (MQ=255)
               gctgctCCGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTc                                                                                                                                >  2:1756229/1‑148 (MQ=255)
               gctgctCCGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTc                                                                                                                                >  2:1834857/1‑148 (MQ=255)
                ctgctCCGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCa                                                                                                                               <  2:885459/148‑1 (MQ=255)
                 tgctCCGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCAc                                                                                                                              <  1:454485/148‑1 (MQ=255)
                  gctCCGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCAcc                                                                                                                             <  2:1377203/148‑1 (MQ=255)
                      cGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCg                                                                                                                         <  1:1117381/148‑1 (MQ=255)
                      cGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCg                                                                                                                         <  2:1852923/148‑1 (MQ=255)
                        gCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGtt                                                                                                                       <  1:2797540/148‑1 (MQ=255)
                           tGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTg                                                                                                                    <  1:1967727/148‑1 (MQ=255)
                            gCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTATGAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTGTTGc                                                                                                                   >  2:1044970/1‑148 (MQ=255)
                               ttttACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTAAGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGt                                                                                                               >  2:2834115/1‑148 (MQ=255)
                                tttACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTg                                                                                                              <  2:1016827/148‑1 (MQ=255)
                                   aCGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGt                                                                                                           <  1:858349/148‑1 (MQ=255)
                                          ggCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAg                                                                                                    <  1:1900431/148‑1 (MQ=255)
                                           gcgcGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGt                                                                                                   <  2:601995/148‑1 (MQ=255)
                                           gcgcGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGt                                                                                                   <  2:1186539/148‑1 (MQ=255)
                                            cgcgTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGtt                                                                                                  <  1:1578801/148‑1 (MQ=255)
                                               gTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGt                                                                                               <  2:1768938/148‑1 (MQ=255)
                                                 aCCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTATCAATGTTGTTAGGGGGGAAGCTCTGGCGGTCAACGGCGGTTTGGAGTGCGTTGGGCAGTTTGTCa                                                                                             >  2:633770/1‑148 (MQ=255)
                                                   cTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGGTGGCCAGTTTGTCcaa                                                                                           >  2:3219773/1‑148 (MQ=255)
                                                   cTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAAATTTTTAAGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCaaa                                                                                           >  2:1131982/1‑148 (MQ=255)
                                                     tCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTGTCAATGTTGTTCGGAGGGAAGCTCATGCTTTCACCGGCGTTTTTGAGTGCGGTTGCCACGTTGTCaaaaa                                                                                         >  2:1005383/1‑147 (MQ=255)
                                                     tCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTATCAATGTTGTAAGGCGGGAAGCTCTGGCTTTAAACGGCGTTTTTCAGTTGGGTTGGCAGTTTTTCaaaaa                                                                                         >  1:2229206/1‑147 (MQ=255)
                                                       cAGTTGTATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCg                                                                                       <  2:1262837/148‑1 (MQ=255)
                                                       cAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCg                                                                                       <  1:2435708/148‑1 (MQ=255)
                                                       cAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCg                                                                                       <  2:1261884/148‑1 (MQ=255)
                                                         gTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGAt                                                                                     <  2:2091448/148‑1 (MQ=255)
                                                         gTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTACAGTGCGTTGGCCAGTTGGTCAAAACCGAt                                                                                     <  2:3448002/148‑1 (MQ=255)
                                                             aaTCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGTGGGAAGCTATGGCTTTAACCGGCGTTTTGAAGTGCGGTGGCCAGTTTGTAAAAACCCATCCaa                                                                                >  2:1187427/1‑147 (MQ=255)
                                                               tctcTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTGTAAATGTTGTACGGAGGGAAGCACTGGCGTTAACCGGCGGTTTGCAGTGCGTTGGGCAGTTTTTCAAAACCGATCCCttt                                                                              >  2:1285804/1‑147 (MQ=255)
                                                               tctcTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCGTTTATGAATGTTGTTAGGAGGGAAGCTCTGGCTTTGACCGGGGGTTTTGGGTGGGGTGGGCAGTTTGTTAAAAACGATACCAtt                                                                               >  1:463932/1‑148 (MQ=255)
                                                                  ctAAATCTTCCTGTCGGAAACCTGCCAGCGCACGGGTAATGCGTTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGTCGTTTTGCTGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACg                                                                           <  1:349232/148‑1 (MQ=255)
                                                                   tAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTGTGAATGTTGTACGGCGGGGAGCTCTGGGGTTGACCGGCGGGTTGGAGTGGGGTGGGCCGGTTGTGAAAAACGGTCCATTGACGc                                                                          >  1:2170634/1‑148 (MQ=255)
                                                                    aaaTCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCa                                                                         <  2:522855/148‑1 (MQ=255)
                                                                     aaTCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCAt                                                                        <  1:2867818/148‑1 (MQ=255)
                                                                       tCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGGCAGTTTGTCAAAACCGATCCATTGACGCATCa                                                                      >  2:3183046/1‑148 (MQ=255)
                                                                         ttCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGt                                                                    <  1:2303588/148‑1 (MQ=255)
                                                                           ccTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTgg                                                                  <  2:949726/148‑1 (MQ=255)
                                                                           ccTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTgg                                                                  <  1:1248521/148‑1 (MQ=255)
                                                                           ccTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTgg                                                                  <  1:431183/148‑1 (MQ=255)
                                                                               aCGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCGTTTCTCAAGGTTGTACGGAGGGAAGCCCAGGCTTTCAACGGCGTTGTTCAGTGCGTTTGCCCGCTTGTCAAAACAGCAACCTTGAAGCATCCGAGGGGAt                                                              >  2:3410593/1‑148 (MQ=255)
                                                                                 ggAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATaa                                                            <  1:543437/148‑1 (MQ=255)
                                                                                  gAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATaaa                                                           <  2:2617915/148‑1 (MQ=255)
                                                                                   aaaCCTGCCAGCGCTAGGGTATTGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTTGCTTTCACCGGAGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAAt                                                          <  1:3451043/148‑1 (MQ=255)
                                                                                        tGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGGTTTCACCGGCGTTTTTCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGGGGGGATAAATCGAAg                                                     >  2:459914/1‑148 (MQ=255)
                                                                                           cAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTa                                                  <  1:1005383/148‑1 (MQ=255)
                                                                                           cAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTa                                                  <  2:3409874/148‑1 (MQ=255)
                                                                                              cgcAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGGCGCATCAGTGGGGATAAATCGAAGTTACGc                                               <  2:3214219/148‑1 (MQ=255)
                                                                                               gcAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCa                                              <  1:2517666/148‑1 (MQ=255)
                                                                                                cAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCAt                                             <  2:815840/148‑1 (MQ=255)
                                                                                                cAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCAt                                             <  1:2039080/148‑1 (MQ=255)
                                                                                                 aaGGGTAATGCGGTAGTGGTTATCGTCGCTTTTATCAAAGTTGTACGGAGGGCAGCTCTGGCTTTCACCGGCGTGTTGCAGTGCGTTGGCCAGTGTGTCAAAACCGATACAATTACGCAACACTGGGGTTAAAAAGAATTTAAGCATa                                            >  2:1805310/1‑148 (MQ=255)
                                                                                                  aGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAg                                           <  2:1602806/148‑1 (MQ=255)
                                                                                                  aGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAg                                           <  2:2096336/148‑1 (MQ=255)
                                                                                                  aGGGTAATGCGGTAGTGGTTATCGTCGCTGTTCTCACTGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAg                                           <  2:2161048/148‑1 (MQ=255)
                                                                                                   gggTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGt                                          <  2:2791892/148‑1 (MQ=255)
                                                                                                   gggTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGGTTGTCAAAACCGATCCAGTGACGCCTCAGGGGGGGTAAATAGAAGTTAAGCATAGt                                          >  2:1087802/1‑148 (MQ=255)
                                                                                                   gggTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCAACGGCGTTTTTCAGTGCGTTGGCCAGTTTGTCAAAACCCATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGt                                          >  2:1087783/1‑148 (MQ=255)
                                                                                                    ggTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTc                                         <  1:1677111/148‑1 (MQ=255)
                                                                                                         tGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTtc                                    <  1:1792075/148‑1 (MQ=255)
                                                                                                           cGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTtctc                                  <  1:2057912/148‑1 (MQ=255)
                                                                                                           cGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCCGCGTTTTGCAGTGCGGTGGCCAGTTTGTCAAAACCGATCCATTGTCGCATCAGTGGGGATAAATAGAAGTTACGCAAAGTCATTtatc                                  >  2:3393695/1‑148 (MQ=255)
                                                                                                            ggTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTcc                                 <  2:1715635/148‑1 (MQ=255)
                                                                                                              tAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCtt                               <  2:1953165/148‑1 (MQ=255)
                                                                                                              tAGTGGTTATCGTCGCTTTTCTAAATGTTTTACGGAGGGAAGCTCAGGCTTTAAACGGCGTTTTTCAGTGCGTTGGCCAGTTTGTCAAAAACGATCCAATTACGCCTCAGTGGGGATAAATCGAAGTTAAGCATAGTCAtttatacat                               >  2:855242/1‑142 (MQ=255)
                                                                                                               aGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCTGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTc                              <  1:1696256/148‑1 (MQ=255)
                                                                                                               aGTGGTTATCGTCGCTTTTCTCAAATTTGTTAGGAGGGACGCTCTGGCTTTCACCCGCGGTTTGCAGTGCGGTGGCCAGATTGTCAAAAACCTTACAATGAAGCCAAAGTGGGGATAAAACGAAGTTAATCAAAGTCATTTATCCtta                              >  2:3042988/1‑147 (MQ=255)
                                                                                                               aGTGGTTATCGTCGCTTTTATCAAATTTTTAAGGAGGGAAACTATGGGTTTAACCGGCGTTTTGCAGTGCGTTTGGCCGTTTGTCAAAACCCATCCAATTTCACAACAATGGGGATAAAAAGAAGTTAACCATAGTCAATTATccgac                              >  2:3376201/1‑145 (MQ=255)
                                                                                                                 tGGTTATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTa                            <  1:2548292/148‑1 (MQ=255)
                                                                                                                    ttATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTaaga                         <  1:1501370/148‑1 (MQ=255)
                                                                                                                     tATCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTaagaa                        <  2:351128/148‑1 (MQ=255)
                                                                                                                      aTCGTCGCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTaagaag                       <  1:3046497/148‑1 (MQ=255)
                                                                                                                       tcgtcgCTTTTCTCAATTTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGGTGGCCAGTTTTTCAAAACCGATCCAATGACGCATCAGTGGGGATAAATAGAAGTTAAGCATAGTCATTTATCCTTCTAATAAGc                      >  2:3222441/1‑148 (MQ=255)
                                                                                                                        cgtcgCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTAAGAAGCg                     <  1:2158646/148‑1 (MQ=255)
                                                                                                                         gtcgCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTAAGAAGCGa                    <  2:2229206/148‑1 (MQ=255)
                                                                                                                          tcgCTTTTCTCATTGTTGTAGGGAGGGAAGCTCTGTCTTTCACCGGCGTTTTGCAGTGCGTTGGACAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTAAGAAGCGag                   <  2:1485521/148‑1 (MQ=255)
                                                                                                                          tcgCTTTTCTAAATGTTGTACGGCGGGAAGCACTGGCTTTCACCGGCGGTTTGCAGTTCGTTGGCCCGTTTGTCAAAAACCATCCAATTACGCAACAGTGGGGGTAAATCGAATTTACGCATAGTAATTTATACATATAAGAAGCCag                   >  2:3002907/1‑148 (MQ=255)
                                                                                                                           cgCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTTTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTAATAAGCGagt                  >  1:1772172/1‑148 (MQ=255)
                                                                                                                           cgCTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTAAGAAGCGagt                  <  2:1668711/148‑1 (MQ=255)
                                                                                                                             cTTTTCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTAAGAAGCGagtaa                <  1:2005499/148‑1 (MQ=255)
                                                                                                                              ttttCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTAAGAAGCGagtaag               >  2:1028604/1‑148 (MQ=255)
                                                                                                                              ttttCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTAAGAAGCGagtaag               <  2:570426/148‑1 (MQ=255)
                                                                                                                              ttttCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCACTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTAAGAAGCGagtaat               >  2:659834/1‑147 (MQ=255)
                                                                                                                                ttCTCAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTAAGAAGCGagtaagta             <  2:506844/148‑1 (MQ=255)
                                                                                                                                    cAATGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTAAGAAGCGAGTAAGTACCTg         >  2:1932002/1‑148 (MQ=255)
                                                                                                                                      aTGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTAAGAAGCGAGTAAGTACCTGCa       >  2:641138/1‑148 (MQ=255)
                                                                                                                                      aTGTTGTACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTAAGAAGCGAGTAAGTACCTGCa       >  1:2737372/1‑148 (MQ=255)
                                                                                                                                       tgttgtACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTAAGAAGCGAGTAAGTACCTGCaa      <  1:1756180/148‑1 (MQ=255)
                                                                                                                                       tgttgtACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTAAGAAGCGAGTAAGTACCTGCaa      <  1:1014964/148‑1 (MQ=255)
                                                                                                                                         ttgtACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTAAGAAGCGAGTAAGTACCTGCAAAt    <  2:1486303/148‑1 (MQ=255)
                                                                                                                                         ttgtACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCACTGCGTTGGCCAGTTTGTAAAAACCGATCAATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCAATTCTCCTTCAAAGAAGAGAGTAAGTACCTGCAAAt    >  1:3323501/1‑148 (MQ=255)
                                                                                                                                           gtACGGAGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTAAGAAGCGAGTAAGTACCTGCAAATcc  <  2:2449544/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
TTTTCTCTTCTTTTGGCTGCTCCGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTTCTCAATGTTGTACGGCGGGAAGCTCTGGCTTTCACCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTAAGAAGCGAGTAAGTACCTGCAAATCC  >  CP009273/3860014‑3860300

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: