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,6530CA99.5% 704.2 / ‑4.3 213E44* (GAA→TAA) ibpAheat shock chaperone
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base A (109/103);  total (109/104)
Fisher's exact test for biased strand distribution p-value = 4.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 213 total aligned reads displayed.

TCAGCGATGCCCTGGTACAGATAGGTGCGCTCTTTTTGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTCGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTG  >  CP009273/3860508‑3860797
                                                                                                                                                 |                                                                                                                                                
tCAGCGATGCCCTGGTACAGATAGGTGCGCTCTTTTTGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGt                                                                                                                                                <  1:1741711/148‑1 (MQ=255)
   gCGATGCCCTGGTACAGATAGGTGCGCTCTTTTTGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTa                                                                                                                                             <  1:2198859/148‑1 (MQ=255)
    cGATGCCCTGGTACAGATAGGTGCGCTCTTTTTGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTAc                                                                                                                                            >  2:1127566/1‑148 (MQ=255)
    cGATGCCCTGGTACAGATAGGTGCGCTCTTTTTGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTAc                                                                                                                                            >  2:1127524/1‑148 (MQ=255)
       tGCCCTGGTACAGATAGGTGCGCTCTTTTTGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAg                                                                                                                                         <  2:628512/148‑1 (MQ=255)
        gCCCTGGTACAGATAGGTGCGCTCTTTTTGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGt                                                                                                                                        <  2:3326799/148‑1 (MQ=255)
                 cAGATAGGTGCGCTCTTTTTGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTta                                                                                                                               <  2:1782234/148‑1 (MQ=255)
                 cAGATAGGTGCGCTCTTTTTGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTta                                                                                                                               <  2:1608718/148‑1 (MQ=255)
                 cAGATAGGTGCGCTCTTTTTGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTta                                                                                                                               <  1:3241547/148‑1 (MQ=255)
                  aGATAGGTGCGCTCTTTTTGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTtat                                                                                                                              <  2:272732/148‑1 (MQ=255)
                   gATAGGTGCGCTCTTTTTGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTtata                                                                                                                             <  1:168545/148‑1 (MQ=255)
                     tAGGTGCGCTCTTTTTGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACg                                                                                                                           >  1:3071083/1‑148 (MQ=255)
                      aGGTGCGCTCTTTTTGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAATGTTATACgg                                                                                                                          <  1:3106462/148‑1 (MQ=255)
                      aGGTGCGCTCTTTTTGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACgg                                                                                                                          <  1:2941216/148‑1 (MQ=255)
                        gTGCGCTCTTTTTGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAg                                                                                                                        <  1:55048/148‑1 (MQ=255)
                        gTGCGCTCTTTTTGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAg                                                                                                                        <  1:3819738/148‑1 (MQ=255)
                        gTGCGCTCTTTTTGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAg                                                                                                                        <  1:2758337/148‑1 (MQ=255)
                              tctTTTTGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAgc                                                                                                                  >  1:1545380/1‑148 (MQ=255)
                              tctTTTTGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAgc                                                                                                                  >  1:2869242/1‑148 (MQ=255)
                               ctTTTTGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAgcc                                                                                                                 >  2:146288/1‑148 (MQ=255)
                                tttttGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAgccg                                                                                                                >  1:774096/1‑148 (MQ=255)
                                 ttttGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAgccgc                                                                                                               <  1:3700721/148‑1 (MQ=255)
                                  tttGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAgccgcc                                                                                                              >  1:1468495/1‑148 (MQ=255)
                                   ttGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCa                                                                                                             <  2:3224384/148‑1 (MQ=255)
                                       tcgtcgGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTAc                                                                                                         >  2:3776444/1‑148 (MQ=255)
                                         gtcgGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATAActc                                                                                                       <  2:165137/148‑1 (MQ=255)
                                          tcgGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTActct                                                                                                      <  2:3083748/148‑1 (MQ=255)
                                           cgGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTActctg                                                                                                     >  2:1292304/1‑148 (MQ=255)
                                              cGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTActctggct                                                                                                  >  1:2988076/1‑148 (MQ=255)
                                                   gCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGgt                                                                                             <  2:1956525/148‑1 (MQ=255)
                                                     aCCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGgttg                                                                                           <  1:244074/148‑1 (MQ=255)
                                                     aCCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGgttg                                                                                           <  2:2388945/148‑1 (MQ=255)
                                                      ccTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGgttgt                                                                                          <  2:1491630/148‑1 (MQ=255)
                                                       cTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGgttgtt                                                                                         >  1:3555707/1‑148 (MQ=255)
                                                         ttCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGtttt                                                                                       <  2:595599/148‑1 (MQ=255)
                                                         ttCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGtttt                                                                                       <  2:3480850/148‑1 (MQ=255)
                                                          tCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTc                                                                                      <  1:3498234/148‑1 (MQ=255)
                                                            accaccAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCtt                                                                                    <  1:2026899/148‑2 (MQ=255)
                                                             ccaccaGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTaa                                                                                   >  1:1663820/1‑148 (MQ=255)
                                                             ccaccaGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTaa                                                                                   >  2:3216204/1‑148 (MQ=255)
                                                               accaGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGt                                                                                 <  2:2802331/148‑1 (MQ=255)
                                                                  agcagATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGt                                                                              >  2:1791943/1‑148 (MQ=255)
                                                                      gATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTaaa                                                                          >  1:2244529/1‑148 (MQ=255)
                                                                         tATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTaaacaa                                                                       >  1:662016/1‑148 (MQ=255)
                                                                              tGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGt                                                                  >  1:1505539/1‑148 (MQ=255)
                                                                              tGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGt                                                                  >  2:2128843/1‑148 (MQ=255)
                                                                                 gCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCaa                                                               <  1:3226156/148‑1 (MQ=255)
                                                                                    gTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATc                                                            >  2:1883645/1‑148 (MQ=255)
                                                                                    gTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATc                                                            >  2:1124445/1‑148 (MQ=255)
                                                                                       aTTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCaa                                                         >  2:1248804/1‑148 (MQ=255)
                                                                                       aTTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCaa                                                         >  1:3086015/1‑148 (MQ=255)
                                                                                        tttCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAAt                                                        <  1:1845609/148‑1 (MQ=255)
                                                                                        tttCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAAt                                                        <  2:2689386/148‑1 (MQ=255)
                                                                                         ttCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATa                                                       >  1:3296351/1‑148 (MQ=255)
                                                                                         ttCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGGTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATa                                                       >  1:2248174/1‑148 (MQ=255)
                                                                                             aGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAg                                                   >  1:2151933/1‑148 (MQ=255)
                                                                                               ttCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGaa                                                 >  1:1632119/1‑148 (MQ=255)
                                                                                               ttCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGaa                                                 >  2:1753529/1‑148 (MQ=255)
                                                                                                     ctcAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTaa                                           >  2:2256777/1‑148 (MQ=255)
                                                                                                     ctcAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTaa                                           <  2:2312374/148‑1 (MQ=255)
                                                                                                     ctcAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTaa                                           >  1:431596/1‑148 (MQ=255)
                                                                                                         gCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCg                                       <  2:2999186/148‑1 (MQ=255)
                                                                                                            aaaCCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCgggg                                    <  1:2694606/148‑1 (MQ=255)
                                                                                                            aaaCCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCgggg                                    <  2:1663820/148‑1 (MQ=255)
                                                                                                             aaCCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGa                                   >  2:3004125/1‑148 (MQ=255)
                                                                                                             aaCCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGa                                   >  1:1337731/1‑148 (MQ=255)
                                                                                                             aaCCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGa                                   >  1:2128599/1‑148 (MQ=255)
                                                                                                               ccAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATa                                 >  1:1425729/1‑148 (MQ=255)
                                                                                                                cAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATaa                                <  2:805838/148‑1 (MQ=255)
                                                                                                                cAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATaa                                <  2:3640554/148‑1 (MQ=255)
                                                                                                                cAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATaa                                >  1:534440/1‑148 (MQ=255)
                                                                                                                cAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATaa                                <  2:2245877/148‑1 (MQ=255)
                                                                                                                cAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATaa                                <  2:1881343/148‑1 (MQ=255)
                                                                                                                 aGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATaaa                               <  2:2499082/148‑1 (MQ=255)
                                                                                                                 aGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATaaa                               <  2:2869242/148‑1 (MQ=255)
                                                                                                                 aGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATaaa                               <  2:1545380/148‑1 (MQ=255)
                                                                                                                  gCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAAt                              <  1:1410247/148‑1 (MQ=255)
                                                                                                                  gCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAAt                              <  2:3107408/148‑1 (MQ=255)
                                                                                                                       aGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAg                         >  2:1595878/1‑148 (MQ=255)
                                                                                                                        gCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGt                        <  1:1143909/148‑1 (MQ=255)
                                                                                                                         cGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGtt                       >  2:2326358/1‑148 (MQ=255)
                                                                                                                            tAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACg                    >  2:1185056/1‑148 (MQ=255)
                                                                                                                            tAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACg                    >  2:3295470/1‑148 (MQ=255)
                                                                                                                            tAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACg                    >  2:1185429/1‑148 (MQ=255)
                                                                                                                             aGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGc                   >  2:704875/1‑148 (MQ=255)
                                                                                                                                aaTGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATa                <  1:176814/148‑1 (MQ=255)
                                                                                                                                aaTGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATa                <  1:683002/148‑1 (MQ=255)
                                                                                                                                  tGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAAt              <  2:388032/148‑1 (MQ=255)
                                                                                                                                    cGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCa            <  2:3452754/148‑1 (MQ=255)
                                                                                                                                      gTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAAt          <  1:3652003/148‑1 (MQ=255)
                                                                                                                                       tAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATa         >  2:156110/1‑148 (MQ=255)
                                                                                                                                          tGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCt      >  1:2238933/1‑148 (MQ=255)
                                                                                                                                          tGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCt      <  1:1563928/148‑1 (MQ=255)
                                                                                                                                          tGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCt      >  1:401278/1‑148 (MQ=255)
                                                                                                                                           ggTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTc     <  1:3353383/148‑1 (MQ=255)
                                                                                                                                           ggTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTc     >  1:1705263/1‑148 (MQ=255)
                                                                                                                                           ggTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTc     >  1:1594444/1‑148 (MQ=255)
                                                                                                                                            gTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTcc    <  2:3595381/148‑1 (MQ=255)
                                                                                                                                             ttttAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCt   <  1:2041728/148‑1 (MQ=255)
                                                                                                                                              tttAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTg  >  2:2167152/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
TCAGCGATGCCCTGGTACAGATAGGTGCGCTCTTTTTGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTCGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTG  >  CP009273/3860508‑3860797

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: