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.

TCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTT‑GGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGC  >  NC_000913/164465‑164762
                                                                                                                                                     |                                                                                                                                                      
 cccAAACATGTCTGGCCGTACGACCCGGCAAATACTGCACCGACGCGACGGACTAAAAATTATTCGTACGGTT‑GGTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGTTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCctc‑t                                                                                                                                                     <  2:4315230/150‑1 (MQ=38)
   cAAACATGTCTGCCCGCACGACCCGGCAAAAACTGCACCGACGCGACGGACGAAAACGTATCCCTACGTTT‑GGTTTTCTCCCCCCCCCTGTGGTTGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑Tcc                                                                                                                                                   <  2:3376841/150‑1 (MQ=255)
     aaCTTGTCGGGCCGGACGACCGGGAAAATAGTGCACCGACGCGACGTACTAAAATTTATTCTTATGTTT‑GTTTTTCTCCCTCTCCCTGTGGGTGTGGGCGGGGTTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCt                                                                                                                                                 <  2:2897542/150‑1 (MQ=255)
     aaCGTGTGTGGCGGGACGACCGGGCAAATCCGGCACGGACGCGACGGACTAAAATGTTTTCTTCTGGTT‑GGTTTTCTCCCTCTCGGTGTGGTTGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCt                                                                                                                                                 <  2:6555474/150‑1 (MQ=255)
      aCATGTCTGGCCGGACGACCCGGAAAATCCTGCACCGCCGCGACGGACGAAAAAGTATTCGTACGGTT‑GGTTTTCCCCCTCTCCCTGGGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTc                                                                                                                                                <  1:8351824/150‑1 (MQ=255)
         tGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTT‑GGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCGCACCCAAACCCGC‑TCCCTCCAg                                                                                                                                             >  1:7705937/1‑150 (MQ=38)
          gTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAATTATTCGTATGGTT‑GGTTTTCTCCCTCTCCCTGTGGGTGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAgg                                                                                                                                            <  2:5873264/150‑1 (MQ=255)
             tGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTT‑GGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTGGACCCGCACCCAAAACCAC‑ACCCACAAggggg                                                                                                                                         >  1:3901739/1‑148 (MQ=37)
              ggCCGGACGCCCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATACGTACGGTT‑GGTTTTCTCCCTCTCCCTGTGGGATAGGGCCGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTTGAGGgaga                                                                                                                                        <  2:3761441/150‑1 (MQ=38)
               gCCGGACGACCCGCCAAATACTGCACCGACGCGACGGACGAAAAATTATTCGTATGGTT‑GTTTTTCTCCCCCTCCCTGTGGTAGTGGGCGGGGTTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGagag                                                                                                                                       <  2:7265758/150‑1 (MQ=255)
               gCCGGACCACCCCAAAAATACTGCACCGACGGGACGGACGAAAATTTTTTCGTATGGTT‑GTTTTCCTCCCTCTCCGTGTGGGTGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGagag                                                                                                                                       <  1:6675425/150‑1 (MQ=255)
                  ggACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTT‑GGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGGGGCCCACCCCCAAAACCAC‑AACCACAAGGGAGAgggg                                                                                                                                    >  1:541593/1‑150 (MQ=25)
                   gccgacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTACGGTT‑GGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGa                                                                                                                                   <  2:2248187/148‑1 (MQ=255)
                      cacCCGCAAAATACCCCCCCGCCGCGCGGGACTAAAATGTATTCGTTCGGTG‑GGTTTTCCCCGTCTCCCTGTGGTTGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTTGAGGGAGAgggttagg                                                                                                                                <  1:5214362/149‑6 (MQ=255)
                         cgggCAAATACGCCCCCGACGCGACGGACTAAAAAGTTTTCGTACGGTT‑GGTTTTCTCCCTCTCCCTGGGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATc                                                                                                                             <  2:3683446/148‑1 (MQ=38)
                             cAAATACGGCACCGCCGCGCCGGCCGAAAATGTATTCGTAAGGTG‑GGTTTTCGCCCTCTCCCGGTGGGAGTGGGCCGGGGTGCGGGCATCAGCGCGCACGGTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGc                                                                                                                         <  2:7636513/150‑1 (MQ=255)
                              aaaTACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTT‑GGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCCCCACCGAAACCCC‑ACCCGCAAGGGGGAGGGGGCCGATCGAGca                                                                                                                        >  2:5307036/1‑150 (MQ=25)
                                   cggCACCGCCGCGCCGGACGAACAAGTATTCGAACGGTT‑GGTTTTCGCCCTCGCCCTGTGGGAGAGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAt                                                                                                                   <  2:7719300/148‑1 (MQ=255)
                                   cTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTT‑GGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGGGGCCCCGCACCCGAACCCAC‑ACCCACCAGGGCGGGGGGGCCCAGCGAGcccccaa                                                                                                                   >  1:4533980/1‑144 (MQ=21)
                                      cACCGACGCGACGGACGAAAATGTTTTCGTCAGGTT‑GGTTTTCCCCCTCTCCCTGTGGGAGCGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtttt                                                                                                                <  1:3793706/150‑1 (MQ=255)
                                       cccgccgCGCCGTCCGAAAAATTTTTTTTATGGTT‑GGTTTTCGCCCTCTGCGGGTGGGTGTGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTg                                                                                                               <  1:9327842/149‑1 (MQ=255)
                                        ccGACGCGCCGGACGACAACTTTTTGGTATGTTT‑GTTTTTCTCCCTCTCCGGGTGGTCGTGGGCGGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                                                              <  2:1864697/150‑1 (MQ=255)
                                        ccGACGCGACGTACTAAAAATTTTTCGTATGGTT‑GTTTTTCCCCCTCTCCGTGTGGTCGCGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                                                              <  1:3087090/150‑1 (MQ=255)
                                           ccgcgacgGACGAAAATGTTTTCTTTTGGTT‑GTTTTTCTCCCTCTCCCTGTGGTAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgagag                                                                                                           <  2:5174296/149‑1 (MQ=38)
                                                   gacgAAAAAGTATTCGTAAGGTT‑GGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGGCCCACACCCCAAACCAC‑ACACACAAGGGGGGGGGGGCCGCTCGAGCCCCAAAATTTAGAGataaaaaa                                                                                                   >  1:4680921/1‑145 (MQ=16)
                                                     cgAAAAAGTTTGCGTATGGTT‑GGTTTTCTCCCTCTCCGGGTGGGAGCGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttctt                                                                                                 <  1:5983761/150‑1 (MQ=255)
                                                       aaaaaGTATTCGTAAGGTT‑GGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGACCCACACCCAAAACCAC‑ACCCACAACGGGGAGGGGACCCATCCAGCACAAAATTTTATAGATAtcttaataa                                                                                               >  2:4679162/1‑145 (MQ=21)
                                                         aatgttttCGTATGTTT‑GGTTTTCTCCCTCTCCCTGTGGGTGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGt                                                                                             <  1:9183158/144‑1 (MQ=38)
                                                            ttATTCTTATGTTT‑GTTTTTCACCCTCTCCCTGTGGTTGTGGGCGGGGGTGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                                                                          <  1:903594/149‑1 (MQ=255)
                                                            gTATTCGTATGGTT‑GGTTTTCTCCCTCTCCGTGTGGGAGAGGGCCGGGTTGAGGGCCTCCGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                                                                          <  1:533878/150‑1 (MQ=255)
                                                                 tcggaGGTTGGGTTTGCCCCCCCCCCCGGGGGGTGGGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTaa                                                                                      <  1:8173030/146‑1 (MQ=255)
                                                                     tgttt‑gtttTACCCCCTGGCCGTGGGGGAGGGGGGGGGGGTGGGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaa                                                                                 <  1:2567434/147‑1 (MQ=255)
                                                                       gtt‑ggttTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGTCACCCTCACCCAAAACCAC‑TACCACAACGGAGAGGGGACCCATCAAGCCCAAATTTTTAGAGATTAATTATTCAGTGTTGTCACaaaaga                                                                               >  2:8252827/1‑150 (MQ=18)
                                                                       gtt‑ggttTTCTCCCTCGCCCTGTGGGAGAGGGCCGGGGGGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaaga                                                                               <  2:1129164/150‑1 (MQ=255)
                                                                        tt‑ggttTTCTCCCTCCCCCTGTGGGAGTGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAgaagaa                                                                              <  2:3901739/150‑1 (MQ=255)
                                                                        tt‑ggttTTCTCCCTCACCCTGTGGGAGAGGGGGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAGGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAgaagaa                                                                              <  1:2979503/150‑1 (MQ=255)
                                                                         t‑ggttTTCTCCCACTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAAc                                                                             <  1:7262536/150‑1 (MQ=255)
                                                                            gttTTCTCCCTCTCCCTGTGGGTGTGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCTCAAATTTTGAGAGTTATCTTCTTCTGTCTTGTTACAGAAGATCAg                                                                           <  1:2308437/150‑1 (MQ=255)
                                                                                tctccctctcccTGTGGGTGAGGGCGGGGTTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGaaaa                                                                       <  1:1545603/150‑1 (MQ=255)
                                                                                tctccccctcccGGTGGGTGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGaaaa                                                                       <  1:977301/150‑1 (MQ=255)
                                                                                      tctccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGGCACCCCCCCCCGAACCCAC‑ACCCACCAGGGGGGGGGGGCCCCTCGCGCCCAAATTTTGTGAGATATATTATTCGGTGTTGTAACACAAGAACAGAAAAAAGGGc                                                                 >  2:8373620/1‑150 (MQ=16)
                                                                                          cctGGGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCtttt                                                             <  2:4534737/150‑1 (MQ=255)
                                                                                                gtggggggCGGGGGTGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCt                                                       <  1:3185960/145‑1 (MQ=255)
                                                                                                ggtgggggCGGGGGGGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCt                                                       <  1:3219750/145‑1 (MQ=255)
                                                                                                ggAGTGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCt                                                       <  1:1896138/150‑1 (MQ=255)
                                                                                                ggAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCt                                                       <  1:1069103/150‑1 (MQ=255)
                                                                                                 gagaGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGAACCCCAACCAAAACCACAACCCCCAAGGG‑GAGGGGGCCCAAAGAGCACAAAAATTTAGAGATATAATAATATGTCGTTTAAAAAAAGAACAAAAAAAAGGGCgttttggccgg                                                      >  2:1401003/1‑143 (MQ=255)
                                                                                                       gCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCAAACCCAA‑TCCCGCAAGGGGGGGGGGACCCATCCAGCCCAAAATTTGAGAGAAAAAATCTTATGTCTTTTAACACAAGAACAGAAAATAGGGCTTTTTCGCCAGAataat                                                >  2:7717035/1‑148 (MQ=255)
                                                                                                         cGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGc                                              <  1:29507/150‑1 (MQ=255)
                                                                                                          ggggggggggCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg                                             <  2:3292240/150‑1 (MQ=255)
                                                                                                          ggggTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg                                             <  1:8362832/150‑1 (MQ=255)
                                                                                                          ggggTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg                                             <  2:694673/150‑1 (MQ=255)
                                                                                                          ggggTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg                                             <  2:6793614/150‑1 (MQ=255)
                                                                                                          ggggTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg                                             <  2:5835836/150‑1 (MQ=255)
                                                                                                            ggTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGga                                           <  1:5239193/150‑1 (MQ=255)
                                                                                                                 gggCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaa                                      <  2:9037056/150‑1 (MQ=255)
                                                                                                                 gggCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaa                                      <  2:5025573/150‑1 (MQ=255)
                                                                                                                  ggactCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaaa                                     <  2:8948418/146‑1 (MQ=255)
                                                                                                                    cctCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGc                                   <  1:3289124/148‑1 (MQ=255)
                                                                                                                     ctCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCa                                  <  2:300753/149‑1 (MQ=255)
                                                                                                                     aTCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCa                                  <  2:3037034/150‑1 (MQ=255)
                                                                                                                      tCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCAt                                 <  2:4474404/150‑1 (MQ=255)
                                                                                                                            cgcACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCgg                           <  1:4130298/150‑1 (MQ=255)
                                                                                                                               aCGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAGGCATGGCCGGGaa                        <  1:6493298/150‑1 (MQ=255)
                                                                                                                                  ttCACCCTCACCCGAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGGATCAGAAGAACAGAAACTCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGa                     <  2:4038731/150‑1 (MQ=255)
                                                                                                                                    caccctcaccctAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGAcc                   <  2:9261873/150‑1 (MQ=255)
                                                                                                                                        ctcaccctAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGa               <  2:6396299/150‑1 (MQ=255)
                                                                                                                                        ctcaccctAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGa               <  2:7448432/150‑1 (MQ=255)
                                                                                                                                        ctcaccctAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGa               <  2:166487/150‑1 (MQ=255)
                                                                                                                                        ctcaccctAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAGTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGa               >  1:7512357/1‑150 (MQ=255)
                                                                                                                                         tcaccctAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAg              >  2:2235569/1‑150 (MQ=255)
                                                                                                                                         tcaccctAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAg              >  2:616667/1‑150 (MQ=255)
                                                                                                                                           accctAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGcc            >  1:610335/1‑150 (MQ=255)
                                                                                                                                           accctAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGcc            >  1:5525963/1‑150 (MQ=255)
                                                                                                                                             cctAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCaa          >  1:3245372/1‑150 (MQ=255)
                                                                                                                                             cctAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCaa          >  1:6745143/1‑150 (MQ=255)
                                                                                                                                              ctAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAAt         <  2:5094194/150‑1 (MQ=255)
                                                                                                                                              ctAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAAt         <  2:4968993/150‑1 (MQ=255)
                                                                                                                                              ctAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCAGGAATGACCGCGAGCCAAt         >  1:8426405/1‑150 (MQ=255)
                                                                                                                                               tAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAAtt        >  2:5423718/1‑150 (MQ=255)
                                                                                                                                                tccccTC‑TCCCTTGAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTg       >  2:5585989/3‑150 (MQ=255)
                                                                                                                                                aaCCCTC‑TCCCTCCAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTg       <  1:1823781/150‑1 (MQ=255)
                                                                                                                                                aaCCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTg       >  2:2940834/1‑150 (MQ=255)
                                                                                                                                                  cccTC‑TCCCTTGAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTgca     >  1:4879972/1‑148 (MQ=255)
                                                                                                                                                  cccTC‑TCCCTTGAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGa     >  2:355900/1‑150 (MQ=255)
                                                                                                                                                  cccTC‑TCCCTTGAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAG‑AAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGAc    >  1:6563137/1‑150 (MQ=255)
                                                                                                                                                  cccTC‑TCCCTTGAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCAACTTATGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCCAATGGa     >  1:4880309/1‑150 (MQ=255)
                                                                                                                                                  cccTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGa     >  1:7960019/1‑150 (MQ=255)
                                                                                                                                                  cccTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCGTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGa     <  1:8919834/150‑1 (MQ=255)
                                                                                                                                                  cccTC‑TACATCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGACTATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGa     <  2:8388119/150‑1 (MQ=255)
                                                                                                                                                   ccTC‑TCCCTTGAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGAc    >  2:8715333/1‑150 (MQ=255)
                                                                                                                                                   ccTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGAc    >  2:6006357/1‑150 (MQ=255)
                                                                                                                                                   ccTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGAc    <  2:8130761/150‑1 (MQ=255)
                                                                                                                                                   ccTC‑TCCCGTGAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGAc    >  1:8267062/1‑150 (MQ=255)
                                                                                                                                                    ctc‑tCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACg   >  1:7123693/1‑150 (MQ=255)
                                                                                                                                                    ctc‑tCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACg   >  1:5795774/1‑150 (MQ=255)
                                                                                                                                                     tc‑tcCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGa     >  2:8643849/1‑147 (MQ=255)
                                                                                                                                                     tc‑tcCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGa     <  1:8643849/147‑1 (MQ=255)
                                                                                                                                                     tc‑tcCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGc  >  1:3276552/1‑150 (MQ=255)
                                                                                                                                                     tc‑tcCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAG‑AAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGa     <  2:5318967/146‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
TCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTT‑GGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGC  >  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: ‑