Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 164,579 G→C 5.1% intergenic (+45/‑151) 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,5790GC 5.1% 820.6 / 2.4 256intergenic (+45/‑151)hrpB/mrcBputative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base G (102/138);  new base C (4/12);  total (106/150)
Fisher's exact test for biased strand distribution p-value = 1.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 264 total aligned reads displayed.

TACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAA  >  NC_000913/164430‑164727
                                                                                                                                                     |                                                                                                                                                    
tACCGTGTGGTGCAAAAAGAGATGAATGGGCGTTAGCCCAAACATGTCGGGCGGGACGACCCGGCAAATAGTGCACCGAGGCGACGGACGAAAAAGTATTCGTACGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAggg                                                                                                                                                      <  2:966357/150‑1 (MQ=38)
  cggTGAGGTGCAAAAAGAGATGAACGGGCGTTATCCCAAACCTGTCTGGCCGGACGACCCGGAAAATACTGCACCGACGCGACGGACGAAAAAGTTTTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCa                                                                                                                                                    <  2:51897/148‑1 (MQ=255)
  ccGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCa                                                                                                                                                    <  2:51806/150‑1 (MQ=255)
  ccGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCa                                                                                                                                                    >  2:8515143/1‑150 (MQ=255)
    gTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATc                                                                                                                                                  >  1:3840441/1‑150 (MQ=255)
    gTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATc                                                                                                                                                  >  1:6543634/1‑150 (MQ=255)
         gTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAgcgc                                                                                                                                             >  1:7764147/1‑150 (MQ=255)
           gCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCgcgcgc                                                                                                                                           >  1:5380533/1‑150 (MQ=255)
                 agagATGAACGGGCTTTATCCCAAACTTGTCTGGCGGGACGACCCGGCAAATACTGCACCGACGCGAGGGACTAAAATGTATTCGTATGTTTGGTTTTCGCCCTCGCCCTGTGGGAGTGGGCCGGGGTGCGGGCATCAGCGCGCACGTtc                                                                                                                                     <  2:8910086/150‑1 (MQ=255)
                   agaTTAAAGGGGGTTACCCCAAACTGGTGGGGCCGGACGACCCGCCAAATACTCCCCGGACGCGACGGACTAAAATGTATTCGTATGGTTGTTTTTCTCCCTCTCCCTGTGGGTGTGGGCCGGGGTGCGGGCATCAGCGCGCACGTtcac                                                                                                                                   <  2:476279/150‑1 (MQ=255)
                    tattaatgggCTTTAACCAAAACCTGTCGGGCGGTACGACCCGGAAATTACTGCCCCGACGCGAGTTAATAAAATTTTTTCTTATGGTTGTTTTTCCCCCTCCCCCTGTGGTTTGGGGCGGGGTTGCGGGCATCAGCGCGCACGTtcacc                                                                                                                                  <  2:3142937/143‑1 (MQ=255)
                    gaTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGtgggcc                                                                                                                                  >  1:9568471/1‑145 (MQ=255)
                       gAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccctc                                                                                                                               >  2:7958573/1‑150 (MQ=255)
                       gAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGtcacccac                                                                                                                               >  2:6443624/1‑148 (MQ=255)
                         atgggCTTTAACCCAAACTGGTCGGGCCGGACGACCCGGAAAATACGGCACCGACGCTACGGACGAAAATTTATTCGTATGGTTGGTTTTCTCCCTCGCCCTGTGGGAGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTtcaccctcac                                                                                                                             <  2:709183/148‑1 (MQ=25)
                         aaGGGCGTTATCCCAACCTGGTCTGGCCGGACGACCCGTCAAATCCTGCACCGACGCGACGGACGAAAATTTATTCGTATGGTTGGTTTTCTCCCTCTCCCTGTGGGTGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTtcaccctcac                                                                                                                             <  2:617059/150‑1 (MQ=37)
                         aaGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccctcac                                                                                                                             >  1:10883788/1‑150 (MQ=255)
                          tgggCGTTAGCCCAAACATGGCTGGCCGGACGACCCGGCAAATCCTGCACCGACGCGACGGACGAAAATGTTTTCTTATGGTTGGTTTTCTCCCTCTCCCTGTGGGTGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccctcacc                                                                                                                            <  2:4574145/149‑1 (MQ=37)
                           gggCGTTATCCCCAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAACAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGGtgccccccaccc                                                                                                                           >  2:10199823/1‑150 (MQ=38)
                            ggCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCAGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTgggcccccacccc                                                                                                                          >  1:3150405/1‑149 (MQ=38)
                                ttATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGACCCCCAACCaaaac                                                                                                                      >  2:240274/1‑148 (MQ=255)
                                ttATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCGCCCCCACCCCaaac                                                                                                                      >  2:7819694/1‑148 (MQ=255)
                                ttATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAGCAGCGCGCACGTGGACCCACACCCCAAcc                                                                                                                      >  1:6369974/1‑150 (MQ=38)
                                ttATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTGTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTgcacccgcccccacacc                                                                                                                      >  1:10745336/1‑145 (MQ=38)
                                ttATCCCAAACATGTCGGGCCGCACGACCCGGAAAATACTGCACCGACGCGACGGACGAAAATGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAAcc                                                                                                                      <  2:1695908/150‑1 (MQ=255)
                                ttATCCCAAACAGGTCTGGCCGGACGACCCGGCAAATACGGCCCCGAGGCTAGTGACTAAAATTTTTTCTTATGTTTGTTTTTCTCCCTACCCCGGTGGGTGCGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAAcc                                                                                                                      <  2:293729/150‑1 (MQ=255)
                                taaccccAAACAGGTGGGGCCGGCCGACCCGGCAAATACTGCACCGCGGCGACGGACTAAAAAGTTTTGGTATGGTTGGTTTTCGCCCTCCCCGTGGGGGAGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAAcc                                                                                                                      <  2:5616464/146‑1 (MQ=255)
                                   ttccggcgcaaaagccctattttctgttcttctgttacaagacagaagaagatatctctcaaaatttgtgctcgatcggtcccctctccctggaGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCctc                                                                                                                   >  1:5770286‑M1/95‑150 (MQ=255)
                                    cccAACCTTGTGGGGCCGGACGCCCCGGCATTTCCTGCCCCGCCGCGCGGGACTAAAATTTTTTCGTATGTTTGGTTTTCTCCCTCTCCGGGGGGGTGTGGGGGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCctct                                                                                                                  <  2:6678617/150‑1 (MQ=255)
                                     ccAATCTTGGCTGGCCGGACCACCGGGAAATTCCGGCCCCTCCGCGCCGGACGAAAATGTATTCTTTGGGTTGGTTTTCGCCCTCTCCCTGGGGGTGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCtctc                                                                                                                 <  1:9849500/150‑1 (MQ=255)
                                     ccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGACCCGCACCCAAAcccacac                                                                                                                 >  1:8093122/1‑146 (MQ=255)
                                      caaccttGTCGGGCGGGACGCCCGGGCAAATCCGGCCCCGCCGCGTGGGCCGAAAATGTTTTCGTTGGGTGGGTTTTCTCCCTCGCCGTGGGGGTGTGGGCGGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCtctcn                                                                                                                <  1:10373179/144‑2 (MQ=255)
                                      cAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGGCCCCCACCCAAAacccacac                                                                                                                >  2:7402179/1‑145 (MQ=37)
                                          cATGTCGGGCCGTCCGCCCCGCAAAAAACGCCCCCGCCGCGACGGACTAAAATTTATTCGTATGTTTGGTTTACCCCCACCCCCGGGGGTAGAGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCa                                                                                                            <  2:2053471/150‑1 (MQ=255)
                                           ctgcgtGGCCGGACGACCCGGCAAATACGGCCCCGACGCGACGGACTAAAAAGTTTTCGTATGGTTGTTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGGGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCaa                                                                                                           <  2:5029656/145‑1 (MQ=37)
                                           aTGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCaa                                                                                                           <  2:2684807/150‑1 (MQ=255)
                                            tgtggggccggccGACCCGCAAAATACTCCCCCGCCGCGCGGGACTAAAAATTTTTCTTTTGGTTGGTTTTCGCCCCCGCCGTGGGGTTGGGGGGGGGGGTGGGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                                          <  1:6223955/145‑1 (MQ=255)
                                                tGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTGGCCCCACCCCCAAACCCCCACCCGCaccgggg                                                                                                      >  1:9826758/1‑144 (MQ=25)
                                                gggccggccGCCCCGGCAAATACGGCCCCGCCGCGCCGGCCGAACAAGTATTCGTACGGTGGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGCGGGCCTCCGCGCGCCCGTTCACCCTCACCCTAACCCTCTCCCTTGAGGgag                                                                                                      <  2:10037604/149‑1 (MQ=25)
                                                 ggccggccTCCCCGCAAAATACTCCCCCGACGCGACTCACTAAAATTTATTGTTTTGGTGGTTTTACCCCCTCGCCGTGTGGGTGCGGGCGGGGGTGCGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGgaga                                                                                                     <  1:2137710/150‑1 (MQ=255)
                                                 ggCCGGACCACCCGGCAAATACTGCACCTACGCGACGGACTAAAATGTTTTCGTATGGTTGTTTTTCTCCAACGCCCTGTGGGTGAGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGgaga                                                                                                     <  2:573528/150‑1 (MQ=255)
                                                   ccGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCACACCCAAACCCACACCCGCCAgggggcgg                                                                                                   >  1:5469430/1‑145 (MQ=37)
                                                       acgacCCGGAAAATACTGCCCCGCGGCGACGGACTAAAATTTATTCGTATGTTTGTTTTTCTCCCTCTCCCTGTGGGTGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGAc                                                                                               <  2:5997348/150‑1 (MQ=255)
                                                         gacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGGCCCGCACCCAAAACCACAACCCCAAGCGGGAGGGGCCCg                                                                                             >  1:8476966/1‑150 (MQ=25)
                                                         gacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCCACCCCAAAACCACAACCACCAGGGGGAGGGGGCCg                                                                                             >  1:1865546/1‑150 (MQ=25)
                                                            ccgtaaaaTCCTGCACCGCCGCGCGGGACCAAAACTTTTTGGTAGGGTTGGTTTGCCCCCCCCCCGTGGGGGTGCGGGCGGGGGGGCGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATc                                                                                          <  2:3860258/145‑1 (MQ=255)
                                                               gCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGGGGGCCCGCCCCCAAAACCACACCCCCCAGCGAGCGGGGCCCGAACGAg                                                                                       >  1:683384/1‑150 (MQ=21)
                                                               gCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCGCACCCACGCCCGCCCGGGAGAGGGGGCCGATCGAg                                                                                       >  2:9814813/1‑150 (MQ=37)
                                                                 aaaTACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCGCACCCAAACCCACACCCACCAGGGCGAGGGGGCCGATCGAGcc                                                                                     >  1:4639582/1‑149 (MQ=37)
                                                                      cTGCCCCGCGGCGACGGACTAAAATGTAGTCGTATGGTTGGTTTTCTCCCTCCCCGGGGGGGCGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAt                                                                                <  1:702286/150‑1 (MQ=255)
                                                                         aagacagaagaagatatctctcaaaatttgtgctcgatcggtcccctctcccttgaGGGAGAGGGTTAGGGTGAGGGCATCAGCGCGCACGTtcacc                                                                                                                                  >  1:10262255‑M1/57‑97 (MQ=255)
                                                                           ccgccgCGACGGACTAAAATGTTTTCGTATGTTTGTTTTTCTCCCTCTCCCTGTGGGAGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                           <  2:8476966/150‑1 (MQ=38)
                                                                           ccGACGCGACGGACTAAAAATTATTCGTATGTTTGTTTTTCTCCCTCTCCGTGGGGTCGCGGGCGGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                           <  1:6056096/150‑1 (MQ=255)
                                                                           ccGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                           <  2:8214071/150‑1 (MQ=255)
                                                                           ccGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                           <  2:2789610/150‑1 (MQ=255)
                                                                             gacgcgacgcACTAAAATTTTTTCGTATGGTTGGTTGTCGCCCTCTCCCTGTGTGAGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgaga                                                                         <  1:6852931/150‑1 (MQ=255)
                                                                              agaagaagatatctctcaaaatttgtgctcgatcggtcccctctcccttgaGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTa                                                                                                                         >  2:1535925‑M1/52‑101 (MQ=255)
                                                                              agaagaagatatctctcaaaatttgtgctcgatcggtcccctctcccttgaGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTa                                                                                                                         <  1:1535925‑M1/50‑1 (MQ=255)
                                                                                     ggacgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCCCCCACACCCCAACCCACAACCACAAGGGACAGGGGCCCGAACGACCCCAAAATTTTGGAGATAtctt                                                                 >  2:10836548/1‑150 (MQ=21)
                                                                                          aaaaaGTATTCGTACGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttcttct                                                            <  2:10046251/150‑1 (MQ=255)
                                                                                          aaaaaGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttcttct                                                            <  1:2686761/150‑1 (MQ=255)
                                                                                             atgTATTCGTATGGTTTGTTTTCTCCCTCTCCCTGTGGGAGCGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTc                                                         <  1:10659855/148‑1 (MQ=255)
                                                                                             aaGTATTCGTATGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTc                                                         <  1:7958573/150‑1 (MQ=255)
                                                                                               gTATTCGTAAGGTTGTTTTTCTCCCACTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                                       <  2:3618661/150‑1 (MQ=255)
                                                                                                 aTTCGTATGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGt                                                     <  1:7361234/150‑1 (MQ=255)
                                                                                                   tCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGTCACCCTCAACCGAACCCACTCCCACAAGGGAGAGGGGACCGATCGAGCACAAAATTTGAGTGATATCTTATTCTGTGTTGTaa                                                   >  1:1001509/1‑150 (MQ=37)
                                                                                                     gtttggttgtttTTCTCCCTCTCCCTGGGGGTGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACa                                                 <  2:683384/146‑1 (MQ=39)
                                                                                                        cggttggttTTCCCCCTCTCCCTGTTGGAGCGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaa                                              <  1:2910042/149‑1 (MQ=255)
                                                                                                        cggttggttTACTCCCTCCCCCTGGGGGAGTGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGTTATCTTCTTCTGTCTTGTAACagaa                                              <  1:921224/149‑1 (MQ=255)
                                                                                                        aGGTTGTTTTTCGCCCTCTCCCTGTGGGTTTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaa                                              <  2:1281411/150‑1 (MQ=255)
                                                                                                         gtttggttTTCGCCCGCTCCGTGGGGGTGTGGGCGGGGGTGAGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaag                                             <  2:187998/148‑1 (MQ=255)
                                                                                                          gttgtttTTCCCCCGCGCCCGGTGGGCGAGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaaga                                            <  2:5155676/150‑1 (MQ=255)
                                                                                                          gttggttTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCCCACCCCAACCCGCCCCCCCAAGGGAGAGGGGGCCGCGCGAGCCCCAAATTTTAGTGATATAATCATCGGTCGTGTAAAagaaga                                            >  2:2573804/1‑150 (MQ=18)
                                                                                                           ttggttTTCTCCCTCGCCGGGGGGGTGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAgaagaa                                           <  1:9375051/150‑1 (MQ=255)
                                                                                                               ttttCTCCCTCTCCCTGTGGGAGAGGGCCGGGGGGAGGGCCTCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGa                                       <  2:5898068/150‑1 (MQ=255)
                                                                                                                      cctctccctGTGGGTGTGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCgg                                <  1:5600894/150‑1 (MQ=255)
                                                                                                                       ctctccctGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGGTCACCCTCACCCGAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCggg                               <  2:9049874/150‑1 (MQ=255)
                                                                                                                       ctcgccgtggGGGAGTGGGGGGGGGTGCGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCggg                               >  1:4144460‑M1/11‑150 (MQ=255)
                                                                                                                          tccctGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTGGACCCGCACCCAAACCCCCGCCCACCAGGGCGGGGGGGCCCCGCGCGCCCAAAGTGTTGGTGGTATAAAGATCTGGGGTGTAACAGAAGAACACGAAATCgggggg                            >  2:7910296/1‑147 (MQ=255)
                                                                                                                             ctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCACACCCGAACCCACACCCTCAAGGGAGAGGGGACCCAACGAGCACAAATTTTGAGAGATATATTCTTCTGTCGTTTAACAGAAGAAAAGAAAATCGGGCGTTTg                         >  2:9101082/1‑150 (MQ=255)
                                                                                                                              tgtgGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGGGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgc                        >  2:7051134/1‑150 (MQ=255)
                                                                                                                              gggGGGAGCGGGCCGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgc                        >  1:1947221‑M1/4‑150 (MQ=255)
                                                                                                                                ggggtgtgggCGGGGGTGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgcgc                      <  1:10086408/143‑1 (MQ=255)
                                                                                                                                gGGGAGTGGGCGGGGGTGAGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgcgc                      >  2:2304225‑M1/2‑150 (MQ=255)
                                                                                                                                 ggtcgcGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGcc                     >  2:9552340‑M1/7‑150 (MQ=255)
                                                                                                                                  gtcgtGGGCGGGGGGGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCt                    >  1:10563457‑M1/6‑150 (MQ=255)
                                                                                                                                  GGAGAGGGCCGGGGTGAGGCCATCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCt                    >  1:2447466‑M1/1‑150 (MQ=255)
                                                                                                                                   GAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTg                   >  2:3578751‑M1/1‑150 (MQ=255)
                                                                                                                                        ggCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACCCCCAAAACCACACCCGCAATGGGGAGGGGGCCGATCGAGCCCAAAATTTGAGATATATATTATTATGTCTTTTAAAAGAAGAAAAGAAAATACGGCTTTTTCGCCAGAaaat              >  2:6401717/1‑150 (MQ=255)
                                                                                                                                         GCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCtctc                                                                                                                 <  2:2441762‑M1/50‑1 (MQ=255)
                                                                                                                                           gggggggCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGc           <  2:3068804/149‑1 (MQ=255)
                                                                                                                                           gggggTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGc           <  1:2324145/149‑1 (MQ=255)
                                                                                                                                            ggggTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg          <  2:3382118/150‑1 (MQ=255)
                                                                                                                                            ggggTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg          <  1:486328/150‑1 (MQ=255)
                                                                                                                                            ggggTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATACCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg          <  1:1308967/150‑1 (MQ=255)
                                                                                                                                            ggggTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg          <  2:8667172/150‑1 (MQ=255)
                                                                                                                                             gggTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCgg         <  1:10426756/150‑1 (MQ=255)
                                                                                                                                              ggTGAGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGga        <  1:355232/150‑1 (MQ=255)
                                                                                                                                                    ggCCTCAGCGCGCACGTTCACCCTCACCCGACCCCTCTCCCTCAAGGGAGAGGGGACCGAGCGAGCACAAATTTTGAGAGATATCTTCTTCGGTCGTGTAACCGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaaa  <  2:10539784/150‑1 (MQ=255)
                                                                                                                                                    ggCATCAGCGCGCACGTTCACTCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGATATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaaa  <  1:515053/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAA  >  NC_000913/164430‑164727

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: