Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 164,621 A→G 5.6% intergenic (+87/‑109) 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,6210AG 5.6% 438.8 / 12.5 147intergenic (+87/‑109)hrpB/mrcBputative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base A (46/85);  new base G (3/5);  total (56/91)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.79e-01

Only 100 of 149 total aligned reads displayed.

ATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCT‑‑TGT‑‑‑‑AACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGA  >  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)
                                                  attttttCGTATGTTTGTTTTTCCCCCTCGCCCGGTGGTTGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTc                                                                                                          <  2:811252/145‑1 (MQ=255)
                                                  aaGTATTCGTACGGTTGTTTTTCTCCCTCTCCCTGTGGGTGAGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTc                                                                                                          <  2:953701/150‑1 (MQ=255)
                                                         gtttcggttggttTACCCCCTCCCCGTGTGGGTGTGGGCGGGGGTGTGGGCCTCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCT‑‑TGT‑‑‑‑AAc                                                                                            <  1:2021792/145‑1 (MQ=255)
                                                          gTAAGGTTGTTTTTCTCCCTCTCCCTGTGGGAGAGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCT‑‑TGT‑‑‑‑AACa                                                                                           <  2:3012568/150‑1 (MQ=255)
                                                           tAAGGTTGGTTTTCTCCCTCTCCCTGTGGGTGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCT‑‑TGT‑‑‑‑AACag                                                                                          <  1:2385231/150‑1 (MQ=255)
                                                            atggttggttTGCTCCCTTTCCGTGTGGGTGTGGGCCGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCT‑‑TGT‑‑‑‑AACaga                                                                                         <  1:864319/148‑1 (MQ=255)
                                                            aaGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCT‑‑TGT‑‑‑‑AACaga                                                                                         <  2:1366492/150‑1 (MQ=255)
                                                            aaGGTTGGTTTTCTCCCTCCCCCTGGGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCT‑‑TGT‑‑‑‑AACaga                                                                                         <  1:2976384/150‑1 (MQ=255)
                                                             aGGGTGGTTTTCTCCCTCTCCCTGGGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCT‑‑TGT‑‑‑‑AACagaa                                                                                        <  2:2409315/150‑1 (MQ=255)
                                                                tttgttTACTCCCCCTCGATGGGGGAGAGGGCGGGGTTGATGGCTTCAGCGCGCACGTTCACCCTCACCCTATCCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCTCAAATTTTGAGAGATATCTTATTCTGTCT‑‑TGT‑‑‑‑AACAgaagaa                                                                                     <  1:1352225/147‑1 (MQ=255)
                                                                 tggttTTCTCCCGCGCCCTGTGGGAGAGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCT‑‑TGT‑‑‑‑AACAGAAGAAc                                                                                    <  1:4119438/150‑1 (MQ=255)
                                                                  ggttTTCTCCCGCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCGAACCCTCTCCCTCAAGGGAGAGGGGACCGAACGAGCACAAATTTTGAGAGATATCTTCTTCTGTCT‑‑AGT‑‑‑‑AACAGAAGAACa                                                                                   <  2:297121/150‑1 (MQ=255)
                                                                     tttCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGTCACCCGCCCCCGAACCCACGCCCACAATGGGGAGGGGCCCGATCCAGCCCAAATTTTGAGAGATATCTTGTTCTGTCT‑‑TGT‑‑‑‑AACAGAAGAACAGaa                                                                                >  2:2059319/1‑150 (MQ=25)
                                                                           cctctccctGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGGTCACCCCACCCCAAACCCAAACCCCCACCGGCGAGGGGACCCGTCGGGCCCAAAATTTGAGAGATATATTGTTCTGTAG‑‑AGT‑‑‑‑AACAGAAGGACCGAAAAGCgg                                                                          >  2:307917/1‑150 (MQ=12)
                                                                           cctctccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCAAAACCTCACCCTCAAGGGGGTGGGGGCCGATCGAGCACAAATTTTGAGGGATATATTCTTCTGTCT‑‑TGT‑‑‑‑GAAAGAAGTAAAGAAAATCgg                                                                          >  2:3094714/1‑150 (MQ=37)
                                                                                ccctGGGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCT‑‑TGT‑‑‑‑AACAGAAGAACAGAAAATCGGGCttt                                                                     <  1:2620588/150‑1 (MQ=255)
                                                                                 cctGGGGGAGTGGGCCGGGGTGAGGGCCTCAGCGCGCCCGTGCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCT‑‑TGT‑‑‑‑AACAGAAGAACAGAAAATCGGGCtttt                                                                    <  1:1562785/150‑1 (MQ=255)
                                                                                    gtgGGTGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCT‑‑TGT‑‑‑‑AACAGAAGAACAGAAAATCGGGCTTTTgcg                                                                 <  1:1557801/150‑1 (MQ=255)
                                                                                     tgGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCACACCCAAAACCACAACCACCAGGGGGAGGGGGCCCATCGAGCCCAAAATTTGAGAGATAAAATATTATGTAT‑‑TGT‑‑‑‑AACAGAAAAACATCAAAAAGGGCGTTTtcgc                                                                >  2:1376930/1‑150 (MQ=255)
                                                                                     ggggtgcgggCGGGGGTGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCT‑‑TGT‑‑‑‑AACAGAAGAACAGAAAATCGGGCTTTTgcgc                                                                <  1:1437939/143‑1 (MQ=255)
                                                                                     ggggggtgggCGGGGGGGGGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCT‑‑TGT‑‑‑‑AACAGAAGAACAGAAAATCGGGCTTTTgcgc                                                                <  2:377491/143‑1 (MQ=255)
                                                                                     ggggAGTGGGCCGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCT‑‑TGT‑‑‑‑AACAGAAGAACAGAAAATCGGGCTTTTgcgc                                                                <  1:3867645/149‑1 (MQ=255)
                                                                                     ggggAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCT‑‑TGT‑‑‑‑AACAGAAGAACAGAAAATCGGGCTTTTgcgc                                                                <  1:3884677/149‑1 (MQ=255)
                                                                                                cGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCT‑‑TGT‑‑‑‑AACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGc                                                     <  1:4108338/150‑1 (MQ=255)
                                                                                                cGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCT‑‑TGT‑‑‑‑AACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGc                                                     <  1:434050/150‑1 (MQ=255)
                                                                                                 ggggTGAAGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCT‑‑TGT‑‑‑‑AACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg                                                    <  2:610582/150‑1 (MQ=255)
                                                                                                   ggTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCT‑‑TGT‑‑‑‑AACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGga                                                  <  1:2711821/150‑1 (MQ=255)
                                                                                                        ggccATCAGCGCCCACGTTCACCCTCACCCTACCCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAATTCGAGAGATATCTTCTTCTGTCT‑‑TGT‑‑‑‑AACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaa                                             <  1:458293/147‑1 (MQ=255)
                                                                                                             tCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCT‑‑TGT‑‑‑‑AACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCAt                                        <  2:3969176/150‑1 (MQ=255)
                                                                                                              cAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCT‑‑TGT‑‑‑‑AACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATg                                       <  1:406393/150‑1 (MQ=255)
                                                                                                                 cgcgcACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCT‑‑TGT‑‑‑‑AACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGcc                                    <  2:2929586/150‑1 (MQ=255)
                                                                                                                       cGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCT‑‑TGT‑‑‑‑AACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAAt                              <  2:484181/150‑1 (MQ=255)
                                                                                                                         ttCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCT‑‑TGT‑‑‑‑AACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGa                            <  2:4072574/150‑1 (MQ=255)
                                                                                                                         ttCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCT‑‑TGT‑‑‑‑AACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGa                            <  2:2805426/150‑1 (MQ=255)
                                                                                                                         ttCACCCTCAACCTAACACAAAACCAAAAGGGAGAGGGGCCCGATCGAGCACAAATTTTGAGAGATAAATTCTTCTGTAT‑‑TGT‑‑‑‑AACAGAAGAACAGAAAATAGGGCTTTTTCGCCGGAATATTTCGGAGAAAAAGCAAGGCCCGGAATGa                            >  2:2806713/1‑150 (MQ=255)
                                                                                                                           caccctcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCT‑‑TGT‑‑‑‑AACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGAcc                          <  2:4051553/150‑1 (MQ=255)
                                                                                                                               ctcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCT‑‑TGT‑‑‑‑AACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGACTGACCGCGa                      >  1:3410238/1‑150 (MQ=255)
                                                                                                                               ctcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCT‑‑TGT‑‑‑‑AACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGa                      >  1:1357841/1‑150 (MQ=255)
                                                                                                                               ctcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCT‑‑TGT‑‑‑‑AACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGa                      >  1:2733391/1‑150 (MQ=255)
                                                                                                                               ctcacccgACCCCTCTCCCTCACGGGAGAGGGGCCCGATCGAGCACAAAGGTTGAGAGATATCTTCTTCTGTCT‑‑TGT‑‑‑‑AACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATGGCGGAGAAAAAGCATGGCCGGGAATGACCGCGa                      <  1:2005657/150‑1 (MQ=255)
                                                                                                                                tcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCT‑‑TGT‑‑‑‑AACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAg                     <  1:1655749/150‑1 (MQ=255)
                                                                                                                                tcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCT‑‑TGT‑‑‑‑AACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAg                     <  2:338824/150‑1 (MQ=255)
                                                                                                                                 caccctACCCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCT‑‑TGT‑‑‑‑AACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGc                    <  2:822990/150‑1 (MQ=255)
                                                                                                                                    cctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCT‑‑TGT‑‑‑‑AACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAACAGCATGGCCGGGAATGACCGCGAGCCaa                 >  2:2263268/1‑150 (MQ=255)
                                                                                                                                    cctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCT‑‑TGT‑‑‑‑AACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCaa                 >  2:2263306/1‑150 (MQ=255)
                                                                                                                                    cctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCT‑‑TGT‑‑‑‑AACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCaa                 >  2:481164/1‑150 (MQ=255)
                                                                                                                                      tAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGATCACAAATTTTGAGAGATATCTTCTTCTGTCT‑‑TGT‑‑‑‑AACAGAAGAACAGAAAATCGGGCTTTTGCGCATGATTATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAAtt               <  2:2014534/150‑1 (MQ=255)
                                                                                                                                      tAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCT‑‑TGT‑‑‑‑AACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGTATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAAtt               <  1:3440302/150‑1 (MQ=255)
                                                                                                                                      tAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCT‑‑TGT‑‑‑‑AACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAAtt               <  2:773930/150‑1 (MQ=255)
                                                                                                                                      tAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCT‑‑TGT‑‑‑‑AACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAAtt               <  1:1260942/150‑1 (MQ=255)
                                                                                                                                      tAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCT‑‑TGT‑‑‑‑AACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAAtt               <  2:48429/150‑1 (MQ=255)
                                                                                                                                      gtccccTCTCCCTTGAGGGAGAGGGTACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCT‑‑TGT‑‑‑‑AACagaa                                                                                        <  1:11191/74‑1 (MQ=37)
                                                                                                                                      gtccccTCTCCCTTGAGGGAGAGGGTACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCT‑‑TGT‑‑‑‑AACagaa                                                                                        >  2:11191/4‑77 (MQ=37)
                                                                                                                                       aaCCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTTTTCTGTTCTTCTGTTACAAGACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGa                      <  1:2036110/150‑1 (MQ=255)
                                                                                                                                       aaCCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCT‑‑TGT‑‑‑‑AACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTg              <  1:1616628/150‑1 (MQ=255)
                                                                                                                                         cccTCTCCCTTGAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCT‑‑TGT‑‑‑‑AACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGa            >  1:853933/1‑150 (MQ=255)
                                                                                                                                         cccTCTCCCTCACGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGAGATCTTCTTCTGTCT‑‑TGT‑‑‑‑AACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGa            <  1:4036839/150‑1 (MQ=255)
                                                                                                                                         cccTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCT‑‑TGT‑‑‑‑AACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGa            <  2:1031077/150‑1 (MQ=255)
                                                                                                                                         cccTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCT‑‑TGT‑‑‑‑AACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGa            >  2:1948620/1‑150 (MQ=255)
                                                                                                                                          ccTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCT‑‑TGT‑‑‑‑AACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGAc           >  1:3230224/1‑150 (MQ=255)
                                                                                                                                           ctctCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCT‑‑TGT‑‑‑‑AACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGACTGACCGCGAGCCAATTGGACg          >  1:2246295/1‑150 (MQ=255)
                                                                                                                                           ctctCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCT‑‑TGT‑‑‑‑AACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACg          >  1:94228/1‑150 (MQ=255)
                                                                                                                                           ctctCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCT‑‑TGT‑‑‑‑AACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACg          >  1:1784146/1‑150 (MQ=255)
                                                                                                                                           ctctCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCT‑‑TGT‑‑‑‑AACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACg          >  1:2172320/1‑150 (MQ=255)
                                                                                                                                             ctcCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCT‑‑TGT‑‑‑‑AACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCa        >  1:2029186/1‑150 (MQ=255)
                                                                                                                                              tcCCTCTAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCT‑‑TGT‑‑‑‑AACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCaa       <  2:2429102/150‑1 (MQ=255)
                                                                                                                                                ccTCACGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCT‑‑TGT‑‑‑‑AACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAg     <  2:3627882/150‑1 (MQ=255)
                                                                                                                                                ccTCAAGGGAGAGGGGACCGATCGAGCACAACTTTTGAGAGATATCTTCTTCTGTCT‑‑TGT‑‑‑‑AACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAg     <  2:2125696/150‑1 (MQ=255)
                                                                                                                                                ccTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCT‑‑TGT‑‑‑‑AACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAg     >  2:2999534/1‑150 (MQ=255)
                                                                                                                                                 cTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCT‑‑TGT‑‑‑‑AACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAgg    >  2:365189/1‑150 (MQ=255)
                                                                                                                                                  tCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCT‑‑TGT‑‑‑‑AACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAggg   >  2:755779/1‑150 (MQ=255)
                                                                                                                                                  tCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCT‑‑TGT‑‑‑‑AACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAggg   >  2:2216185/1‑150 (MQ=255)
                                                                                                                                                  tCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCT‑‑TGT‑‑‑‑AACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAggg   >  1:1456008/1‑150 (MQ=255)
                                                                                                                                                  tCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCT‑‑TGT‑‑‑‑AACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAggg   >  2:1563089/1‑150 (MQ=255)
                                                                                                                                                   cAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCGGTCT‑‑TGT‑‑‑‑AACAGAAGAACAGAAACTCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGa  <  1:1753502/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                            
ATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCT‑‑TGT‑‑‑‑AACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGA  >  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: