Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 164,620 C→T 5.4% intergenic (+86/‑110) 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,6200CT 5.4% 481.5 / 11.1 148intergenic (+86/‑110)hrpB/mrcBputative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base C (52/87);  new base T (3/5);  total (56/92)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 149 total aligned reads displayed.

ATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTG‑TCT‑‑TGT‑‑‑‑A‑ACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGA  >  NC_000913/164473‑164769
                                                                                                                                                   |                                                                                                                                                             
aTGTCGGGCCGGACGACCCGGCAAATCCTGCACCGACGCGACGGACGAAAAAGTTTTCGTACGGTGGGTTTTCCCCCCCGCCGGGTGGGTGAGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCaa                                                                                                                                                             <  2:188782/150‑1 (MQ=255)
    gtGGCCGGACGACCGGGCAAATCCCGCACCGACGCGCCGGACGAAAATGTATGCGTATGGTTGGTTTGCGCCCTCTCCCTGTGGGTGCGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGga                                                                                                                                                         <  1:3962434/149‑1 (MQ=255)
    cgggCCGGACGCCCCGGAAAATACTGCACCGACGCGACGGACTAAAATGTATTCGTATGGTTGGTTTGCCCCCTCTCCGGGGGGGAGTGGGCGGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGga                                                                                                                                                         <  1:3962516/148‑1 (MQ=255)
    cTGGCCGGACGACCCGGCAAATACTGCACCGAGGCGACGGACGAAAACTTTTTGGTTCGGTTGGTTTTCGCCCTCTCCCTGTGGTAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGga                                                                                                                                                         <  2:2802809/150‑1 (MQ=37)
       gCCGGACGCCCCTCAAAAACCTCCACCGCCGCGCCGGCCTAAAATTTTTTCGTAGGTTTGTTTTCCCCCCACGCGGTGGGGGCGTGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGagag                                                                                                                                                      <  1:585107/150‑1 (MQ=255)
         gggACGACCCGGCAAATCCTGCACCGACGCGCCGGACGAAAATTTTTTCGTATGGTTGGTTTCCTCCCTCTCCCTGGGGGTGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAggg                                                                                                                                                    <  1:3054329/149‑1 (MQ=38)
         cgcacgacCGGGCCAATACTGCACCGACGGGAGGGACTTAAAAGTATTCTTACTGTGGGTTTTCTCCAACTCCGTGTGGGAGTGGGCCGGGGGGAGGGCATCAGCGCGCCCGGTCTCCCTCACCCTAACCCTCTCCCTTGAGGGAGAggg                                                                                                                                                    <  2:1686575/147‑1 (MQ=255)
         cGGCCGACCCGGCAAATACTGCCCCGACGCGAGGGACTAAAAAGTTTTCTTACGGTTGGTTTTCGCCCCCCCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAggg                                                                                                                                                    <  2:3254692/150‑1 (MQ=255)
         cGGACGACCCGGCAAATCCTGCACCGCCGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAggg                                                                                                                                                    <  2:592033/150‑1 (MQ=255)
         cGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTGGCCCCACACCCCGACCCCCACCCACCAGGGGGAggg                                                                                                                                                    >  1:1127625/1‑150 (MQ=25)
               acCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGACCCCCACCCACAACCACACCCCCAAAGGAGAGGGGGcccc                                                                                                                                              >  2:2055698/1‑148 (MQ=25)
               acCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGGAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGGCACCCGCACCCACAACCGCACCCGCCAGGGGGGGGGGGCCgc                                                                                                                                              >  2:4089293/1‑149 (MQ=25)
                   ggCAATCACGGCACCGCCGGGCCGGACTAAAATGTATTCGTATGGTTGGTTTCCGCCCTCTCCGTGGGGGTGCGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGa                                                                                                                                          <  1:2749863/150‑1 (MQ=255)
                   ggCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCGGCGCGCACGTTCACCCGCACCCGAAACCAATACCTCAAGGGGGTGGGGAGCggtggg                                                                                                                                          >  2:3312620/1‑145 (MQ=25)
                   ggCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGACCCGCACCCGAACCCACTCCCACCAGGGGGCGGGGGCCgcggga                                                                                                                                          >  2:3312590/1‑145 (MQ=25)
                      aaaaACTCCACCGACGCGACGGACTAAAATGTATTCGTACGGTTGGTTTACCCCCTCTCCGTGTGGGTGTGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGca                                                                                                                                       <  1:3000004/150‑1 (MQ=255)
                      aaaTACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGTGACCCACACCCGAACCCTCACCCACAACGGGGAGGGGGCCCAtacagcc                                                                                                                                       >  2:4050333/1‑144 (MQ=25)
                       aaTACTGCACCGACGCGACGGACGAAAAAGTATTCGTATGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGcac                                                                                                                                      <  2:4119482/150‑1 (MQ=255)
                         tACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCAACCACACCCCAACCCAAACCCACAAGGGAGGGGGGACCGCGCGAGCCCaa                                                                                                                                    >  2:1326277/1‑150 (MQ=25)
                           cTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGCCCCGCACCCCAACCCACACCCACAACGGCGGGGGGCCCCATAGCGCACaaaa                                                                                                                                  >  1:3012772/1‑149 (MQ=21)
                            tGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGGGCGCCCCCCCCCAACCCCCCAACCCCAACGGGGAGGGGGCCCCGCGAGCACcaaat                                                                                                                                 >  2:96788/1‑148 (MQ=16)
                                ccgccgCGCCGGACTAAAACGTATTCGTACGGTTGGTTTTCTCCCTCGCCGTGTGGTTGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                                                                             <  2:3281316/150‑1 (MQ=255)
                                ccGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                                                                             <  2:3391177/150‑1 (MQ=255)
                                        acggactAAAAAGTTTTCGTATGGTTGGTTTTCTCCCTCTCCCTGTGGGTGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtc                                                                                                                     <  2:285867/150‑1 (MQ=255)
                                               aaaaaGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGCCCCACACCCCAAACCGCACCCGCAAGGGCGCGGGGGCCCAGCGAGCCCCAATTTTGTGGGATGtcttgttct                                                                                                              >  2:2025801/1‑150 (MQ=16)
                                                  attttttCGTATGTTTGTTTTTCCCCCTCGCCCGGTGGTTGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTG‑Tc                                                                                                          <  2:811252/145‑1 (MQ=255)
                                                  aaGTATTCGTACGGTTGTTTTTCTCCCTCTCCCTGTGGGTGAGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTG‑Tc                                                                                                          <  2:953701/150‑1 (MQ=255)
                                                         gtttcggttggttTACCCCCTCCCCGTGTGGGTGTGGGCGGGGGTGTGGGCCTCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTG‑TCT‑‑TGT‑‑‑‑A‑Ac                                                                                            <  1:2021792/145‑1 (MQ=255)
                                                          gTAAGGTTGTTTTTCTCCCTCTCCCTGTGGGAGAGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTG‑TCT‑‑TGT‑‑‑‑A‑ACa                                                                                           <  2:3012568/150‑1 (MQ=255)
                                                           tAAGGTTGGTTTTCTCCCTCTCCCTGTGGGTGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTG‑TCT‑‑TGT‑‑‑‑A‑ACag                                                                                          <  1:2385231/150‑1 (MQ=255)
                                                            atggttggttTGCTCCCTTTCCGTGTGGGTGTGGGCCGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTG‑TCT‑‑TGT‑‑‑‑A‑ACaga                                                                                         <  1:864319/148‑1 (MQ=255)
                                                            aaGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTG‑TCT‑‑TGT‑‑‑‑A‑ACaga                                                                                         <  2:1366492/150‑1 (MQ=255)
                                                            aaGGTTGGTTTTCTCCCTCCCCCTGGGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTG‑TCT‑‑TGT‑‑‑‑A‑ACaga                                                                                         <  1:2976384/150‑1 (MQ=255)
                                                             aGGGTGGTTTTCTCCCTCTCCCTGGGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTG‑TCT‑‑TGT‑‑‑‑A‑ACagaa                                                                                        <  2:2409315/150‑1 (MQ=255)
                                                                tttgttTACTCCCCCTCGATGGGGGAGAGGGCGGGGTTGATGGCTTCAGCGCGCACGTTCACCCTCACCCTATCCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCTCAAATTTTGAGAGATATCTTATTCTG‑TCT‑‑TGT‑‑‑‑A‑ACAgaagaa                                                                                     <  1:1352225/147‑1 (MQ=255)
                                                                 tggttTTCTCCCGCGCCCTGTGGGAGAGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTG‑TCT‑‑TGT‑‑‑‑A‑ACAGAAGAAc                                                                                    <  1:4119438/150‑1 (MQ=255)
                                                                  ggttTTCTCCCGCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCGAACCCTCTCCCTCAAGGGAGAGGGGACCGAACGAGCACAAATTTTGAGAGATATCTTCTTCTG‑TCT‑‑AGT‑‑‑‑A‑ACAGAAGAACa                                                                                   <  2:297121/150‑1 (MQ=255)
                                                                     tttCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGTCACCCGCCCCCGAACCCACGCCCACAATGGGGAGGGGCCCGATCCAGCCCAAATTTTGAGAGATATCTTGTTCTG‑TCT‑‑TGT‑‑‑‑A‑ACAGAAGAACAGaa                                                                                >  2:2059319/1‑150 (MQ=25)
                                                                           cctctccctGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGGTCACCCCACCCCAAACCCAAACCCCCACCGGCGAGGGGACCCGTCGGGCCCAAAATTTGAGAGATATATTGTTCTG‑TAG‑‑AGT‑‑‑‑A‑ACAGAAGGACCGAAAAGCgg                                                                          >  2:307917/1‑150 (MQ=12)
                                                                           cctctccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCAAAACCTCACCCTCAAGGGGGTGGGGGCCGATCGAGCACAAATTTTGAGGGATATATTCTTCTG‑TCT‑‑TGT‑‑‑‑G‑AAAGAAGTAAAGAAAATCgg                                                                          >  2:3094714/1‑150 (MQ=37)
                                                                                ccctGGGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTG‑TCT‑‑TGT‑‑‑‑A‑ACAGAAGAACAGAAAATCGGGCttt                                                                     <  1:2620588/150‑1 (MQ=255)
                                                                                 cctGGGGGAGTGGGCCGGGGTGAGGGCCTCAGCGCGCCCGTGCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTG‑TCT‑‑TGT‑‑‑‑A‑ACAGAAGAACAGAAAATCGGGCtttt                                                                    <  1:1562785/150‑1 (MQ=255)
                                                                                    gtgGGTGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTG‑TCT‑‑TGT‑‑‑‑A‑ACAGAAGAACAGAAAATCGGGCTTTTgcg                                                                 <  1:1557801/150‑1 (MQ=255)
                                                                                     tgGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCACACCCAAAACCACAACCACCAGGGGGAGGGGGCCCATCGAGCCCAAAATTTGAGAGATAAAATATTATG‑TAT‑‑TGT‑‑‑‑A‑ACAGAAAAACATCAAAAAGGGCGTTTtcgc                                                                >  2:1376930/1‑150 (MQ=255)
                                                                                     ggggtgcgggCGGGGGTGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTG‑TCT‑‑TGT‑‑‑‑A‑ACAGAAGAACAGAAAATCGGGCTTTTgcgc                                                                <  1:1437939/143‑1 (MQ=255)
                                                                                     ggggggtgggCGGGGGGGGGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTG‑TCT‑‑TGT‑‑‑‑A‑ACAGAAGAACAGAAAATCGGGCTTTTgcgc                                                                <  2:377491/143‑1 (MQ=255)
                                                                                     ggggAGTGGGCCGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTG‑TCT‑‑TGT‑‑‑‑A‑ACAGAAGAACAGAAAATCGGGCTTTTgcgc                                                                <  1:3867645/149‑1 (MQ=255)
                                                                                     ggggAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTG‑TCT‑‑TGT‑‑‑‑A‑ACAGAAGAACAGAAAATCGGGCTTTTgcgc                                                                <  1:3884677/149‑1 (MQ=255)
                                                                                                cGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTG‑TCT‑‑TGT‑‑‑‑A‑ACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGc                                                     <  1:4108338/150‑1 (MQ=255)
                                                                                                cGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTG‑TCT‑‑TGT‑‑‑‑A‑ACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGc                                                     <  1:434050/150‑1 (MQ=255)
                                                                                                 ggggTGAAGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTG‑TCT‑‑TGT‑‑‑‑A‑ACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg                                                    <  2:610582/150‑1 (MQ=255)
                                                                                                   ggTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTG‑TCT‑‑TGT‑‑‑‑A‑ACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGga                                                  <  1:2711821/150‑1 (MQ=255)
                                                                                                        ggccATCAGCGCCCACGTTCACCCTCACCCTACCCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAATTCGAGAGATATCTTCTTCTG‑TCT‑‑TGT‑‑‑‑A‑ACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaa                                             <  1:458293/147‑1 (MQ=255)
                                                                                                             tCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTG‑TCT‑‑TGT‑‑‑‑A‑ACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCAt                                        <  2:3969176/150‑1 (MQ=255)
                                                                                                              cAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTG‑TCT‑‑TGT‑‑‑‑A‑ACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATg                                       <  1:406393/150‑1 (MQ=255)
                                                                                                                 cgcgcACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTG‑TCT‑‑TGT‑‑‑‑A‑ACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGcc                                    <  2:2929586/150‑1 (MQ=255)
                                                                                                                       cGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTG‑TCT‑‑TGT‑‑‑‑A‑ACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAAt                              <  2:484181/150‑1 (MQ=255)
                                                                                                                         ttCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTG‑TCT‑‑TGT‑‑‑‑A‑ACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGa                            <  2:4072574/150‑1 (MQ=255)
                                                                                                                         ttCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTG‑TCT‑‑TGT‑‑‑‑A‑ACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGa                            <  2:2805426/150‑1 (MQ=255)
                                                                                                                         ttCACCCTCAACCTAACACAAAACCAAAAGGGAGAGGGGCCCGATCGAGCACAAATTTTGAGAGATAAATTCTTCTG‑TAT‑‑TGT‑‑‑‑A‑ACAGAAGAACAGAAAATAGGGCTTTTTCGCCGGAATATTTCGGAGAAAAAGCAAGGCCCGGAATGa                            >  2:2806713/1‑150 (MQ=255)
                                                                                                                           caccctcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTG‑TCT‑‑TGT‑‑‑‑A‑ACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGAcc                          <  2:4051553/150‑1 (MQ=255)
                                                                                                                               ctcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTG‑TCT‑‑TGT‑‑‑‑A‑ACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGACTGACCGCGa                      >  1:3410238/1‑150 (MQ=255)
                                                                                                                               ctcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTG‑TCT‑‑TGT‑‑‑‑A‑ACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGa                      >  1:1357841/1‑150 (MQ=255)
                                                                                                                               ctcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTG‑TCT‑‑TGT‑‑‑‑A‑ACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGa                      >  1:2733391/1‑150 (MQ=255)
                                                                                                                               ctcacccgACCCCTCTCCCTCACGGGAGAGGGGCCCGATCGAGCACAAAGGTTGAGAGATATCTTCTTCTG‑TCT‑‑TGT‑‑‑‑A‑ACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATGGCGGAGAAAAAGCATGGCCGGGAATGACCGCGa                      <  1:2005657/150‑1 (MQ=255)
                                                                                                                                tcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTG‑TCT‑‑TGT‑‑‑‑A‑ACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAg                     <  1:1655749/150‑1 (MQ=255)
                                                                                                                                tcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTG‑TCT‑‑TGT‑‑‑‑A‑ACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAg                     <  2:338824/150‑1 (MQ=255)
                                                                                                                                 caccctACCCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTG‑TCT‑‑TGT‑‑‑‑A‑ACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGc                    <  2:822990/150‑1 (MQ=255)
                                                                                                                                    cctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTG‑TCT‑‑TGT‑‑‑‑A‑ACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAACAGCATGGCCGGGAATGACCGCGAGCCaa                 >  2:2263268/1‑150 (MQ=255)
                                                                                                                                    cctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTG‑TCT‑‑TGT‑‑‑‑A‑ACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCaa                 >  2:2263306/1‑150 (MQ=255)
                                                                                                                                    cctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTG‑TCT‑‑TGT‑‑‑‑A‑ACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCaa                 >  2:481164/1‑150 (MQ=255)
                                                                                                                                      tAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGATCACAAATTTTGAGAGATATCTTCTTCTG‑TCT‑‑TGT‑‑‑‑A‑ACAGAAGAACAGAAAATCGGGCTTTTGCGCATGATTATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAAtt               <  2:2014534/150‑1 (MQ=255)
                                                                                                                                      tAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTG‑TCT‑‑TGT‑‑‑‑A‑ACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGTATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAAtt               <  1:3440302/150‑1 (MQ=255)
                                                                                                                                      tAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTG‑TCT‑‑TGT‑‑‑‑A‑ACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAAtt               <  2:773930/150‑1 (MQ=255)
                                                                                                                                      tAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTG‑TCT‑‑TGT‑‑‑‑A‑ACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAAtt               <  1:1260942/150‑1 (MQ=255)
                                                                                                                                      tAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTG‑TCT‑‑TGT‑‑‑‑A‑ACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAAtt               <  2:48429/150‑1 (MQ=255)
                                                                                                                                      gtccccTCTCCCTTGAGGGAGAGGGTACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTG‑TCT‑‑TGT‑‑‑‑A‑ACagaa                                                                                        <  1:11191/74‑1 (MQ=37)
                                                                                                                                      gtccccTCTCCCTTGAGGGAGAGGGTACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTG‑TCT‑‑TGT‑‑‑‑A‑ACagaa                                                                                        >  2:11191/4‑77 (MQ=37)
                                                                                                                                       aaCCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTTTTCTGTTCTTCTGTTACAAGACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGa                      <  1:2036110/150‑1 (MQ=255)
                                                                                                                                       aaCCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTG‑TCT‑‑TGT‑‑‑‑A‑ACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTg              <  1:1616628/150‑1 (MQ=255)
                                                                                                                                         cccTCTCCCTTGAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTG‑TCT‑‑TGT‑‑‑‑A‑ACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGa            >  1:853933/1‑150 (MQ=255)
                                                                                                                                         cccTCTCCCTCACGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGAGATCTTCTTCTG‑TCT‑‑TGT‑‑‑‑A‑ACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGa            <  1:4036839/150‑1 (MQ=255)
                                                                                                                                         cccTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTG‑TCT‑‑TGT‑‑‑‑A‑ACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGa            <  2:1031077/150‑1 (MQ=255)
                                                                                                                                         cccTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTG‑TCT‑‑TGT‑‑‑‑A‑ACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGa            >  2:1948620/1‑150 (MQ=255)
                                                                                                                                          ccTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTG‑TCT‑‑TGT‑‑‑‑A‑ACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGAc           >  1:3230224/1‑150 (MQ=255)
                                                                                                                                           ctctCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTG‑TCT‑‑TGT‑‑‑‑A‑ACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGACTGACCGCGAGCCAATTGGACg          >  1:2246295/1‑150 (MQ=255)
                                                                                                                                           ctctCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTG‑TCT‑‑TGT‑‑‑‑A‑ACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACg          >  1:94228/1‑150 (MQ=255)
                                                                                                                                           ctctCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTG‑TCT‑‑TGT‑‑‑‑A‑ACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACg          >  1:1784146/1‑150 (MQ=255)
                                                                                                                                           ctctCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTG‑TCT‑‑TGT‑‑‑‑A‑ACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACg          >  1:2172320/1‑150 (MQ=255)
                                                                                                                                             ctcCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTG‑TCT‑‑TGT‑‑‑‑A‑ACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCa        >  1:2029186/1‑150 (MQ=255)
                                                                                                                                              tcCCTCTAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTG‑TCT‑‑TGT‑‑‑‑A‑ACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCaa       <  2:2429102/150‑1 (MQ=255)
                                                                                                                                                ccTCACGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTG‑TCT‑‑TGT‑‑‑‑A‑ACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAg     <  2:3627882/150‑1 (MQ=255)
                                                                                                                                                ccTCAAGGGAGAGGGGACCGATCGAGCACAACTTTTGAGAGATATCTTCTTCTG‑TCT‑‑TGT‑‑‑‑A‑ACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAg     <  2:2125696/150‑1 (MQ=255)
                                                                                                                                                ccTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTG‑TCT‑‑TGT‑‑‑‑A‑ACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAg     >  2:2999534/1‑150 (MQ=255)
                                                                                                                                                 cTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTG‑TCT‑‑TGT‑‑‑‑A‑ACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAgg    >  2:365189/1‑150 (MQ=255)
                                                                                                                                                  tCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTG‑TCT‑‑TGT‑‑‑‑A‑ACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAggg   >  2:755779/1‑150 (MQ=255)
                                                                                                                                                  tCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTG‑TCT‑‑TGT‑‑‑‑A‑ACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAggg   >  2:2216185/1‑150 (MQ=255)
                                                                                                                                                  tCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTG‑TCT‑‑TGT‑‑‑‑A‑ACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAggg   >  1:1456008/1‑150 (MQ=255)
                                                                                                                                                  tCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTG‑TCT‑‑TGT‑‑‑‑A‑ACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAggg   >  2:1563089/1‑150 (MQ=255)
                                                                                                                                                   cAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCGG‑TCT‑‑TGT‑‑‑‑A‑ACAGAAGAACAGAAACTCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGa  <  1:1753502/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                             
ATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTG‑TCT‑‑TGT‑‑‑‑A‑ACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGA  >  NC_000913/164473‑164769

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: ‑