Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,129,742 T→C 6.2% P68P (CCA→CCG)  gcvH‑I ← glycine cleavage system protein H

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,129,7420TC 6.2% 182.0 / 2.1 64P68P (CCA→CCG) gcvH‑Iglycine cleavage system protein H
Reads supporting (aligned to +/- strand):  ref base T (33/27);  new base C (2/2);  total (35/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.77e-01

GCGTCCTGGTCGAGCAGAGCGGTTACGGCTTCTGCATCGGCAGGGATGAAGCGGAAGAACCAGCCTTCGCCCATCGGGTCTTCGTTGACCAGTTCAGGGCTGTCTTCGAGCTGGCCGTTGACCTCGACCACTTCACCCGTCAGTGGCATGTACACGCCGCTGGCGGCCTTGACCGATTCGACGGTGGACGCTTCGTTGCCTTTTTCATACTGCTGCAACTCTGGCAGTTGCACGAAAACCACGTCGCCAAGGGCATTCTGGGCGTAGGCGGTGATGCCCACGGTGA  >  NC_002947/1129599‑1129884
                                                                                                                                               |                                                                                                                                              
tcGTCCTGGTCGAGCAGAGCTGTTACGGCAGCTGCATAGGCAGGGATGAAGCGGAAGAACCAGCCTTCGCCCATCGGGTCTTCGTTGACCAGTTCAGGGCTGTCTTCGAGCTGGCCGTTGACCTCGACCACTTCACCCGTCAGTGGCAt                                                                                                                                           <  8:322024/148‑1 (MQ=255)
   tCCTGGTCGAGCAGAGCGGTTACGGCTTCTGCATCGGCAGGGATGAAGCGGAAGAACCAGCCTTCGCCCATCGGGTCTTCGTTGACCAGTTCAGGGCTGTCTTCGAGCTGGCCGTTGACCTCGACCACTTCACCCGTCAGTGGCATGTa                                                                                                                                        <  4:28918/149‑1 (MQ=255)
     cTGGTCGAGCAGAGCGGTTACGGCTTCTGCATCGGCAGGGAAGAAGCGGAAGAACCAGCCTTCGCCCATCGGGTCTTCGTTGACCAGTTCAGGGCTGTCTTCGAGCTGGCCCTTGACCTCGACCACTTCACCCGTCAGTGGCATGTaca                                                                                                                                      >  3:212270/1‑149 (MQ=255)
        gTCGAGCAGAGCGGTTACGGCTTCTGCATCGGCAGGGATGAAGCGGAAGAACCAGCCTTCGCCCATCGGGTCTTCGTTGACGAGTTCAGGGCTGTCTTCGAGCTGGCCGTTGACCTCGACCACTTCACCCGTCAGTGGCATGTACAcg                                                                                                                                    <  1:237220/148‑1 (MQ=255)
                gagCGGTTACGGCTTCTGCATCGGCAGGGATGAAGCGGAAGAACCAGCCTTCGCCCATCGGGTCTTCGTTGACCAGTTCAGGGCTGTCTTCGAGCTGGCCGTTGACCTCGACCACTTCACCCGTCAGTGGCATGTACACGCCGCTggcg                                                                                                                           >  8:74211/1‑149 (MQ=255)
                 tgCGGTTACGGGTTCTGCATCGGCAGGGATGATGCGGAATAACCAGCCTTGGGCCATCGGGTCTTCGTTGACCAGTTCTGGGCTGTCTTCGAGCTGGCCGTTGAGGTCGACCACTTCACCCGTCAGTGGCATGTACACGCCGCTggcgg                                                                                                                          <  4:212270/148‑1 (MQ=255)
                     gTTACGGCTTCTGCATCGGCAGGGATGAAGCGGAAGAACCAGCCTTCGCCCATCGGGTCTTCGTTGACCAGTTCAGGGCTGTCTTCGAGCTGGCCGTTGACCTCGACCACTTCACCCGTCAGCGGCGTGTACATGCCACTgac                                                                                                                            >  1:259983/1‑141 (MQ=255)
                     gTTACGGCTTCTGCATCGGCAGGGATGAAGCGGAAGAACCAGCCTTCGCCCATCGGGTCTTCGTTGACCAGTTCAGGGCTGTCTTCGAGCTGGCCGTTGACCTCGACCACTTCACCCGTCAGCGGCGTGTACATGCCACTgac                                                                                                                            <  2:259983/143‑3 (MQ=255)
                      ttACGGCTTCTGCATCGGCAGGGATGAAGCGGAAGAACCAGCCTTCGCCCATCGGGTCTTCGTTGACCAGTTCAGGGCTGTCTTCGAGCTGGCCGTTGACCTCGACCACTTCACCCGCCAGCGGCGTGTACATGCCACTgacgg                                                                                                                          <  5:231554/144‑1 (MQ=255)
                           gCTTCTGCATCGGCAGGGATGAAGCGGAAGAACCAGCCTTCGACCATCGTGTCTTCGTTGACAAGTTCAGGGCTGTCTTCGAGCTGTCCGTTGACCTCGACCACTTCACACGTCAGTGGCATGTACACGCCGCTGGCGGACttca                                                                                                                    >  8:323189/1‑143 (MQ=255)
                               cTGCATCGGCAGGGATGAAGCGGAAGAACCAGCCTTCGCCCATCGGGTCTTCGTTGACCAGTTCAGGGCTGTCTTCGAGCTGGCCGTTGACCTCGACCACTTCACCCGTCAGTGGCATGTACACGCCGCTGGCGGCCTTGACCGAtt                                                                                                              <  1:87530/147‑1 (MQ=255)
                                     cGGCAGGGATGAAGCGGAAGAACCAGCCTTCGCCCATCGGGTCTTCGTTGACCAGTTCAGGGCTGTCTTCGAGCTGGCCGTTGACCTCGACCACTTCACCCGTCAGTGGCATGTACACGCCGCTGGCGGCCTTGACCGATTCGACGGTg                                                                                                      >  5:361636/1‑149 (MQ=255)
                                           ggATGAAGCGGAAGAACCAGCCTTCGCCCATCGGGTCTTCGTTGACCAGTTCAGGGCTGTCTTCGAGCTGGCCGTTGACCTCGACCACTTCACCCGTCAGTGGCATGTACACGCCGCTGGCGGCCTTGACCGATTCGACGGTGGACGCt                                                                                                >  6:178471/1‑149 (MQ=255)
                                                 aGCGGAAGAACCAGCCTTCGCCCATCGGGTCTTCGTTGACCAGTTCAGGGCTGTCTTCGAGCTGGCCGTTGACCTCGTCCACTTCACCCGTCAGTGGCATGTACACGCCGCTGGCGGCCTTGACCGATTCGACGGTGGACGCTTCGTTg                                                                                          >  8:105086/1‑149 (MQ=255)
                                                 aGCGGAAGAACCAGCCTTCGCCCATCGGGTCTTCGTTGACCAGTTCAGGGCTGTCTTCGAGCTGGCCGTTGACCTCGACCACTTCACCCGTCAGTGGCATGTACACGCCGCTGGCGGCCTTGACCGATTCGACGGTGGACGCTTCGTTg                                                                                          >  8:199024/1‑149 (MQ=255)
                                                  gCGGAAGAACCAGCCTTCGCCCATCGGGTCTTCGTTGACCAGTTCAGGGCTGTCTTCGAGCTGGCCGTTGACCTCGACCACTTCACCCGTCAGTGGCATGTACACGCCGCTGGCGGCCTTCACCGATTCGACGGTGGACGCTTCGTTGc                                                                                         <  1:212737/149‑1 (MQ=255)
                                                     gaagaaCCAGCCTTCGCCCATCGGGTCTTCGTTGACCAGTTCAGGGCTGTCTTCGAGCTGGCCGTTGACCTCGACCACTTCACCCGTCAGTGGCATGTACACGCCGCTGGCGGCCTTGACCGATTCGACGGTGGACGCTTCGTTGCCtt                                                                                      >  7:194022/1‑149 (MQ=255)
                                                       agaaCCAGCCTTCGCCCATCGGGTCTTCGTTGACCAGTTCAGGGCTGTCTTCGAGCTGGCCGTTGACCTCGACCACTTCACCCGTCAGTGGCATGTACACGCCGCTGGCGGCCTTGACCGATTCGACGGTGGACGCTTCGTTGCCtttt                                                                                    <  7:214320/149‑1 (MQ=255)
                                                            cAGCCTTCGCCCATCGGGTCTTCGTTGACCAGTTCAGGGCTGTCTTCGAGCTGGCCGTTGACCTCGACCACTTCACCCGTCAGTGGCATGTACACGCCGCTGGCGGCCTTGACCGATTCGACGGTGGACGCTTCGTTGCCTTTTTCATa                                                                               <  1:37824/149‑1 (MQ=255)
                                                                cTTCGCCCATCGGGTCTTCGTTGACCAGTTCAGGGCTGTCTTCGAGCTGGCCGTTGACCTCGACCACTTCACCCGTCAGTGGCATGTACACGCCGCTGGCGGCCTTGACCGATTCGACGGTGGACGCTTCGTTGCCTTTTTCATActgc                                                                           <  8:221935/149‑1 (MQ=255)
                                                                cTTCGCCCATCGGGTCTTCGTTGACCAGTTCAGGGCTGTCTTCGAGCTGGCCGTTGACCTCGACCACTTCACCCGTCAGTGGCATGTACACGCCGCTGGCGGCCTTGACCGATTCGACGGTGGACGCTTCGTTGCCTTTTTCATActgc                                                                           >  7:38845/1‑149 (MQ=255)
                                                                  tCGCCCATCGGGTATTCGTTTACCAGTTCAGGGCTGTCTTCGAGCTGGCCGTTGACCTCGACCACTTCACCCGTCAGTGGCATGTACACGCCGATGGCGGCCTTGAACGATTCGACGCTGGACTCTTCGTTGCCTTTTTCATActgatg                                                                         >  6:345655/1‑149 (MQ=255)
                                                                    gCCCATCGGGTCTTCGTTGACCAGTTCAGGGCTGTCTTCGAGCTGGCCGTTGACCTCGACCACTTCACCCGTCAGTGGCATGTACACGCCGCTGGCGGCCTTGACCGATTCGACGGTGGACGCTTCGTTGCCTTTTTCATActg                                                                            <  5:234782/144‑1 (MQ=255)
                                                                    gCCCATCGGGTCTTCGTTGACCAGTTCAGGGCTGTCTTCGAGCTGGCCGTTGACCTCGACCACTTCACCCGTCAGTGGCATGTACACGCCGCTGGCGGCCTTGACCGATTCGACGGTGGACGCTTCGTTGCCTTTTTCATActg                                                                            >  6:234782/1‑144 (MQ=255)
                                                                      ccATCGGGTCTTCGTTGACCAGTTCAGGGCTGTCTTCGAGCTGGCCGTTGACCTCGTCCACTTCACCCGTCAGTGGCATGTACACGCCGCTGGCGGCCTTGACCGATTCGACGGTGGACGCTTCGTTGCCTTTTTCATACTGCTGCAAc                                                                     <  7:105086/149‑1 (MQ=255)
                                                                           gggTCTTCGTTGACCAGTTCAGGGCTGTCTTCGAGCTGGGCGTTGACCTCGACCACTTCACCCGTCAGTGGCATGTACACGCCGCTGGCGGGCTTGACCGATTCGACGGTGGACGCTTCGTTGCCTTTTTCATACTGCTGCAACTCTg                                                                 <  8:119495/148‑1 (MQ=255)
                                                                             gTCTTCGTTGACCAGTTCAGGGCTGTCTTCGAGCTGGCCGTTGACCTCGACCACTTCACCCGTCAGTGGCATGTACACGCCGCTGGCGGCCTTGACCGATTCGACGGTGGACGCTTCGTTGCCTTTTTCATACTGCTGCAACTCTgg                                                                >  8:48237/1‑147 (MQ=255)
                                                                               cTTCGTTGACCAGTTCAGGGCTGTCTTCGAGCTGGCCGTTGACCTCGACCACTTCACCCGTCAGTGGCATGTACACGCCGCTGGCGGCCTTGACCGATTCGACGGTGGACGCTTCGTTGCCTTTTTCATACTGCTGCAACTCTGGCAGt                                                            <  8:295103/149‑1 (MQ=255)
                                                                                ttCGTTGACCAGTTCAGGGCTGTCTTCGAGCTGGCCGTTGACCTCGACCACTTCACCCGTCAGTGGCATGTACACGCCGCTGGCGGCCTTGACCGATTCGACGGTGGACGCTTCGTTGCCTTTTTCATACTGCTGCAACTCTGGCAGtt                                                           >  2:219163/1‑149 (MQ=255)
                                                                                  cGTTGACCAGTTCAGGGCTGTCTTCGAGCTGGCCGTTGACCTCGACCACTTCACCCGTCAGTGGCATGTACACGCCGCTGGCGGCCTTTACCGATTCGACGGTGGACGCTTCGTTGCCTTTTTCATACTGCTGCAACTCTGGCAGTTg                                                          >  1:23567/1‑148 (MQ=255)
                                                                                    ttGACCAGTTCAGGGCTGTCTTCGAGCTGGCCGTTGACCTCGACCACTTCACCCGTCAGTGGCATGTACACGCCGCTGGCGGCCTTGACCGATTCGACGGTGGACGCTTCGTTGCCTTTTTCATACTGCTGCAACTCTGGCAGTTGCAc                                                       <  7:323189/149‑1 (MQ=255)
                                                                                       aCCAGTTCAGGGCTGTCTTCGAGCTGGCCGTTGACCTCGACCACTTCACCCGTCAGTGGCATGTACACGCCGCTGGCGGCCTTGACCGATTCGACGGTGGACGCTTCGTTGCCTTTTTCATACTGCTGCAACTCTGGCAGTTGCACGa                                                     >  5:4845/1‑148 (MQ=255)
                                                                                        ccAGTTCAGGGCTGTCTTCGAGCTGGCCGTTGACCTCGACCACTTCACCCGTCAGTGGCATGTACACGCCGCTGGCGGCCTTGACCGATTCGACGGTGGACGCTTCGTTGCCTTTTTCATACTGCTGCAACTCTGGCAGTTGCACGaaa                                                   <  5:178471/149‑1 (MQ=255)
                                                                                         cAGTTCAGGGCTGTCTTCGAGCTGGCCGTTGACCTCGACCACTTCACCCGTCAGTGGCATGTACACGCCGCTGGCGGCCTTGACCGATTCGACGGTGGACGCTTCGTTGCCTTTTTCATACTGCTGCAACTCTGGCAGTTGCACGaaaa                                                  >  8:157564/1‑149 (MQ=255)
                                                                                         cAGTTCAGGGCTGTCTTCGAGCTGGCCGTTGACCTCGACCACTTCACCCGTCAGTGGCATGTACACGCCGCTGGCGGCCTTGACCGATTCGACGGTGGACGCTTCGTTGCCTTTTTCATACTGCTGCAACTCTGGCAGTTGCACGaaaa                                                  >  1:112388/1‑149 (MQ=255)
                                                                                         cAGTTCAGGGCTGTCTTCGAGCTGGCCGTTGACCTCGACCACTTCACCCGTCAGTGGCATGTACACGCCGCTGGCGGCCTTGACCGATTCGACGGTGGACGCTTCGTTGCCTTTTTCATACTGCTGCAACTCTGGCAGTTGCACGaaaa                                                  >  4:134540/1‑149 (MQ=255)
                                                                                            ttCAGGGCTGTCTTCGAGCTGGCCGTTGACCTCGACCACTTCACCCGTCAGTGGCATGTACACGCCGCTGGCGGCCTTGACCGATTCGACGGTGGACGCTTCGTTGCCtttt                                                                                    <  4:353861/112‑1 (MQ=255)
                                                                                            ttCAGGGCTGTCTTCGAGCTGGCCGTTGACCTCGACCACTTCACCCGTCAGTGGCATGTACACGCCGCTGGCGGCCTTGACCGATTCGACGGTGGACGCTTCGTTGCCtttt                                                                                    >  3:353861/1‑112 (MQ=255)
                                                                                               aGGGCTGTCTTCGAGCTGGCCGTAGACCTCGACCACTTCACCAGGCAGTGGCATTTACACGCCGCATTCGGCCTTGACCGCTTCGACGCTGCACGCTTCGTAGCCAATATCATACTGCAGCAACTATGCCAGTAGCACGAAACACACGt                                            >  6:356010/1‑149 (MQ=255)
                                                                                                gggCTGTCTTCGAGCTGGCCGTTGACCTCGACCACTTCACCCGTCAGTGGCATGTACACGCCGCTGGCGGCCTTGACCGATTCGACGGTGGACGCTTCGTTGCCTTTTTCATACTGCTGCAACTCTGGCAGTTGCACGAAAACCACGTc                                           <  7:199024/149‑1 (MQ=255)
                                                                                                 ggCTGTCTTCGAGCTGGCCGTTGACCTCGACCACTTCACCCGTCAGTGGCATGTACACGCCGCTGGCGGCCTTGACCGAtt                                                                                                              <  7:218019/81‑1 (MQ=255)
                                                                                                 ggCTGTCTTCGAGCTGGCCGTTGACCTCGACCACTTCACCCGTCAGTGGCATGTACACGCCGCTGGCGGCCTTGACCGAtt                                                                                                              >  8:218019/1‑81 (MQ=255)
                                                                                                 ggCTGTCTTCGAGCTGGCCGTTGACCTCGACCACTTCACCCGTCAGTGGCATGTACACGCCGCTGGCGGCCTTGACCGATTCGACGGTGGACGCTTCGTTGCCTTTTTCATACTGCTGCAACTCTGGCAGTTGCACGAAAACCACGTCg                                          >  7:341010/1‑149 (MQ=255)
                                                                                                     gTCTTCGAGCTGGCCGTTGACCTCGACCACTTCACCCGTCAGTGGCATGTACACGCCGCTGGCGGCCTTGACCGATTCGACGGTGGACGCTTCGTTGCCTTTTTCATACTGCTGCAACTCTGGCAGTTGCACGAAAACCACGTCGcc                                        >  1:47413/1‑147 (MQ=255)
                                                                                                      tCTTCGAGCTGGCCGTTGACCTCGACCACTTCACCCGTCAGTGGCATGTACACGCCGCTGGCGGCCTTGACCGATTCGACGGTGGACGCTTCGTTGCCTTTTTCATACTGCTGCAACTCTGGCAGTTGCACGAAAACCACGTCGCCAAg                                     >  1:174006/1‑149 (MQ=255)
                                                                                                      tCTTCGAGCTGGCCGTTGACCTCGACCACTTCACCCGTCAGTGGCATGTACACGCCGCTGGCGGCCTTGACCGATTCGACGGTGGACGCTTCGTTGCCTTTTTCATACTGCTGCAACTCTGGCAGTTGCACGAAAACCACGTCGCCAAg                                     <  7:157564/149‑1 (MQ=255)
                                                                                                      tCTTCGAGCTGGCCGTTGACCTCGACCACTTCACCCGTCAGT‑GCATGTACACGCCGCTGGCGGCCTTGACCGATTCGACGGTGGACGCTTCGTTGCCTTTTTCATACTGCTGCAACTCTGGCAGTTGCACGAAAACCACGTCGCCAAgg                                    <  1:207225/149‑1 (MQ=255)
                                                                                                          cGAGCTGGCCGTTGACCTCGACCACTTCACCCGTCAGTGGCATGTACACGCCGCTGGCGGCCTTGACCGATTCGACGGTGGACGCTTCGTTGCCTTTTTCATACTGCTGCAACTCTGGCAGTTGCACGAAAACCACGTCGCCAAGGGCa                                 >  6:132812/1‑149 (MQ=255)
                                                                                                           gAGCTGGCCGTTGACCTCGACCACTTCACCCGTCAGTGGCATGTACACGCCGCTGGCGGCCTTGACCGATTCGACGGTGGACGCTTCGTTGCCTTTTTCATACTGCTGCAACTCTGGCAGTTGCACGAAAACCACGTCGCCAAGGGCAt                                >  4:198762/1‑149 (MQ=255)
                                                                                                                 gCCGTTGACCTCGACCACTTCACCCGTCAGTGGCATGTACACGCCGCTGGCGGCCTTGACCGATTCGACGGTGGACGCTTCGTTGCCTTTTTCATACTGCTGCAACTCTGGCAGTTGCACGAAAACCACGTCGCCAAGGGCATTCTgg                           <  7:48237/148‑1 (MQ=255)
                                                                                                                 gCCGTTGACCTCGACCACTTCACCCGTCAGTGGCATGTACACGCCGCTGGCGGCCTTGACCGATTCGACGGTGGACGCTTCGTTGCCTTTTTCATACTGCTGCAACTCTGGCAGTTGCACGAAAACCACGTCGCCAAGGGCATTCTgg                           <  2:174006/148‑1 (MQ=255)
                                                                                                                          cTCGACCACTTCACCCGTCAGTGGCATGTACACGCCGCTGGCGGCCTTGACCGATTCGACGGTGGACGCTTCGTTGCCTTTTTCATACTGCTGCAACTCTGGCAGTTGCACGAAAACAACGTCGCCAAGGGCATTCTGGGCGTAGGCgg                 >  7:351500/1‑149 (MQ=255)
                                                                                                                           tCGACCACTTCACCCGTCAGTGGCATGTACACGCCGCTGGCGGCCTTGACCGATTCGACGGTGGACGCTTCGTTGCCTTTTTCATACTGCTGCAACTCTGGCAGTTGCACGAAAACCACGTCGCCAAGGGCATTCTGGGCGTAGGCGGt                >  3:64042/1‑149 (MQ=255)
                                                                                                                           tCGACCACTTCACCCGTCAGTGGCATGTACACGCCGCTGGCGGCCTTGACCGATTCGACGGTGGACGCTTCGTTGCCTTTTTCATACTGCTGCAACTCTGGCAGTTGCACGAAAACCACGTCGCCAAGGGCATTCTGGGCGTAGGCGGt                >  1:111963/1‑149 (MQ=255)
                                                                                                                             gACCACTTCACCCGTCAGTGGCATGTACACGCCGCTGGCGGCCTTGACCGATTCGACGGTGGACGCTTCGTTGCCTTTTTCATACTGCTGCAACTCTGGCAGTTGCACGAAAACCACGTCGCCAAGGGCATTCTGGGCGTAGGCGGTGa              <  7:74211/149‑1 (MQ=255)
                                                                                                                             gACCACTTCACCCGTCAGTGGCATGTACACGCCGCTGGCGGCCTTGACCGATTCGACGGTGGACGCTTCGTTGCCTTTTTCATACTGCTGCAACTCTGGCAGTTGCACGAAAACCACGTCGCCAAGGGCATTCTGGGCGTAGGCGGTGa              >  7:144440/1‑149 (MQ=255)
                                                                                                                                cACTTCACCCGTCAGTGGCATGTACACGCCGCTGGCGGCCTTGACCGATTCGACGGTGGACGCTTCGTTGCCTTTTTCATACTGCTGCAACTCTGGCAGTTGCACGAAAACCACGTCGCCAAGGGCATTCTGGGCGTAGGCGGTGATGc           <  8:341010/149‑1 (MQ=255)
                                                                                                                                cACTTCACCCGGCCGTCGCAGGTACACGCCGCGGGCGGCCTTGACCGATTCGCCGGTGGACGCTTCGTTGCCTTTTTCATACTGCTGCAACTCTGGCAGTTGCACGAAAACCACGTCGCCAAGGGCATTCTGGGCGTAGGCGGTGATGc           <  8:272391/149‑1 (MQ=255)
                                                                                                                                   ttCACCCGTCAGTGGCATGTACACGCCGCTGGCGGCCTTGACCGATTCGACGGTGGACGCTTCGTTGCCTTTTTCATACTGCTGCAACTCTGGCAGTTGCACGAAAACCACGTCGCCAAGGGCa                                 >  2:223906/1‑124 (MQ=255)
                                                                                                                                   ttCACCCGTCAGTGGCATGTACACGCCGCTGGCGGCCTTGACCGATTCGACGGTGGACGCTTCGTTGCCTTTTTCATACTGCTGCAACTCTGGCAGTTGCACGAAAACCACGTCGCCAAGGGCa                                 <  1:223906/124‑1 (MQ=255)
                                                                                                                                     cACCCGTCAGTGGCATGTACACGCCGCTGGCGGCCTTGACCGATTCGACGGTGGACGCTTCGTTGCCTTTTTCATACTGCTGCAACTCTGGCAGTTGCACGAAAACCACGTCGCCAAGGGCATTCTGGGCGTAGGCGGTGATGCCCACg      >  4:167551/1‑149 (MQ=255)
                                                                                                                                     aggccGCCAGCGGCGTGTACATGCCACTGACGGCCTTGACCGATTCGACGGTGGACGCTTCGTTGCCTTTTTCATACTGCTGCAACTCTGGCAGTTGCACGAAAACCACGTCGCCAAGGGCATTCTGGGCGTAGGCGGTGATGCCCACg      >  4:167881/4‑149 (MQ=255)
                                                                                                                                      aCCCGTCAGTGGCATGTACACGCCGCTGGCGGCCTTGACCGATTCGACGGTGGACGCTTCGTTGCCTTTTTCATACTGCTGCAACTCTGGCAGTTGCACGAAAACCACGTCGCCAAGGGCATTCTGGGCGTAGGCGGTGATGCCCACg      <  6:4845/148‑1 (MQ=255)
                                                                                                                                         cGTCAGTGGCATGTACACGCCGCTGGCGGCCTTGACCGATTCGACGGTGGACGCTTCGTTGCCTTTTTCATACTGCTGCAACTCTGGCAGTTGCACGAAAACCACGTCGCCAAGGGCATTCTGGGCGTAGGCGGTGATGCCCACGGTGa  >  4:5515/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                              
GCGTCCTGGTCGAGCAGAGCGGTTACGGCTTCTGCATCGGCAGGGATGAAGCGGAAGAACCAGCCTTCGCCCATCGGGTCTTCGTTGACCAGTTCAGGGCTGTCTTCGAGCTGGCCGTTGACCTCGACCACTTCACCCGTCAGTGGCATGTACACGCCGCTGGCGGCCTTGACCGATTCGACGGTGGACGCTTCGTTGCCTTTTTCATACTGCTGCAACTCTGGCAGTTGCACGAAAACCACGTCGCCAAGGGCATTCTGGGCGTAGGCGGTGATGCCCACGGTGA  >  NC_002947/1129599‑1129884

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