Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 161,068 A→C 9.0% intergenic (+47/‑149) hrpB → / → mrcB putative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273161,0680AC 9.0% 372.4 / 7.8 191intergenic (+47/‑149)hrpB/mrcBputative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base A (75/94);  new base C (7/15);  total (82/109)
Fisher's exact test for biased strand distribution p-value = 3.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 193 total aligned reads displayed.

TGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCA  >  CP009273/160922‑161217
                                                                                                                                                  |                                                                                                                                                        
 gtggTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAg                                                                                                                                                      <  1:1052616/148‑1 (MQ=255)
 gAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAg                                                                                                                                                      >  1:5152491/1‑150 (MQ=255)
  aGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCGgc                                                                                                                                                     >  1:1426535/1‑150 (MQ=255)
    gTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAgcgc                                                                                                                                                   >  1:3937381/1‑150 (MQ=255)
      gCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAgcgcgc                                                                                                                                                 >  2:3384213/1‑150 (MQ=255)
        aaaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCAc                                                                                                                                               >  2:616710/1‑150 (MQ=255)
            agagATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACgggg                                                                                                                                           >  2:5029789/1‑147 (MQ=255)
                 ttaaaGGGCTTTATCCCAAAAATGTCGGGCCGGCCGACCGGGCAAATCCTGCACCGACGCGAGGTACTAAAATTTATTCGTATGTTTGTTTTTCCCCCTCTCCGTGTGGGTGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTtcaccct                                                                                                                                      <  1:267760/148‑1 (MQ=255)
                     tgggCGTTATCCCAACCATGTCTGGCCGGACGACCCGGCAAATACTGCACCTACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGGGGGAGAGGGCCGGGGTGTGGGCATCAGCGCGCACGTtcaccctcacc                                                                                                                                  <  2:193693/149‑1 (MQ=255)
                      ggggGTTATCCCAAACATGTCTGGCCGTACGACCCGTAAATTACTGCACCTACGCGACGTACTAAAATTTATTCTTATGTTTTTTTTACTCCATTTCCCTGGGGGTGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTtcaccctcaccc                                                                                                                                 <  1:1751441/150‑1 (MQ=255)
                       ggatttATCCCATCCAGGTCTGGCCGGACGACCCGGCAAATTCTGCTCCGACTCGATGGACGAAAAAGTATTCGTAAGGTTGTTTTTGTCCCTCTCGCTGTGGTAGAGGGCCGGGGTGTGGGCATCAGCGCGCACGTTcacactcaccct                                                                                                                                <  2:1556202/146‑1 (MQ=255)
                       ggCTTTATCCAAAAATTGTCTGGCCGGCCGACCCGGAAAATACTGCACCTACGCAACGTACTAAAATTTATTCTTAAGGTGGTTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGTGGGCATCAGCGCGCACGTTcaccctcaccct                                                                                                                                <  1:1636152/150‑1 (MQ=21)
                       ggCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtgacccgcaccca                                                                                                                                >  1:1553850/1‑149 (MQ=255)
                       ggCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGgctccctcacccg                                                                                                                                >  2:3195280/1‑149 (MQ=255)
                          gTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACAACCaaaa                                                                                                                             >  1:4278211/1‑149 (MQ=255)
                            tATCCAAAACTTTTCTGGCGGGACCCCCCGCAAAATCCTGCACCGACGCGACGGACTAAAATTTATTCGTATGGTTGTTTTTCTCCCTCTCCCTGGGGGAGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAAccc                                                                                                                           <  1:872802/150‑1 (MQ=255)
                                ccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAATGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCAGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTAACCCACACCCAAACCCactc                                                                                                                       >  2:4880300/1‑150 (MQ=38)
                                ccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCGCCCACACCCCAAACCactc                                                                                                                       >  1:2756719/1‑150 (MQ=38)
                                   aaCATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGGCCCGCACCCAAAcccaaacccg                                                                                                                    >  1:2527954/1‑143 (MQ=38)
                                   aaCATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCGGACCCTAACCCACTccca                                                                                                                    >  1:3209367/1‑149 (MQ=255)
                                      aTGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCGCACCCTAACCCACTCCCACaa                                                                                                                 >  1:1230215/1‑150 (MQ=255)
                                          cTGGCCGGACGACCCGGCAAATCCTCCACCGCCGCGACGGACGAAAATGTATTCGTATGGTTGGTTTACACCCTCCCCCTGTGGTTGCGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGga                                                                                                             <  1:3157579/150‑1 (MQ=255)
                                          cTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGACCCGCACCCGAACCCACACCCACAAGGga                                                                                                             >  1:1767924/1‑150 (MQ=38)
                                          cTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCCCACCCACACCCACACCCCCAAaggg                                                                                                             >  1:509332/1‑149 (MQ=37)
                                           tGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTTCACCCTCACCCGAACCCACACCCTCAAggggg                                                                                                            >  2:174463/1‑148 (MQ=255)
                                           tGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCGCCCGCACCCAAACCCACTCCCTCAAGGgag                                                                                                            >  1:1620346/1‑150 (MQ=38)
                                              ccGGACGACCCGGCAAATACTGCACCGACGCGACGTACGAAAAATTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGAGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAgg                                                                                                         <  1:3487768/150‑1 (MQ=255)
                                                 gccgacCCGCAAAATCCTCACCCGACGCGCCGTACTAAAATTTTTTCGTAGGTTTGTTTTACTCCCTCTCCTTGTGGGTGTGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGa                                                                                                      <  1:190722/148‑1 (MQ=255)
                                                 gacgacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCTAACCCACACCCACAAGGGTGAggggt                                                                                                      >  2:4694887/1‑149 (MQ=38)
                                                  actacCCGGCAAATCCGGCACCTCCGCGACGGACTAAAAAGTATTCGTACTGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGTGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGAc                                                                                                     <  1:4680987/147‑1 (MQ=37)
                                                  acgacCCGGCAAATACTGCACCGCCGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGTGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGAc                                                                                                     <  2:2014613/150‑1 (MQ=255)
                                                   cgacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCAACCAAAACCACACCCGCCACGGAGAGGGGGcc                                                                                                    >  1:3057974/1‑150 (MQ=37)
                                                      cccGGCAAATACCGCCCCGACGCGCCGGACGAACATGTATTCGTACGGTTGGTTTTCGCCCTCCCCCGGTGGGAGTGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGAt                                                                                                 <  1:2816234/150‑1 (MQ=255)
                                                       ccGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGAGACCGATc                                                                                                <  2:2996722/150‑1 (MQ=255)
                                                       ccGGCAAATACTCCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATc                                                                                                <  1:4285562/150‑1 (MQ=255)
                                                        cgtcAAACACTGCACCGACGCGACGGACGAAACAGTATTCGTAAGGTTGGTTTTCGCCCTCTCCCTGTGGGAGAGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCg                                                                                               <  2:4942762/147‑1 (MQ=39)
                                                         ggCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGa                                                                                              <  2:4503687/150‑1 (MQ=255)
                                                         ggCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCGCCCGCACCCGCACCCCCACCCCCACGGGGGGGGGGGcccggggc                                                                                              >  1:3092957/1‑144 (MQ=25)
                                                          gCAAATACTGCACCGACGCGCCGGACTAAAAATTATTCTTACGGTTGGTTTTCTCCCTCTCCCTGTGGGAGCGGGCCGGGTTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAg                                                                                             <  2:1828729/150‑1 (MQ=38)
                                                          caaaaTACTGCACCGACGCGACGGACTAAAAATTATTCGTAAGGTTGGTTTTCTCCCTCTCCCGGTGGGAGTGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAg                                                                                             <  2:1828716/148‑1 (MQ=255)
                                                            aaaTACTGCACCGACGCGACGGACTAAAAAGTATTCGAACGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑ca                                                                                        <  2:1753768/150‑1 (MQ=255)
                                                            aaaTACGCCACCGACGCTCCGGACGAAAAATTATTCTTATGTTTGGTTTCCGCCCTCTCCGTGTGGGTGTGGGCCGGGGTGAGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑ca                                                                                        <  1:2397437/150‑1 (MQ=255)
                                                             aaTACTGCACCTACGCGCCGGACTAAAAAGTTTTCGTATGGTTGGTTTACTCCCTATCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑cac                                                                                       <  2:4854313/150‑1 (MQ=38)
                                                             aaTACTGCACCGACGCGACGGACGAAATAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCGCCCCCTAACCCTCTCCCGCAAGGGGGGGGGGACCGAGAGAG‑‑cac                                                                                       >  2:195485/1‑150 (MQ=37)
                                                             aaTACTCCACCGACGCGACGGACTAAAAAGTATTCGTACGGTTGGTTTTATCCCTCTCCCTGGGGGAGAGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑cac                                                                                       <  2:1391896/150‑1 (MQ=255)
                                                             aaTACGGCACCGACGCGCGGGACGAAATAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGGGGGAGAGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑cag                                                                                       <  2:2255796/150‑2 (MQ=255)
                                                                   gCACCCACGCGACGGACGAAACAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAACAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAAttt                                                                                 <  2:1251798/150‑1 (MQ=255)
                                                                     aCCGACGCGACGGACGACAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTg                                                                               <  1:3319068/150‑1 (MQ=255)
                                                                       cgacgcgacgGACTAAATTGTTTTCGTTTGGTTTGTTTTCTCCCTAGCCCTGTGGGTTTGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTgag                                                                             <  1:374483/150‑1 (MQ=255)
                                                                        gacgcgacgGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCTAACCCGATCCCTCAAGGGGGAGGGGACCGAACGAG‑‑CACAACTGTTgagc                                                                            >  2:4341689/1‑149 (MQ=37)
                                                                        gacgcgacgGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCAAACCCACTCCCGCAAGGGAGAGGGGACCCATCGAG‑‑CACAAATTTTgaga                                                                            >  2:4341609/1‑150 (MQ=38)
                                                                         acgctacgTACTAAAATGTTTTCGTAAGGTTGGTTTTCCCCCGCGCCCTGTGGGAGCGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTgagag                                                                           <  1:4152733/150‑1 (MQ=255)
                                                                          cgcgacgGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTgagaga                                                                          <  2:1328034/150‑1 (MQ=255)
                                                                           gcgacgTACTAAAATTTATTGTTATGGTTGGTTTTCTCCCTCTCCCTGTGGTTGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGat                                                                         <  1:1126741/150‑1 (MQ=37)
                                                                               cggacgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGAACGTTCAACCCCACACCAACACACACCCAAAACGGGGAGGGGGCCCATAGAGAACACAATTTTTgagtaaact                                                                       >  2:2130329/1‑144 (MQ=16)
                                                                                ggacgAAAAAGTTTTCGTAAGGTTGGTTTACTCCCTCTCCCGGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATAtctt                                                                    <  1:1750658/150‑1 (MQ=255)
                                                                                 gacgAAAATGTATTCGTAAGGTTGGTTTTCTCCCTCCCCGGGGGGGAGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATAtcttc                                                                   <  1:2163958/150‑1 (MQ=255)
                                                                                 gacgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGACCCACACCACAAACCAAAACCACAAAGGAGAGGGGACCCAAAGAG‑‑CAACAAAATTTAGTGATAtcata                                                                   >  1:721177/1‑147 (MQ=16)
                                                                                 gacgAAAAAGTAGTCGTAAGGGTGGTTTTCTCCCTCTCCCGGTGGGTGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATAtcttc                                                                   <  2:226721/150‑1 (MQ=255)
                                                                                    gAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCGAACCCACACCCACAAAGGAGAGGGGACCCATCGAG‑‑CACAAAATTTGATAGatataatatta                                                                >  2:377091/1‑143 (MQ=25)
                                                                                    gAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAAGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCGCACCCAAACCCACACCCACAAGGGAGAGGGGACCGATCCAC‑‑CACAAATTTTGAGAGATAacttcatc                                                                >  2:1514696/1‑150 (MQ=37)
                                                                                     aaaaaGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCAAACCCGAACCCACTCCCAAAAAGGATAGGGGACCGAACGAG‑‑CACAAATTTTGAGAGATAtattcttaa                                                               >  1:4467205/1‑148 (MQ=25)
                                                                                       aaaGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCAACCAAAACCAAAACCACACCCAAAAAGGAGATGGGACCGATAGAC‑‑AACAAAATTTTATATATAtattattattt                                                             >  2:4438563/1‑146 (MQ=12)
                                                                                        aaGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTc                                                            <  2:2558853/150‑1 (MQ=255)
                                                                                        aaGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACCCCCAAACCCTCTCCCACAAGGGAGAGGGGCCCGATCGAG‑‑CCCAAAGTTTGAGAGATATATTCTTCTGta                                                            >  2:2999250/1‑149 (MQ=37)
                                                                                         aGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGGTGACCCACCCCCGAAACCAATACCTCAAGGGGGAGGGGACCCATAGGGCCCCAAATTTTGGAGATAtattcttagtggg                                                           >  2:4784844/1‑144 (MQ=16)
                                                                                         aGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCt                                                           <  2:2977902/150‑1 (MQ=255)
                                                                                          ttATTCGTAAGGTTGGTTTTCTCCCTCTCCCGGGGGGAGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCtt                                                          <  1:1569504/149‑1 (MQ=255)
                                                                                            aTTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACAACCAAACCCACTCCCACAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGATATATATTAATCTGTCTTGt                                                        >  2:806010/1‑150 (MQ=37)
                                                                                            aTTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCAACCGCACCCAAAACCTCTCCCACAAGGGAGAGGGGACCGATCGAG‑‑CACAAAATTTTAGAGATATATTATTATGTCTTGt                                                        >  2:1523956/1‑150 (MQ=37)
                                                                                               cGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAAc                                                     <  2:3937381/150‑1 (MQ=255)
                                                                                               cGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCAAACCCACTCCCACAAGGGGGAGGGGACCGATCGAG‑‑CACAAATTTTTTGATATATATTATTAAGTATTTTAAc                                                     >  1:2881601/1‑150 (MQ=25)
                                                                                                gTACGGTTGGTTTTCCCCATCTCCCTGGGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACa                                                    <  1:796322/150‑1 (MQ=255)
                                                                                                gTACGGTTGGTTTTCACCCTCTCCGTGGGGGTGTGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACa                                                    <  1:796549/150‑1 (MQ=255)
                                                                                                      ttggttTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCAAAACCTATACCGCAAGGGTGAGGGGACCGATCGAG‑‑CCAAAATTTTTAGAGATATCTTCTTATGTATTTTAAAAgaataa                                              >  2:1287431/1‑150 (MQ=25)
                                                                                                        ggttTACTCCCTCTCCCTGTTGTAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACa                                            <  2:188754/150‑1 (MQ=255)
                                                                                                              ctccctctccctGGGGGTGTGGTCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAAt                                      <  2:77583/150‑1 (MQ=255)
                                                                                                               tccctcaccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATc                                     <  2:1904190/150‑1 (MQ=255)
                                                                                                                ccctctccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCCCGGTCACCCGCACCCAAACCCACACCCTCAAGGGAGTGGGGTCCGATCGCT‑‑CCCAAATGTTGTGAGATATATTAGTCTGTCGTTTAAAAGAAGAACAGaaaaaac                                    >  2:3423446/1‑147 (MQ=19)
                                                                                                                    ctccctGTGGGATCGGGCCGGGGGGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCt                                <  2:2406959/150‑1 (MQ=255)
                                                                                                                      ccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAATAACAGAAAATCGGGCttt                              <  1:1589424/150‑1 (MQ=255)
                                                                                                                      ccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCGCACCCAAAACCACACCCACAAGGGGGAGGGGACCCTGCGCG‑‑CACAAAATTTTCGTGAAAAAATATTATGTCATGTAACAGTAGACAAGAAAATAGGGCgtt                              >  2:4378399/1‑150 (MQ=12)
                                                                                                                       cctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCCCACCCAAAACCACTCCCACAAGGGGGCGGGGCCCCTTCCTG‑‑CACAAAATTTGTGAGATATCATATTATGTCTTGGAACATAAGTACACAAAATCGGGCgttt                             >  2:4049702/1‑150 (MQ=19)
                                                                                                                       cctGGGGGAGAGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCtttt                             <  2:509332/150‑1 (MQ=255)
                                                                                                                         tgtgGGAGTGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAATAACAGAAAATCGGGCTTTTgc                           <  1:1454952/150‑1 (MQ=255)
                                                                                                                           tgGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGACCCACACCCAAAACCTCACCCGCAAGGGGGTGGGGGCCCCTCGGG‑‑CACAAATTTTTTGAGGTATATTATTCTGTGTTGTAAAAGCAGAAAAGCAAATCGGGCTTTTtcgg                         >  2:1905728/1‑149 (MQ=255)
                                                                                                                               agtgggCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGa                     <  1:4906747/147‑1 (MQ=255)
                                                                                                                                 agGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAat                   <  2:2962018/150‑1 (MQ=255)
                                                                                                                                  gggCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCACGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAata                  <  2:4115143/150‑1 (MQ=255)
                                                                                                                                  gggCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAata                  <  2:2066970/150‑1 (MQ=255)
                                                                                                                                   ggCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAatat                 <  1:820246/150‑1 (MQ=255)
                                                                                                                                   ggCCGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCGACCCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAatat                 <  1:3566270/150‑1 (MQ=255)
                                                                                                                                   ggCCGGGGTGAGGGCATCCGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAatat                 <  1:3566043/150‑1 (MQ=255)
                                                                                                                                   ggCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAatat                 <  2:1668328/150‑1 (MQ=255)
                                                                                                                                   ggCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAatat                 <  2:1620346/150‑1 (MQ=255)
                                                                                                                                    gCCGGGGGGCGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATAtt                <  2:3693648/150‑1 (MQ=255)
                                                                                                                                     ccGGGGTGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCGCCCACAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTg               <  1:3959733/150‑1 (MQ=255)
                                                                                                                                       ggggTGAGGGCATCAGCGCGCACGTTCACCCACACCCAAACCCTATCCCACAAGGGAGAGGGGACCCATCGAG‑‑CACAAAATTTGAGAGATAAATTATTATGTCATTTAAAACAAGAACATAAAAACAGGCTTTTTCGCCAGAATATTTCg             >  2:4146565/1‑150 (MQ=255)
                                                                                                                                              gggCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaa      <  2:57680/150‑1 (MQ=255)
                                                                                                                                               gccATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaaa     <  1:1954365/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                        
TGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCA  >  CP009273/160922‑161217

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