Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 161,068 A→C 11.8% 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,0680AC11.8% 234.3 / 3.1 127intergenic (+47/‑149)hrpB/mrcBputative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base A (45/64);  new base C (4/14);  total (49/78)
Fisher's exact test for biased strand distribution p-value = 1.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 129 total aligned reads displayed.

CCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGA‑‑G‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAG  >  CP009273/160919‑161215
                                                                                                                                                     |                                                                                                                                                       
ccGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCa                                                                                                                                                         >  1:747981/1‑150 (MQ=255)
ccGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCa                                                                                                                                                         >  1:1873396/1‑150 (MQ=255)
  gTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATc                                                                                                                                                       >  2:2817096/1‑150 (MQ=255)
   tGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCa                                                                                                                                                      >  2:907760/1‑150 (MQ=255)
     aGGTGCAAAAAGAGATGAAAGGGGGTTATCCCAAACATGTCTGGCCGGACGACCCGGAAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAgc                                                                                                                                                    <  2:1874194/150‑1 (MQ=255)
     aGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAgc                                                                                                                                                    <  2:3125272/150‑1 (MQ=255)
        tGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTACGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAgcgcg                                                                                                                                                 >  2:3225351/1‑150 (MQ=255)
           aaaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCAc                                                                                                                                              >  1:3244126/1‑150 (MQ=255)
           aaaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCAc                                                                                                                                              >  2:3900168/1‑150 (MQ=255)
            aaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACg                                                                                                                                             >  1:3224286/1‑150 (MQ=255)
            aaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACg                                                                                                                                             >  1:366857/1‑150 (MQ=255)
             aaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACgg                                                                                                                                            >  1:2111686/1‑149 (MQ=255)
             aaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACgg                                                                                                                                            >  1:3701404/1‑149 (MQ=255)
             aaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGt                                                                                                                                            >  2:2912428/1‑150 (MQ=255)
              aaGTGATGAAAGGGCGTTATCCAAAACTGGTCTGGCCGGACTACCCGGCAAATACTGCACCGCCGCGAGGGACGAAAAAGTATTCGGAAGGTTGGTTTTCTCCCTCTCCCGGGGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACggt                                                                                                                                           <  2:3078575/150‑3 (MQ=37)
                gagaTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGtgga                                                                                                                                         >  2:3444829/1‑147 (MQ=255)
                 agaTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGtgggg                                                                                                                                        >  1:1059825/1‑146 (MQ=255)
                 agaTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGtgggc                                                                                                                                        >  1:1988465/1‑146 (MQ=255)
                   aTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTggaccc                                                                                                                                      >  1:3535884/1‑150 (MQ=255)
                   aTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTggaccc                                                                                                                                      >  2:1184038/1‑150 (MQ=255)
                    tGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGTGGGAATCAGCGCGCACGTtaaccca                                                                                                                                     >  2:2178521/1‑149 (MQ=255)
                    tGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTgcaaccc                                                                                                                                     >  2:1493000/1‑149 (MQ=255)
                    tGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtgacccg                                                                                                                                     >  1:1022751/1‑149 (MQ=255)
                    tGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtgacccc                                                                                                                                     >  1:1022896/1‑149 (MQ=255)
                      aaaGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGAATTCGTAACGTAGGTTTTCTCCCTCTCCCTGTGGTAGAGGGCCGGGGTGAGGGCATCAGCGAGAACGTtgacccacc                                                                                                                                   >  1:2523865/1‑147 (MQ=38)
                       aaGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCCCGTtcaccccccc                                                                                                                                  >  2:1800085/1‑146 (MQ=255)
                        cgggCGTTTTCCAAAAAAGGTCTGGCCGTACGACCCGGCAAATACTGCCCCTCCGCGCCGGACTAAAAATTTTTCGTATGTTTGGTTTTAGCCCTCTCCCTGTGGTTGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccctcacc                                                                                                                                 <  2:905071/149‑1 (MQ=25)
                           gCGTTATCCCAAACATGTCTGGCGGGACGACCCGGCAAATCCTGCACCGACGCGACGGACGAAAAAGTATTCTTACGTTTGGTTTTCTCCCTCGCCCTGTGGGTGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTa                                                                                                                              <  2:2139288/150‑1 (MQ=38)
                             gTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGGCCCACACCCCAAc                                                                                                                            >  2:403412/1‑150 (MQ=38)
                               tATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATGCGAGCGGTAGGATTTCTCCCTCTCACTCTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCAGGTGCGCCGACCCCCAAAccc                                                                                                                          >  2:291385/1‑150 (MQ=25)
                                acccaacaCATGTCTGTCCGTACGACCCGTCAAATACTGCACCTACGCGACGTACCAAAAAGTATTCTTAAGGTTGTTTTTCTCCCTTTCCCTGTGGTAGTGGGCGGGGGTTAGGGCATCAGCGCGCACGGTCACCCTCACCCTAACCct                                                                                                                         <  2:266391/143‑1 (MQ=255)
                                aTCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGGCCCACACCCCAAcccc                                                                                                                         >  1:2067068/1‑149 (MQ=38)
                                aTCCCAAACATGTCTGGCCGGACGACCCGGAAAATACTGCACCGACGCGACGGACTAAATTGTATTCGTATGTTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCct                                                                                                                         <  2:533306/150‑1 (MQ=38)
                                 tCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGTCACCCACACCCAAAaccac                                                                                                                        >  1:3464790/1‑148 (MQ=255)
                                  cccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGGCCCCCACCCAAAcccaca                                                                                                                       >  1:222949/1‑147 (MQ=38)
                                   ccAAACATGTGGGGCCGGACGACCCGCCAAATACTGCACCGACGCGACGGACGAAAATTTATTCGTATGTTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCtctc                                                                                                                      <  2:3667539/150‑1 (MQ=37)
                                      aaattGTCTGGCCGGACGACCCGGAAAATCCTGCACCGCCGCGAGGCACGAAAAATTTTTGTTATGGTTGTTTTTCGCCCTCGCCCTGGGGGTGTGGGCGGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCt                                                                                                                   <  2:1935858/146‑1 (MQ=255)
                                       acttGTCTGGCCGGACGACCCGGCAAATCCTCAACCGCCGGGACGTACGAAAATTTTTTCGTATGTTTGGTTTTCTCCCTGTCCCGGTGGGATTGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCTCACCCTAACCCTCTCCCTc                                                                                                                  <  2:64665/147‑1 (MQ=25)
                                          tGTCGGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGTACGAAAACTTATTCTTAAGGTTGGTTTTCTCCCTCCCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                                               <  2:2646141/150‑1 (MQ=38)
                                              gggCCGGACGACCGCGCCAATCCTGCCCCGACGCGCCGGACGAAAACGTATGCGTAGGTGTGGTTTTCGCCCTCTCCCTGGGGGATAGGGCCGGGGCGAGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGgag                                                                                                           <  2:4071859/149‑1 (MQ=255)
                                               ggCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAATTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGgaga                                                                                                          <  2:3634138/150‑1 (MQ=255)
                                                 ccGGACGCCCCGCAAAATCCTCCCCCGCCGCGACGGCCGAAAAAGTATTCGTATGTTTGTTTTTCCCCCTCTCCCTGTGGGAGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAgg                                                                                                        <  1:3060765/150‑1 (MQ=255)
                                                   gtacgacCGCGCAACTACTGCACCCACGCGACGCACCAACATTTATTCGTATGCTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCGGGGGTGAGGGAATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAgggg                                                                                                      <  2:4021396/148‑1 (MQ=25)
                                                   ggccgacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTACGGTTGGTTTTCCCCCACTCCCTGTGGGAGAGGGCCGGGGTGCGGGCATCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAgggg                                                                                                      <  2:2375769/147‑1 (MQ=255)
                                                   ggACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCAACCTCAACCTAACCCACAACCACCAGGGAGAgggg                                                                                                      >  2:2249975/1‑150 (MQ=37)
                                                   ggACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCGCACCCAAACCCACACCCGCCAgggggggggg                                                                                                      >  2:480997/1‑150 (MQ=37)
                                                      cgacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTGCACCCGCACCCGAACCCACTACCGCAAGGGCGGGGGGGcc                                                                                                   >  2:2372432/1‑150 (MQ=37)
                                                      cgacCCGGCAAAGACGGCCCCGCCGCGCCGCCCGAAACAGTATTCGTACGGTTGGTTGACTCCCGCGCCCTGGGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGAcc                                                                                                   <  2:3092104/150‑1 (MQ=37)
                                                       gacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGCCCCGAACCCAAAACCACACCCGCAAGGGGGGGGGGGccc                                                                                                  >  1:308239/1‑149 (MQ=25)
                                                          ccGGCCACTCCTGCACCGCCGCCACGGACGAAAAAGTATTCGCAAGGGTGGTGTTCGCCCCCTCCCGGGGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATc                                                                                               <  2:2626686/150‑1 (MQ=255)
                                                          ccGGCAAATACTGCACCGACGCGCCGTACTAAAAAGTATTCTTATGGTTTGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATc                                                                                               <  2:4129016/150‑1 (MQ=38)
                                                          ccGCCAATTCCTGCCCCGCCGCGCCGCACTACAATGTATGCGTCTGGTGGTTTTCCCCCCCCGCCGTGGGGTTGTGGCCCGGGGTGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTANCCCTCTCCCTCAAGGGAGAGGGGACCGATc                                                                                               <  1:3859573/150‑1 (MQ=255)
                                                                 aTACCGCACCGACGCGACGGACGAAATAGTATTCGTAAGGTTGTTTTTCTCCCTCTCCCTGGGGGAGAGGCCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGA‑‑G‑‑caca                                                                                    <  2:3098448/150‑1 (MQ=39)
                                                                    cggCACCGACGCGAGGTACTAAAATTTATTCGTTCGTTTGTTTTACTCCCTCTCCGTGTGGGAGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGA‑‑G‑‑CACAAAt                                                                                 <  1:3719647/148‑1 (MQ=37)
                                                                     tGCACCCACGCGCCGGACGAAAATGTATTCGTATGGTTGTTTTTCTCCCTCGCCCTGTGGGTGAGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGA‑‑G‑‑CACAAAtt                                                                                <  1:4089260/150‑1 (MQ=255)
                                                                     gccACCGCCGCGCCGGACGAAAATTTATTCGTTTGTTTGTTTTTCTCCCTCTCCCTGTGGGTGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGA‑‑G‑‑CACAAAtt                                                                                <  1:1378106/148‑1 (MQ=37)
                                                                     cccACCGACGCGCCGTCCGAAAACGTATTCGTATGGTTGGTTTTCCCCCTCCCCGTGTGGGAGTGGCCCGGGGTTCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGA‑‑G‑‑CACAAAtt                                                                                <  2:3229202/148‑1 (MQ=255)
                                                                             cgcgacgTACGAAAAAGTATTCGTATGGTTGGTTTTCTCCCTATCCCTGTGGGTGTGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGA‑‑G‑‑CACAAATTTTgagaga                                                                        <  2:586526/150‑1 (MQ=255)
                                                                             cgcgacgGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGA‑‑G‑‑CACAAATTTTgagaga                                                                        <  2:3933966/150‑1 (MQ=255)
                                                                             cgcgacgGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGTCCCGCACCCCAACCCACGCCCCCCAAGGAGAGGGGCCCGAGCGA‑‑GCCCACAAAttttttga                                                                          >  1:478624/1‑146 (MQ=21)
                                                                              gcgacgGACGAAAAAGTATTCGTATGGTTGGTTTTCTCCCTCTCCCTGTGGGTGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACACTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGA‑‑G‑‑CACAAATTTTGAGAGat                                                                       <  2:3821151/150‑1 (MQ=255)
                                                                               cgacgGACGAAAAAGTATTGTTAAGGTGGTTTTACCCCCTCTCCCTGGGGGTTTGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGA‑‑G‑‑CACAAATTTTGAGAGata                                                                      <  1:850259/150‑1 (MQ=255)
                                                                                  cggactAAAATGTATTCGTACGGTTGGTTTTCTCCCGCGCCCGGGGGGTGCGGGCCGGGGGGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGA‑‑G‑‑CACAAATTTTGAGAGATAtct                                                                   <  2:1893807/150‑1 (MQ=255)
                                                                                  cggacgAAAATGTATTCGTACGGTTGGTTTTCTCCCTCTCCCTGTGGTATAGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGA‑‑G‑‑CACAAATTTTGAGAGATAtct                                                                   <  2:1893887/150‑1 (MQ=39)
                                                                                   ggacgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCCCCCACACCCAAAACCAAAACCCCAACGGAGAGGGGCCCCCTAGAGCG‑‑CACAAATTGTGAGAAatcat                                                                    >  1:2733135/1‑147 (MQ=16)
                                                                                     acgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCCCACCCAAAACCCCTCCCACCAAGGGGAGGGGACCGCTAGA‑‑C‑‑CCCAAATTTTTAGGGataaaatat                                                                >  2:3037766/1‑144 (MQ=21)
                                                                                         aaaaTTTTTCGTATTGTTGTTTTTCGCCCCCTCCCTGGGGTTGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGA‑‑G‑‑CACAAATTTTGAGAGATATCTTCTTCTg                                                            <  1:3896154/150‑1 (MQ=255)
                                                                                          aatttATTCGTATGGTTGGTTTTCTCCATCTCCCGGTGGTAGTGGGCCGGGGTGAGGTCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGA‑‑G‑‑CACAAATTTTGAGAGATATCTTCTTCTGt                                                           <  1:3645563/146‑1 (MQ=255)
                                                                                          aaaGTTTTCGTTAGGTTGGTTTTCGCCCTCTCCCTGTGGGTGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGA‑‑G‑‑CACAAATTTTGAGAGATATCTTCTTCTGt                                                           <  2:3586527/150‑1 (MQ=255)
                                                                                                 tcgtttggtgggttTTCCCCCTCTCCGGGGGGGTGTGGGCGGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGA‑‑G‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTaa                                                    <  1:1158076/144‑1 (MQ=255)
                                                                                                 tCGTATGGTTGGTTTCCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGA‑‑G‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTaa                                                    <  2:434909/150‑1 (MQ=255)
                                                                                                 tCGTAAGGTTGGTTTTCGCCCTCGCCCTGTGGGAGTGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGA‑‑G‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTaa                                                    <  2:434952/150‑1 (MQ=255)
                                                                                                  cttAAGTTTGTTTTTCTCCCTCTCCCTGTGGGATAGGGCCTGGGTGAGGGCATCAGCGCTCACGTTCACCCTCACCATAACCCTCTCCCTCAAGGGAGAGGGGACCTATCGA‑‑G‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAAc                                                   <  2:225581/148‑1 (MQ=38)
                                                                                                   gTAAGGTTGGTTTACTCCCTCCCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGA‑‑G‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACa                                                  <  2:1038041/150‑1 (MQ=255)
                                                                                                      tggttggttTTCTCCCTCTCCCGGTGGGAGAGGGCCGGGGTTTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGA‑‑G‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaa                                               <  1:4161521/149‑1 (MQ=255)
                                                                                                      aGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCACACCCCAAACCGCACCCGCAAGGGGGAGGGGGCCGATCGG‑‑G‑‑CCCCAAATTTGAGTGATTTATTGTTATGTGTTGTAgaagaa                                               >  2:2269029/1‑150 (MQ=12)
                                                                                                       ggttggttTTCTCCCTCTCCCTGTGGGAGAGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGA‑‑G‑‑CACAAATTTTGAGAGATATCTTCTTCGGTCTTGTAACagaag                                              <  2:1205734/150‑1 (MQ=255)
                                                                                                       ggttggttTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGA‑‑G‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaag                                              <  2:2538124/150‑1 (MQ=255)
                                                                                                       ggttggttTTCTCCCAATCCCTGTGGGTGTGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGA‑‑G‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaag                                              <  2:2346695/150‑1 (MQ=255)
                                                                                                         tgggttTTCTCCCTCTCCCTTTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCATCTACCTCAAGGGAGAGGGGACCGATCGA‑‑G‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAgaagaa                                            <  2:3477434/148‑1 (MQ=255)
                                                                                                           gtttttCTCCCTCTCCCTGGGGGAGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGA‑‑G‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACa                                          <  1:3278955/148‑1 (MQ=255)
                                                                                                            gttGTCCCCCCCCCCGTGGGGGTGTGGGCGGGGGTGAGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGA‑‑G‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAg                                         <  1:2851953/150‑1 (MQ=255)
                                                                                                             ttttCGCCCACTCCCTGTGGGAGCGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGA‑‑G‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGa                                        <  1:170851/150‑1 (MQ=255)
                                                                                                             tgttCCCCCTCTCCCGGGGGGTGTGGGCGGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGA‑‑G‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGa                                        <  1:272498/148‑1 (MQ=255)
                                                                                                               ttCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGAGGGAATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGA‑‑G‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGaaa                                      <  2:636273/150‑1 (MQ=255)
                                                                                                                   ccctctccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCGCACCCCAACCCACACCCACCAAGGGGAGGGGGCCGCAAGA‑‑G‑‑CACAAAAATTTTGAGATATATTATTGTGTGTTGTAAAAAAAGAAAAAAAAAACg                                  >  2:780006/1‑150 (MQ=12)
                                                                                                                         ccctGTGGGCGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGA‑‑G‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCttt                            <  1:2564564/150‑1 (MQ=255)
                                                                                                                         ccctGGGGGTGTGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGA‑‑G‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCttt                            <  1:3008855/150‑1 (MQ=255)
                                                                                                                           cggtgGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGA‑‑G‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCtttta                          <  2:486370/148‑2 (MQ=255)
                                                                                                                                ggAGTGGGCGGGGGGGCGGCCATCAGCGCGCCCGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGA‑‑G‑‑CACCAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCt                     <  2:1126784/150‑1 (MQ=255)
                                                                                                                                     gggCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGA‑‑G‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAata                <  2:1550079/150‑1 (MQ=255)
                                                                                                                                      ggCCGGGGTGCGGGCATCAGCGCGCACGGTCACCCTCACCCGAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGA‑‑G‑‑CACAAATTCTGAGAGATATCTTCTTCTGTATTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAatat               <  2:1650487/150‑1 (MQ=255)
                                                                                                                                      ggCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGA‑‑G‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAatat               <  2:47332/150‑1 (MQ=255)
                                                                                                                                        ccGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGA‑‑G‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTg             <  2:845634/150‑1 (MQ=255)
                                                                                                                                         cGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCGAAACCACACCCACAAGGGGGAGGGGCCCCAAAGA‑‑G‑‑CCCCAAATTTTTGAGATATATTATTATTTATTGTAACAGAAGAAAACAAAAACTGGCGATTGGGCCGGaaaaattc            >  2:3384160/1‑144 (MQ=255)
                                                                                                                                           gggTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGA‑‑G‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCgg          <  1:2478777/150‑1 (MQ=255)
                                                                                                                                                 gggCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGA‑‑G‑‑CACAAATTTTTAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaa    <  2:3681884/150‑1 (MQ=255)
                                                                                                                                                 gggCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGA‑‑G‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaa    <  1:2713590/150‑1 (MQ=255)
                                                                                                                                                  ggCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGA‑‑G‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaaa   <  2:1228842/150‑1 (MQ=255)
                                                                                                                                                   gCATCAGCGCGCACGTTCACCCTCTCCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGA‑‑G‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAg  <  2:2378801/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                       
CCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGA‑‑G‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAG  >  CP009273/160919‑161215

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: