Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 164,581 A→C 20.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,5810AC20.4% 302.0 / 43.5 211intergenic (+47/‑149)hrpB/mrcBputative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base A (69/96);  new base C (4/42);  total (73/138)
Fisher's exact test for biased strand distribution p-value = 1.59e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 213 total aligned reads displayed.

CCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTT‑GGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCA  >  NC_000913/164432‑164730
                                                                                                                                                      |                                                                                                                                                      
ccGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTT‑GGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCa                                                                                                                                                        >  1:5094194/1‑150 (MQ=255)
 cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTT‑GGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAt                                                                                                                                                       >  1:1594318/1‑150 (MQ=255)
 cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTT‑GGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAt                                                                                                                                                       >  1:705044/1‑150 (MQ=255)
 cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTT‑GGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAt                                                                                                                                                       >  2:8362832/1‑150 (MQ=255)
 cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTT‑GGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAt                                                                                                                                                       >  1:7046509/1‑150 (MQ=255)
  gTGAGGTGAAAAATGAGATGAAAGGGCGTTATCCCAACCATGTCTGGCCGGACGAGCCGGAAAATACTGCACCGACGCGACGGACTAAAAAGTATTCTTAAGGTT‑GGTTTTGTCGTTCTCCGTTTGGGAGTGGGCCGGGGTGAGGGCATc                                                                                                                                                      <  2:8632139/150‑1 (MQ=255)
   tGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTT‑GGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCa                                                                                                                                                     >  2:2567434/1‑150 (MQ=255)
   tGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTT‑GGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCa                                                                                                                                                     >  2:9183101/1‑150 (MQ=255)
   tGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTT‑GGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCa                                                                                                                                                     >  1:525791/1‑150 (MQ=255)
   tGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTT‑GGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCa                                                                                                                                                     >  1:4377614/1‑150 (MQ=255)
      ggTGCAAAAAGAGTTGAAGGGGGTTTAGCCAAAACTTGTCTGGCCGGCCGACCCGGAAATTCTTGCACCGACGCGACGGACGAAAATGTTTTGGTTAGTTT‑GTTTTTCTCCGACTCCCTGTGGGTGAGGGCCGGGGTTTGGGCATCAgcg                                                                                                                                                  <  2:7687270/150‑1 (MQ=255)
       gTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTT‑GGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAgcgc                                                                                                                                                 >  2:498844/1‑150 (MQ=255)
        tgaaaaaaGAGATGAACGGGCGTTATCCCAAACCTGTCTGGCCGGACGACCGGGCAAATACTGCACCGACGCGACGGACGAAAATGTATTCGTAAGGTT‑GGTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAgcgcg                                                                                                                                                <  2:1361906/147‑1 (MQ=38)
               agagATGAAAGGGCGTTATCCCAAACTTGTCTGGCCGGACGACCGGGCAAATACTGCACCGACGCGACGGACGAAAATGTTTTCTTAAGGTT‑GGTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTtc                                                                                                                                         <  2:4538650/150‑1 (MQ=38)
               agagATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTGAGGTT‑GGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGCGGGCCTCAGCGCGCACGtgg                                                                                                                                         >  1:6419694/1‑148 (MQ=255)
               agagATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTT‑GGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACgggc                                                                                                                                         >  1:6927201/1‑147 (MQ=255)
               acagaTGAACGGGCGTTATCCCAAACATGTCTGGCCGGACGCCCCGGCAACCACTGCACCGCCGCGACGGACGAAAAATTATTCGTCTGGTT‑GGTTTTCTCCCTCGCCCGGTGGGAGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTtc                                                                                                                                         <  2:7355752/148‑1 (MQ=255)
                  gaTGAACGGGGTTTATCCCAAACATGTGTGGCCGGACCCCCCGGAAAATACTGCCCCGACGCGACGGACGAAAAATTATTCTTAAGGTG‑GTTTTTCTCCCTCTCCCTGTGGGTGTGGGCGGGGGTTAGGGCATCAGCGCGCACGTtcacc                                                                                                                                      <  2:1748666/150‑1 (MQ=25)
                  gaTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTT‑GGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGgcacc                                                                                                                                      >  1:7950076/1‑150 (MQ=255)
                   aTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTT‑GGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTgccccc                                                                                                                                     >  2:4515347/1‑150 (MQ=255)
                      aaaGGGCGTTATCCCAAACAGGTCTGGCCGGACGACCCGGCAAATACGGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTT‑GGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccctca                                                                                                                                  <  2:9342182/150‑1 (MQ=255)
                        tgggCTTTATCCCAAACTTGTCTGGCCGGACGACCCGGAAAATACTGCACCGACGCGCGGTACTAAAATTTTTTCTTATGGTT‑GTTTTTCTCCCTCTCCGTGTGGTTGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTtcaccctcacc                                                                                                                                <  2:2079031/149‑1 (MQ=255)
                        aGGGCGTTATCCCACCCATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTT‑GGTTTTCACCCTCTCCCTGTGGGTGAGGGCGGGGGTGAGGGCATCAGCGCGCACGTtcaccctcacc                                                                                                                                <  2:6224425/150‑1 (MQ=39)
                         gggCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTT‑GGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTggccccacaccc                                                                                                                               >  1:478433/1‑150 (MQ=255)
                             gTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGGATTCGTAAGGTT‑GGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGCGGGGAGGGCATCAGCGCGCACGTggaccccccccccaca                                                                                                                           >  2:7861828/1‑146 (MQ=38)
                                aacccAAACTTGTGGGGCCGGACGACCGGGAAAATACGGCCCCGACGCCACGGACTAAAATGTTTTCGTATGGTT‑GGTTTTCCCCCTCGCCGTGGGGGAGCGGGCCGGGTTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCct                                                                                                                        <  2:8896863/148‑1 (MQ=255)
                                  cccAAACATGTCTGGCCGTACGACCCGGCAAATACTGCACCGACGCGACGGACTAAAAATTATTCGTACGGTT‑GGTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGTTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCctc‑t                                                                                                                     <  2:4315230/150‑1 (MQ=38)
                                  cccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTT‑GGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCGCACCCAAAaccaca                                                                                                                      >  1:3524129/1‑147 (MQ=38)
                                   ccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTT‑GGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCGAAcccacac                                                                                                                     >  1:4984999/1‑146 (MQ=255)
                                   ccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTT‑GGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCACACCCCAAcccacac                                                                                                                     >  1:3357763/1‑146 (MQ=255)
                                    cAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTT‑GGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGGCACCCCCCCCCCCAcccacccc                                                                                                                    >  1:8868720/1‑145 (MQ=37)
                                    cAAACATGTCTGCCCGCACGACCCGGCAAAAACTGCACCGACGCGACGGACGAAAACGTATCCCTACGTTT‑GGTTTTCTCCCCCCCCCTGTGGTTGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑Tcc                                                                                                                   <  2:3376841/150‑1 (MQ=255)
                                         ttttCTGGCCGCACGCCCGGCCAAATACTGCCCCGACGCGCGGGACTAAAATTTTTTCTTATGTTT‑GGTTTTCTCCCTCCCCGTGTGGTCGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCaa                                                                                                              <  1:8320468/147‑1 (MQ=255)
                                         ttgtggggCCGCACGACCCGCAAAATAGTGCCCCGCGGGGCGGGACTAAAATTTTTTCTTTTGTTT‑GTTTTTCCCCCTCTCCGTGTGGGTGTGGGGGGGGTTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCaa                                                                                                              <  2:2529311/144‑1 (MQ=255)
                                          tGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTT‑GGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCGCACCCAAACCCGC‑TCCCTCCAg                                                                                                             >  1:7705937/1‑150 (MQ=38)
                                           gTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAATTATTCGTATGGTT‑GGTTTTCTCCCTCTCCCTGTGGGTGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAgg                                                                                                            <  2:5873264/150‑1 (MQ=255)
                                               ggccggccGCCCCGGCAAATACGGCACCGACGCGCCGGACGAAAAAGTATTCGTACGTTT‑GGTTTTCCCCCACTCCCTGTGGTAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTTGAGGgaga                                                                                                        <  2:5745192/150‑1 (MQ=37)
                                               ggCCGCACTACCGGCAAAATACGGCACCGACGCGAGGTACTAAAATTTTTTCTTACGGTT‑GTTTTACTCCCTCCCCCTTTGGGTTTGGGCCGGGGTGTGGGCCTCAGCGCGCACGTTCACCCCCACCCTAACCCTC‑TCCCTCAAGGgaga                                                                                                        <  1:3122354/150‑1 (MQ=255)
                                                gCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTAGTCGTAAGGTT‑GGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCGCCCCCGAAACCAC‑TCCCTCAAggggggg                                                                                                       >  1:3856213/1‑146 (MQ=37)
                                                  gggACGACCCGGAAAATACTTCACCTACGCGACGTACGAAAAATTATTGTTATGTTT‑GTTTTTAACCCTCTCCCTGTGGGATCGGGCGGGGGGCAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAggg                                                                                                     <  2:2951417/149‑1 (MQ=255)
                                                  cgtacgacCCGGCAAATACTGCACCGAGGCGACGTACTAAAAATTATTCGTATGGTT‑GGTTTTCTCCCTCTCCCTGGGGGTGTGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAggg                                                                                                     <  1:8475060/147‑1 (MQ=255)
                                                   ggACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTT‑GGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGGCCCCCACCCAAAACCAA‑ACCCACAAGGGGGCgggg                                                                                                    >  1:583176/1‑150 (MQ=25)
                                                   ggACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTT‑GGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGGGGCCCACCCCCAAAACCAC‑AACCACAAGGGAGAgggg                                                                                                    >  1:541593/1‑150 (MQ=25)
                                                     ccgacCCGGCAACAACGGCCCCGACGCGCCGGCCGAACAAGTTTTCGTACGGTT‑GGTTTTCCCCCTCTCCCTGGGGGAGCGGGCCGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAgggtta                                                                                                  <  2:8325345/149‑4 (MQ=255)
                                                     ccgacCCGGCAAACACTGCACCGACGCGACGGACGAAAATGTATTCGTACGGTT‑GTTTTTCTCCCTCCCCGTGTGGGAGCGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTTGAGGGAGAgggtta                                                                                                  <  2:4468551/149‑4 (MQ=255)
                                                     accacCCGGAAAACAGTGCACCGACGCGACGGACTAAAAAGTTTTCGTTTGGTT‑GTTTTTCTCCCTCTCCCTGTGGTAGCGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGAc                                                                                                  <  2:6614360/147‑1 (MQ=255)
                                                       cacCCGCAAAATACCCCCCCGCCGCGCGGGACTAAAATGTATTCGTTCGGTG‑GGTTTTCCCCGTCTCCCTGTGGTTGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTTGAGGGAGAgggttagg                                                                                                <  1:5214362/149‑6 (MQ=255)
                                                           cGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTT‑GGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGCCCCCCACCCAAAACCAC‑AACCACAAGGGCGAGGGGGCCGATCg                                                                                            >  2:1545603/1‑150 (MQ=25)
                                                           cGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTT‑GGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGCCCCACACCCAAAACCGC‑ACCCACCAGGGGGGGGGGGCCgcggg                                                                                            >  1:7636546/1‑146 (MQ=25)
                                                                    cTGCCCCGACGCGCCGGACGAAAAAGTATTCGTACGGTT‑GGTTGTCTCCCACTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAt                                                                                   <  2:7718288/150‑1 (MQ=39)
                                                                    cTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTT‑GGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGGGGCCCCGCACCCGAACCCAC‑ACCCACCAGGGCGGGGGGGCCCAGCGAGcccccaa                                                                                   >  1:4533980/1‑144 (MQ=21)
                                                                    cTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTT‑GGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCCCGGTCACCCTCACCCGAAACCAC‑TCCCACAAGGGGGAGGGGACCGATCGAGCACACAt                                                                                   >  2:6896208/1‑150 (MQ=37)
                                                                       ccccGACGCGACGGACTAAAATGTTTTCGTATGGTT‑GTTTTTCTCCCTCTCCGGGTGGTTGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtttt                                                                                <  1:3937985/148‑1 (MQ=255)
                                                                        cccgccgCGCCGTCCGAAAAATTTTTTTTATGGTT‑GGTTTTCGCCCTCTGCGGGTGGGTGTGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTg                                                                               <  1:9327842/149‑1 (MQ=255)
                                                                         ccGACGCGCCGGACGACAACTTTTTGGTATGTTT‑GTTTTTCTCCCTCTCCGGGTGGTCGTGGGCGGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                              <  2:1864697/150‑1 (MQ=255)
                                                                         ccGACGCGACGTACTAAAATTTATTCTTATGTTT‑GGTTTTCTCCCTCTCCCTGTGGGTGTGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                              <  2:5512765/150‑1 (MQ=37)
                                                                          cgacgcgacgGACGAAAAAGTTTTCGTATGGTT‑GTTTTTCTCCCTCTCCGTGTGGGTGAGGGCCGGGTTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgag                                                                             <  2:6527032/150‑1 (MQ=255)
                                                                            ccgcgacgGACGAAAATGTTTTCTTTTGGTT‑GTTTTTCTCCCTCTCCCTGTGGTAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgagag                                                                           <  2:5174296/149‑1 (MQ=38)
                                                                                 acgcacgAAAAAGTATTCGTAAGGTT‑GGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACCTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtc                                                                      <  1:9070010/150‑1 (MQ=255)
                                                                                    gacgAAAAAGTATTCGTAAGGTT‑GGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGGCCCACACCCCAAACCAC‑ACACACAAGGGGGGGGGGGCCGCTCGAGCCCCAAAATTTAGAGataaaaaa                                                                   >  1:4680921/1‑145 (MQ=16)
                                                                                      cgAAAAAGTTTGCGTATGGTT‑GGTTTTCTCCCTCTCCGGGTGGGAGCGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttctt                                                                 <  1:5983761/150‑1 (MQ=255)
                                                                                      cgAAAAAGTATTCGTAAGGTT‑GGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttctt                                                                 <  2:9187416/150‑1 (MQ=255)
                                                                                        aaaaaGTATTCGTAAGGTT‑GGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGACCCACACCCAAAACCAC‑ACCCACAACGGGGAGGGGACCCATCCAGCACAAAATTTTATAGATAtcttaataa                                                               >  2:4679162/1‑145 (MQ=21)
                                                                                          aatgttttCGTATGTTT‑GGTTTTCTCCCTCTCCCTGTGGGTGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGt                                                             <  1:9183158/144‑1 (MQ=38)
                                                                                           acgTATGCGTAGGGTG‑GGTTTTCGCCCGCCCCCTGTGGGAGCGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTc                                                            <  2:7507718/148‑1 (MQ=255)
                                                                                             ttATTCGTATGGTT‑GGTTTTCTCCCTCTCCCTGTGGTTGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                                          <  2:3697747/149‑1 (MQ=255)
                                                                                             gTATTCGTATGGTT‑GGTTTTCTCCCTCTCCGTGTGGGAGAGGGCCGGGTTGAGGGCCTCCGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                                          <  1:533878/150‑1 (MQ=255)
                                                                                             gTATTCGTATGGTG‑GGTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                                          <  1:6661782/150‑1 (MQ=255)
                                                                                                  tcggaGGTTGGGTTTGCCCCCCCCCCCGGGGGGTGGGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTaa                                                      <  1:8173030/146‑1 (MQ=255)
                                                                                                     aaGGTT‑GGTTGTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACaga                                                  <  2:5317733/150‑1 (MQ=255)
                                                                                                      tgttt‑gtttTACCCCCTGGCCGTGGGGGAGGGGGGGGGGGTGGGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaa                                                 <  1:2567434/147‑1 (MQ=255)
                                                                                                       gttt‑ggttTTCTCCCTCTCCGTGTGGGTGAGGGCGGGGTTGAGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaag                                                <  1:4803143/148‑1 (MQ=255)
                                                                                                        gtt‑ggttTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGTCACCCTCACCCAAAACCAC‑TACCACAACGGAGAGGGGACCCATCAAGCCCAAATTTTTAGAGATTAATTATTCAGTGTTGTCACaaaaga                                               >  2:8252827/1‑150 (MQ=18)
                                                                                                        gtt‑ggttTTCTCCCTCGCCCTGTGGGAGAGGGCCGGGGGGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaaga                                               <  2:1129164/150‑1 (MQ=255)
                                                                                                               ttccgccctctccgtGTGGGAGTGGGCGGGGGTGAGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGaa                                         <  2:7302301/145‑1 (MQ=255)
                                                                                                                 tctccctgtccgTGTGGGATAGGGCGGGGTTGGGGGCCTCAGCGCGCACGTTCACCCTCTCCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGaaaa                                       <  2:3582767/150‑1 (MQ=255)
                                                                                                                 tctccctctcccTGTGGGTGAGGGCGGGGTTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGaaaa                                       <  1:1545603/150‑1 (MQ=255)
                                                                                                                 tcgccctctcccgGGGGGAGAGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGaaaa                                       <  1:977176/147‑1 (MQ=255)
                                                                                                                    ccctccccgtGTGGGCGAGGGCGGGGGTGGGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCg                                    <  1:3303700/150‑1 (MQ=255)
                                                                                                                     cgtctccctGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCgg                                   <  2:8498723/148‑1 (MQ=255)
                                                                                                                           cctGGGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCtttt                             <  2:4534737/150‑1 (MQ=255)
                                                                                                                                 ggtgggggCGGGGGGGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCt                       <  1:3219750/145‑1 (MQ=255)
                                                                                                                                 ggAGTGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCt                       <  1:1896138/150‑1 (MQ=255)
                                                                                                                                 ggAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCt                       <  1:1069103/150‑1 (MQ=255)
                                                                                                                                 ggAGAGGGCCGGGGTGAGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCt                       <  2:901485/150‑1 (MQ=255)
                                                                                                                                  gagaGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGAACCCCAACCAAAACCACAACCCCCAAGGG‑GAGGGGGCCCAAAGAGCACAAAAATTTAGAGATATAATAATATGTCGTTTAAAAAAAGAACAAAAAAAAGGGCgttttggccgg                      >  2:1401003/1‑143 (MQ=255)
                                                                                                                                       ggCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCGCACCCAAACCCAC‑ACCCACAACGGAGAGGGGCCCGATAGAGCACAAAATTTTAGGGAAAAATTATTATGTCTTGTTAAAGAAGAACAGAAAAAAGGGCTTTTTCGCCGGAatat                 >  2:6098366/1‑150 (MQ=255)
                                                                                                                                          cGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGc              <  1:29507/150‑1 (MQ=255)
                                                                                                                                           ggggggggggCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg             <  2:3292240/150‑1 (MQ=255)
                                                                                                                                           ggggTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg             <  1:8362832/150‑1 (MQ=255)
                                                                                                                                           ggggTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg             <  2:694673/150‑1 (MQ=255)
                                                                                                                                           ggggTGAGGGCATCAGCGCGCCCGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg             <  2:1041365/150‑1 (MQ=255)
                                                                                                                                           ggggTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg             <  2:5835836/150‑1 (MQ=255)
                                                                                                                                             ggTGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGga           <  2:4028380/150‑1 (MQ=255)
                                                                                                                                                 aGGGCATCAGCGCGCACGTTGACCCACAACCGAAACCAC‑ACCCACAAAGGGGAGGGGACCCAACGAGCACAAAATTTTAGAGATAAAATAATAAGTCTTTTTACAAAAAAAAAGAAAAAACGGCGTTTTCGCCGGAATAATTCGGGGaaa       >  2:3802493/1‑150 (MQ=255)
                                                                                                                                                  gggCCTCAGCGCGCCCGTTCACCCTCACCCTAACCCTC‑TCCCTCAATGGAGAGGGGACCGATCGAGCACATATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaa      <  1:2508510/150‑1 (MQ=255)
                                                                                                                                                   ggCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaaa     <  1:3009625/150‑1 (MQ=255)
                                                                                                                                                    gCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAg    <  2:806771/150‑1 (MQ=255)
                                                                                                                                                    gCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAg    <  1:229463/150‑1 (MQ=255)
                                                                                                                                                      aTCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCa  <  2:3037034/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                      
CCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTT‑GGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCA  >  NC_000913/164432‑164730

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: ‑