Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,365,164 A→C 6.8% intergenic (‑25/‑48) PP_3841 ← / → PP_3842 MerR family transcriptional regulator/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,365,1640AC 6.8% 220.5 / 7.4 74intergenic (‑25/‑48)PP_3841/PP_3842MerR family transcriptional regulator/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (33/36);  new base C (2/3);  total (35/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.94e-01

TTCACTGTAACTGCGATAGTTGTTGGCTTGTCGGTGTGGCTCGATAAGCCCGAGCTGTTCGTAATAGCGCAATGTATGGCGGCTGAGCCCACTCCGGCGCTCCAACTCTCCAATATTCACGTTAAAAATCCTTGACTGTAGAGCAGGCTCTACACCTTAGCCTGCGCCCATCATTCAAACAAGGGTTCAAACATGAGTGTGGAGTGCTTTGTCACCGGCGGCAGCGGCTTTGTGGGGCAACATTTGCTGGCCGGTCTTACTGCGCAAGGGCACAAGACCTGGGTGCTGATGC  >  NC_002947/4365020‑4365311
                                                                                                                                                |                                                                                                                                                    
ttCACTGTAACTGCGATAGTTGTTGGCTTGTCGGTGTGGCTCGATAAGCCCGAGCTGTTCGTAATAGCGCAATGTATGGCGGCTGAGCCCACTCCGGCGCTCCAACTCTCCAATATTCACGTTAAAAATCCTTGACTGTAGAGCAg                                                                                                                                                     <  5:107500/146‑1 (MQ=255)
    cTGTAACTGCGATAGTTGTTGGCTTGTCGGTGTGGCTCGATAAGCCCGAGCTGTTCGTAATAGCGCAATGTATGGCGGCTGAGCCCACTCCGGCGCTCCAACTCTCCAATATTCACGTTAAAAATCCTTGACTGTAGAGCAGGCTCTac                                                                                                                                              <  6:37266/149‑1 (MQ=255)
     tGTAACTGCGATAGTTGTTGGCTTGTCGGTGTGGCTCGATAAGCCCGAGCTGTTAGTAATAGCGCAATGTATGGCGGCTGAGCCCACTCCGGCGCTCCAACTCTCCAATATTCACGTTAAAAATCCTTGACTGTAGAGCAGGCTCTaca                                                                                                                                             >  4:97012/1‑149 (MQ=255)
         aCTGCGATAGTTGTTGGCTTGTCGGTGTGGCTCGATAAGCCCGAGCTGTTCGTAATAGCGCAATGTATGGCGGCTGAGCCCACTCCGGCGCTCCAACTCTCCAATATTCACGTTAAAAATCCTTGACTGTAGAGCCTGCTCTacagtca                                                                                                                                         <  4:277102/149‑5 (MQ=255)
          ccgCGATAGTTGTTGGCTTGTCGGTGTGGCTCGATAAGCCCGAGCTGTTCGTAATAGCGCAATGTATGGCGGCTGAGCCCACTCCGGCGCTCCAACTCTCCAATATTCACGTTAAAAATCCTTGACTGTAGAGCAGGCTCTACACCTTa                                                                                                                                        <  5:97077/147‑1 (MQ=255)
          cTGCGATAGTTGTTGGCTTGTCGGTGTGGCTCGATAAGCCCGAGCTGTTCGTAATAGCGCAATGTATGGCGGCTGAGCCCACTCCGGCGCTCCAACTCTCCAATATTCACGTTAAAAATCCTTGACTGTAGAGCCTGCTCTacagtcaa                                                                                                                                        <  5:41064/149‑6 (MQ=255)
          cTGCGATAGTTGTTGGCTTGTCGGTGTGGCTCGATAAGCCCGAGCTGTTCGTAATAGCGCAATGTATGGCGGCTGAGCCCACTCCGGCGCTCCAACTCTCCAATATTCACGTTAAAAATCCTTGACTGTAGAGCAGGCTCTACACCTTa                                                                                                                                        <  2:267931/149‑1 (MQ=255)
          cTGCGATAGTTGTTGGCTTGGCGGTGTGGCTCGATAAGCCCGAGCTGTTCGTAATAGCGCAATGTATGGCGGCTGAGCCCACTCCGGCGCTCCAACTCTCCAATATTCACGTTAAAAATCCTTGACTGTAGAGCCTGCTCTacagtcaa                                                                                                                                        <  5:101566/149‑6 (MQ=255)
             cGATAGTTGTTGGCTTGTCGGTGTGGCTCGATAAGCCCGAGCTGTTCGTAATAGCGCAATGTATGGCGGCTGAGCCCACTCCGGCGCTCCAACTCTCCAATATTCACGTTAAAAATCCTTGACTGTAGAGCAGGCTCTACACCTTAGc                                                                                                                                      <  5:335831/148‑1 (MQ=255)
                          catGTCGGTGTGGCGCGATAAGCCCGAGCTGTCCGTAATAGCGCAATCTATCCCGCCTTACCCCACTCCGGCGCTCCAACTCTCCAATATTCACCTTAAAAATCCTTGACTGTAGACCAGGCTCTACACCTTAGCCTGCGCCCATCAtt                                                                                                                        <  6:234224/147‑1 (MQ=255)
                              tCGGTGTGGCTCGATAAGCCCGAGCTGTTCGTAATAGCGCAATGTATGGCGGCTGAGCCCACTCCGGCGCTCCAACTCTCCAATATTCACGTTAAAAATCCTTGACTGTAGAGCAGGCTCTACACCTTAGCCTGCGCCCATCATTCaaa                                                                                                                    <  5:179068/149‑1 (MQ=255)
                              tCGGTGTGGCTCGATAAGCCCGAGCTGTTCGTAATAGCGCAATGTATGGCGGCTGAGCCCACTCCGGCGCTCCAACTCTCCAATATTCACGTTAAAAATCCTTGACTGTAGAGCAGGCTCTACACCTTAGCCTGCGCCCATCATTCaaa                                                                                                                    >  2:247387/1‑149 (MQ=255)
                                 gtgtGGCTCGATAAGCCCGAGCTGTTCGTAATAGCGCAATGTATGGCGGCTGAGCCCACTCCGGCGCTCCAACTCTCCAATATTCACGTTAAAAATCCTTGACTGTAGAGCAGGCTCTACACCTTAGCCTGCGCCCATCATTCaaaa                                                                                                                   >  5:175580/1‑146 (MQ=255)
                                 gtgtGGCTCGATAAGCCCGAGCTGTTCGTAATAGCGCAATGTATGGCGGCTGAGCCCACTCCGGCGCTCCAACTCTCCAATATTCACGTTAAAAATCCTTGACTGTAGAGCAGGCTCTACACCTTAGCCTGCGCCCATCATTCAAAc                                                                                                                   >  2:139483/1‑147 (MQ=255)
                                      gCTCGATAAGCCCGAGCTGTTCGTAATAGCGCAATGTATGGCGGCTGAGCCCACTCCGGCGCTCCAACTCTCCAATATTCACGTTAAAAATCCTTGACTGTAGAGCAGGCTCTACACCTTAGCCTGCGCCCATCAtt                                                                                                                        <  7:147717/137‑1 (MQ=255)
                                      gCTCGATAAGCCCGAGCTGTTCGTAATAGCGCAATGTATGGCGGCTGAGCCCACTCCGGCGCTCCAACTCTCCAATATTCACGTTAAAAATCCTTGACTGTAGAGCAGGCTCTACACCTTAGCCTGCGCCCATCAtt                                                                                                                        >  8:147717/1‑137 (MQ=255)
                                      gCTCGATAAGCCCGAGCTGTTCGTAATAGCGCAATGTATGGCGGCTGAGCCCACTCCGGCGCTCCAACTCTCCAATATTCACGTTAAAAATCCTTGACTGTAGAGCAGGCTCTACACCTTAGCCTGCGCCCATCATTCAAACAAg                                                                                                                >  3:238426/1‑145 (MQ=255)
                                      gCTCGATAAGCCCGAGCTGTTCGTAATAGCGCAATGTATGGCGGCTGAGCCCACTCCGGCGCTCCAACTCTCCAATATTCACGTTAAAAATCCTTGACTGTAGAGCAGGCTCTACACCTTAGCCTGCGCCCATCATTCAAACAAg                                                                                                                <  4:238426/145‑1 (MQ=255)
                                      gCTCGATAAGCCCGAGCTGTTCGTAATAGCGCAATGTATGGCGGCTGAGCCCACTCCGGCGCTCCAACTCTCCAATATTCACGTTAAAAATCCTTGACTGTAGAGCAGGCTCTACACCTTAGCCTGCGCCCATCATTCAAACAAg                                                                                                                <  2:194250/145‑1 (MQ=255)
                                      gCTCGATAAGCCCGAGCTGTTCGTAATAGCGCAATGTATGGCGGCTGAGCCCACTCCGGCGCTCCAACTCTCCAATATTCACGTTAAAAATCCTTGACTGTAGAGCAGGCTCTACACCTTAGCCTGCGCCCATCATTCAAACAAg                                                                                                                >  1:194250/1‑145 (MQ=255)
                                      gCTCGATAAGCCCGAGCTGTTCGTAATAGCGCAATGTATGGCGGCTGAGCCCACTCCGGCGCTCCAACTCTCCAATATTCACGTTAAAAATCCTTGACTGTAGAGCAGGCTCTACACCTTAGCCTGCGCCCATCATTCAAACAAGGGtt                                                                                                           >  4:306702/1‑149 (MQ=255)
                                          gATAAGCCCGAGCTGTTCGTAATAGCGCAATGTATGGCGGCTGAGCCCACTCCGGCGCCCCAACTCTCCAATATTCACGTTAAAAATCCTTGACTGTAGAGCAGGCTCTACACCTTAGCCTGCGCCCATCATTCAAACAAGGGTTCaaa                                                                                                       <  1:71517/149‑1 (MQ=255)
                                            tAAGCCCGAGCTGTTCGTAATAGCGCAATGTATGGCGGCTGAGCCCACTCCGGCGCTCCAACTCTCCAATATTCACGTTAAAAATCCTTGACTGTAGAGCAGGCTCTACACCTTAGCCTGCGCCCATCATTCAAACAAGGGTTCAAACa                                                                                                     <  3:144564/149‑1 (MQ=255)
                                              aGCCCGAGCTGTTCGTAATAGCGCAATGTATGGCGGCTGAGCCCACTCCGGCGCTCCAACTCTCCAATATTCACGTTAAAAATCCTTGACTGTAGAGCAGGCTCTACACCTTAGCCTGCGCCCATCATTCAAACAAGGGTTCAAACAt                                                                                                    <  8:7373/148‑1 (MQ=255)
                                              aGCCCGAGCTGTTCGTAATAGCGCAATGTATGGCGGCTGAGCCCACTCCGGCGCTCCAACTCTCCAATATTCACGTTAAAAATCCTTGACTGTAGAGCAGGCTCTACACCTTAGCCTGCGCCCATCATTCAAACAAGGGTTCAAACATg                                                                                                   <  1:139483/149‑1 (MQ=255)
                                               gCCCGAGCTGTTCGTAATAGCGCAATGTATGTCGGCTGAGCCCACTCCGGCGCTCCAACTCTCCAATATTCACGTTAAAAATCCTTGACTGTAGAGCAGGCTCTACACCTTAGCCTGCGCCCATCATTCAAACAAGGGTTCAAACATGa                                                                                                  <  7:270202/149‑1 (MQ=255)
                                               gCCCGAGCTGTTCGTAATAGCGCAATGTATGGCGGCTGAGCCCACTCCGGCGCTCCAACTCTCCAATATTCACGTTAAAAATCCTTGACTGTAGAGCAGGCTCTACACCTTAGCCTGCGCCCATCATTCAAACAAGGGTTCAAACATGa                                                                                                  <  3:306702/149‑1 (MQ=255)
                                                  cGAGCTGTTCGTAATAGCGCAATGTATGGCGGCTGAGCCCACTCCGGCGCTCCAACTCTCCAATATTCACGTTAAAAATCCTTGACTGTAGAGCAGGCTCTACACCTTAGCCTGCGCCCATCATTCAAACAAGGGTTCAAACATGAgtg                                                                                               >  6:69596/1‑149 (MQ=255)
                                                     gCTGTTCGTAATAGCGCAATGTATGGCGGCTGAGCCCACTCCGGCGCTCCAACTCTCCAATATTCACGTTAAAAATCCTTGACTGTAGAGCAGGCTCTACACCTTAGCCTGCGCCCATCATTCAAACAAGGGTTCACACATGAGTGTgg                                                                                            >  8:212979/1‑149 (MQ=255)
                                                         ttCGTATTAGCGCAAAGTATGGCGGCTGAGCCCACTCCGGCGCTCCAACTCTCCAATATTCACGTTAAAAATCCTTGACTGGAGAGCAGGCTCTACACCTTAGCCTGCGCCCATCATTCAAAGAAGGTTTCAAACATGa                                                                                                  <  4:36744/139‑1 (MQ=255)
                                                         ttCGTAATAGCGCAATGTATGGCGGCTGAGCCCACTCCGGCGCTCCAACACTCCAATATTCACGTTAAAAATCCATGACTGTAGCGCAGGCTCTACACCTTAGCCTGCGCCCATCATTCAAACAAGGGTTCAAACATGAgacctcgg                                                                                          >  3:36744/1‑140 (MQ=255)
                                                                tAGCGCAATGTATGGCGGCTGAGCCCACTCCGGCGCTCCAACTCTCCAATATTCACGTTAAAAATCCTTGACTGTAGAGCAGGCTCTACACCTTAGCCTGCGCCCATCATTCAAACAAGGGTTCAAACATGAgt                                                                                                >  1:285540/1‑134 (MQ=255)
                                                                tAGCGCAATGTATGGCGGCTGAGCCCACTCCGGCGCTCCAACTCTCCAATATTCACGTTAAAAATCCTTGACTGTAGAGCAGGCTCTACACCTTAGCCTGCGCCCATCATTCAAACAAGGGTTCAAACATGAgt                                                                                                <  2:285540/134‑1 (MQ=255)
                                                                tAGCGCAATGTATGGCGGCTGAGCCCACTCCGGCGCTCCAACTCTCCAATATTCACGTTAAAAATCCTTGACTGTAGAGCAGGCTCTACACCTTAGCCTGCGCCCATCATTCAAACAAGGGTTCAAACATGAGTGTGGAGt                                                                                         <  8:203918/141‑1 (MQ=255)
                                                                tAGCGCAATGTATGGCGGCTGAGCCCACTCCGGCGCTCCAACTCTCCAATATTCACGTTAAAAATCCTTGACTGTAGAGCAGGCTCTACACCTTAGCCTGCGCCCATCATTCAAACAAGGGTTCAAACATGAGTGTGGAGt                                                                                         >  7:203918/1‑141 (MQ=255)
                                                                tAGCGCAATGAATGGCGGCTGAGCCCACTCCGGCGCTCCAACTCTCCAATATTCACGTTAAAAATCCTTGACTGTAGAGCAGGCTCTACACCTTAGCCTGCGCCCATCATTCAAACAAGGGTTCAAACATGAGTGTGGAGTGCTTTGTc                                                                                 <  8:25472/149‑1 (MQ=255)
                                                                  gcgcAATGTATGGCGGCTGAGCCCACTCCGGCGCTCCAACTCTCCAATATTCACGTTAAAAATCCTTGACTGTAGAGCAGGCTCTACACCTTAGCCTGCGCCCATCATTCAAACAAGGGTTCAAACATGAGTGTGGAGTGCTTTGTCAc                                                                               >  1:163493/1‑149 (MQ=255)
                                                                   cgcAATGTATGGCGGCTGAGCCCACTCCGGCGCTCCAACTCTCCAATATTCACGTTAAAAATCCTTGACTGTAGAGCAGGCTCTACCCCTTAGCCTGCGCCCATCATTCAAACAAGGGTTCAAACATGAGTGTGGAGTGCTTTGTCAcc                                                                              >  7:170145/1‑149 (MQ=255)
                                                                      aaTGTATGGCGGGTTGTCCCAATCCGGCGCTCCAACTCTACAATATTTACGTTAAAAATTCTTGCCTGTTGATGAGGCTCTACACCTTAGCCTGCGCCCATCATTCAAACAAGGGTTCAAACATGAGTGTGGAGTGCTTTGTCACcggc                                                                           <  6:85174/149‑1 (MQ=255)
                                                                      aaTGTATGGCGGCTGAGCCCACTCCGGCGCTCCAACTCTCCAATATTCACGTTAAAAATCCTTGACTGTAGAGCAGGCTCTACACCTTAGCCTGCGCCCATCATTCAAACAAGGGTTCAAACATGAGTGTGGAGTGCTTTGTCACcggc                                                                           >  4:252527/1‑149 (MQ=255)
                                                                            tGGCGGCTGAGCCCACTCCGGCGCTCCAACTCTCCAATATTCACGTTAAAAATCCTTGACTGTAGAGCAGGCTCTACACCTTAGCCTGCGCCCATCATTCAAACAAGGGTTCAAACATGAGTGTGGAGTGCTTTGTCACCGGCGGCAGc                                                                     >  3:293739/1‑149 (MQ=255)
                                                                              gcggcTGAGCCCACTCCGGCGCTCCTACTCTCCAATATTGACGTTAAAAATCCTTGACTGTAGAGCAGGCTCTACACCTTAGCCTGCGACCATCATTCATACAAGGGTTCAAACATGAGTGTGGAGTGCTTTGTCACCGGCGGGAGCgg                                                                   >  6:62187/1‑149 (MQ=255)
                                                                              ccggcTGAGCCCACTCCGGCGCTCCAACTCTCCAATATTCACGTTAAAAATCCTTGACTGTAGAGCAGGCTCTACACCTTAGCCTCCGCCCATCATTCAAACAAGGGTTCAAACATGAGTCTCGAGTCCTTTCTCACCGGAGGCAGCgg                                                                   >  6:66805/2‑149 (MQ=255)
                                                                                  cTGAGCCCACTCCGGCGCTCCAACTCTCCAATATTCACGTTAAAAATCCTTGACTGTAGAGCAGGCTCTACACCTTAGCCTGCGCCCATCATTCAAACAAGGGTTCAAACATGAGTGTGGAGTGCTTTGTCACCGGCGGCAGCGGCttt                                                               <  3:97012/149‑1 (MQ=255)
                                                                                     aGCCCACTCCGGCGCTCCAACTCTCCAATATTCACGTTAAAAATCCTTGACTGTAGAGCAGGCTCTACACCTTAGCCTGCGCCCATCATTCAAACAAGGGTTCAAACATGAGTGTGGAGTGCTTTGTCACCGGCGGCAGCGGCTTtgtg                                                            <  4:293739/149‑1 (MQ=255)
                                                                                     aGCCCACTCCGGCGCTCCAACTCTCCAATATTCACGTTAAAAATCCTTGACTGTAGAGCAGGCTCTACACCTTAGCCTGCGCCCATCATTCAAACAAGGGTTCAAACATGAGTGTGGAGTGCTTTGTCACCGGCGGCAGCGGCTTtgtg                                                            <  1:164843/149‑1 (MQ=255)
                                                                                     aGCCCACTCCGGCGCTCCAACTCTCCAATATTCACGTTAAAAATCCTTGACTGTAGAGCAGGCTCTACACCTTAGCCTGCGCCCATCATTCAAACAAGGGTTCAAACATGAGTGTGGAGTGCTTTGTCACCGGCGGCAGCGGCTTtgtg                                                            <  8:66385/149‑1 (MQ=255)
                                                                                       cccACTCCGGCGCTCCAACTCTCCAATATTCACGTTAAAAATCCTTGACTGTAGAGCAGGCTCTACACCTTAGCCTGCGCCCATCATTCAAACAAGGGCTTCAAACATGAGTGTGGAGTGCTTTGTCACCGGCGGCAGCGGCTTTGTgg                                                           <  1:247387/149‑1 (MQ=255)
                                                                                        ccACTCCGGCGCTCCAACTCTCCAATATTCACGTTAAAAATCCTTGACTGTAGAGCAGCCTCTCCACCGTAGCCTGCGCCCATCATTCAAACAAGGGTTCAAACATGAGTGGGGCGTGCTTTGTCACAGGCGGCAGCGGCTTTGTgggg                                                         >  6:186389/1‑149 (MQ=255)
                                                                                         cACTCCGGCGCTCCAACTCTCCAATATTCACGTTAAAAATCCTTGACTGTAGAGCAGGCTCTACACCTTAGCCTGCGCCCATCATTCAAACAAGGGTTCAAACATGAGTGTGGAGTGCTTTGTCACCGGCGGCAGCGGCTTTGTgggg                                                         >  4:288525/1‑148 (MQ=255)
                                                                                          aCTCCGGCGCTCCAACTCTCCAATATTCACGTTAAAAATCCTTGACTGTAGAGCAGGCTCTACACCTTAGCCTGCGCCCATCATTCAAACAAGGGTTCAAACATGAGTGTGGAGTGCTTTGTCACCGGCGGCAGCGGCTTTGTGGGGc                                                        >  3:18942/1‑148 (MQ=255)
                                                                                           cTCCGGCTCTCAAACTCTCCAATATTCACGTTAAAAATACTTGACTGTAGAGCAGGCTCTACACCTTAGCCTGCGCCCATCATTCAAACAAGGGTTCAAACATGAGTGTGGAGTGCTTTGTCACCGGCGGCAGCGGCTTTGTGGGGCaa                                                      <  4:18942/149‑1 (MQ=255)
                                                                                            tCCGTCGCTCCAACTCTCCAATATTCACGTTAAAAATCCTTGACTGTAGAGCAGGCTCTACACCTTAGCCTGCGCCCATCATTCAAACAAGGGTTCAAACATGa                                                                                                  <  5:276866/104‑1 (MQ=255)
                                                                                            tCCGGCGCTCCAACTCTCCAATATTCAAGTTAAAAATCCTTGACTGTAGAGCAGGCTCTACACCTTAGCCTGCGCCCATCATTCAAACAAGGGTTCAAACATGa                                                                                                  >  6:276866/1‑104 (MQ=255)
                                                                                                  gcTCCAACTCTCCAATATTCACGTTAAAAATCCTTGACTGTAGAGCAGGCTCTACACCTTAGCCTGCGCCCATCATTCAAACAAGGGTTCAAACATGAGTGTGGAGTGCTTTGTCACCGGCGGCAGCGGCTTTGTGGGGCAACATTTGc                                               >  5:235387/1‑149 (MQ=255)
                                                                                                  gcTCCAACTCTCCAATATTCACGTTAAAAATCCTTGACTGTAGAGCAGGCTCTACACCTTAGCCTGCGCCCATCATTCAAACAAGGGTTCAAACATGAGTGTGGAGTGCTTTGTCACCGGCGGCAGCGGCTTTGTGGGGCAACATTTGc                                               >  8:45268/1‑149 (MQ=255)
                                                                                                     ccAACTCTCCAATATTCACGTTAAAAATCCTTGACTGTAGAGCAGGCTCTACACCTTAGCCTGCGCCCATCATTCAAACAAGGGTTCAAACATGAGTGTGGAGTGCTTTGTCACCGGCGGCAGCGGCTTTGTGGGGCAACATTTGCTgg                                            >  4:220568/1‑149 (MQ=255)
                                                                                                       aaCTCTCCAATATTCACGTTAAAAATCCTTGACTGTAGAGCAGGCTCTACACCTTAGCCTGCGCCCATCATTCAAACAAGGGTTCAAACATGAGTGTGGAGTGCTTTGTCACCGGCGGCAGCGGCTTTGTGGGGCAACATTTGCt                                              >  4:209865/1‑145 (MQ=255)
                                                                                                       aaCTCTCCAATATTCACGTTAAAAATCCTTGACTGTAGAGCAGGCTCTACACCTTAGCCTGCGCCCATCATTCAAACAAGGGTTCAAACATGAGTGTGGAGTGCTTTGTCACCGGCGGCAGCGGCTTTGTGGGGCAACATTTGCt                                              <  3:209865/145‑1 (MQ=255)
                                                                                                                   ttCACGTTAAAAATCCTTGACTGTAGAGCAGGCTCTACACCTTAGCCTGCGCCCATCATTCAAACAAGGGTTCACACATGAGTGTGGAGTGCTTTGTCACCGGCGGCAGCGGCTTTGTGGGGCAACATTTGCTGGCCGGTCTTACTgc                               <  7:212979/148‑1 (MQ=255)
                                                                                                                     cACGTTAAAAATCCTTGACTGTAGAGCAGGCTCTACACCTTAGCCTGCGCCCATCATTCAAACAAGGGTTCAAACATGAGTGTGGAGTGCTTTGTCACCGGCGGCAGCGGCTTTGTGGGGCAACATTTGCTGGCCAGTCTTACTgcgc                             <  6:235387/148‑1 (MQ=255)
                                                                                                                                 ccTTGACTGTAGAGCAGGCTCTACAGCTTAGCCTGCGCCCTTCTTTCTTACAAGGGTTGATACATGTGTTTGGTGTGGTTTGTCACGGGCGGCAGGGGCTTTGTGGGGCAACTTTTGCTGGCCGGTCTTAGTGCGCATGGGCTCAAGa                 <  4:108381/148‑1 (MQ=255)
                                                                                                                                  ctaaggtGTAGAGCCTGCTCTACACCTTAGCCTGCGCCCATCATTCAAACAAGGGTTCAAACATGAGTGTGGAGTGCTTTGTCACCGGCGGCAGCGGCTTTGTGGGGCAACATTTGCTGGCCGGTCTTACTGCGCAAGGGCACAAGAcc               >  4:224057/7‑149 (MQ=255)
                                                                                                                                   ttGACTGTAGAGCAGGCTCTACACCTTAGCCTGCGCCCATCATTCAAACAAGGGTTCAAACATGAGTGTGGAGTGCTTTGTCACCGGCGGCAGCGGCTTTGTGGGGCAACATTTGCTGGCCGGTCTTACTGCGCAAGGGCACAAGACCt              >  2:266182/1‑149 (MQ=255)
                                                                                                                                      aCTGTAGAGCAGGCTCTACACCTTAGCCTGCGCCCATCATTCAAACAAGGGTTCAAACATGAGTGTGGAGTGCTTTGTCACCGGCGGCAGCGGCTTTGTGGGGCAACATTTGCTGGCCGGTCTTACTGCGCAAGGGCACAAGACCTgg            <  5:69596/148‑1 (MQ=255)
                                                                                                                                        tGTAGAGCCTGCTCTACACCTTAGCCTGCGCCCATAATTCAAACAAGGGTTCAAACATGAGTGTGGAGTGCTTTGTCACCGGCGGCAGCGGCTTTGTGGGGCAACATTTTCTGTCCGGTCTTACTGCGCAAGGGCACAAGACCTGGGTg         >  4:40444/1‑149 (MQ=255)
                                                                                                                                        tGTAGAGCAGGCTCTACACCTTAGCCTGCGCCCATCATTCAAACAAGGGTTCAAACATGAGTGTGGAGTGCTTTGTCACCGGCGGCAGCGGCTTTGTGGGGCAAc                                                     >  2:112032/1‑105 (MQ=255)
                                                                                                                                        tGTAGAGCAGGCTCTACACCTTAGCCTGCGCCCATCATTCAAACAAGGGTTCAAACATGAGTGTGGAGTGCTTTGTCACCGGCGGCAGCGGCTTTGTGGGGCAAc                                                     <  1:112032/105‑1 (MQ=255)
                                                                                                                                          tAGCGCAGGCTCTACACCTTAGCCTGCGCCCATCATTCAAACAAGGGTTCAAACATGAGTGTGGAGTGCTTTGTCACCGGCGGCAGCGGCTTTGTGGGGCAACATATGCTGGCCGGTCTTACTGCGGAAGGGCACAAGACCTGTGTGCt       >  2:303420/1‑149 (MQ=255)
                                                                                                                                            gagCAGGCTCTACACCTTAGCCTGCGCCCATCATTCAAACAAGGGTTCAAACATGAGTGTGGAGTGCTTTGTCACCGGCGGCAGCGGCTTTGTGGGGCAACAATTGCTGGCCGGTCTTACTGCGCAAGGGCACAAGACCTGGGTGCTGa     >  4:54419/1‑149 (MQ=255)
                                                                                                                                             agCAGGCTCTACACCTTAGCCTGCGCCCATCATTCAAACAAGGGTTCAAACATGAGTGTGGAGTGCTTTGTCACCGGCGGCAGCGGCTTTGTGGGGCAACATTTGCTGGCCGGTCTTACTGCGCAAGGGCACAAGACCTGGGTGCTGAt    <  8:170145/149‑1 (MQ=255)
                                                                                                                                              gCAGGCTCTACACCTTAGCCTGCGCCCATCATTCAAACAAGGGTTCAAACATGAGTGTGGAGTGCTTTGTCACCGGCGGCAGCGGCTTTGTGGGGCAACATTTGCTGGCCGGTCTTACTGCGCAAGGGCACAAGACCTGGGTGCTGATg   >  4:237820/1‑149 (MQ=255)
                                                                                                                                               cAGGCTTTACAACTTAGCCTGCGCCCATCATTCAAACAAGGGTTCAAACATGAGTGTGGAGTGCTTTGTCACCGGCGGCAGCGGCTTTGTGGGGCAACATTTGCTGGCCGGTCTTACTGCGCAAGGGCACAAGACCTGGGTGCTGATgc  >  5:301725/1‑149 (MQ=255)
                                                                                                                                               cAGGCTCTACACCTTAGCCTGCGCCCATCATTCAAACAAGGGTTCAAACATGAGTGTGGAGTGCTTTGTCACCGGCGGCAGCGGCTTTGTGGGGCAACATTTGCTGGCCGGTCTTACTGCGCAAGGGCACAAGACCTGGGTGCTGATgc  <  3:220568/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
TTCACTGTAACTGCGATAGTTGTTGGCTTGTCGGTGTGGCTCGATAAGCCCGAGCTGTTCGTAATAGCGCAATGTATGGCGGCTGAGCCCACTCCGGCGCTCCAACTCTCCAATATTCACGTTAAAAATCCTTGACTGTAGAGCAGGCTCTACACCTTAGCCTGCGCCCATCATTCAAACAAGGGTTCAAACATGAGTGTGGAGTGCTTTGTCACCGGCGGCAGCGGCTTTGTGGGGCAACATTTGCTGGCCGGTCTTACTGCGCAAGGGCACAAGACCTGGGTGCTGATGC  >  NC_002947/4365020‑4365311

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