Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 164,581 A→C 17.6% 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,5810AC17.6% 377.4 / 29.7 210intergenic (+47/‑149)hrpB/mrcBputative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base A (62/107);  new base C (9/32);  total (71/139)
Fisher's exact test for biased strand distribution p-value = 9.71e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 214 total aligned reads displayed.

CCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCA  >  NC_000913/164432‑164730
                                                                                                                                                     |                                                                                                                                                     
 cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAt                                                                                                                                                      >  2:2390542/1‑150 (MQ=255)
  gTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATc                                                                                                                                                     >  2:133194/1‑150 (MQ=255)
     aGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAgc                                                                                                                                                  >  1:641440/1‑150 (MQ=255)
       gTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAgcgc                                                                                                                                                >  1:1854850/1‑150 (MQ=255)
           aaaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCAc                                                                                                                                            >  2:726837/1‑150 (MQ=255)
                 agaTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAACAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcac                                                                                                                                      >  2:4249508/1‑150 (MQ=255)
                 agaTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGtggac                                                                                                                                      >  1:101452/1‑146 (MQ=255)
                  gaTGAAAGGGCGTTATCCCAAACATGTGTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACTAAAATTTATTCTTATGGTTGGTTTTCTCCCTCGCCCTGTGGGTGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcacc                                                                                                                                     <  2:1922758/150‑1 (MQ=38)
                     gAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTgcccccgc                                                                                                                                  >  2:7427658/1‑148 (MQ=255)
                       aaGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTgcacccgcac                                                                                                                                >  2:6759403/1‑150 (MQ=255)
                        tggggtttAACCCAAACTTGTCTGGCCGGACGACCGGGCTAATACGGCACCGACGCGACGGACTAAAAAGTATTCGTACGTTTGTTTTTCCCCATCTCCGTGGGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccctcacc                                                                                                                               <  2:1780686/144‑1 (MQ=25)
                        aGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTggacccacacc                                                                                                                               >  1:1984513/1‑150 (MQ=255)
                        aGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGgcccccgcacc                                                                                                                               >  1:3367676/1‑150 (MQ=255)
                          ggCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTgcacccacaccca                                                                                                                             >  2:1976687/1‑149 (MQ=255)
                            ctttATCCCAAACTTGTGGGGCGGGACGACCCGCAAAACACTGCACCGCCGCGACGCACTAAAATTTATGCGTACGTTTGGTTTTCTCCCTCTCCCGGTGGGTGCGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTaa                                                                                                                           <  1:6803842/148‑1 (MQ=255)
                            ctttAACCCAAACCTGTGGGGCCGGCCTACCCGGAAAATACGGCACCGACGCGTGGGACTAAAATTTATTTTTTTGGTTGGTTTTCGCCCTGTCCGTGGGGGTGTGGGCGGGGGTGGGGCCCTCAGCGCGCACGTTCACCCTCACCCTaa                                                                                                                           <  1:2057548/148‑1 (MQ=255)
                            cGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAGCAGCGCGCACGTgccccctcaaacaaa                                                                                                                           >  1:6144211/1‑144 (MQ=38)
                            cGTAATCCCAAACACCACCGCCCGCACCACCCGGCAAAAACCCCACCCACCCCACGCACCAAAAACAATACGAACGCTGGCTTTACACCCAACCCCTGTGGTACAGGCCCGGGGTCAGGCCATCAGCGCGCACGTTCACCCTCACCCTaa                                                                                                                           <  2:6696829/150‑1 (MQ=255)
                             gTTATCCCAAACTTGTCTGGCCGGACGACCCGGAAAATACTGCACCGACGCGAGGGACTAAAAATTTTCCGTACGGTTGTTTTTCTCCCTCTCCGTGTGGGTGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAAc                                                                                                                          <  2:3039212/150‑1 (MQ=37)
                              ttATCCCAAACTTGTGTGGCGGGAGGACCCGTAAAATACTGCACCTAGGCGACGGACTAAAAATTTTTGTTATGGTTGTTTTTCGCCCTCTCCCTGTGTTTGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAAcc                                                                                                                         <  2:852999/150‑1 (MQ=255)
                              ttATCCCAAACATGTGTGGCCGGACGACCCGGCAAATACTGCCCCGAGGCGAGGGACTAAAAAGTATTGGTACGGTTGTTTTTCTCCCTCTCCCGGTGGGTGCGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAAcc                                                                                                                         <  2:432799/150‑1 (MQ=255)
                              ttATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGACCCACACCCaaaac                                                                                                                         >  2:318098/1‑148 (MQ=255)
                              ttATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACCCCCAAAcc                                                                                                                         >  2:6056337/1‑150 (MQ=255)
                               tATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGACCCACACCCAAAacc                                                                                                                        >  2:1095112/1‑150 (MQ=255)
                                 tCCCAAACATGTCCGGCCGTACGCCCCGGCAAACTCTCCACCGACGCTACGGACGAAAAAGAATCCGTAAGGTTGGTTTTCTCCGTCTCCCTGTGGGTGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCctc                                                                                                                      <  2:4987875/150‑1 (MQ=37)
                                   ccAACCTTGTCGGGCCGGACGACCCGGCAATTACTGCACCGACGCGACGGACTAAAAAGTATTGGTATGGTTGGTTTTCTCCCTCTCCGTGTGGGTGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCtctc                                                                                                                    <  2:7101491/150‑1 (MQ=37)
                                   ccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGACCCGCACCCAAAaccaaaa                                                                                                                    >  1:7162814/1‑146 (MQ=38)
                                    caaccctGTCTGGCCGGACGACCCGGAAAATACTGCACCGACGCGACGGACCAAAATTTATTCGTATGGTTGGTTTTCTCCATCTCCGTGTGGGTGAGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTcc                                                                                                                   <  2:4892450/144‑1 (MQ=37)
                                       aCATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATGCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTTCAACCACACCCGTAACCTCTCCCTc                                                                                                                >  2:3424619/1‑150 (MQ=38)
                                          ggTCTGGCCGGACCACCCGCCAAATACTCCACCGACGCGACGCACGAAAAAGTATTCGAAAGGTTGGTTTACCCCCACCCCCTGGGGCAGAGGGCCGGGGTGAGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                                             <  2:3352407/149‑1 (MQ=255)
                                              tGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCCCCCCCCCCCCAAACCCCACCCACAAcggcg                                                                                                         >  2:7392474/1‑148 (MQ=25)
                                              tGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGGGCCCCACACCCCAACCCCCACCCGCCAggggg                                                                                                         >  2:797835/1‑148 (MQ=25)
                                               ggCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCCCACCCGAACCCACACCCACAAGGgagc                                                                                                        >  2:4733411/1‑149 (MQ=38)
                                                gCCGGACGCCCCGGCAACTCCTCCCCCGCCGCGCCGGACGAAAAAGTATTCGTATGGTTGGTTTTCCCCCGCCCCCTGGGGGCGTGGGCCGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGagag                                                                                                       <  2:6366982/150‑1 (MQ=255)
                                                  ggcacgccCCGCCAACTACGGCACCGACGCGCCGTCCGAAAAATTATTCTTACGGTTGTTTTACTCCCTCCCCCTGGGGGAGTGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAggg                                                                                                     <  1:628594/147‑1 (MQ=255)
                                                  cGGACGCCCCGGCAAATCGGGCACCGCCGCGCGGGACTAAAATTTTTTCGTCTGGTGGGTTTTCGCCCTCGCCCGGGGGGTGCGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCACGGGAGAggg                                                                                                     <  1:5997945/150‑1 (MQ=255)
                                                  cGGACGACCCGGCAAATAGTGCACCGACGCGACGGACTAAAATGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAggg                                                                                                     <  2:8129009/150‑1 (MQ=255)
                                                       gacCCGGCAAATACTGCCCCGACGCGCCGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAgggttagg                                                                                                <  2:5349464/150‑6 (MQ=255)
                                                       gacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAACAGCGCGCACGTTCACCCACCCCCCAACCCACGCCCCCAACGGAGGGGGGGccc                                                                                                >  2:8154710/1‑149 (MQ=37)
                                                        acCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCTCCCGCACCCAAAACCACACCCACAAGGGGGGGGGGCCCGa                                                                                               >  1:5933713/1‑150 (MQ=25)
                                                             gCAAATCCTGCACCGACGCGACGGACGAAACAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAg                                                                                          <  2:6212994/150‑1 (MQ=255)
                                                             gCAAATACTGCACCGACGCGACGGACGAAAGAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTGCCCCCCCACCCAAACCCACTCCCCAAAGGGGGCGGGGGCCgggagcg                                                                                          >  2:7205950/1‑144 (MQ=25)
                                                             gCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTGACCCACACCCAAAACCACGACCACAAGGGCGAGGGGAccccgagag                                                                                          >  1:5227524/1‑143 (MQ=25)
                                                             gCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGGAAGCTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGGGGGGACCGCTCGAg                                                                                          >  2:5910109/1‑150 (MQ=38)
                                                                      gCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCCACCACACCCGAACCCACAACCACCAGGGTGAGGGGACCGAGCGATCCCAAAatt                                                                                 >  2:638784/1‑150 (MQ=25)
                                                                       ccccgccgCGACGGACGAAAATGTATTCGACCGGTGGGTTTTCCCCCACGCCCTGTGGGCGTGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtttt                                                                                <  1:7269541/148‑1 (MQ=255)
                                                                        aCCGACGCGACGGACGAAAAAGTATGCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAGCCGCGCGCACGGGCACCCGCACCCGCACCCCCTCCCGCCCGGGCGGGGGGCCCGCGCGAGCACaacggtgg                                                                               >  2:3381857/1‑144 (MQ=21)
                                                                         ccgccgCGACGGACTAAAATGTATTCTTATGTTTGTTTTCCCCCCTCTCCGTGGGGGTGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                              <  1:3478796/150‑1 (MQ=255)
                                                                         ccGACGCTCGGGACGAAAATGTTTTCGTACGGTGGTTTTTCCCCCTCCCCCTGTGGGTGTGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                              <  2:7804406/150‑1 (MQ=255)
                                                                         ccGACGCGTGGGAGTAAAATGTTTTGGTTTGGTTGGTTTTCTCCCTCTCGGTGGGGGTGTGGGCGGGGGTGGGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                              <  1:6776729/150‑1 (MQ=255)
                                                                          cgccgcgacgGACTAAAAAGTTTTCGTTCGGTTGGTTTTCTCCCTCGCCCGGGGGGAGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgag                                                                             <  1:1976687/147‑1 (MQ=255)
                                                                               cgacgGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCCCCAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGata                                                                        >  2:7825376/1‑150 (MQ=255)
                                                                                gacggacgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCCCCCACACCCCAAACCACAACCACACAGGGGGGGGGGCCGAACGGGCCCAAAATATTAGAGatat                                                                       >  1:7987441/1‑150 (MQ=16)
                                                                                   gtacgAAAAAGTATTCGTATGGTTGGTTTTCTCGCTCTCGCTGTGGGTGAGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtctt                                                                    <  1:4393630/148‑1 (MQ=255)
                                                                                   ggacgAAAATGTATTGGTACGGTTGGTTTTCTCCCTCGCCGTGTGGGTGTGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtctt                                                                    <  1:424062/150‑1 (MQ=255)
                                                                                        aaaaaGTATTCGTATGGTTGTTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttcttct                                                               <  2:2540934/150‑1 (MQ=255)
                                                                                        aaaaaGTATTCGTACGGTTGGTTTTCTCCCTCTCCCTGTGGGTGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttcttct                                                               <  2:8138737/150‑1 (MQ=255)
                                                                                          aatgTATTCGAATGGTTGGTTTTCGCCCTCGCCCTGGGTTTGATGGCCGGGGAGAGTGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGt                                                             <  2:158769/147‑1 (MQ=255)
                                                                                           attttttCTTAAGTTTGGTTTTCCCCCCCCCCGGGGGGGTGCGGGCGGGGGTGAGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTc                                                            <  2:3815173/145‑1 (MQ=255)
                                                                                           aattATTCGTAAGGTTGGTTTTCTCCCTCTCCCGGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTc                                                            <  2:5205346/147‑1 (MQ=255)
                                                                                             tttttCGTACGGTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                                          <  2:4574825/147‑1 (MQ=255)
                                                                                             gTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGGGGGAGAGGGCGGGGGTGTGGGCATCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                                          <  1:7490340/150‑1 (MQ=255)
                                                                                              ttttCGTATGGTTGGTTTTCCCCCTCTCCGTGTGGGTGTGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTg                                                         <  1:3059655/148‑1 (MQ=255)
                                                                                              tATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTg                                                         <  2:6192405/150‑1 (MQ=255)
                                                                                               aTTCGTAAGGTTGGTTTTCTCCCTCGCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGt                                                        <  1:2036318/150‑1 (MQ=255)
                                                                                                 tCGTACGGTTGGTTTACGCCCACCCCCTGGGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTaa                                                      <  2:6953983/150‑1 (MQ=255)
                                                                                                  cGTAGGGTTGGTTTTCTCCCACACCCTGGGGGAGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAAc                                                     <  2:2286866/150‑1 (MQ=255)
                                                                                                   gTATGGTTTGTTTTCCCCCTCGCCGGGGGGTTGTGGGCGGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACa                                                    <  1:1674931/150‑1 (MQ=255)
                                                                                                   gTATGGTTGTTTTTCCCCCCCGCCGGGGGGGTGTGGGCGGGGGTGCGGCCCTCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACa                                                    <  1:2390542/150‑1 (MQ=255)
                                                                                                    gaaGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTATGTCTTGTAACag                                                   <  2:6628168/149‑1 (MQ=255)
                                                                                                     ttggttgtttTTCCCCCCCTCCGGGGGGTTGTGGGGGGGGGTGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACaga                                                  <  1:77789/148‑1 (MQ=255)
                                                                                                     atggttggttTTCTCCCACCCCGGGGGGGAGTGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACaga                                                  <  2:2925456/148‑1 (MQ=255)
                                                                                                     atggttggttTCCTCCCTCTCCCTGTGGGTGTGGGCCGGGGTGTGGGCATCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACaga                                                  <  1:3230715/148‑1 (MQ=255)
                                                                                                      aGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaa                                                 <  2:1883205/150‑1 (MQ=255)
                                                                                                         ttgtttTTCCCCCCCCCCGGGTGGGTGTGGGCGGGGTTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAgaagaa                                              <  1:1095112/150‑1 (MQ=255)
                                                                                                         ttggttTTCTCCCTCGCCCTGTGGGAGTGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAgaagaa                                              <  1:4035410/150‑1 (MQ=255)
                                                                                                           ggttTTCGCCCTCGCCCTGTGGGTGTGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACa                                            <  2:5812381/150‑1 (MQ=255)
                                                                                                                   ccctctccctGTGGTTGAGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCg                                    <  2:3211383/150‑1 (MQ=255)
                                                                                                                          cctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCAAAACCACACCCGCAAGGGGGGGGGGGCCCATCGAGCACAAAATTTGTGTGTTATCATCTTCTGTCTGGGTACAGAAGAACAGCAAAACGGGcggtt                             >  2:7445306/1‑146 (MQ=19)
                                                                                                                           ctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCGCCCCCACCCAAACCCACACCCACAACGGAGAGGGGGCCGCTCGAGCACAAAAATTTATAGATATATTATTATTTATTGTAACAGAAGAACAGAAAAACGGGCttgtt                            >  1:4726506/1‑147 (MQ=255)
                                                                                                                              tgGGAGCGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgcgc                         <  1:1756884/150‑1 (MQ=255)
                                                                                                                               gggTGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGcc                        <  2:5933713/150‑1 (MQ=255)
                                                                                                                               gggCGAGGGCCGGGGTGCGGGCATCCGCGCGCCCGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAACATCGGGCTTTTGCGcc                        <  1:2608847/150‑1 (MQ=255)
                                                                                                                               gggAGAGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGcc                        <  2:6788097/150‑1 (MQ=255)
                                                                                                                               gggAGAGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGcc                        <  2:2005203/150‑1 (MQ=255)
                                                                                                                                ggAGGGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCt                       <  1:7392421/150‑1 (MQ=255)
                                                                                                                                     gggCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAata                  <  2:2720307/150‑1 (MQ=255)
                                                                                                                                      ggCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAGGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAatat                 >  1:971339/1‑150 (MQ=255)
                                                                                                                                      ggCCGGGGTGAGGGCATCAGCGCGCACGTGCAACCACACCCCAACCCACAACCACAAAGGAGATGGGACCCCTCCCGCACAAATTTTGTGCGATAAATTAATATTTCTTGTACACCAAGGACAGCAAAACCGGCTTTTGCGccggcaaaa                 >  2:3852220/1‑143 (MQ=255)
                                                                                                                                       gcgggggTGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATAtt                <  1:3037414/147‑1 (MQ=255)
                                                                                                                                        ccGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTg               <  1:4099535/150‑1 (MQ=255)
                                                                                                                                          ggggTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg             <  2:7751998/150‑1 (MQ=255)
                                                                                                                                          ggggTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGAGATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg             <  2:6598783/150‑1 (MQ=255)
                                                                                                                                          ggggTGAGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg             <  2:3619300/150‑1 (MQ=255)
                                                                                                                                            ggTGAGGGCCTCAGCGCGGACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGga           <  2:6316760/150‑1 (MQ=255)
                                                                                                                                            ggTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGga           <  2:2245899/150‑1 (MQ=255)
                                                                                                                                              tgagcccATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGgaga         <  2:2094491/144‑1 (MQ=255)
                                                                                                                                                  ggCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaaa     <  2:6486117/150‑1 (MQ=255)
                                                                                                                                                   ccATCAGCGCGCACCCTCACCCACACCCTAACCCTCACCCTCAACGGAGAGGGGACCGATCGAGCACAAATATAGAGAGATATCTACTTCTGTCTTGTAACACAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAg    <  2:7217270/149‑1 (MQ=255)
                                                                                                                                                     aTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCTAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCa  <  1:726800/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: