Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,380,163 G→T 5.9% 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 5.9% 217.2 / 3.6 68intergenic (+74/+34)recN/furDNA repair/recombination protein/ferric iron uptake global transcriptional repressor Fur
Reads supporting (aligned to +/- strand):  ref base G (30/34);  new base T (2/2);  total (32/36)
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.89e-01

GGCGGGATCGACCTGACCAAGGAATCCCTGGCCCATGCACGCAAGATGGTGGTGAGCGGCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGCTCGGCTACGATTTCGCGCTGGCGCTTCTCGATTTCGGCATCCATGAATTCGATGACT  >  NC_002947/5380021‑5380311
                                                                                                                                              |                                                                                                                                                    
ggcggGATCGACCTGACCAAGGAATCCCTGGCCCATGCACGCAAGATGGTGGTGAGCGGCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTc                                                                                                                                                <  7:52695/149‑1 (MQ=255)
         gACCTGACCAAGGAATCCCTGGCCCATGCACGCAAGATGGTGGTGAGCGGCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCTTCGGGTCGCCTTTgt                                                                                                                                        <  4:1503/148‑2 (MQ=255)
         gACCTGACCAAGGAATCCCTGGCCCATGCACGCAAGATGGTGGTGAGCGGCAAGGCCTGAGCCTTGCGTTTCCAGGACGATTCGCGGGCACAGACAACACATCGTTTAGACACAGAAAGCACAAAGGAGACCCGCAGGTCGCCatttc                                                                                                                                        >  7:140053/1‑146 (MQ=255)
         gACCTGACCAAGGAATCCCTGGCCCATGCACCCAACAGGGTGGTGAGCGGCAAGGCCTGAGCCTTGCTTATCCAGGACGATTCGCGGGCACAGCCAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGCCCCGAAGGACGCCTTTGc                                                                                                                                        >  7:145065/1‑148 (MQ=255)
          aCCTGACCAAGGAATCCCTGGCCCATGCACGCAAGATGGTGGTGAGCGGCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTa                                                                                                                                      >  4:330415/1‑149 (MQ=255)
            cTGACCAAGGAATCCCTGGCCCATGCACGCAAGATGGTGGTGAGCGGCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATc                                                                                                                                    <  5:5059/149‑1 (MQ=255)
                 cAAGGAATCCCTGGCCCATGCACGCTAGATGGTGGTTAGCGGCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATaa                                                                                                                               <  4:283432/149‑1 (MQ=255)
                 cAAGGAATCCCTGGCCCATGCACGCAAGATGGTGGTTAGCGGCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATaa                                                                                                                               <  5:245645/149‑1 (MQ=255)
                 cAAGGAATCCCTGGCCCATGCACGCAAGATGGTGGTGAGCGGCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTg                                                                                                                                         <  6:278804/139‑1 (MQ=255)
                 cAAGGAATCCCTGGCCCATGCACGCAAGATGGTGGTGAGCGGCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTg                                                                                                                                         >  5:278804/1‑139 (MQ=255)
                 cAAGGAATCCCTGGCCCATGCACGCAAGATGGTGGTGAGCGGCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATaa                                                                                                                               <  5:123176/149‑1 (MQ=255)
                 cAAGGAATCCCTGGCCCATGCACGCAAGATGGTGGTGAGCGGCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATaa                                                                                                                               >  2:183457/1‑149 (MQ=255)
                  aaGGAATCCCTGGCCCATGCACGCAAGATGGTGGTGAGCGGCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAAc                                                                                                                              >  7:330260/1‑149 (MQ=255)
                      aaTCCCTGGCCCATGCACGCAAGATGGTGGTGAGCGGCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGaaa                                                                                                                          <  8:151245/149‑1 (MQ=255)
                      aaTCCCTGGCCCATGCACGCAAGATGGTGGTGAGCGGCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACAGCGTTTAGCCCCAGAAAGCACAAAGGCGCCCGGAAGGGCGCGGGTGGTAGCGATAACGaa                                                                                                                           >  2:331427/1‑148 (MQ=255)
                              gCCCATGCACGCAAGATGGTGGTGAGCGGCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGt                                                                                                                    <  3:317578/147‑1 (MQ=255)
                              gCCCATGCACGCAAGATGGTGGTGAGCGGCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTa                                                                                                                  >  2:323218/1‑149 (MQ=255)
                                ccATGCACGCAAGATGGTCGTGAGCGGCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTAct                                                                                                                >  7:320098/1‑149 (MQ=255)
                                    gCACGCAAGATGGTGGTGAGCGGCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGCCAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTActtctt                                                                                                            >  6:260440/1‑149 (MQ=255)
                                     cACGCAAGATGGTGGTGAGCGGCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTActtcttc                                                                                                           <  1:13470/149‑1 (MQ=255)
                                       cGCAAGATGGTGGTGAGCGGCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTActtcttctt                                                                                                         <  7:217020/149‑1 (MQ=255)
                                            gATGGTGGTGAGCGGCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCa                                                                                                    <  2:21680/149‑1 (MQ=255)
                                                 tggtgAGCGGCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTAc                                                                                               <  6:168599/149‑1 (MQ=255)
                                                    tgAGCGGCAAGGCCTGAGACTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTAGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAg                                                                                             >  5:364469/1‑148 (MQ=255)
                                                      aGCGGCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGa                                                                                           <  8:320098/148‑1 (MQ=255)
                                                        cGGCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCa                                                                                        <  6:5505/149‑1 (MQ=255)
                                                            aaGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGAtt                                                                                    >  4:187678/1‑149 (MQ=255)
                                                               tccTGAGCCTTGTTTTTCCAGGACGCTTCGCGGGCACAGACAACACATCTTTTAGCCCCATCAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAACAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATtgtg                                                                                 <  8:140053/148‑1 (MQ=255)
                                                               gCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATtgtg                                                                                 <  1:183457/149‑1 (MQ=255)
                                                               gCCTGAGCATTGCTTTTCCAGGACGATGCGCGGGCACAGACAACACATCTTTTAGCCCAAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATtgtg                                                                                 <  8:145065/149‑1 (MQ=255)
                                                                ctcaggCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTgg                                                                                <  8:111520/144‑1 (MQ=255)
                                                                 cTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCTTCGGGTCGCCTTTgt                                                                                                                                        >  7:275004/1‑91 (MQ=38)
                                                                 cTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCTTCGGGTCGCCTTTgt                                                                                                                                        <  8:275004/92‑2 (MQ=38)
                                                                      ccTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCACAATCTGGTcctgt                                                                          <  7:75706/149‑4 (MQ=255)
                                                                            ttttCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTc                                                                     >  4:184740/1‑148 (MQ=255)
                                                                            ttttCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCg                                                                    >  7:133154/1‑149 (MQ=255)
                                                                                 cAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGc                                                                >  7:200887/1‑148 (MQ=255)
                                                                                  aGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCa                                                              >  1:186404/1‑149 (MQ=255)
                                                                                   ggACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGAGCAGATTGTGGTCCACCAGCTCGAAGCCAt                                                             >  5:97447/1‑149 (MQ=255)
                                                                                    gACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCAt                                                             >  6:44448/1‑148 (MQ=255)
                                                                                     aCGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGc                                                           <  1:331427/149‑1 (MQ=255)
                                                                                     aCGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGc                                                           <  6:38550/149‑1 (MQ=255)
                                                                                     aCGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGc                                                           <  3:330415/149‑1 (MQ=255)
                                                                                                   caGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGCTCGGCTACGATTTc                                             <  8:330260/149‑1 (MQ=255)
                                                                                                        aacaCATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGCTCGGCTACGATTTCGCGCt                                        >  3:108901/1‑149 (MQ=255)
                                                                                                        aacaCATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGCTCGGCTACGATTTCGCGCt                                        >  1:214587/1‑149 (MQ=255)
                                                                                                        aacaCATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGCTCGGCTACGATTTCGCGCt                                        >  5:277588/1‑149 (MQ=255)
                                                                                                                  ttAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGCTCGGCTACGATTTCGCGCTGGCGCtt                                 <  6:97447/146‑1 (MQ=255)
                                                                                                                   tAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGCTCGGCTACGATTTCGCGCTGGCGCTtctc                              <  2:214587/148‑1 (MQ=255)
                                                                                                                       cccAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGCTCGGCTACGATTTCGCGCTGGCGCTTCTCGAttt                         <  6:277588/149‑1 (MQ=255)
                                                                                                                        ccAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGCTCGGCTACGATTTCGCGCTGGCGCTTCTCGATTTc                        <  1:148371/149‑1 (MQ=255)
                                                                                                                        ccAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGCTCGGCTACGATTTCGCGCTGGCGCTTCTCGATTTc                        >  6:172650/1‑149 (MQ=255)
                                                                                                                           gAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGCTCGGCTACGATTTCGCGCTGGCGCTTCTCGATTTCGGc                     <  3:187678/149‑1 (MQ=255)
                                                                                                                                 gcaAAGGCGACCTTCGGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGCTCGGCTACGATTTCGCGCTGGCGCTTCTCGATTTCGGCATCCAt               >  5:34089/2‑149 (MQ=255)
                                                                                                                                 acaAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGCTCGGCTACGATTTCGCGCTGGCGCTTCTCGATTTCGGCATCCAt               >  5:379136/1‑149 (MQ=255)
                                                                                                                                 acaAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGCTCGGCTACGATTTCGCGCTGGCGCTTCTCGATTTCGGCATCCAt               <  1:323218/149‑1 (MQ=255)
                                                                                                                                 acaAAGGCGACACGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGCTCGGCTACGATTTCGCGCTGGAGCTTCTCGATTTCGGCATCCAt               >  2:131485/1‑149 (MQ=255)
                                                                                                                                   aaaGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGCTCGGCTACGATTTCGCGCTGGCGCTTCTCGATTTCGGCATCCATg              >  4:247075/1‑148 (MQ=255)
                                                                                                                                      ggCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGCTCGGCTACGATTTCGCGCTGGCGCTTCTCGATTTCGGCATCCATGAAtt          >  5:227113/1‑149 (MQ=255)
                                                                                                                                       gCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGCTCGGCTACGATTTCGCGCTGGCGCTTCTCGATTTCGGCATCCATGAATTc         <  5:260440/149‑1 (MQ=255)
                                                                                                                                        cGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGCTCGGCTACGATTTCGCGCTGGCGCTTCTCGATTTCGGCATCCATGAATTc         >  6:322232/1‑148 (MQ=255)
                                                                                                                                          aCCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGCTCGGCTACGATTTCGCGCTGGCGCTTCTCGATTTCGGCATCCATGAATTCGAt      <  5:44448/149‑1 (MQ=255)
                                                                                                                                          aCCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGCTCGGCTACGATTTCGCGCTGGCGCTTCTCGATTTCGGCATCCATGAATTCGAt      <  6:227113/149‑1 (MQ=255)
                                                                                                                                           cccGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGCTCGGCTACGATTTCGCGCTGGCGCTTCTCGATTTCGGCATCCATGAATTCGATg     >  8:246966/1‑149 (MQ=255)
                                                                                                                                            ccGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGCTCGGCTACGATTTCGCGCTGGCGCTTCTCGATTTCGGCATCCATGAATTCGATg     >  7:189840/1‑148 (MQ=255)
                                                                                                                                            ccGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGCTCGGCTACGATTTCGCGCTGGCGCTTCTCGATTTCGGCATCCATGAATTCGATGa    <  3:184740/149‑1 (MQ=255)
                                                                                                                                            ccGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGCTCGGCTACGATTTCGCGCTGGCGCTTCTCGATTTCGGCATCCATGAATTCGATGa    <  2:300719/149‑1 (MQ=255)
                                                                                                                                             cGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGCTCGGCTACGATTTCGCGCTGGCGCTTCTCGATTTCGGCATCCATGAATTCGATGAc   >  8:56506/1‑149 (MQ=255)
                                                                                                                                              gAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGCTCGGCTACGATTTCGCGCTGGCGCTTCTCGATTTCGGCATCCATGAATTCGATGACt  >  7:335863/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
GGCGGGATCGACCTGACCAAGGAATCCCTGGCCCATGCACGCAAGATGGTGGTGAGCGGCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGCTCGGCTACGATTTCGCGCTGGCGCTTCTCGATTTCGGCATCCATGAATTCGATGACT  >  NC_002947/5380021‑5380311

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: