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.

CCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCA  >  NC_000913/164432‑164730
                                                                                                                                                      |                                                                                                                                                      
ccGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCa                                                                                                                                                        >  1:5094194/1‑150 (MQ=255)
 cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAt                                                                                                                                                       >  1:1594318/1‑150 (MQ=255)
 cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAt                                                                                                                                                       >  1:705044/1‑150 (MQ=255)
 cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAt                                                                                                                                                       >  2:8362832/1‑150 (MQ=255)
 cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAt                                                                                                                                                       >  1:7046509/1‑150 (MQ=255)
  gTGAGGTGAAAAATGAGATGAAAGGGCGTTATCCCAACCATGTCTGGCCGGACGAGCCGGAAAATACTGCACCGACGCGACGGACTAAAAAGTATTCTTAAGGTTGGTTTTGTCGTTCTCCGTTTGGGAGTGGGCCGGGGTGAGGGCATc                                                                                                                                                      <  2:8632139/150‑1 (MQ=255)
   tGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCa                                                                                                                                                     >  2:2567434/1‑150 (MQ=255)
   tGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCa                                                                                                                                                     >  2:9183101/1‑150 (MQ=255)
   tGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCa                                                                                                                                                     >  1:525791/1‑150 (MQ=255)
   tGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCa                                                                                                                                                     >  1:4377614/1‑150 (MQ=255)
      ggTGCAAAAAGAGTTGAAGGGGGTTTAGCCAAAACTTGTCTGGCCGGCCGACCCGGAAATTCTTGCACCGACGCGACGGACGAAAATGTTTTGGTTAGTTTGTTTTTCTCCGACTCCCTGTGGGTGAGGGCCGGGGTTTGGGCATCAgcg                                                                                                                                                  <  2:7687270/150‑1 (MQ=255)
       gTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAgcgc                                                                                                                                                 >  2:498844/1‑150 (MQ=255)
        tgaaaaaaGAGATGAACGGGCGTTATCCCAAACCTGTCTGGCCGGACGACCGGGCAAATACTGCACCGACGCGACGGACGAAAATGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAgcgcg                                                                                                                                                <  2:1361906/147‑1 (MQ=38)
               agagATGAAAGGGCGTTATCCCAAACTTGTCTGGCCGGACGACCGGGCAAATACTGCACCGACGCGACGGACGAAAATGTTTTCTTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTtc                                                                                                                                         <  2:4538650/150‑1 (MQ=38)
               agagATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTGAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGCGGGCCTCAGCGCGCACGtgg                                                                                                                                         >  1:6419694/1‑148 (MQ=255)
               agagATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACgggc                                                                                                                                         >  1:6927201/1‑147 (MQ=255)
               acagaTGAACGGGCGTTATCCCAAACATGTCTGGCCGGACGCCCCGGCAACCACTGCACCGCCGCGACGGACGAAAAATTATTCGTCTGGTTGGTTTTCTCCCTCGCCCGGTGGGAGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTtc                                                                                                                                         <  2:7355752/148‑1 (MQ=255)
                  gaTGAACGGGGTTTATCCCAAACATGTGTGGCCGGACCCCCCGGAAAATACTGCCCCGACGCGACGGACGAAAAATTATTCTTAAGGTGGTTTTTCTCCCTCTCCCTGTGGGTGTGGGCGGGGGTTAGGGCATCAGCGCGCACGTtcacc                                                                                                                                      <  2:1748666/150‑1 (MQ=25)
                  gaTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGgcacc                                                                                                                                      >  1:7950076/1‑150 (MQ=255)
                   aTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTgccccc                                                                                                                                     >  2:4515347/1‑150 (MQ=255)
                      aaaGGGCGTTATCCCAAACAGGTCTGGCCGGACGACCCGGCAAATACGGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccctca                                                                                                                                  <  2:9342182/150‑1 (MQ=255)
                        tgggCTTTATCCCAAACTTGTCTGGCCGGACGACCCGGAAAATACTGCACCGACGCGCGGTACTAAAATTTTTTCTTATGGTTGTTTTTCTCCCTCTCCGTGTGGTTGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTtcaccctcacc                                                                                                                                <  2:2079031/149‑1 (MQ=255)
                        aGGGCGTTATCCCACCCATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCACCCTCTCCCTGTGGGTGAGGGCGGGGGTGAGGGCATCAGCGCGCACGTtcaccctcacc                                                                                                                                <  2:6224425/150‑1 (MQ=39)
                         gggCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTggccccacaccc                                                                                                                               >  1:478433/1‑150 (MQ=255)
                             gTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGGATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGCGGGGAGGGCATCAGCGCGCACGTggaccccccccccaca                                                                                                                           >  2:7861828/1‑146 (MQ=38)
                                aacccAAACTTGTGGGGCCGGACGACCGGGAAAATACGGCCCCGACGCCACGGACTAAAATGTTTTCGTATGGTTGGTTTTCCCCCTCGCCGTGGGGGAGCGGGCCGGGTTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCct                                                                                                                        <  2:8896863/148‑1 (MQ=255)
                                  cccAAACATGTCTGGCCGTACGACCCGGCAAATACTGCACCGACGCGACGGACTAAAAATTATTCGTACGGTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGTTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCctct                                                                                                                     <  2:4315230/150‑1 (MQ=38)
                                  cccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCGCACCCAAAaccaca                                                                                                                      >  1:3524129/1‑147 (MQ=38)
                                   ccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCGAAcccacac                                                                                                                     >  1:4984999/1‑146 (MQ=255)
                                   ccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCACACCCCAAcccacac                                                                                                                     >  1:3357763/1‑146 (MQ=255)
                                    cAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGGCACCCCCCCCCCCAcccacccc                                                                                                                    >  1:8868720/1‑145 (MQ=37)
                                    cAAACATGTCTGCCCGCACGACCCGGCAAAAACTGCACCGACGCGACGGACGAAAACGTATCCCTACGTTTGGTTTTCTCCCCCCCCCTGTGGTTGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTcc                                                                                                                   <  2:3376841/150‑1 (MQ=255)
                                         ttttCTGGCCGCACGCCCGGCCAAATACTGCCCCGACGCGCGGGACTAAAATTTTTTCTTATGTTTGGTTTTCTCCCTCCCCGTGTGGTCGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCaa                                                                                                              <  1:8320468/147‑1 (MQ=255)
                                         ttgtggggCCGCACGACCCGCAAAATAGTGCCCCGCGGGGCGGGACTAAAATTTTTTCTTTTGTTTGTTTTTCCCCCTCTCCGTGTGGGTGTGGGGGGGGTTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCaa                                                                                                              <  2:2529311/144‑1 (MQ=255)
                                          tGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCGCACCCAAACCCGCTCCCTCCAg                                                                                                             >  1:7705937/1‑150 (MQ=38)
                                           gTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAATTATTCGTATGGTTGGTTTTCTCCCTCTCCCTGTGGGTGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAgg                                                                                                            <  2:5873264/150‑1 (MQ=255)
                                               ggccggccGCCCCGGCAAATACGGCACCGACGCGCCGGACGAAAAAGTATTCGTACGTTTGGTTTTCCCCCACTCCCTGTGGTAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGgaga                                                                                                        <  2:5745192/150‑1 (MQ=37)
                                               ggCCGCACTACCGGCAAAATACGGCACCGACGCGAGGTACTAAAATTTTTTCTTACGGTTGTTTTACTCCCTCCCCCTTTGGGTTTGGGCCGGGGTGTGGGCCTCAGCGCGCACGTTCACCCCCACCCTAACCCTCTCCCTCAAGGgaga                                                                                                        <  1:3122354/150‑1 (MQ=255)
                                                gCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTAGTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCGCCCCCGAAACCACTCCCTCAAggggggg                                                                                                       >  1:3856213/1‑146 (MQ=37)
                                                  gggACGACCCGGAAAATACTTCACCTACGCGACGTACGAAAAATTATTGTTATGTTTGTTTTTAACCCTCTCCCTGTGGGATCGGGCGGGGGGCAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAggg                                                                                                     <  2:2951417/149‑1 (MQ=255)
                                                  cgtacgacCCGGCAAATACTGCACCGAGGCGACGTACTAAAAATTATTCGTATGGTTGGTTTTCTCCCTCTCCCTGGGGGTGTGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAggg                                                                                                     <  1:8475060/147‑1 (MQ=255)
                                                   ggACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGGCCCCCACCCAAAACCAAACCCACAAGGGGGCgggg                                                                                                    >  1:583176/1‑150 (MQ=25)
                                                   ggACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGGGGCCCACCCCCAAAACCACAACCACAAGGGAGAgggg                                                                                                    >  1:541593/1‑150 (MQ=25)
                                                     ccgacCCGGCAACAACGGCCCCGACGCGCCGGCCGAACAAGTTTTCGTACGGTTGGTTTTCCCCCTCTCCCTGGGGGAGCGGGCCGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAgggtta                                                                                                  <  2:8325345/149‑4 (MQ=255)
                                                     ccgacCCGGCAAACACTGCACCGACGCGACGGACGAAAATGTATTCGTACGGTTGTTTTTCTCCCTCCCCGTGTGGGAGCGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGAGAgggtta                                                                                                  <  2:4468551/149‑4 (MQ=255)
                                                     accacCCGGAAAACAGTGCACCGACGCGACGGACTAAAAAGTTTTCGTTTGGTTGTTTTTCTCCCTCTCCCTGTGGTAGCGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGAc                                                                                                  <  2:6614360/147‑1 (MQ=255)
                                                       cacCCGCAAAATACCCCCCCGCCGCGCGGGACTAAAATGTATTCGTTCGGTGGGTTTTCCCCGTCTCCCTGTGGTTGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGAGAgggttagg                                                                                                <  1:5214362/149‑6 (MQ=255)
                                                           cGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGCCCCCCACCCAAAACCACAACCACAAGGGCGAGGGGGCCGATCg                                                                                            >  2:1545603/1‑150 (MQ=25)
                                                           cGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGCCCCACACCCAAAACCGCACCCACCAGGGGGGGGGGGCCgcggg                                                                                            >  1:7636546/1‑146 (MQ=25)
                                                                    cTGCCCCGACGCGCCGGACGAAAAAGTATTCGTACGGTTGGTTGTCTCCCACTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAt                                                                                   <  2:7718288/150‑1 (MQ=39)
                                                                    cTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGGGGCCCCGCACCCGAACCCACACCCACCAGGGCGGGGGGGCCCAGCGAGcccccaa                                                                                   >  1:4533980/1‑144 (MQ=21)
                                                                    cTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCCCGGTCACCCTCACCCGAAACCACTCCCACAAGGGGGAGGGGACCGATCGAGCACACAt                                                                                   >  2:6896208/1‑150 (MQ=37)
                                                                       ccccGACGCGACGGACTAAAATGTTTTCGTATGGTTGTTTTTCTCCCTCTCCGGGTGGTTGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtttt                                                                                <  1:3937985/148‑1 (MQ=255)
                                                                        cccgccgCGCCGTCCGAAAAATTTTTTTTATGGTTGGTTTTCGCCCTCTGCGGGTGGGTGTGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTg                                                                               <  1:9327842/149‑1 (MQ=255)
                                                                         ccGACGCGCCGGACGACAACTTTTTGGTATGTTTGTTTTTCTCCCTCTCCGGGTGGTCGTGGGCGGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                              <  2:1864697/150‑1 (MQ=255)
                                                                         ccGACGCGACGTACTAAAATTTATTCTTATGTTTGGTTTTCTCCCTCTCCCTGTGGGTGTGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                              <  2:5512765/150‑1 (MQ=37)
                                                                          cgacgcgacgGACGAAAAAGTTTTCGTATGGTTGTTTTTCTCCCTCTCCGTGTGGGTGAGGGCCGGGTTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgag                                                                             <  2:6527032/150‑1 (MQ=255)
                                                                            ccgcgacgGACGAAAATGTTTTCTTTTGGTTGTTTTTCTCCCTCTCCCTGTGGTAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgagag                                                                           <  2:5174296/149‑1 (MQ=38)
                                                                                 acgcacgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACCTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtc                                                                      <  1:9070010/150‑1 (MQ=255)
                                                                                    gacgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGGCCCACACCCCAAACCACACACACAAGGGGGGGGGGGCCGCTCGAGCCCCAAAATTTAGAGataaaaaa                                                                   >  1:4680921/1‑145 (MQ=16)
                                                                                      cgAAAAAGTTTGCGTATGGTTGGTTTTCTCCCTCTCCGGGTGGGAGCGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttctt                                                                 <  1:5983761/150‑1 (MQ=255)
                                                                                      cgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttctt                                                                 <  2:9187416/150‑1 (MQ=255)
                                                                                        aaaaaGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGACCCACACCCAAAACCACACCCACAACGGGGAGGGGACCCATCCAGCACAAAATTTTATAGATAtcttaataa                                                               >  2:4679162/1‑145 (MQ=21)
                                                                                          aatgttttCGTATGTTTGGTTTTCTCCCTCTCCCTGTGGGTGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGt                                                             <  1:9183158/144‑1 (MQ=38)
                                                                                           acgTATGCGTAGGGTGGGTTTTCGCCCGCCCCCTGTGGGAGCGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTc                                                            <  2:7507718/148‑1 (MQ=255)
                                                                                             ttATTCGTATGGTTGGTTTTCTCCCTCTCCCTGTGGTTGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                                          <  2:3697747/149‑1 (MQ=255)
                                                                                             gTATTCGTATGGTTGGTTTTCTCCCTCTCCGTGTGGGAGAGGGCCGGGTTGAGGGCCTCCGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                                          <  1:533878/150‑1 (MQ=255)
                                                                                             gTATTCGTATGGTGGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                                          <  1:6661782/150‑1 (MQ=255)
                                                                                                  tcggaGGTTGGGTTTGCCCCCCCCCCCGGGGGGTGGGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTaa                                                      <  1:8173030/146‑1 (MQ=255)
                                                                                                     aaGGTTGGTTGTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACaga                                                  <  2:5317733/150‑1 (MQ=255)
                                                                                                      tgtttgtttTACCCCCTGGCCGTGGGGGAGGGGGGGGGGGTGGGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaa                                                 <  1:2567434/147‑1 (MQ=255)
                                                                                                       gtttggttTTCTCCCTCTCCGTGTGGGTGAGGGCGGGGTTGAGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaag                                                <  1:4803143/148‑1 (MQ=255)
                                                                                                        gttggttTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGTCACCCTCACCCAAAACCACTACCACAACGGAGAGGGGACCCATCAAGCCCAAATTTTTAGAGATTAATTATTCAGTGTTGTCACaaaaga                                               >  2:8252827/1‑150 (MQ=18)
                                                                                                        gttggttTTCTCCCTCGCCCTGTGGGAGAGGGCCGGGGGGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaaga                                               <  2:1129164/150‑1 (MQ=255)
                                                                                                               ttccgccctctccgtGTGGGAGTGGGCGGGGGTGAGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGaa                                         <  2:7302301/145‑1 (MQ=255)
                                                                                                                 tctccctgtccgTGTGGGATAGGGCGGGGTTGGGGGCCTCAGCGCGCACGTTCACCCTCTCCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGaaaa                                       <  2:3582767/150‑1 (MQ=255)
                                                                                                                 tctccctctcccTGTGGGTGAGGGCGGGGTTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGaaaa                                       <  1:1545603/150‑1 (MQ=255)
                                                                                                                 tcgccctctcccgGGGGGAGAGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGaaaa                                       <  1:977176/147‑1 (MQ=255)
                                                                                                                    ccctccccgtGTGGGCGAGGGCGGGGGTGGGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCg                                    <  1:3303700/150‑1 (MQ=255)
                                                                                                                     cgtctccctGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCgg                                   <  2:8498723/148‑1 (MQ=255)
                                                                                                                           cctGGGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCtttt                             <  2:4534737/150‑1 (MQ=255)
                                                                                                                                 ggtgggggCGGGGGGGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCt                       <  1:3219750/145‑1 (MQ=255)
                                                                                                                                 ggAGTGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCt                       <  1:1896138/150‑1 (MQ=255)
                                                                                                                                 ggAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCt                       <  1:1069103/150‑1 (MQ=255)
                                                                                                                                 ggAGAGGGCCGGGGTGAGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCt                       <  2:901485/150‑1 (MQ=255)
                                                                                                                                  gagaGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGAACCCCAACCAAAACCACAACCCCCAAGGGGAGGGGGCCCAAAGAGCACAAAAATTTAGAGATATAATAATATGTCGTTTAAAAAAAGAACAAAAAAAAGGGCgttttggccgg                      >  2:1401003/1‑143 (MQ=255)
                                                                                                                                       ggCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCGCACCCAAACCCACACCCACAACGGAGAGGGGCCCGATAGAGCACAAAATTTTAGGGAAAAATTATTATGTCTTGTTAAAGAAGAACAGAAAAAAGGGCTTTTTCGCCGGAatat                 >  2:6098366/1‑150 (MQ=255)
                                                                                                                                          cGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGc              <  1:29507/150‑1 (MQ=255)
                                                                                                                                           ggggggggggCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg             <  2:3292240/150‑1 (MQ=255)
                                                                                                                                           ggggTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg             <  1:8362832/150‑1 (MQ=255)
                                                                                                                                           ggggTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg             <  2:694673/150‑1 (MQ=255)
                                                                                                                                           ggggTGAGGGCATCAGCGCGCCCGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg             <  2:1041365/150‑1 (MQ=255)
                                                                                                                                           ggggTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg             <  2:5835836/150‑1 (MQ=255)
                                                                                                                                             ggTGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGga           <  2:4028380/150‑1 (MQ=255)
                                                                                                                                                 aGGGCATCAGCGCGCACGTTGACCCACAACCGAAACCACACCCACAAAGGGGAGGGGACCCAACGAGCACAAAATTTTAGAGATAAAATAATAAGTCTTTTTACAAAAAAAAAGAAAAAACGGCGTTTTCGCCGGAATAATTCGGGGaaa       >  2:3802493/1‑150 (MQ=255)
                                                                                                                                                  gggCCTCAGCGCGCCCGTTCACCCTCACCCTAACCCTCTCCCTCAATGGAGAGGGGACCGATCGAGCACATATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaa      <  1:2508510/150‑1 (MQ=255)
                                                                                                                                                   ggCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaaa     <  1:3009625/150‑1 (MQ=255)
                                                                                                                                                    gCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAg    <  2:806771/150‑1 (MQ=255)
                                                                                                                                                    gCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAg    <  1:229463/150‑1 (MQ=255)
                                                                                                                                                      aTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCa  <  2:3037034/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                      
CCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCA  >  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: