Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,093,724 T→C 6.2% K254R (AAG→AGG)  hisP ← histidine lysine/arginine/ornithine ABC transporter ATP binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,093,7240TC 6.2% 181.9 / 2.0 64K254R (AAG→AGG) hisPhistidine lysine/arginine/ornithine ABC transporter ATP binding protein
Reads supporting (aligned to +/- strand):  ref base T (22/38);  new base C (2/2);  total (24/40)
Fisher's exact test for biased strand distribution p-value = 6.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.03e-01

GTGCTGGGCCAGATGAGAAAACCGATTCGCTGGGTGGTCATAGGGTACGGTCCGAAACCTGAGGTCGTTCGA‑‑GTCTGTGCGCCGGGTCAGGCTGCGCTCGCCGCGCAGCATGCCCTATTCTGGTGCAAAAATGCAGCTTTTACTTCAGGCTGCCGGAGAGGAACTGCTGCAGGCGCTCCGACTGCGGGTTGGCCAGCACTTCACGTGGGCAGCCGGTTTCTTCCACCAGGCCTTTGTGCAGGAACACCAACTGGTTGGACACCTCGCGGGCAAAGCCCATTTCGTGGGTCAC  >  NC_002947/5093581‑5093872
                                                                                                                                                   |                                                                                                                                                     
gTGCTGGGCCAGATGAGAAAACCGATTCGCTGGGTGGTCATAGGGTACGGTCCGAAACCTGAGGTCGTTCGA‑‑GTCTGTGCGCCGGGTCAGGCTGCGCTCGCCGCGCAGCATGCCCTATTCTGGTGCAAAAATGCAGCTTTTACTTCa                                                                                                                                                    <  7:291520/147‑1 (MQ=255)
      ggCCAGATGAGAAAACCGATTCGCTGGGTGGTCATAGGGTACGGTCCGAAACCTGAGGTCGTTCGA‑‑GTCTGTGCGCCGGGTCAGGCTGCGCTCGCCGCGCAGCATGTCCTATTCTGGTGCAAAAATGCAGCTTTTACTTCAGGCTGCCg                                                                                                                                            <  7:91013/149‑1 (MQ=255)
      ggCCAGATGAGAAAACCGATTCGCTGGGTGGTCATAGGGTACGGTCCGAAACCTGAGGTCGTTCGA‑‑GTCTGTGCGCCGGGTCAGGCTGCGCTCGCCGCGCAGCATGCCCTATTCTGGTGCAAAAATGCAGCTTTTACTTCAGGCTGc                                                                                                                                              <  1:356556/147‑1 (MQ=255)
          agttGAGAAAACCGATTCGCTGGGTGGTCATAGGGTACGGTCCGAAACCTGAGGTCGTTCGA‑‑GTCTGTGCGCCGGGTCAGGCTGCGCTCGCCGCGCAGCATGCCCTATTCTGGTGCAAAAATGCAGCTTTTACTTCAGGCTGCCGgaga                                                                                                                                        <  7:151296/146‑1 (MQ=255)
          aGATGAGAAAACCGATTCGCTGGGTGGTCATAGGGTACGGTCCGAAACCTGAGGTCGTTCGA‑‑GTCTGTGCGCCGGGTCAGGCTGCGCTCGCCGCGCAGCATGCCCTATTCTGGTGCAAAAATGCAGCTTTTACTTCAGGCTGCCGgaga                                                                                                                                        <  7:75294/149‑1 (MQ=255)
          aGATGAGAAAACCGATTCGCTGGGTGGTCATAGGGTACGGTCCGAAACCTGAGGTCGTTCGA‑‑GTCTGTGCGCCGGGTCAGGCTGCGCTCGCCGCGCAGCATGCCCTATTCTGGTGCAAAAATGCAGCTTTTACTTCAGGCTGCCGgaga                                                                                                                                        <  2:139331/149‑1 (MQ=255)
          aGATGAGAAAACCGATTCGCTGGGTGGTCATAGGGTACGGTCCGAAACCTGAGGTCGTTCGA‑‑GTCTGTGCGCCGGGTCAGGCAGCGCTCGCCGCGCAGCATGCCCTATTCTGGTGCAAAAATGCAGCTTTTACTTCAGGCTGCCGga                                                                                                                                          >  1:381757/1‑147 (MQ=255)
            aTGAGAAAACCGATTCGCTGGGTGGTCATAGGGTACGGTCCGAAACCTGAGGTCGTTCGA‑‑GTCTGTGCGCCGGGTCAGGCTGCGCTCGCCGCGCATCATGCCCTATTCTGGTGCAAAAATGCAGCTTTTACTTCAGGCTGCCGGAGAgg                                                                                                                                      <  5:357938/149‑1 (MQ=255)
             tGAGAAAACCGATTCGCTGTGTGGTCATAGTGTACTGTCCGAAACCTGAGGTCGTTACGACTTTGTGTGCGCCCCCGTCAGGCTGCGCTCGCCGCGCAGCATGCCCTATTCTGGTGCAAAAATGCAGCTTTTACTTCAGGCTGCCGgag                                                                                                                                         <  1:320222/149‑1 (MQ=255)
               cgaAAACCGATTCGCGGGGTGGGCATAGGGTACGGTCCGAAACCTGAGGTCTTTCGA‑‑GTCTGTGCGCCGGGTCAGGCTGCGCTCGCCGCGCAGCATGCCCTATTCTGGTGCAAAAATGCAGCTTTTACTTCAGGCTGCCGGAGAGGaa                                                                                                                                    <  7:187425/147‑1 (MQ=255)
               agaAAACCGATTCGCTGGGTGGTCATAGGGTACGGTCCGAAACCTGAGGTCGTTCGA‑‑GTCTGTGCGCCGGGTCAGGCTGCGCTCGCCGCGCAGCATGCCCTATTCTGGTGCAAAAATGCAGCTTTTACTTCAGGCTGCCGGAGAGGaa                                                                                                                                    <  7:271388/148‑1 (MQ=255)
                    cccGATTCGCTGGGTGGGCCTAGGGTACGGTCCGAAACCTGAGGTCGTTCGA‑‑GTCTGTGCGCCGGGTCAGGCTGCGCTCGCCGCGCAGCATGCCCTATTCTGGTGCAAAAATGCAGCTTTTACTTCAGGCTGCCGGAGAGGAActgctg                                                                                                                              <  8:57801/148‑1 (MQ=255)
                    aCCGATTCGCTGGGTGGTCATAGGGTACGGTCCGAAACCTGAGGTCGTTCGA‑‑GTCTGTGCGCCGGGTCAGGCTGCGCTCGCCGCGCAGCATGCCCTATTCTGGTGCAAAAATGCAGCTTTTACTTCAGGCTGCCGGAGAGGAActgctg                                                                                                                              <  5:196484/149‑1 (MQ=255)
                      cGATTCGCTGGGTGGTCATAGGGTACGGTCCGAAACCTGAGGTCGTTCGA‑‑GTCTGTGCGCCGGGTCAGGCTGCGCTCGCCGCGCAGCATGCCCTATTCTGGTGCAAAAATGCAGCTTTTACTTCAGGCTGCCGGAGAGGAActgctg                                                                                                                              >  3:361310/1‑147 (MQ=255)
                       gATTCGCTGGGTGGTCATAGTGTACGGTCCGAAACCTGAGGTCGTTCGA‑‑GTCTGTGCGCCGGGTCAGGCAGCGCTCGCCGCGCAGCATGCCCTATTCTGGTGCAAAAATGCAGCTTTTACTTCAGGCTGCCGGAGAGGAACTGCTGCAg                                                                                                                           <  2:381757/149‑1 (MQ=255)
                            gCTGGGTGGTCATAGGGTACGGTCCGAAACCTGAGGTCGTTCGA‑‑GTCTGTGCGCCGGGTCAGGCTGCGCTCGCCGCGCAGCATGCCCTATTCTGGTGCAAAAATGCAGCTTTTACTTCAGGCTGCCGGAGAGGAACTGCTGCAGGCGCt                                                                                                                      >  2:137724/1‑149 (MQ=255)
                                    gtCATAGGGTACGGTCCGAAACCTGAGGTCGTTCGA‑‑GTCTGTGCGCCGGGTCAGGCTGCGCTCGCCGCGCAGCATGCCCTATTCTGGTGCAAAAATGCAGCTTTTACTTCAGGCTGCCGGAGAGGAACTGCTGCAGGCGCTCCGACTGc                                                                                                              >  3:143794/1‑149 (MQ=255)
                                    gtCATAGGGTACGGTCCGAAACCTGAGGTCGTTCGA‑‑GTCTGTGCGCCGGGTCAGGCTGCGCTCGCCGCGCAGCATGCCCTATTCTGGTGCAAAAATGCAGCTTTTACTTCAGGCTGCCGGAGAGGAACTGCTGCAGGCGCTCCGACTGc                                                                                                              <  5:185542/149‑1 (MQ=255)
                                     tCATAGGGTACGGTCCGAAACCTGAGGTCGTTCGA‑‑GTCTGTGCGCCGGGTCAGGCTGCGCTCGCCGCGCAGCATGCCCTATTCTGGTGCAAAAATGCAGCTTTTACTTCAGGCTGCCGGAGAGGAACTGCTGCAGGCGCTCCGACTGCg                                                                                                             >  5:231702/1‑149 (MQ=255)
                                     tCATAGGGTACGGTCCGAAACCTGAGGTCGTTCGA‑‑GTCTGTGCGCCGGGTCAGGCTGCGCTCGCCGCGCAGCATGCCCTATTCTGGTGCAAAAATGCAGCTTTTACTTCAGGCTGCCGGAGAGGAACTGCTGCAGGCGCTCCGACTGCg                                                                                                             <  4:143794/149‑1 (MQ=255)
                                       aTAGGGTACGGTCCGAAACCTGAGGTCGTTCGA‑‑GTCTGTGCGCCGGGTCAGGCTGCGCTCGCCGCGCAGCATGCCCTATTCTGGTGCAAAAATGCAGCTTTTACTTCAGGCTGCCGGAGAGGAACtgctgc                                                                                                                             >  7:112520/1‑131 (MQ=255)
                                       aTAGGGTACGGTCCGAAACCTGAGGGCGTTCGA‑‑GTCTGTGCGCCGGGTCAGGCTGCGCTCGCCGCGCAGCATGCCCTATTCTGGTGCAAAAATGCAGCTTTTACTTCAGGCTGCCGGAGAGGAACtgctgc                                                                                                                             <  8:112520/131‑1 (MQ=255)
                                          gggTACGGTCCGAAACCTGAGGTCGTTCGA‑‑GTCTGTGCGCCGGGTCAGGCTGCGCTCGCCGCGCAGCATGCCCTATTCTGGTGCAAAAATGCAGCTTTTACTTCAGGCTGCCGGAGAGGAACTGCTGCAGGCGCTCCGACTGCGGGtt                                                                                                        <  8:324694/148‑1 (MQ=255)
                                            gTACGGTCCGAAACCTGAGGTCGTTCGA‑‑GTCTGTGCGCCGGGTCAGGCTGCGCTCGCCGCGCAGCATGCCCTATTCTGGTGCAAAAATGCAGCTTTTACTTCAGGCTGCCGGAGAGGAACTGCTGCAGGCGCTCCGACTGCGGGTTGGc                                                                                                     >  8:180594/1‑149 (MQ=255)
                                            gTACGGTCCGAAACCTGAGGTCGTTCGA‑‑GTCTGTGCGCCGGGTCAGGCTGCGCTCGCCGCGCAGCATGCCCTATTCTGGTGCAAAAATGCAGCTTTTACTTCAGGCTGCCGGAGAGGAACTGCTGCAGGCGCTCCGACTGCGGGTTGGc                                                                                                     <  4:50884/149‑1 (MQ=255)
                                             tACGGTCCGAAACCTGAGGTCGTTCGA‑‑GTCTGTGCGCCGGGTCAGGCTGCGCTCGCCGCGCAGCATGCCCTATTCTGGTGCAAAAATGCAGCTTTTACTTCAGGCTGCCGGAGAGGAACTGCTGCAGGCGCTCCGACTGCGGGTTGGcc                                                                                                    <  6:231702/149‑1 (MQ=255)
                                                 gTCCGAAACCTGAGGTCGTTCGA‑‑GTCTGTGCGCCGGGTCAGGCTGCGCTCGCCGCGCAGCATGCCCTATTCTGGTGCAAAAATGCAGCTTTTACTTCAGGCTGCCGGAGAGGAACTGCTGCAGGCGCTCCGACTGCGGGTTg                                                                                                       >  7:85508/1‑142 (MQ=255)
                                                 gTCCGAAACCTGAGGTCGTTCGA‑‑GTCTGTGCGCCGGGTCAGGCTGCGCTCGCCGCGCAGCATGCCCTATTCTGGTGCAAAAATGCAGCTTTTACTTCAGGCTGCCGGAGAGGAACTGCTGCAGGCGCTCCGACTGCGGGTTg                                                                                                       <  8:85508/142‑1 (MQ=255)
                                                 gTCCGAAACCTGAGGTCGTTCGA‑‑GTCTGTGCGCCGGGTCAGGCTGCGCTCGCCGCGCAGCATGCCCTATTCTGGTGCAAAAATGCAGCTTTTACTTCAGGCTGCCGGAGAGGAACTGCTGCAGGCGCTCCGACTGCGGGTTGGCCAGCa                                                                                                <  1:264577/149‑1 (MQ=255)
                                                  tCCGAAACCTGAGGTCGTTCGA‑‑GTCTGTGCGCCGGGTCAGGCTGCGCTCGCCGCGCAGCATGCCCTATTCTGGTGCAAAAATGCAGCTTTTACTTCAGGCTGCCGGAGAGGAACTGCTGCAGGCGCTCCGACTGCGGGTTGGCCAGCAc                                                                                               >  3:211052/1‑149 (MQ=255)
                                                    cGAAACCTGAGGTCGTTCGA‑‑GTCTGTGCGCCGTGTCAGGATGCGCTCGCCGCGCAGCATGCCCTATTCTGGTGCAAAAATGAAGCTTTTAATTCAGCATCAAGGAGAGGAAATGCTCCAGCAGCTCAGACTGCGGGTTGGCCAGCACt                                                                                              >  8:351631/1‑148 (MQ=255)
                                                     gAAACCTGAGGTCGTTCGA‑‑GTCTGTGCGCCGGGTCAGGCTGCGCTCGCCGCGCAGCATGCCCTATTCTGGTGCAAAAATGCAGCTTTTACTTCAGGCTGCCGGAGAGGAACTGCTGCAGGCGCTCCGACTGCGGGTTGGCCAGCACtt                                                                                             >  7:350970/1‑148 (MQ=255)
                                                         ccTGAGGTCGTTCGA‑‑GTCTGTGCGCCGGGTCAGGCTGCGCTCGCCGCGCAGCATGCCCTATTCTGGTGCAAAAATGCAGCTTTTACTTCAGGCTGCCGGAGAGGAACTGCTGCAGGCGCTCCGACTGCGGGTTGGCCAGCACTTCACGt                                                                                        <  1:137724/149‑1 (MQ=255)
                                                             aGGTCGTTCGA‑‑GTCTGTGCGCCGGGTCAGGCTGCGCTCGCCGCGCAGCATGCCCTATTCTGGTGCAAAAATGCAGCTTTTACTTCAGGCTGCCGGAGAGGAACTGCTGCAGGCGCTCCGACTGCGGGTTGGCCAGCACTTCACGTggg                                                                                     >  4:381229/1‑148 (MQ=255)
                                                             aGGTCGTTCGA‑‑GTCTGTGCGCCGGGTCAGGCTGCGCTCGCCGCGCAGCATGCCCTATTCTGGTGCAAAAATGCAGCTTTTACTTCAGGCTGCCGGAGAGGAACTGCTGCAGGCGCTCCGACTGCGGGTTGGCCAGCACTTCACGTgg                                                                                      <  8:39226/147‑1 (MQ=255)
                                                                tCGTTCGA‑‑GTCTGTGCGCCGGGTCAGGCTGCGCTCGCCGCGCAGCATGCCCTATTCTGGTGCAAAAATGCAGCTTTTACCTCAGGCTGCCGGAGAGGAACTGCTGCAGGCGCTCCGACTGCGGGTTGGCCAGCACTTCACGTGGGCAGc                                                                                 >  7:220723/1‑149 (MQ=255)
                                                                   ttCGA‑‑GTCTGTGCGCCGGGTCAGGCTGCGCTCGCCGCGCAGCATGCCCTATTCTGGTGCAAAAATGCAGCTTTTACTTCAGGCTGCCGGAGAGGAACTGCTGCAGGCGCTCCGACTGCGGGTTGGCCAGCACTTCACGTGGGCAGCCgg                                                                              >  8:74105/1‑149 (MQ=255)
                                                                           gTCTGTGCGCCGGGTCAGGCTGCGCTCGCCGCGCAGCATGCCCTATTCTGGTGCAAAAATGCAGCTTTTACTTCAGGCTGCCGGAGAGGAACTGCTGCAGGCGCTCCGACTGCGGGTTGGCCAGCACTTCACGTGGGCAGCCGGTttct                                                                         >  2:236571/1‑149 (MQ=255)
                                                                              tgtgCGCCGGGTCAGGCTGCGCTCGCCGCGCAGCATGCCCTATTCTGGTGCAAAAATGCAGCTTTTACTTCAGGCTGCCGGAGAGGAACTGCTGCAGGCGCTCCGACTGCGGGTTGGCCAGCACTTCACGTGGGCAGCCGGTttcttc                                                                       >  2:127132/1‑148 (MQ=255)
                                                                                         gTCAGGCTGCGCTCGCCGCGCAGCATGCCCTATTCTGGTGCAAAAATGCAGCTTTTACTTCAGGCTGCCGGAGAGGAACTGCTGCAGGCGCTCCGACTGCGGGTTGGCAAGCACTTAACGTGGGCAGCCGGTTTCTTCCACCAGCCCtt                                                            >  8:135302/1‑149 (MQ=255)
                                                                                          tCAGGCTGCGCTCGCCGCGCAGCATGCCCTATTCTGGTGCAAAAATGCAGCTTTTACTTCAGGCTGCCGGAGAGGAACTGCTGCAGGCGCTCCGACTGCGGGTTGGCCAGCACTTCACGTGGGCAGCCGGTTTCTTCCACCAGGCCttt                                                           >  1:305827/1‑149 (MQ=255)
                                                                                               cTGCGCTCGCCGCGCAGCATGCCCTATTCTGGTGCAAAAATGCAGCTTTTACTTCAGGCTGCCGGAGAGGAACTGCTGCAGGCGCTCCGACTGCGGGTTGGCCAGCACTTCACGTGGGCAGCCGGTTTCTTCCACCAGGCCTTTGTGCa                                                      <  8:127848/149‑1 (MQ=255)
                                                                                                  cgcTCGCCGCGCAGCATGCCCTATTCTGGTGCAAAAATGCAGCTTTTACTTCAGGCTGCCGGAGAGGAACTGCTGCAGGCGCTCCGACTGCGGGTTGGCCAGCACTTCACGTGGGCAGCCGGTTTCTTCCACCAGGCCTTTGTGCAg                                                     <  7:74105/147‑1 (MQ=255)
                                                                                                  cgcTCGCCGCGCAGCATGCCCTATTCTGGTGCAAAAATGCAGCTTTTACTTCAGGCTGCCGGAGAGGAACTGCTGCAGGCGCTCCGACTGCGGGTTGGCCAGCACTTCACGTGGGCAGCCGGTTTCTTCCACCAGGCCTTTGTGCAg                                                     <  6:143598/147‑1 (MQ=255)
                                                                                                    cTCGCCGCGCAGCATGCCCTATTCTTGTGCAAAAATGCAGCTTTTACTTCAGGCTGCCGGAGAGGAACTGCTGCAGGCGCTCCTACTGCGGGTTGGCCAGCACTTCACGTGGGCAGCCGGTTTCTTCCACCAGGCCTTTGTGCAGGaa                                                  <  7:351631/148‑1 (MQ=255)
                                                                                                       gccgcGCAGCATGCCCTATTCTGGTGCAAAAATGCAGCTTTTACTTCAGGCTGCCGGAGAGGAACTGCTGCAGGCGCTCCGACTGCGGGTTGGCCAGCACTTCACGTGGGCAGCCGGTTTCTTCCACCAGGCCTTTGTGCAGGAacac                                               <  7:135302/148‑1 (MQ=255)
                                                                                                          gcgcAGCATGCCCTATTCTGGTGCAAAAATGCAGCTTTTACTTCAGGCTGCCGGAGAGGAACTGCTGCAGGCGCTCCGACTGCGGGTTGGCCAGCACTTCACGTGGGCAGCCGGTTTCTTCCACCAGGCCTTTGTGCAGGAACACCAAc                                           >  2:76095/1‑149 (MQ=255)
                                                                                                          gcgcAGCATGCCCTATTCTGGTGCAAAAATGCAGCTTTTACCTCAGGCTGCCGGAGAGGAACTGCTGCAGGCGCTCCGACTGCGGGTTGGCCAGCACTTCACGTGGGCAGCCGGTTTCTTCCACCAGGCCTTTGTGCAGGAACACCAAc                                           <  8:220723/149‑1 (MQ=255)
                                                                                                             cagcaTGCCCTAGTCTGGTGCAAAAATGCAGCTTTTACTTCAGGCTGCCGGAGAGGAACTGCTGCAGGCGCTCCGACTGCGGGTTGGCCAGCACTTCACGTGGGCAGCCGGTTTCTTCCACCAGGCCTTTGTGCAGGAACACCAACTgg                                        <  2:45786/149‑1 (MQ=255)
                                                                                                               gcaTGCCCTATTCTGGTGCAAAAATGCAGCTTTTACTTCAGGCTGCCGGAGAGGAACTGCTGCAGGCGCTCCGACTGCGGGTTGGCCAGCACTTCACGTGGGCAGCCGGTTTCTTCCACCAGGCCTTTGTGCAGGAACACCAACTGGtt                                      <  4:21307/149‑1 (MQ=255)
                                                                                                                cttGCCCTATTCTGGTGCAAAAATGCAGCTTTTACTTCAGGCTGCCGGAGAGTAACTGCTGCAGGGGCTCCGACTGCGGGTTGGCCAGCACTTCACGTGGGCTGCCGGTTTCTTCCACCAGGCCTTTGTGCAGGAACACCAACTGGTTg                                     <  1:157031/147‑1 (MQ=255)
                                                                                                                cctGCCCTATTCTGGTGCAAAAATGCAGCTTTTACTTCAGGCTGCCGGAGAGGAACTGCTGCAGGCGCTCCGACTGCGGGTTGGCCAGCACTTCACGTGGGCAGCCGGTTTCTTCCACCAGGCCTTTGTGCAGGAACACCAACTGGTTg                                     <  4:180949/147‑1 (MQ=255)
                                                                                                                   gCCCTATACTGGTGCAAAAATGCAGCTTATACTTCAGGCTGCCGGAGAGGAACTGCTGCAGGCGCTCCGACTGCGGGTTGGCCAGCACTTCACGTGGGCAGCCGGTTTCTTCCACCAGGCCTTTGTGCAGGAACACCAACTGGTTGGac                                  >  2:255004/1‑149 (MQ=255)
                                                                                                                      cTATTCTGGTGCAAAAATGCAGCTTTTACTTCAGGCTGCCGGAGAGGAACTGCTGCAGGCGCTCCGACTGCGGGTTGGCCAGCACTTCACGTGGGCAGCCGGTTTCTTCCACCAGGCCTTTGTGCAGGAACACCAACTGGTTGGACAcc                               >  4:8025/1‑149 (MQ=255)
                                                                                                                        aTTCTGGTGCAAAAATGCAGCTTTTACCTCAGGCTGCCGGAGAGGAACTGCTGCAGGCGCTCCGACTGCGGGTTGGCCAGCACTTCACGTGGGCAGCCGGTTTCTTCCACCAGGCCTTTGTGCAGGAacac                                               >  4:110951/1‑131 (MQ=255)
                                                                                                                        aTTCTGGTGCAAAAATGCAGCTTTTACCTCAGGCTGCCGGAGAGGAACTGCTGCAGGCGCTCCGACTGCGGGTTGGCCAGCACTTCACGTGGGCAGCCGGTTTCTTCCACCAGGCCTTTGTGCAGGAacac                                               <  3:110951/131‑1 (MQ=255)
                                                                                                                              gTGCAAAAATGCAGCTTTTACTTCAGGCTGCCGGAGAGGAACTGCTGCAGGCGCTCCGACTGCGGGTTGGCCAGCACTTCACGTGGGCAGCCGGTTTCTTCCACCAGGCCTTTGTGCAGGAACACCAACTGGTTGGACACCTCGCgg                         <  1:236571/147‑1 (MQ=255)
                                                                                                                               tGCAAAAATGCAGCTTTTACTTCAGGCTGCCGGAGAGGAACTGCTGCAGGCGCTCCGACTGCGGGTTGGCCAGCACTTCACGTGGGCAGCCGGTTTCTTCCACCAGGCCTTTGTGCAGGAACACCAACTGGTTGGACACCTCGCGGGc                       >  2:282483/1‑148 (MQ=255)
                                                                                                                                            cTTTTACTTCAGGCTGCCGGAGAGGAACTGCTGCAGGGGCTCCGACTGCGGGTTGGCCAGCACTTCACGTGGGCAGCCGGTTTCTTCCACCAGGCCTTTGTGCAGGAACACCAACTGGTTGGACACCTCGCGGGCAAAGCCCATTTCGt         <  8:350970/149‑1 (MQ=255)
                                                                                                                                            cTTTTACTTCAGGCTGCCGGAGAGGAACTGCTGCAGGCGCTCCGACTGCGGGTTGGCCAGCACTTCACGTGGGCAGCCGGTTTCTTCCACCAGGCCTTTGTGCAGGAACACCAACTGGTTGGACACCTCGCGGGCAAAGCCCATTTCGt         <  5:13402/149‑1 (MQ=255)
                                                                                                                                             ttttACTTCAGGCTGCCGGAGAGGAACTGCTGCAGGCGCTCCGACTGCGGGTTGGCCAGCACTTCACGTGGGCAGCCGGTTTCTTCCACCAGGCCTTTGTGCAGGAACACCAACTGGTTGGACACCTCGCGGGCAAAGCCCATTTCGTg        >  6:133900/1‑149 (MQ=255)
                                                                                                                                                tACTTCAGGCTGCCTGAGAGGAACTGCTGCAGGCGCTCCGACTGCGGGTTTATCTAGCACTTCTCGTGGGCAGCCGGTTTCTTCCACCAGGCCTTTGTGCAGGAACACCAACTGGTTGGACACCTCGCGGGCAAAGCCCATTTCGTgg       <  1:127132/148‑1 (MQ=255)
                                                                                                                                                 aCTTCAGGCTGCCGGAGAGGAACTGCTGCAGGCGCTCCGACTGCGGGTTGGCCAGCACTTCACGTGGGCAGCCGGTTTCTTCCACCAGGCCTTTGTGCAGGAACACCAACTGGTTGGACACCTCGCGGGCAAAGCCCATTTCGTGGGt     <  3:8025/148‑1 (MQ=255)
                                                                                                                                                 aCTTCAGGCTGCCGGAGAGGAACTGCTGCAGGCGCTCCGACTGCGGGTTGGCCAGCACTTCACGTGGGCAGCCGGTTTCTTCCACCAGGCCTTTGTGCAGGAACACCAACTGGTTGGACACCTCGCGGGCAAAGCCCATTTCGTGGGt     <  4:361310/148‑1 (MQ=255)
                                                                                                                                                   ttCAGGCTGCCGGAGAGGAACTGCTGCAGGCGCTCCGACTGCGGGTTGGCCAGCACTTCACGTGGGCAGCCGGTTTCTTCCACCAGGCCTTTGTGCAGGAACACCAACTGGTTGGACACCTCGCGGGCAAAGCCCATTTCGTGGGTcac  >  8:123492/1‑149 (MQ=255)
                                                                                                                                                   ttCAGGCTGCCGGAGAGGAACTGCTGCAGGCGCTCCGACTGCGGGTTGGCCAGCACTTCACGTGGGCAGCAGGTTTCTTCCACCAGGCGTTTGTGCAGGAACACCAACTGCTTGGACACCTCGCGGGCAAAGCCCATTTCGTGGGacac  >  5:125869/1‑149 (MQ=255)
                                                                                                                                                   ttCAGGCTGCCGGAGAGGAACTGCTGCAGGCGCTCCGACTGCGGGTACGCCAGCACTTCACGTGGGCAGCCTGTTTCTTCCACCAGGCCTTTGTGCAGGAACACCAACAGCTTGGACAGCTCGCGGGCAAAGACCATTTCGTGGGTcac  >  7:162173/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GTGCTGGGCCAGATGAGAAAACCGATTCGCTGGGTGGTCATAGGGTACGGTCCGAAACCTGAGGTCGTTCGA‑‑GTCTGTGCGCCGGGTCAGGCTGCGCTCGCCGCGCAGCATGCCCTATTCTGGTGCAAAAATGCAGCTTTTACTTCAGGCTGCCGGAGAGGAACTGCTGCAGGCGCTCCGACTGCGGGTTGGCCAGCACTTCACGTGGGCAGCCGGTTTCTTCCACCAGGCCTTTGTGCAGGAACACCAACTGGTTGGACACCTCGCGGGCAAAGCCCATTTCGTGGGTCAC  >  NC_002947/5093581‑5093872

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