Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,860,653 C→A 100% E44* (GAA→TAA)  ibpA ← heat shock chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,860,6530CA96.8% 571.5 / 9.0 186E44* (GAA→TAA) ibpAheat shock chaperone
Reads supporting (aligned to +/- strand):  ref base C (0/6);  new base A (75/105);  total (75/111)
Fisher's exact test for biased strand distribution p-value = 8.26e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01
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 186 total aligned reads displayed.

AGCGATGCCCTGGTACAGATAGGTGCGCTCTTTTTGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTCGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCT  >  CP009273/3860510‑3860796
                                                                                                                                               |                                                                                                                                               
aGCGATGCCCTGGTACAGATAGGTGCGCTCTTTTTGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCt                                                                                                                                             <  2:734175/148‑1 (MQ=255)
 gCGATGCCCTGGTACAGATAGGTGCGCTCTTTTTGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTa                                                                                                                                            <  1:165401/148‑1 (MQ=255)
 gCGATGCCCTGGTACAGATAGGTGCGCTCTTTTTGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTa                                                                                                                                            <  2:1768313/148‑1 (MQ=255)
  cGATGCCCTGGTACAGATAGGTGCGCTCTTTTTGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTAc                                                                                                                                           >  1:166316/1‑148 (MQ=255)
         cTGGTACAGATAGGTGCGCTCTTTTTGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCa                                                                                                                                    <  1:1115942/148‑1 (MQ=255)
          tGGTACAGATAGGTGCGCTCTTTTTGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCaa                                                                                                                                   <  2:3027371/148‑1 (MQ=255)
          tGGTACAGATAGGTGCGCTCTTTTTGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCaa                                                                                                                                   >  1:1405856/1‑148 (MQ=255)
            gTACAGATAGGTGCGCTCTTTTTGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACg                                                                                                                                 <  2:1638150/148‑1 (MQ=255)
            gTACAGATAGGTGCGCTCTTTTTGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACg                                                                                                                                 <  1:750497/148‑1 (MQ=255)
            gTACAGATAGGTGCGCTCTTTTTGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACg                                                                                                                                 <  1:750626/148‑1 (MQ=255)
            gTACAGATAGGTGCGCTCTTTTTGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACg                                                                                                                                 <  1:2666708/148‑1 (MQ=255)
                tgATTGGTGAGCTCTTTTTGTTCGTCGGCGTGATCTCCTTTCACCAGGAGCAGATTATGCTGGGAGGTAATTTACAGTTCTCTGTCAGCAAATCCAGCCTCAGCGATAGCAATGCTGTAATGGTTTTAGTCTACCTGTTCAACGTtat                                                                                                                             <  1:472214/147‑1 (MQ=255)
                aGATAGGTGCGCTCTTTTTGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTtat                                                                                                                             <  1:1539394/148‑1 (MQ=255)
                aGATAGGTGCGCTCTTTTTGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTtat                                                                                                                             <  1:472210/148‑1 (MQ=255)
                    aGGTGCGCTCTTTTTGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACgg                                                                                                                         >  1:2075587/1‑148 (MQ=255)
                       tGCGCTCTTTTTGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATATCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAgg                                                                                                                      <  2:578605/148‑1 (MQ=255)
                       tGCGCTCTTTTTGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAgg                                                                                                                      <  2:2935784/148‑1 (MQ=255)
                         cgcTCTTTTTGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGt                                                                                                                    <  2:32209/148‑1 (MQ=255)
                           ctctTTTTGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAg                                                                                                                  >  1:1179400/1‑148 (MQ=255)
                            tctTTTTGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAgc                                                                                                                 <  1:1851652/148‑1 (MQ=255)
                                 ttGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCa                                                                                                            >  2:2609267/1‑148 (MQ=255)
                                  tGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCAt                                                                                                           <  1:3097264/148‑1 (MQ=255)
                                     tcgtcgGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTAc                                                                                                        >  2:2032408/1‑148 (MQ=255)
                                      cgtcgGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTAct                                                                                                       <  2:390751/148‑1 (MQ=255)
                                      cgtcgGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTAct                                                                                                       <  2:2181850/148‑1 (MQ=255)
                                        tcgGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTActct                                                                                                     <  2:1483912/148‑1 (MQ=255)
                                         cgGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTActctg                                                                                                    <  2:1410409/148‑1 (MQ=255)
                                            cGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTActctggct                                                                                                 >  1:2991569/1‑148 (MQ=255)
                                             gTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTActctggctc                                                                                                >  1:1943019/1‑148 (MQ=255)
                                               gAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTActctggctctg                                                                                              >  2:1555015/1‑148 (MQ=255)
                                                aGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTActctggctctgg                                                                                             >  1:501689/1‑148 (MQ=255)
                                                  cACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGgtt                                                                                           >  2:2277698/1‑148 (MQ=255)
                                                   aCCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTCGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGgttg                                                                                          <  1:2442710/148‑1 (MQ=255)
                                                      tttCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGttt                                                                                       >  2:934871/1‑148 (MQ=255)
                                                       ttCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGtttt                                                                                      >  2:2647710/1‑148 (MQ=255)
                                                         caccacCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCt                                                                                    >  1:882731/1‑148 (MQ=255)
                                                         caccacCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCt                                                                                    <  1:531469/148‑1 (MQ=255)
                                                           ccaccaGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTaa                                                                                  >  1:1058973/1‑148 (MQ=255)
                                                             accaGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGt                                                                                <  2:3333088/148‑1 (MQ=255)
                                                              ccaGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTg                                                                               >  1:1685667/1‑148 (MQ=255)
                                                               cagcagATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTgg                                                                              <  1:1335419/148‑1 (MQ=255)
                                                               cagcagATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTgg                                                                              <  2:2058999/148‑1 (MQ=255)
                                                                       tATCCTGGGCGGTAATTTACAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTaaacaa                                                                      <  2:2311802/148‑1 (MQ=255)
                                                                          ccTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACg                                                                   >  1:2076634/1‑148 (MQ=255)
                                                                              ggCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAGTGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCa                                                               >  1:212561/1‑148 (MQ=255)
                                                                              ggCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCa                                                               >  1:2723818/1‑148 (MQ=255)
                                                                              ggCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCa                                                               >  1:1808399/1‑148 (MQ=255)
                                                                                cGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTCGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCaaa                                                             <  2:1177917/148‑1 (MQ=255)
                                                                                cGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTCGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCaaa                                                             <  2:650350/148‑1 (MQ=255)
                                                                                cGGTAAGTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTCGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCaaa                                                             <  2:2435735/148‑1 (MQ=255)
                                                                                  gTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATc                                                           <  2:1738679/148‑1 (MQ=255)
                                                                                         ccAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGc                                                    <  2:1179400/148‑1 (MQ=255)
                                                                                           aGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAg                                                  <  1:2370295/148‑1 (MQ=255)
                                                                                             ttCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGaa                                                >  1:2859446/1‑148 (MQ=255)
                                                                                               cGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTAAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACg                                              <  1:1611903/148‑1 (MQ=255)
                                                                                                gCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACgg                                             <  1:3408822/148‑1 (MQ=255)
                                                                                                gCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACgg                                             >  1:2431101/1‑148 (MQ=255)
                                                                                                gCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACgg                                             <  2:2637173/148‑1 (MQ=255)
                                                                                                 ctctCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGt                                            <  2:206750/148‑1 (MQ=255)
                                                                                                 ctctCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGt                                            <  2:694886/148‑1 (MQ=255)
                                                                                                 ctctCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGt                                            >  2:988801/1‑148 (MQ=255)
                                                                                                  tctcAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTa                                           <  1:3187008/148‑1 (MQ=255)
                                                                                                   ctcAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTaa                                          <  2:1058973/148‑1 (MQ=255)
                                                                                                   ctcAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTaa                                          >  2:460207/1‑148 (MQ=255)
                                                                                                    tcAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTACCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTaaa                                         <  1:2636153/148‑1 (MQ=255)
                                                                                                     cAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAg                                        >  1:272165/1‑148 (MQ=255)
                                                                                                     cAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAg                                        <  2:1829066/148‑1 (MQ=255)
                                                                                                       gCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCg                                      >  1:1002947/1‑148 (MQ=255)
                                                                                                         aaaaCCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCggg                                    <  1:117568/148‑1 (MQ=255)
                                                                                                         aaaaCCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCggg                                    <  2:1777616/148‑1 (MQ=255)
                                                                                                          aaaCCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCgggg                                   <  2:493160/148‑1 (MQ=255)
                                                                                                           aaCCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGa                                  >  1:1137923/1‑148 (MQ=255)
                                                                                                           aaCCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGa                                  <  1:1358067/148‑1 (MQ=255)
                                                                                                           aaCCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGa                                  >  1:1383440/1‑148 (MQ=255)
                                                                                                           aaCCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGa                                  >  2:1756381/1‑148 (MQ=255)
                                                                                                           aaCCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGa                                  <  1:2032071/148‑1 (MQ=255)
                                                                                                            aCCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGAt                                 <  2:157770/148‑1 (MQ=255)
                                                                                                                gCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAAt                             <  1:3394853/148‑1 (MQ=255)
                                                                                                                gCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAAt                             <  2:796630/148‑1 (MQ=255)
                                                                                                                 ccACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATc                            >  2:2745751/1‑148 (MQ=255)
                                                                                                                    caGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCaaa                         <  2:585935/148‑1 (MQ=255)
                                                                                                                     aGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAg                        <  2:2437259/148‑1 (MQ=255)
                                                                                                                     aGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAg                        <  2:842026/148‑1 (MQ=255)
                                                                                                                     aGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAg                        <  2:3429882/148‑1 (MQ=255)
                                                                                                                      gCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGt                       <  1:3054359/148‑1 (MQ=255)
                                                                                                                         aTAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTAc                    >  2:2874146/1‑148 (MQ=255)
                                                                                                                          tAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACg                   <  1:2568490/148‑1 (MQ=255)
                                                                                                                            gCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCa                 >  1:2203038/1‑148 (MQ=255)
                                                                                                                            gCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCa                 >  1:2919860/1‑148 (MQ=255)
                                                                                                                              aaTGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATa               >  2:1057001/1‑148 (MQ=255)
                                                                                                                              aaTGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATa               <  1:398181/148‑1 (MQ=255)
                                                                                                                                    gTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAAt         >  2:2266285/1‑148 (MQ=255)
                                                                                                                                    gTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAAt         <  1:1019587/148‑1 (MQ=255)
                                                                                                                                      aaTGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCa           <  1:1492254/144‑1 (MQ=255)
                                                                                                                                        tGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCt     >  1:2097334/1‑148 (MQ=255)
                                                                                                                                        tGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCt     <  2:558323/148‑1 (MQ=255)
                                                                                                                                         ggTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTc    <  2:2674667/148‑1 (MQ=255)
                                                                                                                                         ggTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTc    >  1:2815414/1‑148 (MQ=255)
                                                                                                                                           ttttAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCt  >  1:2680918/1‑148 (MQ=255)
                                                                                                                                           ttttAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCt  >  1:3457472/1‑148 (MQ=255)
                                                                                                                                               |                                                                                                                                               
AGCGATGCCCTGGTACAGATAGGTGCGCTCTTTTTGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTCGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCT  >  CP009273/3860510‑3860796

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