Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,092,669 A→G 6.2% L318S (TTG→TCG)  PP_4482 ← AraC family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,092,6690AG 6.2% 181.8 / 2.6 64L318S (TTG→TCG) PP_4482AraC family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (28/32);  new base G (2/2);  total (30/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.28e-01

GTACGGCGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTCCCGCGGCGTTGCACCAAAGAAATTACGGTAGGCACTGGAAAAATGCGGCCCCGAGGAGAAGCCACAGGACAGCCCGATCTGGATGATCGACTTGCTGGTCTGCATC  >  NC_002947/5092532‑5092815
                                                                                                                                         |                                                                                                                                                  
gTACGGCGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACggg                                                                                                                                         >  1:338825/1‑149 (MQ=255)
 tACGGCGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGgc                                                                                                                                        <  2:2014/149‑1 (MQ=255)
      cGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGCACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGCGCTCAATTCAACCGGgctgctg                                                                                                                                   >  8:369619/1‑149 (MQ=255)
         tGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGgctgctgct                                                                                                                                 <  2:378718/148‑1 (MQ=255)
           ttGATGTGTTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCg                                                                                                                              <  2:357892/149‑1 (MQ=255)
            tGATGGGTTGACGGGGAGATCGAGCGCCTCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCg                                                                                                                              <  5:32125/148‑1 (MQ=255)
                 tgtTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCt                                                                                                                         <  3:363828/148‑1 (MQ=255)
                  gtTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGa                                                                                                                       >  7:231221/1‑149 (MQ=255)
                   ttGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCACCCCCAAGGCATACCCCCCGCACCCCCCCGAAATCTCAGGCCTTTTCAGCCGGCGCCGAGCTGAATTCAAAGGGGCTGCTGCTGCTACGCTGAt                                                                                                                      >  4:78513/1‑149 (MQ=255)
                   ttGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCGATTCAAACGGGCTGCTGCTGCGACGCTGAt                                                                                                                      >  6:107533/1‑149 (MQ=255)
                   ttGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCGATTCAAACGGGCTGCTGCTGCGACGCTGAt                                                                                                                      >  6:101487/1‑149 (MQ=255)
                   ttGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGAt                                                                                                                      <  8:231221/149‑1 (MQ=255)
                         gggAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGt                                                                                                                <  6:131194/149‑1 (MQ=255)
                         gggAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGt                                                                                                                <  8:329241/149‑1 (MQ=255)
                                cGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTCCcg                                                                                                         >  7:5509/1‑149 (MQ=255)
                                    cgccgcGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTCCCgcggc                                                                                                     >  7:74155/1‑149 (MQ=255)
                                    cgccgcGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTCCCgcggc                                                                                                     >  6:412455/1‑149 (MQ=255)
                                      ccgcGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTCCCGCGGCGt                                                                                                   >  1:294123/1‑149 (MQ=255)
                                          gTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTCCCGCGGCGTTGCa                                                                                               >  1:182385/1‑149 (MQ=255)
                                              gcggcgCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCCCCCGAAATCTCAGCCCTTTTCCGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTCCCGCGGCGTTGCACCa                                                                                            >  2:74956/1‑148 (MQ=255)
                                               cggcgCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTCCCGCGGCGTTGCACCaaa                                                                                          <  8:133305/149‑1 (MQ=255)
                                               cggcgCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTCCCGCGGCGTTGCACCaaa                                                                                          <  8:74155/149‑1 (MQ=255)
                                                     tCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTCCCGCGGCGTTGCACCAAAGAAAt                                                                                     >  6:166161/1‑148 (MQ=255)
                                                                  cAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTCCCGCGGCGTTGCAc                                                                                              >  6:209230/1‑126 (MQ=255)
                                                                  cAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTCCCGCGGCGTTGCAc                                                                                              <  5:209230/126‑1 (MQ=255)
                                                                  cAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTCCCGCGGCGTTGCACCAAAGAAATTACGGTAGGCACTg                                                                       >  4:310039/1‑149 (MQ=255)
                                                                      acCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTCCCgcggc                                                                                                     <  5:20248/115‑1 (MQ=255)
                                                                      acCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTCCCgcggc                                                                                                     >  6:20248/1‑115 (MQ=255)
                                                                        cAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTCCCGCGGCGTTGCACCAAAGAAATTACGGTAGGCACTGGaaaaa                                                                 >  6:179664/1‑149 (MQ=255)
                                                                         aaCCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCGATTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTCCCGCGGCGTTGCACCAAAGAAATTACGGTAGGCACTGGAAAAAt                                                                <  5:107533/149‑1 (MQ=255)
                                                                         aaCCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCGATTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTCCCGCGGCGTTGCACCAAAGAAATTACGGTAGGCACTGGAAAAAt                                                                <  5:101487/149‑1 (MQ=255)
                                                                          aCCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTCCCGCGGCGTTGCACCAAAGAAATTACGGTAGGCACTGGAAAAATg                                                               <  3:89451/149‑1 (MQ=255)
                                                                            cccAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTCCCGCGGCGTTGCACCAAAGAAATTACGGTAGGCACTGGAAAAATGCg                                                             <  6:121060/149‑1 (MQ=255)
                                                                             ccAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTCCCGCGGCGTTGCACCAAAGAAATTACGGTAGGCACTGGAAAAATGCgg                                                            <  2:182385/149‑1 (MQ=255)
                                                                              cAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTCCCGCGGCGTTGCACCAAAGAAATTACGGTAGGCACTGGAAAAATGCGGc                                                           >  1:252638/1‑149 (MQ=255)
                                                                              cAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTCCCGCGGCGTTGCACCAAAGAAATTACGGTAGGCACTGGAAAAATGCGGc                                                           >  1:259244/1‑149 (MQ=255)
                                                                              cAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTCCCGCGGCGTTGCACCAAAGAAATTACGGTAGGCACTGGAAAAATGCGGc                                                           >  6:165703/1‑149 (MQ=255)
                                                                                aGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTCCCGCGGCGTTGCACCAAAGAAATTACGGTAGGCACTGGAAAAATGCGGccc                                                         <  2:252638/149‑1 (MQ=255)
                                                                                aGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTCCCGCGGCGTTGCACCAAAGAAATTACGGTAGGCACTGGAAAAATGCGGccc                                                         <  2:259244/149‑1 (MQ=255)
                                                                                       cacCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTCGCGCGGCGTTGCACCAAAGAAATTACGGTAGGCACTGGAAAAATGc                                                              <  1:307144/137‑1 (MQ=255)
                                                                                        acCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTCCCGCGGCGTTGCACCAAAGAAATTACGGTAGGCACTGGAAAAATGCGGCCCCgaggag                                                  >  7:398175/1‑148 (MQ=255)
                                                                                          cccGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTCCCGCGGCGTTGCACCAAAGAAATTACGGTAGGCACTGGAAAAATGCGGCCCCGAGGAGAAg                                               >  8:364437/1‑149 (MQ=255)
                                                                                           ccGACCCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTCCCGCGGCGTTGCACCAAAGAAATTACGGTAGGCACTGGAAAAATGCGGCCCCGAGGAGAAGc                                              <  2:294123/149‑1 (MQ=255)
                                                                                             gAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTCCCGCGGCGTTGCACCAAAGAAATTACGGTAGGCACTGGAAAAATGCGGCCCCGAGGAGAAGCca                                            <  3:310039/149‑1 (MQ=255)
                                                                                                  cacccGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTCCCGCGGCGTTGCACCAAAGAAATTACGGTAGGCACTGGAAAAATGCGGCCCCGAGGAGAAGCCACAGGa                                       <  7:364437/149‑1 (MQ=255)
                                                                                                       gAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTCCCGCGGCGTTGCACCAAAGAAATTACGGTAGGCACTGGAAAAATGCGGCCCCGAGGAGAAGCCACAGGACAGcc                                  <  5:179664/149‑1 (MQ=255)
                                                                                                        aaaTCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTCCCGCGGCGTTGCACCAAAGAAATTACGGTAGGCACTGGAAAAATGc                                                              <  4:77726/120‑1 (MQ=255)
                                                                                                        aaaTCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTCCCGCGGCGTTGCACCAAAGAAATTACGGTAGGCACTGGAAAAATGc                                                              >  3:77726/1‑120 (MQ=255)
                                                                                                         aaTCTGAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAATAGGGCTGCTTCTGCGTCGCTGATTGCGGTCTTCCCGCGGCGTTGCACCAAAGATATTACGGTAGGCACTGGAAAAATGCGGCGCCGAGGAGAATCCACAGGACAGCCCg                                <  8:5509/149‑1 (MQ=255)
                                                                                                         aaTCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTCCCGCGGCGTTGCACCAAAGAAATTACGGTAGGCACTGGAAAAATGCGGCCCCGAGGAGAAGCCACAGGACAGCCCg                                <  5:412455/149‑1 (MQ=255)
                                                                                                                gCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTACTGCGACGCTGATTGCGGTCTTCCCGCGGCGTTGCACCAAAGAAATTACGGTAGGCACTGGAAAAATGCGGCCCCGAGGAGAAGCCACAGGACAGCCCGATCTGga                         <  2:338825/149‑1 (MQ=255)
                                                                                                                    ttttCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTCCCGCGGCGTTGCACCAAAGAAATTACGGTAGGCACTGGAAAAATGCGGCCCCGAGGAGAAGCCACAGGACAGCCCGATCt                            <  7:214702/142‑1 (MQ=255)
                                                                                                                    ttttCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTCCCGCGGCGTTGCACCAAAGAAATTACGGTAGGCACTGGAAAAATGCGGCCCCGAGGAGAAGCCACAGGACAGCCCGATCt                            >  8:214702/1‑142 (MQ=255)
                                                                                                                    ttttCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTCCCGCGGCGTTGCACCAAAGAAATTACGGTAGGCACTGGAAAAATGCGGCCCCGAGGAGAAGCCACAGGACAGCCCGATCTGgatgat                     >  6:221848/1‑149 (MQ=255)
                                                                                                                      ttCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTCCCGCGGCGTTGCACCAAAGAAATTACGGTAGGCACTGGAAAAATGCGGCCCCGAGGAGAAGCCACAGGACAGCCCGATCTGGATGATc                    <  5:165703/148‑1 (MQ=255)
                                                                                                                         aGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTCCCGCGGCGTTGCACCAAAGAAATTACGGTAGGCACTGGAAAAATGCGGCCCCGAGGAGAAGCCACAGGACAGCCCGATCTGGATGATCGACt                >  1:312876/1‑149 (MQ=255)
                                                                                                                            cGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTCCCGCGGCGTTGCACCAAAGAAATTACGGTAGGCACTGGAAAAATGCGGCCCCGAGGAGAAGCCACAGGACAGCCCGATCTGGATGATCGACTTGc             >  7:254349/1‑149 (MQ=255)
                                                                                                                                gcCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTCCCGCGGCGTTGCACCAAAGAAATTACGGTAGGCACTGGAAAAATGCGGCCCCGAGGAGAAGCCACAGGACAGCCCGATCTGGATGATc                    >  7:262021/1‑138 (MQ=255)
                                                                                                                                gcCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTCCCGCGGCGTTGCACCAAAGAAATTACGGTAGGCACTGGAAAAATGCGGCCCCGAGGAGAAGCCACAGGACAGCCCGATCTGGATGATc                    <  8:262021/138‑1 (MQ=255)
                                                                                                                                    aGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTCCCGCGGCGTTGCACCAAAGAAATTACGGTAGGCACTGGAAAAATGCGGCCCCGAGGAGAAGCCACAGGACAGCCCGATCTGGATGATCGACTTGCTGGTCTGc     <  2:312876/149‑1 (MQ=255)
                                                                                                                                     gCTCAATTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTCCCGCGGCGTTGCACCAAAGAAATTACGGTAGGCACTGGAAAAATGCGGCCCCGAGGAGAAGCCACAGGACAGCCCGATCTGGATGATCGACTTGCTGGTCTGCa    <  3:193952/149‑1 (MQ=255)
                                                                                                                                      cTCAATTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTCCCGCGGCGTTGCACCAAAGAAATTACGGTAGGCACTGGAAAAATGCGGCCCCGAGGAGAAGCCACAGGACAGCCCGATCTGGATGATCGACTTGCTGGTCTGCAt   >  2:131573/1‑149 (MQ=255)
                                                                                                                                        cAATTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTCCCGCGGCGTTGCACCAAAGAAATTACGGTAGGCACTGGAAAAATGCGGCCCCGAGGAGAAGCCACAGGACAGCCCGATCTGGATGATCGACTTGCTg           >  6:303734/1‑139 (MQ=255)
                                                                                                                                        cAATTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTCCCGCGGCGTTGCACCAAAGAAATTACGGTAGGCACTGGAAAAATGCGGCCCCGAGGAGAAGCCACAGGACAGCCCGATCTGGATGATCGACTTGCTg           <  5:303734/139‑1 (MQ=255)
                                                                                                                                         aaTTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTCCCGCGGCGTTGCACCAAAGAAATTACGGTAGGCACTGGAAAAATGCGGCCCCGAGGAGAAGCCACAGGACAGCCCGATCTGga                         >  6:151989/1‑124 (MQ=255)
                                                                                                                                         aaTTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTCCCGCGGCGTTGCACCAAAGAAATTACGGTAGGCACTGGAAAAATGCGGCCCCGAGGAGAAGCCACAGGACAGCCCGATCTGga                         <  5:151989/124‑1 (MQ=255)
                                                                                                                                         aaTTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTCCCGCGGCGTTGCACCAAAGAAATTACGGTAGGCACTGGAAAAATGCGGCCCCGAGGAGAAGCCACAGGACAGCCCGATCTGGATGATCGACTTGCTGGTCTGCATc  <  8:398175/147‑1 (MQ=255)
                                                                                                                                         aaTTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTCCCGCGCCGTTGCACCAAAGAAATTACGGTAGGCACTGGAAAAATGCGGCCCCGAGGAGAAGCCACAGGACAGCCCCAACTGGATGATCGACTTGCTGGTCTGCAt   >  8:202267/1‑146 (MQ=255)
                                                                                                                                         |                                                                                                                                                  
GTACGGCGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTCCCGCGGCGTTGCACCAAAGAAATTACGGTAGGCACTGGAAAAATGCGGCCCCGAGGAGAAGCCACAGGACAGCCCGATCTGGATGATCGACTTGCTGGTCTGCATC  >  NC_002947/5092532‑5092815

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: ‑