Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,380,163 G→T 8.6% intergenic (+74/+34) recN → / ← fur DNA repair/recombination protein/ferric iron uptake global transcriptional repressor Fur

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,380,1630GT 8.6% 214.6 / 8.6 70intergenic (+74/+34)recN/furDNA repair/recombination protein/ferric iron uptake global transcriptional repressor Fur
Reads supporting (aligned to +/- strand):  ref base G (33/31);  new base T (2/4);  total (35/35)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.47e-01

GGATCGACCTGACCAAGGAATCCCTGGCCCATGCACGCAAGATGGTGGTGAGCGGCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGCTCGGCTACGATTTCGCGCTGGCGCTTCTCGATTTCGGCATCCATGAATTCGATGACT  >  NC_002947/5380025‑5380311
                                                                                                                                          |                                                                                                                                                     
ggATCGACCTGACCAAGGAATCCCTGGCCCATGCACGCAAGATGGTGGTGAGCGGCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCTTCGGGTCGCCt                                                                                                                                             <  2:22708/149‑1 (MQ=255)
  aTCGACCTGACCAAGGAATCCCTGGCCCATGCACGCAAGATGGTGGTGAGCGGCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCTTCGGGTCGCCttt                                                                                                                                           <  6:334469/149‑1 (MQ=255)
   tCGACCTGACCAAGGAATCCCTGGCCCATGCACGCAAGATGGTGGTGAGCGGCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCttt                                                                                                                                           <  4:176822/148‑1 (MQ=255)
     gACCTGACCAAGGAATCCCTGGCCCATGCACGCAAGATGGTGGTGAGCGGCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCTTCGGGTCGCCTTTgtg                                                                                                                                        <  7:38028/149‑3 (MQ=255)
      aCCTGACCAAGGAATCCCTGGCCCATGCACGCAAGATGGTGGTGAGCGGCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCt                                                                                                                                        >  8:197424/1‑148 (MQ=255)
      aCCTGACCAAGGAATCCCTGGCCCATGCACGCAAGATGGTGGTGAGCGGCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTa                                                                                                                                       >  7:109254/1‑149 (MQ=255)
       ccTGACCAAGGAATCCCTGGCCCATGCACGCAAGATGGTGGTGAGCGGCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTAt                                                                                                                                      >  4:32498/1‑149 (MQ=255)
           aCCAAGGAATCCCTGGCCCATGCACGCAAGATGGTGGTGAGCGGCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATcc                                                                                                                                    <  3:264025/147‑1 (MQ=255)
             cAAGGAATCCCTGGCCCATGCACGCAAGATGGTGGTGAGCGGCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATaa                                                                                                                                >  5:122191/1‑149 (MQ=255)
             cAAGGAATCCCTGGCCCATGCACGCAAGATGGTGGTGAGCGGCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATaa                                                                                                                                >  6:227018/1‑149 (MQ=255)
               aGGAATCCCTGGCCCATGCACGCAAGATGGTGGTGAGCGGCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAAc                                                                                                                               >  6:219475/1‑148 (MQ=255)
                      ccTGGCCCATGCACGCAAGATGGTGGTGAGCGGCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAAcc                                                                                                                       <  7:197424/149‑1 (MQ=255)
                       cTGGCCCATGCACGCAAGATGGTGGTGAGCGGCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCg                                                                                                                      >  7:252778/1‑149 (MQ=255)
                       cTGGCCCATGCACGCAAGATGGTGGTGAGCGGCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCg                                                                                                                      >  3:131221/1‑149 (MQ=255)
                           cccATGCACGCAAGATGGTGGTGAGCGGCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTAc                                                                                                                  <  1:64146/149‑1 (MQ=255)
                                 cACGCAAGATGGTGGTGAGCGGCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTActtcttc                                                                                                            <  5:104940/149‑1 (MQ=255)
                                        gATGGTGGTGAGCGGCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCa                                                                                                     >  2:194393/1‑149 (MQ=255)
                                          tggtggTGAGCGGCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACg                                                                                                   <  8:252778/149‑1 (MQ=255)
                                               gtgAGCGGCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGATACTTCTTCTTACGCACGTACAg                                                                                              >  6:319362/1‑149 (MQ=255)
                                                 gAGCGGCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGa                                                                                            >  3:167162/1‑149 (MQ=255)
                                                  aGCTGCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATACCGAAAAACCGTTACTTCTTCTTACGCACGTACAGGa                                                                                            <  6:324821/148‑1 (MQ=255)
                                                  aGCGGCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGAc                                                                                           <  1:194393/149‑1 (MQ=255)
                                                      gCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGa                                                                                       <  6:122191/149‑1 (MQ=255)
                                                       cAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGAt                                                                                      >  2:143977/1‑149 (MQ=255)
                                                            ccTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATtgtg                                                                                  <  7:75280/148‑1 (MQ=255)
                                                             cTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTActtcttc                                                                                                            <  2:321459/121‑1 (MQ=255)
                                                             cTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTActtcttc                                                                                                            >  1:321459/1‑121 (MQ=255)
                                                              tGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTc                                                                               >  3:124308/1‑149 (MQ=255)
                                                                 gCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTccac                                                                            <  4:131221/149‑1 (MQ=255)
                                                                   cTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTccac                                                                            >  4:50328/1‑147 (MQ=255)
                                                                   cTTGCTTTTCCAGGACGATTAGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATACATAACGAAAAACCCTTACTTCTTCTTACGCACGTACACGACCAGATTGTGGTccacca                                                                          >  7:225161/1‑149 (MQ=255)
                                                                    ttGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCGATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCCCCAg                                                                         >  4:168603/1‑149 (MQ=255)
                                                                      gCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTACCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATACCGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCACAATCTGGTcctgtac                                                                         <  6:95065/147‑6 (MQ=255)
                                                                       cTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCt                                                                       <  5:172625/148‑1 (MQ=255)
                                                                           tCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGaa                                                                   <  5:219475/148‑1 (MQ=255)
                                                                             cAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGc                                                                 >  1:164734/1‑148 (MQ=255)
                                                                              aGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCa                                                              >  7:243728/1‑149 (MQ=255)
                                                                                gACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCAc                                                                                                    >  7:75943/1‑110 (MQ=255)
                                                                                gACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCAc                                                                                                    <  8:75943/110‑1 (MQ=255)
                                                                                gACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATg                                                            >  3:273569/1‑149 (MQ=255)
                                                                                     ttCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGCTCg                                                        <  4:167162/148‑1 (MQ=255)
                                                                                        gcgGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGCTCGGCTa                                                    <  3:32498/149‑1 (MQ=255)
                                                                                        gcgGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCATTTGCTAACCAATACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCACGAAGCCATGCTCGGCTa                                                    >  5:228971/1‑149 (MQ=255)
                                                                                                  acaacaCATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGCTCGGCTACGATTTcgcg                                          >  5:291535/1‑149 (MQ=255)
                                                                                                   caacaCATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCAt                                                             >  7:292417/1‑129 (MQ=255)
                                                                                                   caacaCATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCAt                                                             <  8:292417/129‑1 (MQ=255)
                                                                                                      cacaTCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGCTCGGCTACGATTTCGCGCTgg                                      >  3:69779/1‑149 (MQ=255)
                                                                                                       acaTCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGCTCGGCTACGATTTCGCGCTGgc                                     >  1:166940/1‑149 (MQ=255)
                                                                                                       acaTCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGCTCGGCTACGATTTCGCGCTGgc                                     >  6:51202/1‑149 (MQ=255)
                                                                                                         aTCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCACAATCTGGTcctgtac                                                                         >  2:289493/1‑107 (MQ=255)
                                                                                                              ttAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGCTCGGCTACGATTTCGCGCTGGCGCTtc                                <  5:227018/147‑1 (MQ=255)
                                                                                                                 gCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGCTCGGCTACGATTTCGCGCTGGCGCTTCTCGAt                           <  2:166940/149‑1 (MQ=255)
                                                                                                                    ccAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGCTCGGCTACGATTTCGCGCTGGCGCTTCTCGATTTc                        >  5:299354/1‑149 (MQ=255)
                                                                                                                     cAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGCTCGGCTACGATTTCGCGCTGGCGCTTCTCGATTTCg                       <  3:50328/149‑1 (MQ=255)
                                                                                                                        aaaGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCa                                                                                                                                   <  8:171263/39‑1 (MQ=255)
                                                                                                                        aaaGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCa                                                                                                                                   >  7:171263/1‑39 (MQ=255)
                                                                                                                        aaaGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGCTCGGCTACGATTTCGCGCTGGCGCTTCTCGATTTCGGCa                    <  4:273569/149‑1 (MQ=255)
                                                                                                                           tagcaAAGGCGACCTTCGGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGCTCGGCTACGATTTCGCGCTGGCGCTTCTCGATTTCGGCATc                  >  7:11537/4‑148 (MQ=255)
                                                                                                                              caAAGGCGACCTTCGGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGCTCGGCTACGATTTcgcg                                          <  8:133623/121‑1 (MQ=255)
                                                                                                                              caAAGGCGACCTTCGGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGCTCGGCTACGATTTcgcg                                          >  7:133623/1‑121 (MQ=255)
                                                                                                                               aaaGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGCTCGGCTACGATTTCGCGCTGGCGCTTCTCGATTTCGGCATCCATg              >  6:133788/1‑148 (MQ=255)
                                                                                                                                aaGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGCTCGGCTACGATTTCGCGCTGGCGCTtc                                >  2:40528/1‑129 (MQ=255)
                                                                                                                                aaGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGCTCGGCTACGATTTCGCGCTGGCGCTTCTCGATTTCGGCATCCATGaa            <  4:124308/149‑1 (MQ=255)
                                                                                                                                aaGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCAACAGCTCGAAGCACATGCTCGGCTACGATTTCGCGCTGGCGCTtc                                <  1:40528/130‑1 (MQ=255)
                                                                                                                                  ggCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGCTCGGCTACGATTTCGCGCTGGCGCTTCTCGATTTCGGCATCCATGAAtt          <  6:299354/149‑1 (MQ=255)
                                                                                                                                   gCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGCTCGGCTACGATTTCGCGCTGGCGCTTCTCGATTTCGGCATCCATGAATTc         >  2:111356/1‑149 (MQ=255)
                                                                                                                                     gACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGCTCGGCTACGATTTCGCGCTGGCGCTTCTCGATTTCGGCATCCATGAATTCGa       <  4:69779/149‑1 (MQ=255)
                                                                                                                                     gACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGCTCGGCTACGATTTCGCGCTGGCGCTTCTCGATTTCGGCATCCATGAATTCGa       <  5:133788/149‑1 (MQ=255)
                                                                                                                                     gACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGCTCGGCTACGATTTCGCGCTGGCGCTTCTCGATTTCGGCATCCATGAATTCGa       <  1:111356/149‑1 (MQ=255)
                                                                                                                                         cGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGCTCGGCTACGATTTCGCGCTGGCGCTTCTCGATTTCGGCATCCATGAATTCGATGAc   <  5:51202/149‑1 (MQ=255)
                                                                                                                                          gAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGCTCGGCTACGATTTCGCGCTGGCGCTTCTCGATTTCGGCATCCATGAATTCGATGACt  <  2:164734/149‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                                     
GGATCGACCTGACCAAGGAATCCCTGGCCCATGCACGCAAGATGGTGGTGAGCGGCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGCTCGGCTACGATTTCGCGCTGGCGCTTCTCGATTTCGGCATCCATGAATTCGATGACT  >  NC_002947/5380025‑5380311

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