Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 164,581 A→C 15.3% 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,5810AC15.3% 195.3 / 10.8 116intergenic (+47/‑149)hrpB/mrcBputative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base A (27/69);  new base C (2/18);  total (29/87)
Fisher's exact test for biased strand distribution p-value = 1.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 120 total aligned reads displayed.

CCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAAT‑ACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCA  >  NC_000913/164432‑164730
                                                                                                                                                      |                                                                                                                                                     
ccGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAAT‑ACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCa                                                                                                                                                       <  2:2719185/150‑1 (MQ=255)
 cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAAT‑ACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAt                                                                                                                                                      <  2:4142420/150‑1 (MQ=255)
       gTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAAT‑ACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAgcgc                                                                                                                                                >  2:3782485/1‑150 (MQ=255)
       gTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAAT‑ACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAgcgc                                                                                                                                                >  2:1984728/1‑150 (MQ=255)
        tGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAAT‑ACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCcgcgcg                                                                                                                                               >  1:2490508/1‑150 (MQ=255)
              atgtgaTGAAAGGGCGTTATCCCAAACTTGTCTGGCGGGACGACCCGGAAAAT‑ACTGCACCGACGAGACGGACGAAAAAGTATTCGTATGGTTGGTTTTCTCCCTCTCCCTGTGGGTGTGGGCCGGGGTGAGGGCATCAGCGCGCACGtt                                                                                                                                         <  2:2700444/146‑1 (MQ=38)
              aaGAGATGATTGGGCTTTATCCCAAACTTGTGGGGCCGTACGACCGGGAAAAT‑ACTGCACCGACGCGACGGACTAAAATGTTTTCTTATGTTTGTTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGTTGTGGGCATCAGCGCGCACGtt                                                                                                                                         <  2:3866992/150‑1 (MQ=21)
              aaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAAT‑ACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCCCggg                                                                                                                                         >  1:9558/1‑148 (MQ=255)
              aaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAAT‑ACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGtg                                                                                                                                         >  2:3912156/1‑149 (MQ=255)
                 agaTGACCGGGCGTTATCCCCAACATGTCTGGCGGGACGCCCCGGCACAT‑ACTGCACCGAGGCGACGGACGAAAAAGTATTCGTACGGTTGGTTTTCTCCCTCTCCCTGTGGGTGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcac                                                                                                                                      <  2:3055815/150‑1 (MQ=37)
                 agaTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAAT‑ACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGtggcc                                                                                                                                      >  1:2003930/1‑146 (MQ=255)
                  tattaatgggCGTTATCACAAACTTTTGTGGCCGGACGACCCGGCAAAA‑ACTGCACCGAGGCGCCGTACTAAATTTTTTTCTTTTGGTTGTTTTTCCCCCTCCCCGTGTGGGAGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTtcacc                                                                                                                                     <  2:3785614/143‑1 (MQ=255)
                  gatgatTGGGCGTTATCCCAAACTTGTCTGGCCGGACGACCCGGCAAAT‑ACTGCACCGACGCGAGGGACGAAAAAGTTTTCGTACGGTTGGTTTTCGCCCTGTCCCTGTGGGTGTGGGCCGGGGTGTGGGCATCAGCGCGCACGTtcacc                                                                                                                                     <  2:893818/150‑1 (MQ=255)
                  gaTGCACGGGCTTTATCCCAAACTGGTCTGGCCGGACGACCCGGAAAAT‑CCTGCACCGCCGCGACGGACGAAAAAGTTTTCGTATGGTTGGTTTTCTCCCTCTCCGTGTGGGAGTGGGCCGGGGTGTGGGCATCAGCGCGCACGTtcacc                                                                                                                                     <  2:5310542/150‑1 (MQ=255)
                    tGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAAT‑ACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGtgggcccg                                                                                                                                   >  1:1911418/1‑143 (MQ=255)
                      aaaGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAAT‑ACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccctca                                                                                                                                 >  2:3968180/1‑150 (MQ=255)
                       acgggCGTTATCCCAATCTTGTCGGGCCGGACGCCCCGGCAAAT‑ACTGCCCCGACGCGAGGGACTAAAAAGTATTCGTATGGTTGGTTTTCTCCCTCTCCCTGGGGTTGTGGGCCGGGGTGCGGGCATCAGCGCGCACGTtcaccctcac                                                                                                                                <  2:4180123/148‑1 (MQ=255)
                        tgggCTTTAACCAAAACTTGTCTGGCCGGACGACCCGGAAAAT‑ACTGCACCGACGCGACGGACTAAAATTTATTGTTAGGGTTGGTTTACCCCCACGCCCTGTGGGTTTGGGCGGGGGTGTGGGCATCAGCGCGCACGTtcaccctcacc                                                                                                                               <  2:2102292/149‑1 (MQ=255)
                        cgggCTTTACCCCAAACAGGTCTGGCCGGACGACCCGCCAAAT‑ACTGCCCCGACGCGCCGGACGAAAATTTTTTGTTTTGGTTGGTTTTCTCCCTCTCCGTGTGGGTGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTtcaccctcacc                                                                                                                               <  2:2109126/149‑1 (MQ=25)
                             gTTACCCCAATCTTGTCGGGCGGGCCGA‑CCGGTAAATGCCGGCCCCGACGGGACGTACTAAAATGTATTGGTAGGGTTGGTTTTCTCCCCCTCCCTGTGGGAGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTatc                                                                                                                          <  2:4210093/150‑3 (MQ=255)
                              ttATCCCAAACATGTCTGGCCGGACGACCCGGCAAAT‑ACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGTCCCACACCCaaaac                                                                                                                         >  1:965308/1‑148 (MQ=38)
                                aacccAAACATCACTGCCCGCACGACCCCCAAAAA‑ACTCCACCCACCCCACGCACGAAAAACTATTCGTAAGGTTGCTTTACACCCACACCCTCTGGGACAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCct                                                                                                                       <  2:3207817/148‑1 (MQ=255)
                                    cAAACATGTCTGGCCGGACGACCCGGAAAAT‑ACGGCACCGACGCGACGGACGAAAAAGTATACTTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTcc                                                                                                                   <  2:2774557/150‑1 (MQ=38)
                                         ttttCTGGCCGGCCCCCCCGCAAAAT‑AGTGCACCGCCGCGAGGTACCAAAATTTATTCTTACTTTTGTTTTCCCCCCCCCCCGTGGGGGTGCGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCaa                                                                                                              <  1:1720012/147‑1 (MQ=255)
                                         ttGTCGGGCCGGACGACCCGGAAATT‑ACTGCACCGAGGCGACGTACTAAAAAGTATTCGTATGGTTGGTTTTCTCCCTCCCCCTGGGGTTGAGGGCCGGGTTTAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCaa                                                                                                              <  2:1773130/149‑1 (MQ=255)
                                         atgtggggCCGGACGACCCGGCAAAT‑ACGGCACGGACGCGACGGACTAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCaa                                                                                                              <  2:4376797/144‑1 (MQ=38)
                                         aTGTCGGGCCGGACGACCCGCCAAAT‑ACTGCCCCGACGCGCCGTACGAAAATGTATTCGTACGGTTGGTTTTCTCCCTCTCCCTGGGGGCTCGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCaa                                                                                                              <  2:4376748/150‑1 (MQ=255)
                                          ttgggggCCGGACGCCCCGGAAAAT‑CCTGCACCGCGGCGAGGGACGAAAATGTTTTCGTTGGGTTGGTTTTCGCCCTCTCCCTGTGGGTGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                                             <  1:2648209/145‑1 (MQ=25)
                                          tGTCGGGCCGGACGACCCGGAAAAT‑ACTCCACCTACGCGACGGACTAAAATGTATTCGTATGGTTGGTTTTCCCCCTCTCCCTTTGGGATTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                                             <  2:3893837/150‑1 (MQ=37)
                                            tCTGGCCGGACGACCCGGCAAAT‑ACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGGCCCCCACCCCAACCCACAccccccacgg                                                                                                           >  2:641800/1‑143 (MQ=37)
                                              tGGCCGGACGACCCGGCAAAT‑ACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTGGACCCGCACCCGAAACCACGCCcacaccggcg                                                                                                         >  1:3001512/1‑144 (MQ=37)
                                              tGGCCGGACGACCCGGCAAAT‑ACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCGCACCCGAACCCTCTCCCGCCAGGgag                                                                                                         >  2:792226/1‑150 (MQ=255)
                                              tGGCCGGACGACCCGGCAAAT‑ACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGCCCCACACCCAAACCCGCTCCCCCAAgggtg                                                                                                         >  1:4397130/1‑148 (MQ=37)
                                              tGGCCGGACGACCCGGCAAAT‑ACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGAGCCCACACCCAAACCCACTCACACAAgggtg                                                                                                         >  1:3001511/1‑148 (MQ=37)
                                                gCCGGACGACCCGGAAAAT‑ACTGCACCGCCGCGACGGACGAAAATTTTTTCGTATGGTTGGTTTCCCCCCCCTCCCTGGGGGAGCGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGagag                                                                                                       <  1:4899699/150‑1 (MQ=255)
                                                 ccGGACGACCCGGCAAAT‑ACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGAGAgg                                                                                                      <  2:5427290/150‑1 (MQ=255)
                                                 ccGGACGACCCGGCAAAT‑ACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGAGAgg                                                                                                      <  2:897620/150‑1 (MQ=255)
                                                 ccGGACGACCCGCAAAAT‑ACGCCACGGACGCGCGGGACTAAAATTTTTTCGTCTGGTTGTTTTTACCCCGCGCCCTGTGGGAGCGGGGGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAgg                                                                                                      <  1:5226318/150‑1 (MQ=255)
                                                  cGGCCGCCCCGGCAAAT‑CCGCCCCCGCCCCGCCGGAGGAAAACGTCTTCGTACGGTGGGTTTCCTCCCTCGCCCGGGGGGCGTGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAggg                                                                                                     <  2:2922165/150‑1 (MQ=255)
                                                  cGGACGACCCGGCAAAT‑ACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCCCACCCAAACCCTCTCCCGCAAGGGAGgggg                                                                                                     >  2:4135861/1‑150 (MQ=38)
                                                  cGGACGACCCGGCAAAT‑ACTGCACCGACGCGACGGACGAAAAAGTATTCGGACGGGTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGGGCGCACGGTCACCCTCACCCGCACCCTCTCCCGCCCGGGAGAggg                                                                                                     >  2:4135392/1‑150 (MQ=37)
                                                    gccgccccGGAAAAA‑AGTGCACCTAGGCGCCGGCCTAAAATTTTTTCTTATGGTTGTTTTTCTCCCTCTCCGGGGGGGTGTGGGCGGGGTTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGa                                                                                                   <  1:1139393/145‑1 (MQ=255)
                                                          ccGGCAAAT‑ACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCCCCCCCACCCAAACCCCCCCCCCAAAGGGGGTGGGGGcccaac                                                                                             >  2:1710160/1‑146 (MQ=25)
                                                             gCAAAG‑ACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAg                                                                                          <  2:5042921/150‑1 (MQ=255)
                                                               aaaT‑ACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGGGGCCCACCCCCAAACCCCCGCCCACACGGGGGGGGGGGCCCCTCGAGcc                                                                                        >  2:2603040/1‑149 (MQ=21)
                                                                aaT‑ACGGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGGTTTCTCCCTCTCCCGGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGTGAGGGGACCGATCGAGcac                                                                                       <  2:3530845/150‑1 (MQ=255)
                                                                     cTGCACCGACGCGACGGACGAACAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAt                                                                                   <  2:823151/150‑1 (MQ=255)
                                                                      tGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGAGCCACACCCCAAACCACACCCACCACGGGGAGGGGCCCGCTCGGGCCCaaaat                                                                                  >  1:809734/1‑148 (MQ=21)
                                                                      tGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGGCCCACAACCAAACCCACTACCACCAGGGAGGGGGGCCCGGGAGAGCACAAAta                                                                                  >  2:266111/1‑149 (MQ=21)
                                                                        cACCGACGCGACGGACTAAAATGTTTTCGTAAGGTTGGTTTTCTCCCTCTCCCGGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtttt                                                                                <  2:543414/150‑1 (MQ=255)
                                                                         aCCGACGCGACGTACCAAAATGTTTTCGTATGTTTGTTTTCCCCCCTCCCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTg                                                                               <  2:5149604/150‑1 (MQ=38)
                                                                          ccGACGCGACGGACTAAAAAGTATTCGTACGGTTGGTTTTCTCCCTCTCCCTGTGGGTGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                              <  1:3968700/150‑1 (MQ=255)
                                                                          ccGACGCGACGGACGAAAAATTTTTCGTACGGTTGGTTTTCTCCCTCTCCCTGTGGGTGTGGGCCGGGGTGCGGGCATCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                              <  1:4555868/150‑1 (MQ=38)
                                                                          ccGACGCGACGGACGAAAAAGTATTCGTAAGGTGGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                              <  2:4455943/150‑1 (MQ=255)
                                                                           cgacgcgaggGACGAAAATGTTTTCGTACGGTTGGTGTTCTCCCTCCCCCGGTGGGTGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgag                                                                             <  1:1238378/150‑1 (MQ=38)
                                                                            gaggcgacgGACGAAAATGTTTTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgaga                                                                            <  2:2662313/147‑1 (MQ=255)
                                                                             acgcgccgGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgagag                                                                           <  2:3929264/150‑1 (MQ=255)
                                                                                            aattATTCGTAAGGTTGGTTTTCTCCCCCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTc                                                            <  1:5466108/147‑1 (MQ=255)
                                                                                              gTATTCGTATGGTTGGTTTTCTCCCTCTCCGTGGGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                                          <  2:2209148/150‑1 (MQ=255)
                                                                                                  tCGTATGTTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTaa                                                      <  2:3026337/150‑1 (MQ=255)
                                                                                                      atgtttggttTCCTCCCTCCCCGTGGGGTGGTGGGCCGGGTTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACaga                                                  <  1:2165767/146‑1 (MQ=255)
                                                                                                      aatgttggttTTCCCCCTCTCCCTGTGGGAGAGGGCCGTCGTGAGGGCATCACCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGCTGACCGATCGAGCACAAATTTAGAGAGATATATTCTTCTGTCTTGTAACaca                                                  <  1:1830041/147‑3 (MQ=38)
                                                                                                       tggttggttTTCCCCCCCTCCGGGGGGGTGCGGGCGGGGTGGTGGGACTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaa                                                 <  2:9558/149‑1 (MQ=255)
                                                                                                       aGGTTGGTTTTCACCCTCACCCTGTGGGTGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaa                                                 <  1:4788362/150‑1 (MQ=255)
                                                                                                           tggttTTCCCCCTCGCCGTGTGGGTGCGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAAc                                             <  1:304733/150‑1 (MQ=255)
                                                                                                                  ctccctctccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCCCCCGAGCCCAAACCCACACCCAACAAGGAGAGGGGCCCGCTCGAGCACAAATGTTGAGAGAAATCATATTCAGTGTTGTAACAGAACAACAGaaaaa                                      >  2:2285069/1‑149 (MQ=16)
                                                                                                                   tccctctccgtGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATc                                     <  2:1760805/150‑1 (MQ=255)
                                                                                                                   tcccccccccgGTGGGAGAGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATc                                     <  2:2030071/150‑1 (MQ=255)
                                                                                                                   cccccccccctGCGGCAGAGGCCCGGGGAGCGCCCCTCAGCGCCCACGTTCACCCACACCCTAACCCTCACCCACAACGGAGAGGGGACCGATCGAGCACAAAATTTGAGAGATATCTTCTTCTGACTTGTAACAGAAGAACAGAAAATc                                     <  2:1208798/143‑1 (MQ=255)
                                                                                                                    ccctctccgtGTGGGAGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCg                                    <  1:2252226/150‑1 (MQ=255)
                                                                                                                    ccctctccctGTGGGTGAGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCg                                    <  1:2252437/150‑1 (MQ=255)
                                                                                                                      ctctccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCggg                                  <  2:3479033/150‑1 (MQ=255)
                                                                                                                       tcgcccgGTGGGTGTGGGCGGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGc                                 <  1:3949615/147‑1 (MQ=255)
                                                                                                                           cctGGGGGTGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCtttt                             <  1:1216789/150‑1 (MQ=255)
                                                                                                                           ccacagggAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGCCCCACACCCAAACCCACACCCCCAAGGGGGCGGGGCCCCAAAGAGCCCAAAAGTTTAGAGATAAAATAGTCTGTCTTGTAAAACAATAACAGCAAAAAGGGCtttt                             >  2:1800096/6‑150 (MQ=255)
                                                                                                                             tgtgGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgc                           <  2:3501762/150‑1 (MQ=255)
                                                                                                                             tgggggAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgc                           <  2:5175163/147‑1 (MQ=255)
                                                                                                                               tgGGTGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgcgc                         <  1:3444659/150‑1 (MQ=255)
                                                                                                                               ggggAGAGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgcgc                         <  2:3931816/149‑1 (MQ=255)
                                                                                                                                gggggggggCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGcc                        <  2:924817/144‑1 (MQ=255)
                                                                                                                                 ggtgtgggCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCt                       <  2:1029852/145‑1 (MQ=255)
                                                                                                                                 ggcgagGGCGGGGGTGTGGGCCTCATCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACAGATCGAGCACAAATTTTGAGAGATATCTTCGTCTGTCTTGTAACAGCAGAACAGAAAATCGGGCTTTTGCGCCt                       <  1:496026/147‑1 (MQ=255)
                                                                                                                                 ggAGAGGGCCGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCt                       <  2:129068/150‑1 (MQ=255)
                                                                                                                                 ggAGAGGGCCGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCt                       <  1:496067/150‑1 (MQ=255)
                                                                                                                                  gagcGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTg                      <  2:4274172/150‑1 (MQ=255)
                                                                                                                                   agagGGCCGGGGGGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGa                     <  2:5325084/150‑1 (MQ=255)
                                                                                                                                    gtgggCCGGGTGGCGGGCATCAGCGCGCACGTTCACCCTCACCCGAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGaa                    <  2:4831338/148‑1 (MQ=255)
                                                                                                                                     cgggCCGGGGGGAGGGCCTCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAat                   <  1:1324637/149‑1 (MQ=255)
                                                                                                                                       ggCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCAAAACCACACCCACAAGGGAGTGGGGGCCGCTAGTGCCCAAAAATTTAGTGATATAATATTAAGTGTTGTTAAAGAAGTACACAAAATCGGGCGTTTTCGCCGGAatat                 >  2:2746333/1‑150 (MQ=255)
                                                                                                                                          gggggTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGc              <  1:4808749/149‑1 (MQ=255)
                                                                                                                                          cGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACCGAAAATCGGGCTTTTGCGCCTGAATATTGc              <  2:3770657/150‑1 (MQ=255)
                                                                                                                                           ggggTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg             <  2:4376681/150‑1 (MQ=255)
                                                                                                                                           ggggTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg             <  1:2603040/150‑1 (MQ=255)
                                                                                                                                             ggTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGga           <  1:196616/150‑1 (MQ=255)
                                                                                                                                             ggTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCAATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGga           <  2:1523559/150‑1 (MQ=255)
                                                                                                                                                   ggCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaaa     <  1:3491620/150‑1 (MQ=255)
                                                                                                                                                    gCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAg    <  1:437586/150‑1 (MQ=255)
                                                                                                                                                      aTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCa  <  2:3324774/150‑1 (MQ=255)
                                                                                                                                                      aTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCa  <  2:446817/150‑1 (MQ=255)
                                                                                                                                                      aTCAGCGCGCACGTTCACCCACAACCAAAACCAAAACCACAAGGGAGAGGGGACCCCTAGAGCCCAAAATTTTAGAGAAATAATATTATGTCTTGTAACAGAAGAACACAAAAACGGGCATTTGCGCCAGAAAATTTCTGATAAAAAGCa  >  2:2284516/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
CCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAAT‑ACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCA  >  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: ‑