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

Only 100 of 109 total aligned reads displayed.

CATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAG‑AAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGG  >  NC_000913/164472‑164767
                                                                                                                                                    |                                                                                                                                                    
cATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTc                                                                                                                                                      <  8:338200/149‑1 (MQ=255)
  tGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCaa                                                                                                                                                    <  2:195930/149‑1 (MQ=255)
  tGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCaa                                                                                                                                                    <  4:89045/149‑1 (MQ=255)
  tGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTGCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCaa                                                                                                                                                    <  1:374404/149‑1 (MQ=255)
    tCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                                                                                   <  4:238538/148‑1 (MQ=255)
    tCTGGACTGACGACCCGGCAAATACTGCACCGACGCGACGTACTTAAAAGTATTCTTTAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACACTCACCCTAACCCACTCCCTCCAgg                                                                                                                                                  <  6:165092/149‑1 (MQ=37)
          cGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCCCCCTAACCCTCTCCCTCAAGGGAGAgg                                                                                                                                            >  1:116872/1‑149 (MQ=255)
          cGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGAGAgg                                                                                                                                            <  6:294384/149‑1 (MQ=255)
          cGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAgg                                                                                                                                            <  3:323606/149‑1 (MQ=255)
             acgacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGGGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGa                                                                                                                                         <  7:132135/149‑1 (MQ=255)
               gacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGAcc                                                                                                                                       >  5:73320/1‑149 (MQ=255)
               gacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTAGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTAACCCACCCCATAACCCTCTCCCGCAAGGGAGAGGGGccc                                                                                                                                       >  7:11060/1‑149 (MQ=37)
                  ccGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGAt                                                                                                                                    <  6:160193/149‑1 (MQ=255)
                  ccGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGGGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGAt                                                                                                                                    <  4:307348/149‑1 (MQ=255)
                   cGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATc                                                                                                                                   >  7:322089/1‑149 (MQ=255)
                    ggCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTTCACCCGCACCCTAACCCACCACCGCCAGGGAGAGGGGGCCgcacg                                                                                                                                  >  7:224492/1‑145 (MQ=37)
                             tGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAt                                                                                                                         >  2:155875/1‑149 (MQ=255)
                                  cgacgcgacgGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                                                                    <  6:6241/149‑1 (MQ=255)
                                   gacgcgacgGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgag                                                                                                                   <  6:73320/149‑1 (MQ=255)
                                    acgcgacgGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgaga                                                                                                                  <  1:245464/149‑1 (MQ=255)
                                       cgacgGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGat                                                                                                               <  7:320693/149‑1 (MQ=255)
                                        gacggacgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAgg                                                                                                                                                  >  4:245035/1‑113 (MQ=255)
                                        gacggacgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAgg                                                                                                                                                  <  3:245035/113‑1 (MQ=255)
                                        gacggacgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGat                                                                                                               >  5:133625/1‑148 (MQ=255)
                                         acggacgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGGGAGata                                                                                                              >  4:117430/1‑148 (MQ=255)
                                         acggacgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGata                                                                                                              >  2:267977/1‑148 (MQ=255)
                                           ggacgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtct                                                                                                           <  8:322089/149‑1 (MQ=255)
                                               gAATAAGTTTTCGTATGGTTGGTTTTCTCCCTCTCGCTGTGGGTGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttctt                                                                                                       <  4:171119/149‑1 (MQ=38)
                                               gAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGTGCCGGGGTGAGGGCATCAGAGCGCACGTTCCCCCTCCCCCTAACCCTCTCCCGGAAGGGAGGGGGGAACGCTCGCCCACAAATTTTGAGAGATAtcttctt                                                                                                       >  7:311030/1‑149 (MQ=37)
                                                aaaaaGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttcttc                                                                                                      >  7:142423/1‑149 (MQ=255)
                                                      tATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                                                                                <  5:281327/149‑1 (MQ=255)
                                                        ttCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTg                                                                                               <  1:97093/148‑1 (MQ=255)
                                                           gTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCACCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                                                                    >  5:284862/1‑124 (MQ=255)
                                                           gTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCACCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                                                                    <  6:284862/124‑1 (MQ=255)
                                                               ggcgggttTTCTCCCTCCCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaa                                                                                       <  8:151720/145‑1 (MQ=255)
                                                                gttggttTTCGCCCGCTCCATGGGGTTGTGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaa                                                                                       <  8:224492/148‑1 (MQ=255)
                                                                   ggttTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAgaagaa                                                                                    <  1:267977/148‑1 (MQ=255)
                                                                   ggttTTCTCCCGCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAgaagaa                                                                                    <  5:259648/148‑1 (MQ=255)
                                                                          tccctctccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAG‑AAAAt                                                                           >  2:21882/1‑149 (MQ=255)
                                                                           ccctctccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCCCCCTCACCCTAACCCTCTCCCGCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAG‑AAAATc                                                                          >  5:410720/1‑149 (MQ=255)
                                                                           ccctctccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGCGCGGGGCCCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAG‑AAAATc                                                                          >  5:134811/1‑149 (MQ=255)
                                                                           ccctctccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                                                                    <  7:100306/108‑1 (MQ=255)
                                                                            cctctccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAG‑AAAATCg                                                                         >  4:66395/1‑149 (MQ=255)
                                                                             ctctccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCCCCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGCTCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAG‑ATAATCgg                                                                        >  1:162873/1‑149 (MQ=255)
                                                                             ctctccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAG‑AAAATCgg                                                                        >  1:312642/1‑149 (MQ=255)
                                                                             ctctccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAGCTTCTTCTGTCTTGTAACAGAAGAACAG‑AAAATCgg                                                                        <  2:116872/149‑1 (MQ=255)
                                                                              tctccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAG‑AAAATCggg                                                                       >  2:82747/1‑149 (MQ=255)
                                                                                  cctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGTGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAG‑AAAATCGGGCttt                                                                   >  7:116350/1‑149 (MQ=255)
                                                                                  cctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAG‑AAAATCGGGCttt                                                                   <  6:133625/149‑1 (MQ=255)
                                                                                        ggAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAG‑AAAATCGGGCtt                                                                    <  7:17539/142‑1 (MQ=255)
                                                                                        ggAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCAACCGAACCCTATGCCGCAAGGGCGAGGGGGCCCAAACAGCACAAAATTTGAGAGCTATCTTCATTTGTGTTGGAAAAGAATAACAC‑AAAATAGGGGtttatatag                                                             >  8:17539/1‑143 (MQ=255)
                                                                                        ggAGAGGGCCGGGGTGAGGGCATAAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGCCCGATCGGGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAAATCGGGCTTTTccgc                                                              >  1:27534/1‑149 (MQ=255)
                                                                                            agGGCCGGGGTGAGGGCATCAGCGCGCACGTTCCCCCGCACCCAAACCCACTCCCGCAAGGGAGAGGGGGCCGATCGCGCCCAAAATTTGAGAGATATATTATTCTGTCTTGTAACAGAAGAACAG‑AAAATAGGGGGTTTTCGCCGGa                                                          >  6:257763/1‑148 (MQ=255)
                                                                                             gggCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGCGCACACATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAG‑AAAATCGGGCTTTTGCGCCTGaa                                                         >  5:91620/1‑148 (MQ=255)
                                                                                              ggCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAG‑AAAATCGGTCTTTTGCGCCTGAata                                                       >  2:347851/1‑149 (MQ=255)
                                                                                               gCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGCCCGCTCGCGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAG‑AAAATCGGGCTTTTTCGCCTGAatat                                                      >  5:15441/1‑149 (MQ=255)
                                                                                                     gTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGGGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAG‑AAAATCGGGCTTTTGCGCCTGAATATTGCGga                                                <  3:117430/149‑1 (MQ=255)
                                                                                                     gTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAG‑AAAATCGGGCTTTTGCGCCTGAATATTGCGga                                                <  5:257763/149‑1 (MQ=255)
                                                                                                           gCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAG‑AAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaaa                                          <  1:155875/149‑1 (MQ=255)
                                                                                                             aTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAG‑AAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAg                                         >  1:76992/1‑148 (MQ=255)
                                                                                                              tCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAG‑AAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCa                                       >  2:217385/1‑149 (MQ=255)
                                                                                                              tCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAG‑AAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCa                                       <  1:71325/149‑1 (MQ=255)
                                                                                                              tCAGCGCGCACGTTCACCCTCACCCTAACCCTATCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAG‑AAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGc                                        >  4:385442/1‑148 (MQ=255)
                                                                                                              tCAGCGCGCACGTTCACCCTCACCCTAACCCGCTCCCTCAAGGGAGAGGGGCCCGATCGAGCACAAATTTTGAGAGCTATCTTCTTCTGTCTTGTAACAGAAGAACAG‑AAAATCGGGCTTTTGCGCCGGAATATTGCGGAGAAAAAGCa                                       >  8:326581/1‑149 (MQ=255)
                                                                                                              tCAGCGCGCACGTTCACCCGCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAACTTCTTCGGTCTTGTAAAAGAAGAAAAG‑AAAATACGGCTATTGCGCCTGAATATTGCGGAGAAAAAGCa                                       >  2:71325/1‑149 (MQ=255)
                                                                                                                aGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGTAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAG‑AAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCAt                                      <  1:303047/148‑1 (MQ=255)
                                                                                                                               cctcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAG‑AAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACcgc                      <  8:142423/149‑1 (MQ=255)
                                                                                                                                  caccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAG‑AAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAg                   >  5:292722/1‑149 (MQ=255)
                                                                                                                                   accctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAG‑AAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGa                    <  1:21882/147‑1 (MQ=255)
                                                                                                                                     cctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAG‑AAAATCGGGCTTTTGCGCCTGAatat                                                      >  1:274714/1‑111 (MQ=255)
                                                                                                                                     cctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAG‑AAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCAt                                      >  2:314692/1‑127 (MQ=255)
                                                                                                                                     cctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAG‑AAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGcc                 >  2:142772/1‑148 (MQ=255)
                                                                                                                                       gtccccTCTCCCTTGAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAG‑AAAATCGGGCTTTTGCGCCTGAATATTGCGga                                                <  8:121288/112‑1 (MQ=255)
                                                                                                                                       gtccccTCTCCCTTGAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAG‑AAAATCGGGCTTTTGCGCCTGAATATTGCGga                                                >  7:121288/4‑115 (MQ=255)
                                                                                                                                       gtccccTCTCCCTTGAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAG‑AAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaaa                                          <  5:270322/118‑1 (MQ=255)
                                                                                                                                       gtccccTCTCCCTTGAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAG‑AAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaaa                                          >  6:270322/4‑121 (MQ=255)
                                                                                                                                       gtccccTCTCACTTGAGGGAGAGGGGACCGATAGAGCACAAATTTTGAGAGATATCTTATTCTGTCTTGTAACAGAAGAACAG‑AAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAAt              >  6:50293/4‑149 (MQ=255)
                                                                                                                                         aCCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAG‑AAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAAtt             <  6:15441/148‑1 (MQ=255)
                                                                                                                                         aCCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAG‑AAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTg            <  6:292722/149‑1 (MQ=255)
                                                                                                                                         aCCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAG‑AAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTg            <  6:134811/149‑1 (MQ=255)
                                                                                                                                          cccTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAG‑AAAATCGGGCTTTTGCGCCTGAata                                                       <  6:309037/105‑1 (MQ=255)
                                                                                                                                          cccTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAG‑AAAATCGGGCTTTTGCGCCTGAata                                                       >  5:309037/1‑105 (MQ=255)
                                                                                                                                          cccTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAG‑AAAATCGGGCTTTTGCGCCTGAATATTGc                                                   >  3:139244/1‑109 (MQ=255)
                                                                                                                                          cccTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAG‑AAAATCGGGCTTTTGCGCCTGAATATTGc                                                   <  4:139244/109‑1 (MQ=255)
                                                                                                                                          cccTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAG‑AAAATCGGGCTTTTGCGCCTGAATATTGCGTAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTg            >  6:66300/1‑148 (MQ=255)
                                                                                                                                           ccTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAG‑AAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGa          >  4:116114/1‑149 (MQ=255)
                                                                                                                                            ctctCCCTCAAGGGAGAGTGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAG‑AAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGAc         <  8:116350/149‑1 (MQ=255)
                                                                                                                                             tctcCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAG‑AAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACg        >  3:178917/1‑149 (MQ=255)
                                                                                                                                               tcCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAG‑AAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaa                                           <  8:419316/112‑1 (MQ=255)
                                                                                                                                               tcCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAG‑AAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaa                                           >  7:419316/1‑112 (MQ=255)
                                                                                                                                               tcCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAG‑AAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGCCCGCGAGCCAATTGGACGc       >  6:361775/1‑148 (MQ=255)
                                                                                                                                               tcCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAG‑AAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACg        <  2:82795/147‑1 (MQ=255)
                                                                                                                                                 ccTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAG‑AAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCa      >  4:150440/1‑147 (MQ=255)
                                                                                                                                                 ccTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAG‑AAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCa      >  6:315929/1‑147 (MQ=255)
                                                                                                                                                 ccTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAG‑AAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCa      >  1:229252/1‑147 (MQ=255)
                                                                                                                                                  cTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAG‑AAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCAt                                      >  5:303790/1‑114 (MQ=255)
                                                                                                                                                  cTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAG‑AAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCAt                                      <  6:303790/114‑1 (MQ=255)
                                                                                                                                                   tCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAG‑AAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAgg  >  6:159524/1‑149 (MQ=255)
                                                                                                                                                   tCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAG‑AAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAgg  >  2:237228/1‑149 (MQ=255)
                                                                                                                                                    cAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGTTATCTTCTTCTGTCTTGTAACAGTAGAACAGAAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAgg  <  2:27534/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAG‑AAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGG  >  NC_000913/164472‑164767

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑