Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 513,469 C→A 7.9% intergenic (‑116/‑75) crp ← / → yhfA DNA‑binding transcriptional dual regulator/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947513,4690CA 7.9% 217.5 / 4.9 70intergenic (‑116/‑75)crp/yhfADNA‑binding transcriptional dual regulator/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (31/33);  new base A (2/4);  total (33/37)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

TTTATCTTGGCGGGTAGGGCGGAGGCAACCATCACGAAATCCTGTTCGATATGACGCATAGGCTTTTTTAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGTAGATGAAGGCACGTATCCAGTGGGCCGGTGAGGCAATGTTCCTCGGCGAATCGGGGAGCGGCCATGTCGTCG  >  NC_002947/513322‑513613
                                                                                                                                                    |                                                                                                                                                
tttATCTTGGCGGGTAGGGCGGAGGCAACCATCACGAAATCCTGTTCGATATGACGCATAGGCTTTTTTAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGaca                                                                                                                                                 >  1:143523/1‑149 (MQ=255)
      ttGGCGGGTAGGGCGGAGGCAACCATCACGAAATCCTGTTCGATATGACGCATAGGCTTTTTTAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGac                                                                                                                                                  >  3:234887/1‑142 (MQ=255)
      ttGGCGGGTAGGGCGGAGGCAACCATCACGAAATCCTGTTCGATATGACGCATAGGCTTTTTTAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGac                                                                                                                                                  <  4:234887/142‑1 (MQ=255)
      ttGGCGGGTAGGGCGGAGGCAACCATCACGAAATCCTGTTCGATATGACGCATAGGCTTTTTTAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACAtctc                                                                                                                                           >  5:269331/1‑149 (MQ=255)
         gCGGGTAGGGCGGAGGCAACCATCACGAAATCCTGTTCGATATGACGCATAGGCTTTTTTAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTGTCGGTTTAGTACCGAAACATCTCGta                                                                                                                                        <  6:428311/149‑2 (MQ=38)
         gCGGGTAGGGCGGAGGCAACCATCACGAAATCCTGTTCGATATGACGCATAGGCTTTTTTAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTGTCGGTTTAGTACCGAAACATCTCGta                                                                                                                                        <  1:31714/149‑2 (MQ=38)
                gggCGGAGGCAACCATCACGAAATCCTGTTCGATATGACGCATAGGCTTTTTTAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTcc                                                                                                                                 <  3:92079/149‑1 (MQ=255)
                gggCGGAGGCAACCATCACGAAATCCTGTTCGATATGACGCATAGGCTTTTTTAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTc                                                                                                                                  <  4:381936/148‑1 (MQ=255)
                gggCGGAGGCAACCATCACGAAATCCTGTTCGATATGACGCATAGGCTTTTTTAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTc                                                                                                                                  <  1:387709/148‑1 (MQ=255)
                 ggCGGAGGCAACCATCACGAAATCCTGTTCGATATGACGCATAGGCTTTTTTAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTcc                                                                                                                                 >  5:332616/1‑148 (MQ=255)
                  gCGGAGGCAACCATCACGAAATCCTGTTCGATATGACGCATAGGCTTTTTTAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCca                                                                                                                               >  2:65608/1‑149 (MQ=255)
                     gAGGCAACCATCACGAAATCCTGTTCGATATGACGCATAGGCTTTTTTAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAg                                                                                                                            <  1:105736/149‑1 (MQ=255)
                          aaCCATCACGAAATCCTGTTCGATATGACGCATAGGCTTTTTTAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGca                                                                                                                       >  6:381581/1‑149 (MQ=255)
                             cATCACGAAATCCTGTTCGATATGACGCATAGGCTTTTTTAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACaa                                                                                                                    <  7:348853/149‑1 (MQ=255)
                             cATCACGAAATCCTGTTCGATATGACGCATAGGCTTTTTTAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACaa                                                                                                                    <  3:288439/149‑1 (MQ=255)
                              aTCACGAAATCCTGTTCGATATGACGCATAGGCTTTTTTAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAg                                                                                                                   <  6:332616/149‑1 (MQ=255)
                               tCACGAAATCCTGTTCGATATGACGCATAGGCTTTTTTAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAgg                                                                                                                  >  6:420289/1‑149 (MQ=255)
                                 aCGAAATCCTTTTCGATATGACGCATAGGCTTTTTTAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGAACACATCTCGTGCTTTTCCCACAGCCGCACAAGGt                                                                                                                 <  1:65608/149‑1 (MQ=255)
                                 aCGAAATCCTGTTCGATATTACGCATAGGCTTTTTTAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGt                                                                                                                 <  4:252841/148‑1 (MQ=255)
                                   gAAATCCTGTTCGATATGACGCATAGGCTTTTTTAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAg                                                                                                                   >  3:93292/1‑144 (MQ=255)
                                   gAAATCCTGTTCGATATGACGCATAGGCTTTTTTAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAg                                                                                                                   <  4:93292/144‑1 (MQ=255)
                                   gAAATCCTGTTCGATATGACGCATAGGCTTTTTTAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCt                                                                                                               >  4:44264/1‑148 (MQ=255)
                                     aaTCCTGTTCGATATGACGCATAGGCTTTTTTAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCt                                                                                                               >  4:11203/1‑146 (MQ=255)
                                     aaTCCTGTTCGATATGACGCATAGGCTTTTTTAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAg                                                                                                            <  1:422435/149‑1 (MQ=255)
                                       tCCTGTTCGATATGACGCATAGGCTTTTTTAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCt                                                                                                          <  4:210430/149‑1 (MQ=255)
                                       tCCTGTTCGATATGACGCATAGGCTTTTTTAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCt                                                                                                          <  4:324580/149‑1 (MQ=255)
                                         cTGTTCGATATGACGCATAGGCTTTTTTAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTg                                                                                                         <  3:11203/148‑1 (MQ=255)
                                         cTGTTCGATATGACGCATAGGCTTTTTTAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGa                                                                                                        <  7:165367/149‑1 (MQ=255)
                                          tGTTCGATATGACGCATAGGCTTTTTTAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGa                                                                                                        >  5:209123/1‑148 (MQ=255)
                                          tGTTCGATATGACGCATAGGCTTTTTTAAAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATCTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGa                                                                                                        >  1:25540/1‑148 (MQ=255)
                                                   tGACGCATAGGCTTTTTTAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGACACCCtttt                                                                                               <  6:166477/148‑1 (MQ=255)
                                                    gACGCATAGGCTTTTTTAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTcc                                                                                             >  7:108450/1‑149 (MQ=255)
                                                            ggCTTTTTTAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCa                                                                                     <  3:274604/149‑1 (MQ=255)
                                                             gCTTTTTTAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAg                                                                                    <  3:44264/149‑1 (MQ=255)
                                                               ttttttAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAg                                                                                    >  3:327683/1‑147 (MQ=255)
                                                                   ttAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTcc                                                                              >  7:22614/1‑149 (MQ=255)
                                                                      gagGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCggg                                                                           >  8:64342/1‑149 (MQ=255)
                                                                         gTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAACTCAGGAGTCCgcg                                                                           >  2:215073/1‑144 (MQ=255)
                                                                                   ggAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGTAGATGAAGGCAc                                                              >  7:123584/1‑149 (MQ=255)
                                                                                     agagCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTACGAAAATCAGGAGTCAGGGTAGATGAAGGCACg                                                             >  7:283574/1‑148 (MQ=255)
                                                                                      gagCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGTAGATGAAGGCACGTa                                                           >  6:265698/1‑149 (MQ=255)
                                                                                          gcCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGTAGATGAAGGCACGTATCCa                                                       <  1:155698/149‑1 (MQ=255)
                                                                                          gcCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAAACGAAACATATCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGACACACTTTTCCGAAAATCAGGAGTCAGGGTAGATGAAGGCACGTAta                                                         >  6:264590/1‑146 (MQ=255)
                                                                                           ccAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGTAGATGAAGGCACGTATCCAg                                                      >  2:48199/1‑149 (MQ=255)
                                                                                             aGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGTAGATGAAGGCACGTATCCAGTg                                                    <  8:22614/149‑1 (MQ=255)
                                                                                             aGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGTAGATGAAGGCACGTATCCAGTg                                                    >  3:188683/1‑149 (MQ=255)
                                                                                                  tttGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGTAGATGAAGGCACGTATCCAGTGGGCCg                                               <  5:381581/149‑1 (MQ=255)
                                                                                                     ggCGCCAGCTTATCAGACAATACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGTAGATGAAGGCACGTATCCAGTGGGCCGGt                                             <  4:204804/148‑1 (MQ=255)
                                                                                                         ccAGCTTATCCGACACTACGCGATGTGTCGGTACTAAACCTACACATCTCGTGCTGTTCCCACAGCCGCACAAGGTCTAAGCTTACACGCTTTTCCGAAAATCAGGAGTCCGGGTAGATGAAGGCACGTATCCAGTGGGCCGGTGAGGc                                        <  4:286313/149‑1 (MQ=255)
                                                                                                             cTTATCAGACACTACGAGATGTTTCGGTACTAAACCGAAACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCTGGTAGATGAAGGCACGTATCCAGTGGGCCGGTGAGGCAATg                                    <  3:128167/149‑1 (MQ=255)
                                                                                                                aTCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGTAGATGAAGGCACGTATCCAGTGGGCCGGTGAGGCAATGt                                   <  2:143523/147‑1 (MQ=255)
                                                                                                                 tCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGTAGATGAAGGCACGTATCCAGTGGGCCGGTGAGGCAATGTTc                                 <  2:294136/148‑1 (MQ=255)
                                                                                                                   aGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGTAGATGAAGGCACGTATCCAGTGGGCCGGTGAGGCAATGTTCCTc                              <  8:123584/149‑1 (MQ=255)
                                                                                                                       agcaCGAGATGTGTCGGTTTAGTACCGAAACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGTAGAt                                                                      >  2:104896/4‑105 (MQ=25)
                                                                                                                       agcaCGAGATGTGTCGGTTTAGTACCGAAACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGTAGAt                                                                      <  1:104896/102‑1 (MQ=25)
                                                                                                                         tACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCAGAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGTAGATGAAGGCACGTATCAAGTGGGCAGGTGAGGCAATGTTCCTCGGCGaa                        >  6:305123/1‑149 (MQ=255)
                                                                                                                         tACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGTAGATGAAGGCACGTATCCAGTGGGCCGGTGAGGCAATGTTCCTCGGCGac                        >  4:273527/1‑148 (MQ=255)
                                                                                                                         tACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGTAGATGAAGGCACGTATCCAGTGGGCCGGTGAGGCAATGTTCCTCGGCGaa                        >  2:164244/1‑149 (MQ=255)
                                                                                                                          aCGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAg                                                                                    <  2:254330/88‑1 (MQ=255)
                                                                                                                          aCGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAg                                                                                    >  1:254330/1‑88 (MQ=255)
                                                                                                                             agaTGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGTAGATGAAGGCACGTATCCAGTGGGCCGGTGAGGCAATGTTCCTCGGCGAATCgg                    <  1:164244/149‑1 (MQ=255)
                                                                                                                              taTGTTTCGGTACTAAAGCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGTAGATGAAGGCACGTATCCAGTGGGCCGGTGAGGCAATGTTCCTCGGCGAATCggg                   <  5:265698/148‑1 (MQ=255)
                                                                                                                                   ttCGGTACTAACCCGACACATCTCGTGATTTTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGCAGTCCGGGTACATGAAGGCACGTATCCATTGGGCCGATGATGCAATGTTACTCGGCGAATCGGGGatag              >  7:312023/1‑146 (MQ=255)
                                                                                                                                      ggTACTAAACCTACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCACGAGACCGGGTAGATCAAGGAACGTATCCAGTGGGCCGGTGAGGCAATGTTCCTCGTCCAATCGGGCAGCGGc            >  2:147923/1‑148 (MQ=255)
                                                                                                                                         aCTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGGAGATGAAGGCACGGATCCAGTGGGCCGGTGAGGCAATGTTCCTCGGCGAATCGGGGAGCGGCCATg        >  8:274745/1‑149 (MQ=255)
                                                                                                                                          cTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGTAGATGAAGGCACGTATCCAGTGGGCCGGTGAGGCAATGTTCCTCGGCGAATCGGGGAGCGGCCATg        <  3:215125/148‑1 (MQ=255)
                                                                                                                                           tAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGTAGATGAAGGCACGTATCCAGTGGGCCGGTGAGGCAATGTTc                                 <  6:305761/122‑1 (MQ=255)
                                                                                                                                           tAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGTAGATGAAGGCACGTATCCAGTGGGCCGGTGAGGCAATGTTc                                 >  5:305761/1‑122 (MQ=255)
                                                                                                                                            aaaCCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGTAGATGAAGGCACGTATCCAGTGGGCCGGTGAGGCAATGTTCCTCGGCGAATCGGGGAGCGGCCATgtcg     >  7:172239/1‑149 (MQ=255)
                                                                                                                                             aaCCGACACATGTCGTGCTTTTCCCACAGTCGCAGAATGTCTCTGCTGACACCGTTTTCCGAAAATCAGGAGTCCGGGTAGATGAAGGCACGTATCCAGTGGGCCGGTGAGGCAATGTTCCTCGGCGAATCGGGGAGCGGCCATgtcg     <  1:147923/148‑1 (MQ=255)
                                                                                                                                              aCCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGTAGATGAAGGCACGTATCCAGTGGGCCGGTGAGGCAATGTTCCTCGGCGAAACGGGGAGCGGCCATgtcgtc   <  5:420289/149‑1 (MQ=255)
                                                                                                                                                cGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGTAGATGAAGGCACGTATCCAGTGGGCCGGTGAGGCAATGTTCCTCGGCGAATCGGGGAGCGGCCATGtcgtcg  >  5:211212/1‑148 (MQ=255)
                                                                                                                                                cGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGTAGATGAAGGCACGTATCCAGTGGGCCGGTGAGGCAATGTTCCTCGGAGCATCGGGGAGCGGCCATGtcgtcg  >  8:231599/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
TTTATCTTGGCGGGTAGGGCGGAGGCAACCATCACGAAATCCTGTTCGATATGACGCATAGGCTTTTTTAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGTAGATGAAGGCACGTATCCAGTGGGCCGGTGAGGCAATGTTCCTCGGCGAATCGGGGAGCGGCCATGTCGTCG  >  NC_002947/513322‑513613

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: