Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 164,581 A→C 21.4% 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
*NC_000913164,5810AC21.4% 90.1 / 5.5 63intergenic (+47/‑149)hrpB/mrcBputative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base A (22/26);  new base C (2/13);  total (24/39)
Fisher's exact test for biased strand distribution p-value = 3.27e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAG  >  NC_000913/164432‑164728
                                                                                                                                                     |                                                                                                                                                   
ccGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCa                                                                                                                                                     >  1:1357729/1‑150 (MQ=255)
   tGAGGTGCAAAACGAGATGAACGGGCGTTATCCCAACCATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGCCGGCCGAAAACGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGGGGGAGAGGGCCGGGGTGAGGGCATCa                                                                                                                                                  <  2:1848540/150‑1 (MQ=39)
   tGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCa                                                                                                                                                  <  2:1567570/150‑1 (MQ=255)
        tGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAgcgcg                                                                                                                                             >  2:2609018/1‑150 (MQ=255)
          cAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCa                                                                                                                                           >  1:1611577/1‑150 (MQ=255)
                   aTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGttgcccc                                                                                                                                  >  1:1206168/1‑145 (MQ=255)
                   aTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGtggcccc                                                                                                                                  >  2:416863/1‑144 (MQ=255)
                   aTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGtgggccc                                                                                                                                  >  1:496226/1‑144 (MQ=255)
                   aTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAGCAGCGCGCACGTtcaccc                                                                                                                                  >  1:1206008/1‑150 (MQ=255)
                     gAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGttgaacccc                                                                                                                                >  2:1431542/1‑143 (MQ=255)
                     gAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTgcacccgc                                                                                                                                >  2:1848521/1‑148 (MQ=255)
                          ggCGTTATCCCAAACATGTCTGGCCGGACGACGCGGCAAATCCTGCCCCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCGGGGGTGTGGGCATCAGCGCGCACGTTcaccctcaccct                                                                                                                           <  2:1712867/150‑1 (MQ=38)
                           gCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTgctcccccccccca                                                                                                                          >  2:510034/1‑148 (MQ=255)
                           gCGGTATCCCAAACATGTCTGGCCGGACGACCCGGCACATACTGCACCGACGCGACGGACGAAAAAGTATTCGTCAGGTTGGTGTTCTCCCTCTCCCTGTGGGAGCGGGCCGGGGTGCGGGCATCAGCGCGCACGTgcacccacaccccc                                                                                                                          >  2:1504435/1‑148 (MQ=38)
                             gTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTGCCCCCCCACCCAAAc                                                                                                                        >  2:1512792/1‑150 (MQ=38)
                                 tCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGGGCCCCCCCCCCCCAAcccac                                                                                                                    >  2:1531007/1‑148 (MQ=38)
                                 tCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCGCACCCaaaaacac                                                                                                                    >  2:1148740/1‑145 (MQ=255)
                                 tCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCCCACCCCAAaccac                                                                                                                    >  2:1148733/1‑148 (MQ=38)
                                   ccAAACTTGTCTGGCGGGCCGCCCGGGATAATAGTGCCCGGACGGGAGGGACTAAAATGTTTTCTTAAGGTTGGTTTTCTCCCTCGCCGTGTGGGTGTGGGCGGGGGTGCGGGCATCAGCGCACACGTTCACCCTCACCCTAACCCtctc                                                                                                                  <  2:2303885/150‑1 (MQ=255)
                                   ccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCtctc                                                                                                                  <  2:2164553/150‑1 (MQ=255)
                                     aaaCATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTGGACCCACCCCCCAACCCACAccc                                                                                                                >  1:2255238/1‑150 (MQ=37)
                                       aCATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTc                                                                                                              <  2:2173310/150‑1 (MQ=255)
                                         aTGTCTGGCCGGACGACCCGGAAAATACTGCACCGACGCGACGTACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCCCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCaa                                                                                                            <  2:474287/150‑1 (MQ=255)
                                             cTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTCAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGCCCCTCCCCCAAACCCGCACCCGCAAgggg                                                                                                        >  2:969642/1‑149 (MQ=37)
                                               ggCGGGACGCCCCGGAAAATAGTGCCCCGCCGCGCCGGACTAAAATGTTTTCGTCCGGTTGGTTTTCCCCCGCTCCCTGGGGGAGTGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGgaga                                                                                                      <  2:2518580/150‑1 (MQ=255)
                                               ggCCGGACGCCCCGGCAACTCCTGCCCCGCCGCGCCGGACGAAAAAGTTTTCGTAGGGTTGTTTTTCCCCCCCGCCCTGTGGGTGTGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGgaga                                                                                                      <  2:2576064/150‑1 (MQ=255)
                                               ggCCGGACGACCCGCCACCAACTCCCCCGACGCCACGGACGAAAAAGTATGCGCAAGGTTGTTATACTCCCTCCCCCTGTGGTAGAGGGCCGGGGCGCGGGCATCCGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAATGgaga                                                                                                      <  2:1771432/150‑1 (MQ=255)
                                                  cGGACGACCCGTAAAATACGGCACCGACGCGACGGACTAAAATGTTTTCGTCTGGTTGGTTTTCTCCCCCGCCGTGTGGGTGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAggg                                                                                                   <  1:2831984/150‑1 (MQ=255)
                                                    gacgacCCGCCAACCCCCGCCCCCCCGCGCCGGCCGAACCAGGCTTCGAAAGGGTGTATTTCGCCCCCTCCCTGGGGGAGTGGCCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCCCCCTAACCCTCTCCCTCAAGGGAGAGGGGa                                                                                                 <  2:791025/150‑1 (MQ=255)
                                                          ccGGAAAATACTGCACCGACGCGCCGGACGAAAAAGTTTTCGTACGTTTGGTGTCCGCCCACTCCCTGTGGGTGTGGGCCGGGGGGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATc                                                                                           <  2:2259917/150‑1 (MQ=255)
                                                                 atacgccACCGCCGCGAGGGACTAAAAAGTTGTCGTATGTTTGTTTTCCACCCCCTCCGTGTGGGAGTGGGCGGGGGTGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGcaca                                                                                    <  1:163709/144‑1 (MQ=255)
                                                                       ccccGACGCGCCGGACGAAAAATTTTTCGTATGGTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtttt                                                                              <  2:2714678/148‑1 (MQ=38)
                                                                       cACCGACGCGCGGGACCAAAACTTATTCTTATGTTTGTTTTTCTCCCTCGCCGTGGGGTAGCGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtttt                                                                              <  1:1993780/150‑1 (MQ=255)
                                                                         ccGACGCGACGGCCTAAACATTTTGCTTACGGTTGGTTTTCTCCCTCTCCTTGGGGGTGTGGGCGGGGTTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGATGGGACCGATCGAGCACAAATTTTga                                                                            <  1:2494669/150‑1 (MQ=255)
                                                                                   ggacgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCCCCCCCAAAACCACACCCACAAGGGGGAGGGGGCCCCGCGCGCACAAATTTTGTGGGATAtatt                                                                  >  1:2138982/1‑150 (MQ=21)
                                                                                    gacgacAACGTTTTCGTATGGTTGTTTTTCTCCCTCTCCCTGTGGGTGCGGGCCGGGGTGAGGGCCTCCGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttc                                                                 <  2:1631330/150‑1 (MQ=38)
                                                                                      cgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCACACCCAAAACCACACCCACAACGGGGTGGGGGCCGATAGAGCAAAAAATTTTTGAGATAtattaat                                                               >  1:969340/1‑147 (MQ=18)
                                                                                             gTATTGTTATGGTTTTTTTCCCCCCTCCCCCTGTGGTAGAGGGCGGGGGTGCGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                                        <  1:741114/150‑1 (MQ=255)
                                                                                             gTATTCGTTAGGTTGGTTTTCTCCCTCTCCCTGGGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                                        <  1:2182702/150‑1 (MQ=255)
                                                                                             gTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                                        <  1:2182659/150‑1 (MQ=255)
                                                                                                  ggtatgtttggttTTCGCCCTCTCCCTGTGGGAGAGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAAc                                                   <  1:1356986/143‑1 (MQ=39)
                                                                                                   gTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTTCCCCCACCCCCGAAACCACGCCCGCAAGGGAGGGGGGCCCCATCGCGCCCCCAATTTTTGAGCTATCATATTCTGTGTTGTAAcc                                                  >  2:1896068/1‑149 (MQ=12)
                                                                                                      cggttggtgGCCTCCGGCTCCCGGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaa                                               <  2:2567176/149‑1 (MQ=255)
                                                                                                        gttgtttTTCTCCCTCTCCGTGTGGGTGTGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaaga                                             <  2:2774913/150‑1 (MQ=255)
                                                                                                        gttgggtTTCGCCCCCGCCGTGGGGGTGCGGGCGGGGGTGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACACATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaaga                                             <  2:276453/150‑1 (MQ=255)
                                                                                                           ggttTTCTCCCTCTCCCTGTGGGAGTGGGCGGGGTTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACa                                          <  1:1448773/150‑1 (MQ=255)
                                                                                                           ggttTACGCCCTCGCCCTGTGGGAGAGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACACATTTTGAGAGATATCTTCTTCTGTCGTGTAACAGAAGAACa                                          <  2:681990/150‑1 (MQ=255)
                                                                                                              tttCTCCCTCTCCGTGTGGGTGTGGGCGGGGGTGTGGGCATCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAgaagaagatat                                       <  1:1512792/150‑6 (MQ=255)
                                                                                                                  tcactctccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCACAACCTAAAACGCTCCCGAAAAGGTGATGGGTGCCAGCGAGCACAAATTTTTAGAGATATCATATTATGTCGTTTTAAAGGAGAAAACaaaaaa                                   >  2:1182614/4‑148 (MQ=12)
                                                                                                                   ccctctccctGTGGGAGAGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGCGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCg                                  <  2:2896765/150‑1 (MQ=255)
                                                                                                                   ccctcgcccgGGGGGAGTGGGCGGGGGTGAGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCg                                  <  1:2805529/150‑1 (MQ=255)
                                                                                                                               gggAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGcc                      <  2:1869463/150‑1 (MQ=255)
                                                                                                                                ggAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCGCCCGCACCCAAACCCACAACCACCACGGGGAGGGGGCCCAAAGAGCACAAAAATTTTGAGATATATTATAATTTGTTTTTACAGAATAAAACAAAAAAGGGCTGTTTCGccg                     >  2:2044409/1‑149 (MQ=255)
                                                                                                                                       gCCGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATAtt              <  2:782092/150‑1 (MQ=255)
                                                                                                                                          ggggTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg           <  2:2528429/150‑1 (MQ=255)
                                                                                                                                            ggTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGga         <  1:1593413/150‑1 (MQ=255)
                                                                                                                                            ggTGAGGGCATCAGCGCGCACGTTAACCCACACCCAAAACCCAACCCACAAGGGTGAGGGGACCCCTAGACCACCAAATTTTAGATATATAATATTATGTATTTTAACACAAAAACACCAAATAGGGCGTTTTCGCCGGAATAtttatgg         >  2:1762374/1‑145 (MQ=255)
                                                                                                                                            ggTGAGGGCATCAGCGCGCACGGTGACCCACACCCAAAACCACACCCACAAAGGCCTGGGGGCCGATCGAGCCCAAAATTTGTGAGAAATATAATAAGGTCTTGTTACCAAACAACAGGAAAAACGGCGTTTTCGCCGGAAAATTggggg         >  2:1762365/1‑149 (MQ=255)
                                                                                                                                                 gggCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaa    <  2:2584643/150‑1 (MQ=255)
                                                                                                                                                  ggCATCAGCGCGCACGTTCACCCTCACCCTAACCCTTTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaaa   >  2:42102/1‑150 (MQ=255)
                                                                                                                                                  ggCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCGTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaaa   <  2:1765087/150‑1 (MQ=255)
                                                                                                                                                  ggCATCAGCGCGCACGTTCACCCTCACCCGAACCCTCTCCCTCAAGGGAGAGGGGACCGATAGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaaa   <  2:933210/150‑1 (MQ=255)
                                                                                                                                                  ggCATCAGCGCGCACGTACACCCACCCCCTACCCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTACCAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaaa   <  1:1682256/150‑1 (MQ=255)
                                                                                                                                                   gCATCAGCGCGCCCGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAg  <  1:1401189/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAG  >  NC_000913/164432‑164728

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: