Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,614,794 T→C 6.8% F389L (TTC→CTC)  PP_4085 → rhs family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,614,7940TC 6.8% 165.8 / 2.6 59F389L (TTC→CTC) PP_4085rhs family protein
Reads supporting (aligned to +/- strand):  ref base T (29/26);  new base C (2/2);  total (31/28)
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.40e-01

TCCCCTATGAAGTGGAGCTGGCCCGCGTCCCACATCCTCAAGGCGGTACGCTGTGGATCTACATCGACAACGACGGTAACCGACTGGAACTGGGCCGCCTCAGTGCCGGTGATGCCTTCATCAGTGCCGTGGACGGCTTGGCCTTCTTCCAGCTGGAGGACGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTATCGGGTCTTCGATACCGACCCGCTGAATCCGAAACGCTCACGGCTGATCAGGCTCGGCGACCGCAATCTCAATTGCCTCGACCTGCTGTACGAC  >  NC_002947/4614648‑4614940
                                                                                                                                                  |                                                                                                                                                  
tCCCCTATGAAGTGGAGCTGGCCCGCGTCCCACATCCTCAAGGCGGTACGCTGTGGATCTACATCGACAACGACGGTAACCGACTGGAACTGGGCCGCCTCAGTGCCGGTGATGCCTTCATCAGTGCCGTGGACGGCTTGGCCttcttc                                                                                                                                                  <  3:301897/149‑1 (MQ=255)
  cccTATGAAGTGGAGCTGGCCCGCGTCCCACATCCTCAAGGCGGTACGCTGTGGATCTACATCGACAACGACGGTAACCGACTGGAACTGGGCCGCCTCAGTGCCGGTGATGCCTTCATCAGTGCCGTGGACGGCTTGGCCTTCTTcc                                                                                                                                                 <  2:210873/148‑1 (MQ=255)
   ccTATGAAGTGGAGCTGGCCCGCGTCCCACATCCTCAAGGCGGTACGCTGTGGATCTACATCGACAACGACGGTAACCGACTGGAACTGGGCCGCCTCAGTGCTGGTGATGCCTTCATCAGTGCCGTGGACGGCTTGGCCTTCTTCCAg                                                                                                                                               >  5:363232/1‑149 (MQ=255)
   ccTATGAAGTGGAGCTGGCCCGCGTCCCACATCCTCAAGGCGGTACGCTGTGGATCTACATCGACAACGACGGTAACCGACTGGAACTGGGCCGCCTCAGTGCCGGTGATGCCTTCATCAGTGCCGTGGACGGCTTGGCCTTCTTCCAg                                                                                                                                               >  2:242423/1‑149 (MQ=255)
         aaGTGGAGCTGGCCCGCGTCCCACATCCTCAAGGCGGTACGCTGTGGATCTACATCGACAACGACGGTAACCGACTGGAACTGGGCCGCCTCAGTGCCGGTGATGCCTTCATCAGTGCCGTGGACGGCTTGGCCTTCTTCCAGCTgga                                                                                                                                          <  4:147825/148‑1 (MQ=255)
          aGTGGAGCTGGCCCGCGTCCCACATCCTCAAGGCGGTACGCTGTGGATCTACATCGACAACGACGGTAAACGACTGGAACTGGGCCGCCTCAGTGCCGGTGATGCCTTCATCAGTGCCGTGGACGGCTTGGCCTTCTTCCAGCTggag                                                                                                                                         >  6:93019/1‑148 (MQ=255)
              gAGCTGGCCCGCGTCCCACATCCTCAAGGCGGTACGCTGTGGATCTACATCGACAACGACGGTAACCGACTGGAACTGGGCCGCCTCAGTGCCGGTGATGCCTTCATCAGTGCCGTGGACGGCTTGGCCTTCTTCCAGCTGGAGGACg                                                                                                                                     <  7:222609/148‑1 (MQ=255)
              gAGCTGGCCCGCGTCCCACATCCTCAAGGCGGTACGCTGTGGATCTACATCGACAACGACGGTAACCGACTGGAACTGGGCCGCCTCAGTGCCGGTGATGCCTTCATCAGTGCCGTGGACGGCTTGGCCTTCTTCCAGCTGGAGGACg                                                                                                                                     <  4:131824/148‑1 (MQ=255)
              gAGCTGGCCCGCGTCCCACATCCTCAAGGCGGTACGCTGTGGATCTACATCGACAACGACGGTAACCGACTGGAACTGGGCCGCCTCAGTGCCGGTGATGCCTTCATCAGTGCCGTGGACGGCTTGGCCTTCTTCCAGCTGGAGGACg                                                                                                                                     >  3:131824/1‑148 (MQ=255)
               aGCTGGCCCGCGTCCCACATCCTCAAGGCGGTACGCTGTGGATCTACATCGACAACGACGGTAACCGACTGGAACTGGGCCGCCTCAGTGCTGGTGATGCCTTCATCAGTGCCGTGGACGGCTTGGCCTTCTTCCAGCTGGAGGACGGc                                                                                                                                   <  6:363232/149‑1 (MQ=255)
                gCTGGCCCGCGTCCCACATCCTCAAGGCGGTACGCTGTGGATCTACATCGACAACGACGGTAACCGCCTGGAACTGGGCCGCCTCAGTGCCGGTGATGCCTTCATCACTGCCGTGGACGGCTTGGCCTTCTTCCAGCTGGAGGACGGc                                                                                                                                   >  7:279111/1‑148 (MQ=255)
                 cTGGCCCGCGTCCCACATCCTCAAGGCGGTACGCTGTGGATCTACATCGACAACGACGGTAACCGACTGGAACTGGGCCGCCTCAGTGCCGGTGATGCCTTCATCAGTGCCGTGGACGGCTTGGCCTTCTTCCAGCTGGAGGACGGCCa                                                                                                                                 <  1:425497/149‑1 (MQ=255)
                            cccACATCCTCAAGGCGGTACGCTGTGGATCTACATCGACAACGACGGTAACCGACTGGAACTGGGCCGCCTCAGTGCCGGTGATGCCTTCATCAGTGCCGTGGACGGCTTGGCCTTCTTCCAGCTGGAGGACGGCCAGAc                                                                                                                              >  5:392247/1‑141 (MQ=255)
                            cccACATCCTCAAGGCGGTACGCTGTGGATCTACATCGACAACGACGGTAACCGACTGGAACTGGGCCGCCTCAGTGCCGGTGATGCCTTCATCAGTGCCGTGGACGGCTTGGCCTTCTTCCAGCTGGAGGACGGCCAGAc                                                                                                                              <  6:392247/141‑1 (MQ=255)
                              cacaTCCTCAAGGCGGTACGCTGTGGATCTACATCGACAACGACGGTAACCGACTGGAACTGGGCCGCCTCAGTGCCGGTGATGCCTTCATCAGTGCCGTGGACGGCTTGGCCTTCTTCCAGCTGGAGGACGGCCAGACCGTTGTCg                                                                                                                      >  6:111934/1‑147 (MQ=255)
                               acaTCCTCAAGGCGGTACGCTGTGGATCTACATCGACAACGACGGTAACCGACTGGAACTGGGCCGCCTCAGTCCCGGTGATGCCTTCATCAGTGCCGTTGACGGCTTGGCCTTCTTCCAGCTGGAGGACGGCCAGACCGTTGTCGAAg                                                                                                                   <  3:391720/149‑1 (MQ=255)
                                        aGGCGGTACGCTGTGGATCTACATCGACAACGACGGTAACCGACTGGAACTGGGCCGCCTCAGTGCCGGTGATGCCTTCATCAGTGCCGTGGACGGCTTGGCCTTCTTCCAGCTGGAGGACGGCCAGACCGTTGTCGAAGATATCAAt                                                                                                           >  2:264400/1‑148 (MQ=255)
                                           cGGTACGCTGTGGATCTACATCGACAACGACGGTAACCGACTGGAACTGGGCCGCCTCAGTGCCGGTGATGCCTTCATCAGTGCCGTGGACGGCTTGGCCTTCTTCCAGCTGGAGGACgg                                                                                                                                    <  7:363774/120‑1 (MQ=255)
                                           cGGTACGCTGTGGATCTACATCGACAACGACGGTAACCGACTGGAACTGGGCCGCCTCAGTGCCGGTGATGCCTTCATCAGTGCCGTGGACGGCTTGGCCTTCTTCCAGCTGGAGGACgg                                                                                                                                    >  8:363774/1‑120 (MQ=255)
                                             gTACGCTGTGGATCTACATCGACAACGACGGTAACCGACTGGAACTGGGCCGCCTCAGTGCCGGTGATGCCTTCATCAGTGCCGTGGACGGCTTGGCCTTCTTCCAGCTGGAGGACGGCCAGACCGTTGTCGAAGATATCAATGAAggg                                                                                                     <  4:50254/149‑1 (MQ=255)
                                             gTACGCTGTGGATCTACATCGACAACGACGGTAACCGACTGGAACTGGGCCGCCTCAGTGCCGGTGATGCCTTCATCAGTGCCGTGGACGGCTTGGCCTTCTTCCAGCTGGAGGACGGCCAGACCGTTGTCGAAGATATCAATGAAggg                                                                                                     <  2:162113/149‑1 (MQ=255)
                                                   tgtgGATCTACATCGACAACGACGGTAACCGACTGGAACTGGGCCGCCTCAGTGCCGGTGATGCCTTCATCAGTGCCGTGGACGGCTTGGCCTTCTTCCAGCTGGAGGACGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTAt                                                                                               >  4:198939/1‑149 (MQ=255)
                                                    gtgGATCTACATCGACAACGACGGTAACCGACTGGAACTGGGCCGCCTCAGTGCCGGTGATGCCTTCATCAGTGCCGTGGACGGCTTGGCCTTCTTCCAGCTGGAGGACGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTATc                                                                                              >  8:298247/1‑149 (MQ=255)
                                                     tgGATCTACATCGACAACGACGGTAACCGACTGGAACTGGGCCGCCTCAGTGCCGGTGATGCCTTCATCAGTGCCGTGGACGGCTTGGCCTTCTTCCAGCTGGAGGACGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTATc                                                                                              <  7:167833/148‑1 (MQ=255)
                                                           tACATCGACAACGACGGTAACCGACTGGAACTGGGCCGCCTCAGTGCCGGTGATGCCTTCATCAGTGCCGTGGACGGCTTGGCCTTCTTCCAGCTGGAGGACGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTATCGGGTCtt                                                                                       >  3:195631/1‑149 (MQ=255)
                                                              aTCGACAACGACGGTAACCGACTGGAACTGGGCCGCCTCAGTGCCGGTGATGCCTTCATCAGTGCCGTGGACGGCTTGGCCTTCTTCCAGCTGGAGGACGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTATCGGGTCTTCGa                                                                                    >  5:48333/1‑149 (MQ=255)
                                                                cGACAACGACGGTAACCGACTGGAACTGGGCCGCCTCAGTGCCGGTGATGCCTTCATCAGTGCCGTGGACGGCTTGGCCTTCTTCCAGCTGGAGGACGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTATCGGGTCTTCGATa                                                                                  >  7:142933/1‑149 (MQ=255)
                                                                     acgacgGTAACCGACTGGAACTGGGCCGCCTCAGTGCCGGTGATGCCTTCATCAGTGCCGTGGACGGCTTGGCCTTCTTCCAGCTGGAGGACGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTATCGGGTCTTCGATACCGAc                                                                             <  7:298247/149‑1 (MQ=255)
                                                                     acgacgGTAACCGACTGGAACTGGGCCGCCTCAGTGCCGGTGATGCCTTCATCAGTGCCGTGGACGGCTTGGCCTTCTTCCAGCTGGAGGACGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTATCGGGTCTTCGATACCGAc                                                                             >  8:56936/1‑149 (MQ=255)
                                                                            tAACCGACTGGAACTGGGCCGCCTCAGTGCCGGTGATGCCTTCATCAGTGCCGTGGACGGCTTGGCCTTCTTCCAGCTGGAGGACGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTATCGGGTCTTCGATACCGACCCGCTGa                                                                      <  5:93019/149‑1 (MQ=255)
                                                                              aCCGACTGGAACTGGGCCGCCTCAGTGCCGGTGATGCCTTCATCAGTGCCGTGGACGGCTTGGCCTTCTTCCAGCTGGAGGACGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTATCGGGTCTTCGATACCGACCCGCTGAAt                                                                    >  7:383212/1‑149 (MQ=255)
                                                                               ccGACTGGAACTGGGCCGCCTCAGTGCCGGTGATGCCTTCATCAGTGCCGTGGACGGCTTGGCCTTCTTCCAGCTGGAGGACGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTATCGGGTCTTCGATACCGACCCGCTGAATc                                                                   <  1:242423/149‑1 (MQ=255)
                                                                                cGACTGGAACTGGGCCGCCTCAGTGCCGGTGATGCCTTCATCAGTGCCGTGTACGGCTTGGCCTTCTTCCAGCTGGAGGACGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTATCGGGTCTTCGATACCGACCCGCTGAATcc                                                                  <  3:248629/149‑1 (MQ=255)
                                                                                cGACTGGAACTGGGCCGCCTCAGTGCCGGTGATGCCTTCATCAGTGCCGTGGACGGCTTGGCCTTCTTCCAGCTGGAGGACGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTATCGGGTCTTCGATACCGACCCGCTGAATcc                                                                  >  3:277099/1‑149 (MQ=255)
                                                                                     ggAACTGGGCCGCCTCAGTGCCGGTGATGCCTTCATCAGTGCCGTGGACGGCTTGGCCTTCTTCCAGCTGGAGGACGGCCAGACCGtt                                                                                                                          <  1:384096/88‑1 (MQ=255)
                                                                                     ggAACTGGGCCGCCTCAGTGCCGGTGATGCCTTCATCAGTGCCGTGGACGGCTTGGCCTTCTTCCAGCTGGAGGACGGCCAGACCGtt                                                                                                                          >  2:384096/1‑88 (MQ=255)
                                                                                      gAACTGGGCCGCCTCAGTGCCGGTGATGCCTTCATCAGTGCCGTGGACGGCTTGGCCTTCTTCCAGCTGGAGGACGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTATCGGGTCTTCGATACCGACCCGCTGAATCCGAAACg                                                            >  5:3302/1‑149 (MQ=255)
                                                                                             gccgccTCAGTGCCGGTGATGCCTTCATCAGTGCCGTGGACGGCTTGGCCTTCTTCCAGCTGGAGGACGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTATCGGGTCTTCGATACCGACCCGCTGAATCCGAAACGCTCACgg                                                     <  1:231730/149‑1 (MQ=255)
                                                                                             gccgccTCAGTGCCGGTGATGCCTTCATCAGTGCCGTGGACGGCTTGGCCTTCTTCCAGCTGGAGGACGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTATCGGGTCTTCGATACCGACCCGCTGAATCCGAAACGCTCACgg                                                     >  1:304188/1‑149 (MQ=255)
                                                                                             gccgccTCAGTGCCGGTGATGCCTCCATCAGTGCCGTGGACGGCTTGGCCTTCTTCCAGCTGGAGGACGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTATCGGGTCTTCGATACCGACCCGCTGAATCCGAAACGCTCACgg                                                     >  8:37416/1‑149 (MQ=255)
                                                                                                   tCAGTGCCGGTGATGCCTTCATCAGTGCCGTGGACGGCTTGGCCTTCTTCCAGCTGGAGGACGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTATCGGGTCTTCGATACCGACCCGCTGAATCCGAAACGCTCACGGCTGAt                                                >  2:142350/1‑148 (MQ=255)
                                                                                                   tCAGTGCCGGAGATGCCTTCATCAGTGCCGTGGACGGCTTGGCCTTCTTCCAGCTGGAGGACGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTATCGGGTCTTCGATACCGACCCGCTGAATCCGAAACGCTCACGGCTGAt                                                >  8:139612/1‑148 (MQ=255)
                                                                                                          cGGTGATGCCTTCATCAGTGCCGTGGACGGCTTGGCCTTCTTCCAGCTGGAGGACGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTATCGGGTCTTCGATACCGACCCGCTGAATCCGAAACGCTCACGGCTGATCAGGCTCg                                        <  1:264400/149‑1 (MQ=255)
                                                                                                                tGCCTTCATCAGTGCCGTGGACGGCTTGGCCTTCTTCCAGCTGGAGGACGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTATCGGGTCTTCGATACCGACCCGCTGAATCCGAAACGCTCACGGCTGATCAGGCTCGGCGAcc                                  <  7:56936/149‑1 (MQ=255)
                                                                                                                  ccTTCATCAGTGCCGTGGACGGCTTGGCCTTCTTCCAGCTGGAGGACGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTATCGGGTCTTCGATACCGACCCGCTGAATCCGAAACGCTCACGGCTGATCAGGCTCGGCGACCg                                 <  8:142933/148‑1 (MQ=255)
                                                                                                                     tcatcaGTGCCGTGGACGGCTTGGCCTTCTTCCAGCTGGAGGCCGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTATCGGGTCTTCGATACCGACCCGCTGAATCCGAAACGCTCACGGCTGATCAGGCTCGGCGCCGGCact                             >  8:324787/1‑147 (MQ=255)
                                                                                                                         caGTGCCGTGGACGGCTTGGCCTTCTTCCAGCTGGAGGACGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTATCGGGTCTTCGATACCGACCCGCTGAATCCGAAACGCTCACGGCTGATCAGGCTCGGCGACCGCAAtctc                          >  8:60894/1‑148 (MQ=255)
                                                                                                                         caGTGCCGTGGACGGCTTGGCCTTCTTCCAGCTGGAGGACGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTATCGGGTCTTCGATACCGACCCGCTGAATCCGAAACGCTCACGGCTGATCAGGCTCGGCGACCGCAAtctc                          >  3:397727/1‑148 (MQ=255)
                                                                                                                         caGTGCCGTGGACGGCTTGGCCTTCTTCCAGCTGGAGGACGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTATCGGGTCTTCGATACCGACCCGCTGAATCCGAAACGCTCACGGCTGATCAGGCTCGGCGACCGCAAactc                          >  7:118562/1‑148 (MQ=255)
                                                                                                                           gTGCCGTGGACGGCGTGGCCTTCTTCCAGCTCGAGGACGGCCAGCCCGTTGTCGAAGATATCAATGAAGGGCTTTATCGCTTCTTCGATACCGACCCGCTGAATCCGAAACGCTCACCGCTGATCAGGCTCGGCGACCGCAATCTCAAt                       <  8:279111/149‑1 (MQ=255)
                                                                                                                                    aCGGCTTGGCCTTCTTCCAGCTGGAGGACGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTATCGGGTCTTCGATACCGACCCGCTGAATCCGAAACGCTCACGGCTGATCAGGCTCGGCGACCGCAATCTCAATTGCCTCGAc              <  8:383212/149‑1 (MQ=255)
                                                                                                                                     cggtctGGCCGTCCTCCAGCTGGAAGAAGGCCAAGCCGTTGTCGAAGATATCAATGAAGGGCTTTATCGGGTCTTCGATACCGACCCGCTGAATCCGAAACGCTCACGGCTGATCAGGCTCGGCGACCGCAATCTCAATTGCCTc                 >  5:156001/6‑145 (MQ=38)
                                                                                                                                     cggtctGGCCGTCCTCCAGCTGGAAGAAGGCCAAGCCGTTGTCGAAGATATCAATGAAGGGCTTTATCGGGTCTTCGATACCGACCCGCTGAATCCGAAACGCTCACGGCTGATCAGGCTCGGCGACCGCAATCTCAATTGCCTc                 <  6:156001/140‑1 (MQ=38)
                                                                                                                                        tctGGCCGTCCTCCAGCTGGAAGAAGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTATCGGGTCTTCGATACCGACCCGCTGAATCCGAAACGCTCACGGCTGATCAGGCTCGGCGACCGCAATCTCAATTGCCTCGACctgc          <  4:65895/147‑1 (MQ=255)
                                                                                                                                        tatGGCCGTCCTCCAGCTGGAAGAAGGCCAGACCGTTGTCTAAGATATCAATGAAGGGCTTTATCGGGTCTTCGATACCGACCCGCTGAATCCGAAACGCTCACGGCTGATCAGGCTCGGCGACCGCAATCTCAATTGCCTCGACctgc          >  8:250617/3‑149 (MQ=255)
                                                                                                                                        cTTGGCCTTCTTCCAGCTGGAGGACGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTATCGGGTCTTCGATACCGACCCGCTGAATCCGAAACGCTCACGGCTGATCAGGCTCGGCGACCGCAATCTCAATTGCCTctccatgc          >  4:307958/1‑142 (MQ=255)
                                                                                                                                         ttGGCCTTCTTCCAGCTGGAGGACGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTATCGGGTCTTCg                                                                                     <  3:118642/73‑1 (MQ=255)
                                                                                                                                         ttGGCCTTCTTCCAGCTGGAGGACGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTATCGGGTCTTCg                                                                                     >  4:118642/1‑73 (MQ=255)
                                                                                                                                           ggCCTTCTTCCAGCTGGAGGACGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTATCGGGTCTTCGATACCGACCCGCTGAATCCGAAACGCTCACGGCTGATCAGGCTCGGCGACCGCAATCTCAATTGCCTCGACctgctg        <  3:198939/148‑1 (MQ=255)
                                                                                                                                            gCCTTCTTCCAGCTGGAGGACGGCCAGACCGTTGTCGAAGATATAAATGAAGGGCTTTATCGGGTCTTCGATACCGACCCGCTGAATCCGAAACGCTCACGGCTGATCAGGCTCGGCGACCGCAATCTCAATTGCCTCGACCTGCTGt       <  8:118562/148‑1 (MQ=255)
                                                                                                                                                 cttcCAGCTGGAGGACGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTATCGTGTCTTCGATACCGACCCGCTGAATCCGAAACGCTCACGGCTGATCAGGCTCGGCGACCGCAATCTCAATTGCCTCGACCTGCTGTacgac  <  5:111934/148‑1 (MQ=255)
                                                                                                                                                 cttcCAGCTGGAGGACGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTATCGGGTCTTCGATACCGACCCGCTGAATCCGAAACGCTCACGGCTGATCAGGCTCGGCGACCGCAATCTCAATTGCCTCGACCTGCTGTacga   <  4:397727/147‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
TCCCCTATGAAGTGGAGCTGGCCCGCGTCCCACATCCTCAAGGCGGTACGCTGTGGATCTACATCGACAACGACGGTAACCGACTGGAACTGGGCCGCCTCAGTGCCGGTGATGCCTTCATCAGTGCCGTGGACGGCTTGGCCTTCTTCCAGCTGGAGGACGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTATCGGGTCTTCGATACCGACCCGCTGAATCCGAAACGCTCACGGCTGATCAGGCTCGGCGACCGCAATCTCAATTGCCTCGACCTGCTGTACGAC  >  NC_002947/4614648‑4614940

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