Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 161,068 A→C 39.5% 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,0680AC39.5% 22.5 / 22.9 60intergenic (+47/‑149)hrpB/mrcBputative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base A (21/15);  new base C (5/19);  total (26/34)
Fisher's exact test for biased strand distribution p-value = 7.24e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAA  >  CP009273/160920‑161214
                                                                                                                                                     |                                                                                                                                                  
cGTGAGGTGCAAAAAGAGATGAAAGGGCTTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACTAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCAt                                                                                                                                                   <  2:2757225/150‑1 (MQ=255)
cGTGAGGTGAAAAAATAGATGAATGGGCTTTATCCAAAACTTGTCGGGCCGGACGACCCGGCAAATACGGCCCCGACGCGACGGACTAAAATGTTTTCGTATGGTTGGTTTTCTCCCTCGCCCTGTGGGTGCGGGCCGGGGTGAGGGCAt                                                                                                                                                   <  2:1653529/150‑1 (MQ=255)
  tgtggTGCAAAAAGAGTTGAATGGGGGTTAACCCAAACCTGTGGGGCGGTAGGACCCGGCAAATACTGCACCGACGCGCGGGACTAAAACTTTTTGGTATGGTTGGTTTTCTCCCTCGCCCTGTGGTAGTGGGCCGGGGTTTGGGCATCa                                                                                                                                                 <  2:3780006/147‑1 (MQ=255)
  tGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTAGTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCa                                                                                                                                                 >  1:4610820/1‑150 (MQ=255)
   gAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAg                                                                                                                                                >  2:1974777/1‑150 (MQ=255)
    aGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAACAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCCgc                                                                                                                                               >  2:800085/1‑150 (MQ=255)
    aGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAgc                                                                                                                                               >  2:4467304/1‑150 (MQ=255)
    aGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCgc                                                                                                                                               >  2:799961/1‑150 (MQ=255)
    aGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAgc                                                                                                                                               >  2:4467182/1‑150 (MQ=255)
      gTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAgcgc                                                                                                                                             >  1:3360215/1‑150 (MQ=255)
       tGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAgcgcg                                                                                                                                            >  1:3629772/1‑150 (MQ=255)
          aaaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCAc                                                                                                                                         >  2:2956680/1‑150 (MQ=255)
           aaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACg                                                                                                                                        >  2:4047044/1‑150 (MQ=255)
           aaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGAAAGGTTGGTTTTCTCCCTCTCCCAGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACg                                                                                                                                        >  1:308181/1‑150 (MQ=255)
            aaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAGCAGCGCGCACGt                                                                                                                                       >  2:1279634/1‑150 (MQ=255)
             aaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGtg                                                                                                                                      >  2:4195646/1‑149 (MQ=255)
             aaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGAGTGAGGGCATCCGCGCGCACGtg                                                                                                                                      >  2:618923/1‑149 (MQ=255)
              agagATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCCCgggc                                                                                                                                     >  2:4229400/1‑147 (MQ=255)
              agagATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACgggg                                                                                                                                     >  1:3288930/1‑147 (MQ=255)
                  tttactgggCGTTAACCAAAACTTGTGGGGCCGGACTACCCGTCAAATACTGCACCTCCGCTAGGGACTAAAATGTTTTGTTATGTTTGTTTTTCTCCCTCTCCCTGTGGTTGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTtcaccc                                                                                                                                 <  2:3125178/144‑1 (MQ=16)
                  aTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACTCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGTCCTCCGCGCGCACGtgggccc                                                                                                                                 >  2:1469626/1‑144 (MQ=38)
                      atgggttttATCCCAACCTGGTGGGGCGGGACGCCCGGCCAAATACTGCACCGACGCGCGGGACTAAAATGTTTTCTTATGTTTGTTTTTCCCCCTCGCCCGGTGGGTGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTtcaccctcac                                                                                                                             <  2:1916248/143‑1 (MQ=255)
                      aaGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTggacccacac                                                                                                                             >  2:2735070/1‑150 (MQ=255)
                         ggCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTggacccacacccc                                                                                                                          >  2:3126775/1‑149 (MQ=255)
                         ggCGTTAACCCACACATGTCGGGCCGGCCGCCCCGGCAAATACTGCCCCGACGCGACGGACGAAAACGTATGCGTACGGTTGGTTTTCTCCCCCTCCCGGTGGGTGTGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTcaccctcaccct                                                                                                                          <  2:4695262/150‑1 (MQ=255)
                           cGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTgcacccacaaccgca                                                                                                                        >  2:4722790/1‑147 (MQ=255)
                                   cAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCCCGGGCCCCCGCACCCCAAcccaccac                                                                                                                >  1:632554/1‑145 (MQ=37)
                                    aaaCATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGCCCCGCACCCCCAACCACTccc                                                                                                               >  2:4219366/1‑150 (MQ=37)
                                    aaaCATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGCCCCACCCCCAAACCCCCAccc                                                                                                               >  2:2955652/1‑150 (MQ=37)
                                     aaCATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTTGCCCCCCACCCCAAaccaaacccg                                                                                                              >  1:1976949/1‑143 (MQ=37)
                                     aaCATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAGCGGCGCGCACGTGGCCCCGCCCCCAAACCCACcccca                                                                                                              >  1:1976957/1‑149 (MQ=37)
                                       ctggTCGGGCGGGACGCCCGGCAAAATACTGCCCCGCCGCGACGGCCTAAAATTTTTTGTTATGTTTGGTTTTCCCCAGCCCCGTGTGGTAGTGGGCCGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCa                                                                                                            <  2:4684574/147‑1 (MQ=255)
                                       cATGTCTGGCCGCACGACCCGCCAAATACGGCCCCGACGCGACGGACTAAAAAGTATTCTTATGGTTGTTTTTCCCCCTCCCCCTGTGGGTGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCa                                                                                                            <  2:1978753/150‑1 (MQ=255)
                                              ggccggccGCCCCGGAAAATACGCCCCCGCCGCGCCGGCCTAAAATGTTTTCGTTTGGGTGGTTTGCGCCCTCGCCGTGGGGGTGTGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGgaga                                                                                                     <  2:2422347/150‑1 (MQ=255)
                                                        cccGGCAAATCGTCCACGGCCGCGCCGGACTAAAATTTATTCTTATGGTTGTTTTTCCCCCTCGCCCTGTGGGAGTGGGCGGGGGTTCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGAt                                                                                           <  2:707486/150‑1 (MQ=255)
                                                                              gcgccgGACGACAACGTAGTCGTATGGGTGGTTTTCTCCCTCGCCCTGTGGGTGGGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGat                                                                      <  2:3141675/150‑1 (MQ=38)
                                                                              gcgacgGACGACAAATTATTCGTATGGTGGGTTTTCTCCCTCTCCCTGTGGGTGTGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGat                                                                      <  2:3141159/150‑1 (MQ=38)
                                                                                   ggcccaaaatGTATTGTTATGTTTGGTTTTCCCCCACCCCCTGGGGTTCTGGGCGGGGGTGTGGGCATCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtctt                                                                 <  2:294427/145‑1 (MQ=255)
                                                                                                   ttatggttggttTACACCCCCCCCGTGGGGGTGTGGCCCGGGGGTCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACa                                                 <  2:1859139/146‑1 (MQ=255)
                                                                                                   gTATGGTGGGTTTTCGCCCTCGCCCTGGGGGAGAGGGCCGGGGGGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACa                                                 <  2:4610820/150‑1 (MQ=255)
                                                                                                           gtttttCGCCCGCTCCGTGGGGGTGTGGGCGGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACa                                         <  2:1540790/148‑1 (MQ=255)
                                                                                                           ggttTTCTCCCTCTCCGTGTGGGATTGGGCGGGGGTGAGGGCCTCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACa                                         <  2:3172484/150‑1 (MQ=255)
                                                                                                           ggttTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGCCCCGCAACCCAACCCACACCCACCCAGGAGGGGGGGCCGCTCGAGCCAAAATTTTGTGAGATATCTTATTATGTCGTGTTAAacaacaaca                                         >  2:3554635/1‑150 (MQ=12)
                                                                                                               gagcccctccccctGGGGGTTGGGGCGGGGTTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGaaa                                     <  1:3488751/146‑1 (MQ=255)
                                                                                                                             gggggAGTGGGCGGGGGGGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCGCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAACTCGGGCTTTTgcg                       <  2:1256951/148‑1 (MQ=255)
                                                                                                                               gggtgtgggCGGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGcc                     <  2:2689099/144‑1 (MQ=255)
                                                                                                                               gggtgtgggCGGGGGTGAGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGcc                     <  2:2757500/144‑1 (MQ=255)
                                                                                                                                 gtgtgggCCGGGGTGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAACATCGGGCTTTTGCGCCTg                   <  2:4292552/146‑1 (MQ=255)
                                                                                                                                 gagtGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCGCACCCGAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTg                   <  2:4291729/150‑1 (MQ=255)
                                                                                                                                  tgggggCGGGGGTGGGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGa                  <  1:799961/147‑1 (MQ=255)
                                                                                                                                  cgtggcccGGGGGGCGGCCCTCAGCGCGCACGTGCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTgc                  <  1:800085/144‑2 (MQ=255)
                                                                                                                                   gggggggggggTGGGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGaa                 <  1:396993/143‑1 (MQ=255)
                                                                                                                                     gggCGGGGGTGGGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAata               <  2:4207710/150‑1 (MQ=255)
                                                                                                                                      ggCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAatat              <  2:1586896/150‑1 (MQ=255)
                                                                                                                                      gcccGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAATTTGAGCGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAatat              <  2:3076107/148‑1 (MQ=255)
                                                                                                                                       gcgggggggAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATAtt             <  2:238293/147‑1 (MQ=255)
                                                                                                                                       gCCGGGGTGAGGGCCTCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATAtt             <  2:238286/150‑1 (MQ=255)
                                                                                                                                             gTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGgag       <  2:3068617/150‑1 (MQ=255)
                                                                                                                                                  ggCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaaa  <  2:3939485/150‑1 (MQ=255)
                                                                                                                                                  gccATCCGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaaa  <  1:1376594/148‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
CGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAA  >  CP009273/160920‑161214

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: