Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 164,581 A→C 20.7% 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.7% 192.1 / 15.3 123intergenic (+47/‑149)hrpB/mrcBputative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base A (32/63);  new base C (6/22);  total (38/85)
Fisher's exact test for biased strand distribution p-value = 2.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 125 total aligned reads displayed.

CCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCA  >  NC_000913/164432‑164730
                                                                                                                                                     |                                                                                                                                                     
 cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAt                                                                                                                                                      >  1:451544/1‑150 (MQ=255)
 cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAt                                                                                                                                                      >  2:1200040/1‑150 (MQ=255)
 cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAt                                                                                                                                                      >  2:3088161/1‑150 (MQ=255)
  gTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATc                                                                                                                                                     >  1:480864/1‑150 (MQ=255)
          cAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGTCGGGGTGAGGGCATCAGCGCGCa                                                                                                                                             >  2:4836622/1‑150 (MQ=255)
           aaaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCAc                                                                                                                                            >  2:3332490/1‑150 (MQ=255)
            aaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACg                                                                                                                                           >  2:874066/1‑150 (MQ=255)
             aaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGt                                                                                                                                          >  2:4038042/1‑150 (MQ=255)
                 cgaTGAACGGGCGTTATCCCAACCATGTCTGGCCGGACGCCCCGGCAAATACGGCACCGACGCGACGGACGAAAAAGTTTGCTTAAGGTTGTTTTTCTCCCTCTCCCTGTGGTTGAGGGCCGGGTTGAGGGCATCAGCGCGCACGTtcac                                                                                                                                      <  2:4649209/149‑1 (MQ=255)
                  taTGAACGGGCGTTACCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGCCGGACGAAAATGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcacc                                                                                                                                     <  2:4954851/149‑1 (MQ=38)
                   aTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccc                                                                                                                                    >  1:4637082/1‑150 (MQ=255)
                     gAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTtcaccctc                                                                                                                                  >  1:3843570/1‑150 (MQ=255)
                     gAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccctc                                                                                                                                  >  1:4512819/1‑150 (MQ=255)
                     gAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcacccac                                                                                                                                  >  1:4283554/1‑148 (MQ=255)
                       aaGGGCGTTACCCCAACCAGGCCGGGCCGGACGACCCGGCAAATCCTGCCCCGACGCGACGGACGAAAAAGTATTCGTACGGTTGGTTTTCTCCCTCGCCCTGTGGGAGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTtcaccctcac                                                                                                                                <  2:1174453/150‑1 (MQ=37)
                         gggCGTTATCCCAAACATGTCTGGCCGGACGACCCGGAAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTtcaccctcaccc                                                                                                                              <  2:1187937/150‑1 (MQ=255)
                         gggCGTTATCCCAAACATGTCTGGCCGGACGACCCGGAAAATACTGAACCGACGCGACGGAGGAAAATGTATTCGTATGGTTGGTTTTCTCCCTCTCCCTGTGGGTGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccctcaccc                                                                                                                              <  2:1187996/150‑1 (MQ=38)
                          ggCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTgcacccacacccc                                                                                                                             >  2:3030252/1‑149 (MQ=255)
                            cGTTATCCCAACCATGTCTGGCCGGACGACCCGGCAAACACTGCACCGACGCGCCGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGCGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTaa                                                                                                                           <  2:3288066/150‑1 (MQ=255)
                             gTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTTTTCGTAAGGTTGGTTTTCTCCCTCGCCCTGTGGGAGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAAc                                                                                                                          <  2:1975511/150‑1 (MQ=255)
                             gTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTGGCCCCACACCCAAAc                                                                                                                          >  2:3993118/1‑150 (MQ=38)
                             gTTATCCCAAACATGTCGGGCCGGACGACCCGGAAAATACGGCACCGACGCGACGGACTAAAAAGTATTCTTAAGTTTGGTTTGCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAAc                                                                                                                          <  2:711541/150‑1 (MQ=38)
                             gTTAGCCCAACCAGGTCGGGCCGGACGACCCGGCAAATCCTGCACCGACGCGCCGGACGAAAAAGTATTCGTAAGGGTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAAc                                                                                                                          <  2:711507/150‑1 (MQ=37)
                              ttttCCCAAACTTTTGTGGCCGGACGACCCGGCAAAGACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAAcc                                                                                                                         <  2:834657/147‑1 (MQ=38)
                              ttttCCCAAACTTTTCTGGCCGGACGACCCGGAAAATACGGCACCGCCGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCGTGTGGGTGAGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAAcc                                                                                                                         <  2:834817/147‑1 (MQ=38)
                              ttATCCCACAGCTGTGGGGCCGGACGCCCCGGCAAATACTGCACCGACGCGCCGGAGGAAAATGTTTGGGTATGGTTGGTTTTCGCCCTCGCCGGGGGGGTGAGGGCGGGGGGGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAAcc                                                                                                                         <  2:4534264/150‑1 (MQ=255)
                              ttATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGCCGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAAcc                                                                                                                         <  2:834651/150‑1 (MQ=255)
                              ttATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCCGCGCGCAGGTTCACCCGCACCCGAAcc                                                                                                                         >  2:3269308/1‑150 (MQ=39)
                              ttATCCCAAACATGTCGGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGTTGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAAcc                                                                                                                         <  2:1237082/150‑1 (MQ=255)
                                  cccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATGCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTGCGCCCGCACCCAAGcccaca                                                                                                                     >  1:3248782/1‑147 (MQ=37)
                                   ccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTGCACCCGCCCCCTAAaccaaaa                                                                                                                    >  1:2030030/1‑146 (MQ=38)
                                   ccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCGCAACCAAAcccacac                                                                                                                    >  1:3695009/1‑146 (MQ=38)
                                     aaaCATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCCCCCGCACCCCAACCCCCGccc                                                                                                                  >  2:4890369/1‑150 (MQ=38)
                                      aaCATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCCGCGCGCACGGGCCCCCCCCCCCCAACCCCCAcccc                                                                                                                 >  1:4156677/1‑149 (MQ=37)
                                         atttgtGGCCGGACGACCCGGAAAATACTGCACCTACGCGACGTACTAAAATTTATTCTTATGGTGGGTTTTCTCCCTCCCCGTGTGGGTGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCaa                                                                                                              <  2:1103105/145‑1 (MQ=255)
                                         aggTCGGGCCGGACGACCCGGCAAATCCTCCCCCGACGCGCCGGACTAAACCGTATTCGTATGGTGGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCaa                                                                                                              <  2:3300421/148‑1 (MQ=25)
                                         aTGTCTGGCGGTACGCCCCGGCAAATACTGCACCGAGGCGACGGACTAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGTGAGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCaa                                                                                                              <  2:1699713/150‑1 (MQ=38)
                                          tgacTGGCCGCACGACCCGGAAAATACTCCACCGACGCGACGCACCAAAAACTATACGTAAGGTTGTTTTTCTCCCTCACCCCGTGGGAGAGGCCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                                             <  2:2545958/147‑1 (MQ=37)
                                          tGTCTGGCCGGACGACCCGGCAAATCCGGCACCCCCGCGCCGGACCAAAATTTTTTCGAACGGTTGGTTTTCCCCCTCCCCCGTTGGGAGCGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                                             <  2:3506668/150‑1 (MQ=255)
                                           gTCGGGCCGGACCACCCGCCAAACACGCCCCCGCCGCGCCGGCCGAAAACGTTTTCTTATGGTTTGTTTTCCCCCTCGCCGTGTGGGAGAGGGCCGGGGTTTGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAgg                                                                                                            <  2:3034367/150‑1 (MQ=255)
                                                 ccGCACGACCCGCCAAATCCGGCACCGACGCGACGGACGAAAATTTATTCTTACGGTTGGTTTGCCCCCCCCCCCGGTGGGAGTGGGCCGGGGTGAGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGAGAgg                                                                                                      <  2:2399012/150‑1 (MQ=255)
                                                   ggACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCCGCGCGCACGTGGCCCCACACCCGAACCCACCCCCCCACgggggggggg                                                                                                    >  2:4931599/1‑150 (MQ=25)
                                                   ggACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCACACCCGAAAACACTCCCACAAaggcgggtgg                                                                                                    >  2:4931555/1‑143 (MQ=25)
                                                    gacgccCCGGCAAATACGGCACCGACGCGACGGACGACAAAGTATTCGTAAGGTTGGTTTTCTCCCTCGCCCTGTGGGAGAGGGCCGGGGTGAGGGCCACAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGa                                                                                                   <  2:3351156/150‑1 (MQ=39)
                                                    gacgacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGa                                                                                                   <  2:5081377/150‑1 (MQ=255)
                                                    gacgacCCGGCAAACACGCCCCCCCCGCGACGTACGAAAATGTATTCTTATGGTGGGTTTTCACCCCCTCCCTGTGGTTGCGGGCCGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGa                                                                                                   <  1:4932994/150‑1 (MQ=255)
                                                    ccccccccGGAAAATACTGCACCGCCCCGACGGACTAAAAATTATTCTTACGGTTGTTTTTCTCCCTCTCCCGGTGGGAGAGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGa                                                                                                   <  1:3355653/145‑1 (MQ=255)
                                                     acgacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCCCACCCCGACCCGCGCCCCCACGGGGGAgggggg                                                                                                  >  1:989918/1‑148 (MQ=25)
                                                       gacCCGGCAAATACTGCACCGACGCGACGTACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTTAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCg                                                                                                <  2:2906863/150‑1 (MQ=255)
                                                       gacCCGCAAAACACTGCACCGCCGCGCCGGAAGAAAATGTTTTCGTATGGTTGTTTTCCTCCCTCTCCGTGTGGGAGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCg                                                                                                <  1:3214391/150‑1 (MQ=255)
                                                          ccGCCAAATACGGCACCGACGCGCCGGACTAAAAAGTATTCGTAAGGTTGGTTTCCTCCCACTCCCTGTGGGAGAGGGCGGGGGGGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATc                                                                                             <  2:122378/150‑1 (MQ=255)
                                                           cGGCAAATACTGCACCGACGCGACGGACTAAAATTTTTTCGTACGGTTGGTTTTCTCCCTCTCCCTGTGGGTGTGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCg                                                                                            <  2:4501213/150‑1 (MQ=38)
                                                           cGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCACAGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCCCACCCAAACCCTCGCCCGCAAGGGCGAGGGGGCCgcgcg                                                                                            >  1:4236200/1‑146 (MQ=37)
                                                             gtaaaTACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCCCCCACACCCAAACCCACACCCACCACGGGGAGGGGGCCCATCGac                                                                                          >  2:1907412/3‑149 (MQ=25)
                                                                   aCTGCACCTACGCGACGGACTAAAAAGTATTCGTAAGTTTGGTTTTCTCCCACTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACaaa                                                                                    <  2:3760663/150‑1 (MQ=255)
                                                                       cACCGCGGCGCCGGACTAAAATTTATTCGTTTGGTTGTTTTACCCCCACTCCGTGGGGGAGTGGGCGGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtttt                                                                                <  1:2555015/150‑1 (MQ=255)
                                                                         cgtacgcgactGACTAAAATTTATTCTTATGGTTGTTTTTCTCCCTCTCCCTTTGGTTGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                              <  1:796438/147‑1 (MQ=255)
                                                                         ccgccgCGCCGGCCGAAACTGTATGCGTATGGTGGGATTCCGCCCACTCCCTGTGGGCGAGGGCCGGGGTGTGCGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGTGGACCGATCGAGCACAAATTTTga                                                                              <  2:3091221/150‑1 (MQ=255)
                                                                         ccgccgCGACGGACTAAAAAGTATTGGTACGGTTGTTTTTCTCCCTCTCCCGGGGGGTGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                              <  2:3685889/150‑1 (MQ=38)
                                                                         ccGACGCGACTTACTAAATTTTATTGGTTTGTTTGGTTGTCGCCCTCCCCGTGTGGTTGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                              <  1:2023905/150‑1 (MQ=255)
                                                                         ccGACGCGACGGACGAAAATGTATTGGTATGGTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                              <  2:2591858/150‑1 (MQ=255)
                                                                         ccGACGCGACGGACGAAAAAGTATTCGTATGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                              <  2:2734210/150‑1 (MQ=255)
                                                                           gccgcgacgGCCGAAAATGTTTTGGTATGGTTGTTTTGCCCCCTCGCCGGGGGGGTGTGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgaga                                                                            <  2:525303/148‑1 (MQ=255)
                                                                           gacgcgacgGACGACCACGTATTCGTATGGTTGGTTTTCCCCCACTCCCTGTGGGTGCGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgaga                                                                            <  2:2712852/150‑1 (MQ=255)
                                                                           gacgcgacgGACGAAAATGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgaga                                                                            <  2:2713088/150‑1 (MQ=255)
                                                                           gacgcgacgGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGTGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgaga                                                                            <  2:4250193/150‑1 (MQ=255)
                                                                                 acggactAAAATTTATTCTTTTGTTTGGTTTTCTCCCTCTCCGTGTGGGTGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTATCCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtc                                                                      <  1:2985512/150‑1 (MQ=37)
                                                                                    gacgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttc                                                                   <  2:4700384/150‑1 (MQ=255)
                                                                                      cgAAAAAGTTTTCGTATGGTTGGTTTTCTCCCTCTCCCTGTGGGTGAGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttctt                                                                 <  2:2669266/150‑1 (MQ=255)
                                                                                      cgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCCAACCAACTCCCACAAGGGTGAGGGGAGCGATAGAGCCCAAATATTTGGAGATAtattaat                                                                 >  2:3232253/1‑147 (MQ=25)
                                                                                       gAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCACACCCAAAACCACACCCACAAAGGAGAGGGGACCGAACGAGCACAAAATTTGATAGAAAtattaatc                                                                >  2:3410738/1‑146 (MQ=25)
                                                                                        aaaaaGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttcttct                                                               <  2:5104183/150‑1 (MQ=255)
                                                                                             gTATTCGTATGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                                          <  1:3823951/150‑1 (MQ=255)
                                                                                             gTATTCGTAAGGTTGGTTTTCTCCCTCGCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                                          <  1:3824204/150‑1 (MQ=255)
                                                                                             gTATTCGTAAGGTGGGTTTTCCCCCTCGCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                                          <  2:4199474/150‑1 (MQ=255)
                                                                                               tttCGTACGGTTGTTTTTCTCCCTCTCCGTGTGGGAGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGt                                                        <  1:5139795/149‑1 (MQ=39)
                                                                                                   gTATGGTTGGTTGCCTCCCTCGCCGTGTGGGTGCGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGACTTGTAACa                                                    <  1:2538153/150‑1 (MQ=255)
                                                                                                      aGGTGGGTTTCCCCCCTCTCCCTGTGGGTGCGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaa                                                 <  2:4970512/150‑1 (MQ=255)
                                                                                                       gtttggttTTCTCCCTCTCCCTGTGGGTGAGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaag                                                <  1:2938907/148‑1 (MQ=255)
                                                                                                        gttggttTTCTCCCTCTCCCTGTGGGAGAGGGTGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaaga                                               <  1:4836622/150‑1 (MQ=255)
                                                                                                           gtttttCCCCCTCCCCGGGTGGGTGTGGGCGGGGGTGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACa                                            <  1:3935136/148‑1 (MQ=255)
                                                                                                               ccggccctctcccaCAGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCGCACCCCAACCCGCTCCCGCAAGGGAGTGGGGGCCGATCGAGCCCAAAATTTGGGAGATATATTATTATGTCTTGTAACAGGCGAACAGaaa                                        >  2:823657/5‑150 (MQ=255)
                                                                                                                  tccctctccctGTGGGAGTGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATc                                     <  1:929950/150‑1 (MQ=255)
                                                                                                                   ccctctccctGGGGGCGTGGGCCGGGGTGAGGGCATCCGCGCGCCCGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAAGCg                                    <  2:4555377/150‑1 (MQ=255)
                                                                                                                   cccacgccctGTGGGAGAGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCg                                    <  2:1242200/144‑1 (MQ=255)
                                                                                                                         ccctGTGGTTGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCttt                              <  1:1854698/150‑1 (MQ=255)
                                                                                                                             gggggTGAGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgcg                          <  1:1945056/148‑1 (MQ=255)
                                                                                                                                ggtgagGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCt                       <  1:448302/147‑1 (MQ=255)
                                                                                                                                 gagaGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTg                      <  2:3843570/150‑1 (MQ=255)
                                                                                                                                  agagGGCCGGGGTGAGGGCATCAGCGCGCACGTGCCCCCACAACCGAACCCACACCCGCCACGGCGCGGGGACCGCACGCGCACAAATATTGAGAGATATAATCTTCTGTCGTGTTACACAAGTCCAGGAAAAAGGGCGTTTGCGccggt                     >  2:2134344/1‑147 (MQ=255)
                                                                                                                                      ggCCGGGGTGAGGGCATCAGCGCGCACGTTGACCCCCACCCAAAACCACACCCACACGGGGGAGGGGCCCCAACGAGCACAAAAATTGAGAGGTATCATATTATGTCTTGTAACAGAAGAAAAGAAAAAACGGCTTTTTCGCCGGAaaat                 >  2:597083/1‑150 (MQ=255)
                                                                                                                                       gCCGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATAtt                <  2:2135134/150‑1 (MQ=255)
                                                                                                                                          ggggTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg             <  1:1365235/150‑1 (MQ=255)
                                                                                                                                          ggggTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg             <  1:2944150/150‑1 (MQ=255)
                                                                                                                                          ggggTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg             <  1:1141669/150‑1 (MQ=255)
                                                                                                                                                 gggCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaa      <  2:5094941/150‑1 (MQ=255)
                                                                                                                                                  ggCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaaa     <  1:1200040/150‑1 (MQ=255)
                                                                                                                                                  gccATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaaa     <  2:488085/148‑1 (MQ=255)
                                                                                                                                                    cATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGc   <  1:128893/150‑1 (MQ=255)
                                                                                                                                                     aTCAGCGCGCACGTTCACCCGAAACCCAACACACTCCCACAAGGGAGAGGGGACCGCTAGAGCACAAATTTTGATAGATATTTTATTCTGTATTTTAACAGAAGAACAGCAAAAACGGCTTTTGCGCCAGAAAAATGCGGAGAAAAAGCa  >  2:3868344/1‑150 (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: