Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 164,613 T→A 11.2% intergenic (+79/‑117) 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,6130TA11.2% 714.6 / 16.1 234intergenic (+79/‑117)hrpB/mrcBputative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base T (38/161);  new base A (29/2);  total (71/163)
Fisher's exact test for biased strand distribution p-value = 4.42e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 253 total aligned reads displayed.

TCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGC  >  NC_000913/164465‑164762
                                                                                                                                                     |                                                                                                                                                      
 cccAAACATGTCTGGCCGTACGACCCGGCAAATACTGCACCGACGCGACGGACTAAAAATTATTCGTACGGTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGTTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCctct                                                                                                                                                     <  2:4315230/150‑1 (MQ=38)
   cAAACATGTCTGCCCGCACGACCCGGCAAAAACTGCACCGACGCGACGGACGAAAACGTATCCCTACGTTTGGTTTTCTCCCCCCCCCTGTGGTTGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTcc                                                                                                                                                   <  2:3376841/150‑1 (MQ=255)
     aaCTTGTCGGGCCGGACGACCGGGAAAATAGTGCACCGACGCGACGTACTAAAATTTATTCTTATGTTTGTTTTTCTCCCTCTCCCTGTGGGTGTGGGCGGGGTTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCt                                                                                                                                                 <  2:2897542/150‑1 (MQ=255)
     aaCGTGTGTGGCGGGACGACCGGGCAAATCCGGCACGGACGCGACGGACTAAAATGTTTTCTTCTGGTTGGTTTTCTCCCTCTCGGTGTGGTTGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCt                                                                                                                                                 <  2:6555474/150‑1 (MQ=255)
      aCATGTCTGGCCGGACGACCCGGAAAATCCTGCACCGCCGCGACGGACGAAAAAGTATTCGTACGGTTGGTTTTCCCCCTCTCCCTGGGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTc                                                                                                                                                <  1:8351824/150‑1 (MQ=255)
         tGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCGCACCCAAACCCGCTCCCTCCAg                                                                                                                                             >  1:7705937/1‑150 (MQ=38)
          gTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAATTATTCGTATGGTTGGTTTTCTCCCTCTCCCTGTGGGTGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAgg                                                                                                                                            <  2:5873264/150‑1 (MQ=255)
             tGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTGGACCCGCACCCAAAACCACACCCACAAggggg                                                                                                                                         >  1:3901739/1‑148 (MQ=37)
              ggCCGGACGCCCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATACGTACGGTTGGTTTTCTCCCTCTCCCTGTGGGATAGGGCCGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGgaga                                                                                                                                        <  2:3761441/150‑1 (MQ=38)
               gCCGGACGACCCGCCAAATACTGCACCGACGCGACGGACGAAAAATTATTCGTATGGTTGTTTTTCTCCCCCTCCCTGTGGTAGTGGGCGGGGTTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGagag                                                                                                                                       <  2:7265758/150‑1 (MQ=255)
               gCCGGACCACCCCAAAAATACTGCACCGACGGGACGGACGAAAATTTTTTCGTATGGTTGTTTTCCTCCCTCTCCGTGTGGGTGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGagag                                                                                                                                       <  1:6675425/150‑1 (MQ=255)
                  ggACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGGGGCCCACCCCCAAAACCACAACCACAAGGGAGAgggg                                                                                                                                    >  1:541593/1‑150 (MQ=25)
                   gccgacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTACGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGa                                                                                                                                   <  2:2248187/148‑1 (MQ=255)
                      cacCCGCAAAATACCCCCCCGCCGCGCGGGACTAAAATGTATTCGTTCGGTGGGTTTTCCCCGTCTCCCTGTGGTTGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGAGAgggttagg                                                                                                                                <  1:5214362/149‑6 (MQ=255)
                         cgggCAAATACGCCCCCGACGCGACGGACTAAAAAGTTTTCGTACGGTTGGTTTTCTCCCTCTCCCTGGGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATc                                                                                                                             <  2:3683446/148‑1 (MQ=38)
                             cAAATACGGCACCGCCGCGCCGGCCGAAAATGTATTCGTAAGGTGGGTTTTCGCCCTCTCCCGGTGGGAGTGGGCCGGGGTGCGGGCATCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGc                                                                                                                         <  2:7636513/150‑1 (MQ=255)
                              aaaTACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCCCCACCGAAACCCCACCCGCAAGGGGGAGGGGGCCGATCGAGca                                                                                                                        >  2:5307036/1‑150 (MQ=25)
                                   cggCACCGCCGCGCCGGACGAACAAGTATTCGAACGGTTGGTTTTCGCCCTCGCCCTGTGGGAGAGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAt                                                                                                                   <  2:7719300/148‑1 (MQ=255)
                                   cTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGGGGCCCCGCACCCGAACCCACACCCACCAGGGCGGGGGGGCCCAGCGAGcccccaa                                                                                                                   >  1:4533980/1‑144 (MQ=21)
                                      cACCGACGCGACGGACGAAAATGTTTTCGTCAGGTTGGTTTTCCCCCTCTCCCTGTGGGAGCGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtttt                                                                                                                <  1:3793706/150‑1 (MQ=255)
                                       cccgccgCGCCGTCCGAAAAATTTTTTTTATGGTTGGTTTTCGCCCTCTGCGGGTGGGTGTGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTg                                                                                                               <  1:9327842/149‑1 (MQ=255)
                                        ccGACGCGCCGGACGACAACTTTTTGGTATGTTTGTTTTTCTCCCTCTCCGGGTGGTCGTGGGCGGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                                                              <  2:1864697/150‑1 (MQ=255)
                                        ccGACGCGACGTACTAAAAATTTTTCGTATGGTTGTTTTTCCCCCTCTCCGTGTGGTCGCGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                                                              <  1:3087090/150‑1 (MQ=255)
                                           ccgcgacgGACGAAAATGTTTTCTTTTGGTTGTTTTTCTCCCTCTCCCTGTGGTAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgagag                                                                                                           <  2:5174296/149‑1 (MQ=38)
                                                   gacgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGGCCCACACCCCAAACCACACACACAAGGGGGGGGGGGCCGCTCGAGCCCCAAAATTTAGAGataaaaaa                                                                                                   >  1:4680921/1‑145 (MQ=16)
                                                     cgAAAAAGTTTGCGTATGGTTGGTTTTCTCCCTCTCCGGGTGGGAGCGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttctt                                                                                                 <  1:5983761/150‑1 (MQ=255)
                                                       aaaaaGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGACCCACACCCAAAACCACACCCACAACGGGGAGGGGACCCATCCAGCACAAAATTTTATAGATAtcttaataa                                                                                               >  2:4679162/1‑145 (MQ=21)
                                                         aatgttttCGTATGTTTGGTTTTCTCCCTCTCCCTGTGGGTGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGt                                                                                             <  1:9183158/144‑1 (MQ=38)
                                                            ttATTCTTATGTTTGTTTTTCACCCTCTCCCTGTGGTTGTGGGCGGGGGTGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                                                                          <  1:903594/149‑1 (MQ=255)
                                                            gTATTCGTATGGTTGGTTTTCTCCCTCTCCGTGTGGGAGAGGGCCGGGTTGAGGGCCTCCGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                                                                          <  1:533878/150‑1 (MQ=255)
                                                                 tcggaGGTTGGGTTTGCCCCCCCCCCCGGGGGGTGGGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTaa                                                                                      <  1:8173030/146‑1 (MQ=255)
                                                                     tgtttgtttTACCCCCTGGCCGTGGGGGAGGGGGGGGGGGTGGGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaa                                                                                 <  1:2567434/147‑1 (MQ=255)
                                                                       gttggttTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGTCACCCTCACCCAAAACCACTACCACAACGGAGAGGGGACCCATCAAGCCCAAATTTTTAGAGATTAATTATTCAGTGTTGTCACaaaaga                                                                               >  2:8252827/1‑150 (MQ=18)
                                                                       gttggttTTCTCCCTCGCCCTGTGGGAGAGGGCCGGGGGGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaaga                                                                               <  2:1129164/150‑1 (MQ=255)
                                                                        ttggttTTCTCCCTCCCCCTGTGGGAGTGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAgaagaa                                                                              <  2:3901739/150‑1 (MQ=255)
                                                                        ttggttTTCTCCCTCACCCTGTGGGAGAGGGGGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAGGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAgaagaa                                                                              <  1:2979503/150‑1 (MQ=255)
                                                                         tggttTTCTCCCACTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAAc                                                                             <  1:7262536/150‑1 (MQ=255)
                                                                            gttTTCTCCCTCTCCCTGTGGGTGTGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCTCAAATTTTGAGAGTTATCTTCTTCTGTCTTGTTACAGAAGATCAg                                                                           <  1:2308437/150‑1 (MQ=255)
                                                                                tctccctctcccTGTGGGTGAGGGCGGGGTTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGaaaa                                                                       <  1:1545603/150‑1 (MQ=255)
                                                                                tctccccctcccGGTGGGTGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGaaaa                                                                       <  1:977301/150‑1 (MQ=255)
                                                                                      tctccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGGCACCCCCCCCCGAACCCACACCCACCAGGGGGGGGGGGCCCCTCGCGCCCAAATTTTGTGAGATATATTATTCGGTGTTGTAACACAAGAACAGAAAAAAGGGc                                                                 >  2:8373620/1‑150 (MQ=16)
                                                                                          cctGGGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCtttt                                                             <  2:4534737/150‑1 (MQ=255)
                                                                                                gtggggggCGGGGGTGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCt                                                       <  1:3185960/145‑1 (MQ=255)
                                                                                                ggtgggggCGGGGGGGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCt                                                       <  1:3219750/145‑1 (MQ=255)
                                                                                                ggAGTGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCt                                                       <  1:1896138/150‑1 (MQ=255)
                                                                                                ggAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCt                                                       <  1:1069103/150‑1 (MQ=255)
                                                                                                 gagaGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGAACCCCAACCAAAACCACAACCCCCAAGGGGAGGGGGCCCAAAGAGCACAAAAATTTAGAGATATAATAATATGTCGTTTAAAAAAAGAACAAAAAAAAGGGCgttttggccgg                                                      >  2:1401003/1‑143 (MQ=255)
                                                                                                       gCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCAAACCCAATCCCGCAAGGGGGGGGGGACCCATCCAGCCCAAAATTTGAGAGAAAAAATCTTATGTCTTTTAACACAAGAACAGAAAATAGGGCTTTTTCGCCAGAataat                                                >  2:7717035/1‑148 (MQ=255)
                                                                                                         cGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGc                                              <  1:29507/150‑1 (MQ=255)
                                                                                                          ggggggggggCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg                                             <  2:3292240/150‑1 (MQ=255)
                                                                                                          ggggTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg                                             <  1:8362832/150‑1 (MQ=255)
                                                                                                          ggggTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg                                             <  2:694673/150‑1 (MQ=255)
                                                                                                          ggggTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg                                             <  2:6793614/150‑1 (MQ=255)
                                                                                                          ggggTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg                                             <  2:5835836/150‑1 (MQ=255)
                                                                                                            ggTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGga                                           <  1:5239193/150‑1 (MQ=255)
                                                                                                                 gggCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaa                                      <  2:9037056/150‑1 (MQ=255)
                                                                                                                 gggCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaa                                      <  2:5025573/150‑1 (MQ=255)
                                                                                                                  ggactCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaaa                                     <  2:8948418/146‑1 (MQ=255)
                                                                                                                    cctCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGc                                   <  1:3289124/148‑1 (MQ=255)
                                                                                                                     ctCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCa                                  <  2:300753/149‑1 (MQ=255)
                                                                                                                     aTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCa                                  <  2:3037034/150‑1 (MQ=255)
                                                                                                                      tCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCAt                                 <  2:4474404/150‑1 (MQ=255)
                                                                                                                            cgcACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCgg                           <  1:4130298/150‑1 (MQ=255)
                                                                                                                               aCGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAGGCATGGCCGGGaa                        <  1:6493298/150‑1 (MQ=255)
                                                                                                                                  ttCACCCTCACCCGAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGGATCAGAAGAACAGAAACTCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGa                     <  2:4038731/150‑1 (MQ=255)
                                                                                                                                    caccctcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGAcc                   <  2:9261873/150‑1 (MQ=255)
                                                                                                                                        ctcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGa               <  2:6396299/150‑1 (MQ=255)
                                                                                                                                        ctcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGa               <  2:7448432/150‑1 (MQ=255)
                                                                                                                                        ctcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGa               <  2:166487/150‑1 (MQ=255)
                                                                                                                                        ctcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAGTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGa               >  1:7512357/1‑150 (MQ=255)
                                                                                                                                         tcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAg              >  2:2235569/1‑150 (MQ=255)
                                                                                                                                         tcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAg              >  2:616667/1‑150 (MQ=255)
                                                                                                                                           accctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGcc            >  1:610335/1‑150 (MQ=255)
                                                                                                                                           accctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGcc            >  1:5525963/1‑150 (MQ=255)
                                                                                                                                             cctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCaa          >  1:3245372/1‑150 (MQ=255)
                                                                                                                                             cctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCaa          >  1:6745143/1‑150 (MQ=255)
                                                                                                                                              ctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAAt         <  2:5094194/150‑1 (MQ=255)
                                                                                                                                              ctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAAt         <  2:4968993/150‑1 (MQ=255)
                                                                                                                                              ctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCAGGAATGACCGCGAGCCAAt         >  1:8426405/1‑150 (MQ=255)
                                                                                                                                               tAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAAtt        >  2:5423718/1‑150 (MQ=255)
                                                                                                                                                tccccTCTCCCTTGAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTg       >  2:5585989/3‑150 (MQ=255)
                                                                                                                                                aaCCCTCTCCCTCCAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTg       <  1:1823781/150‑1 (MQ=255)
                                                                                                                                                aaCCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTg       >  2:2940834/1‑150 (MQ=255)
                                                                                                                                                  cccTCTCCCTTGAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTgca     >  1:4879972/1‑148 (MQ=255)
                                                                                                                                                  cccTCTCCCTTGAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGa     >  2:355900/1‑150 (MQ=255)
                                                                                                                                                  cccTCTCCCTTGAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGAc    >  1:6563137/1‑150 (MQ=255)
                                                                                                                                                  cccTCTCCCTTGAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCAACTTATGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCCAATGGa     >  1:4880309/1‑150 (MQ=255)
                                                                                                                                                  cccTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGa     >  1:7960019/1‑150 (MQ=255)
                                                                                                                                                  cccTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCGTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGa     <  1:8919834/150‑1 (MQ=255)
                                                                                                                                                  cccTCTACATCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGACTATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGa     <  2:8388119/150‑1 (MQ=255)
                                                                                                                                                   ccTCTCCCTTGAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGAc    >  2:8715333/1‑150 (MQ=255)
                                                                                                                                                   ccTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGAc    >  2:6006357/1‑150 (MQ=255)
                                                                                                                                                   ccTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGAc    <  2:8130761/150‑1 (MQ=255)
                                                                                                                                                   ccTCTCCCGTGAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGAc    >  1:8267062/1‑150 (MQ=255)
                                                                                                                                                    ctctCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACg   >  1:7123693/1‑150 (MQ=255)
                                                                                                                                                    ctctCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACg   >  1:5795774/1‑150 (MQ=255)
                                                                                                                                                     tctcCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGa     >  2:8643849/1‑147 (MQ=255)
                                                                                                                                                     tctcCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGa     <  1:8643849/147‑1 (MQ=255)
                                                                                                                                                     tctcCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGc  >  1:3276552/1‑150 (MQ=255)
                                                                                                                                                     tctcCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGa     <  2:5318967/146‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
TCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGC  >  NC_000913/164465‑164762

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: