Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,129,753 C→T 6.4% V65M (GTG→ATG)  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,7530CT 6.4% 180.8 / 2.1 62V65M (GTG→ATG) gcvH‑Iglycine cleavage system protein H
Reads supporting (aligned to +/- strand):  ref base C (33/25);  new base T (2/2);  total (35/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.86e-01

GTCGAGCAGAGCGGTTACGGCTTCTGCATCGGCAGGGATGAAGCGGAAGAACCAGCCTTCGCCCATCGGGTCTTCGTTGACCAGTTCAGGGCTGTCTTCGAGCTGGCCGTTGACCTCGACCACTTCACCCGTCAGTGGCATGTACACGCCGCTGGCGGCCTTGACCGATTCGACGGTGGACGCTTCGTTGCCTTTTTCATACTGCTGCAACTCTGGCAGTTGCACGAAAACCACGTCGCCAAGGGCATTCTGGGCGTAGGCGGTGATGCCCACGGTGACGCTGCCATCGGCTT  >  NC_002947/1129607‑1129899
                                                                                                                                                  |                                                                                                                                                  
gTCGAGCAGAGCGGTTACGGCTTCTGCATCGGCAGGGATGAAGCGGAAGAACCAGCCTTCGCCCATCGGGTCTTCGTTGACGAGTTCAGGGCTGTCTTCGAGCTGGCCGTTGACCTCGACCACTTCACCCGTCAGTGGCATGTACAcg                                                                                                                                                   <  1:237220/148‑1 (MQ=255)
        gagCGGTTACGGCTTCTGCATCGGCAGGGATGAAGCGGAAGAACCAGCCTTCGCCCATCGGGTCTTCGTTGACCAGTTCAGGGCTGTCTTCGAGCTGGCCGTTGACCTCGACCACTTCACCCGTCAGTGGCATGTACACGCCGCTggcg                                                                                                                                          >  8:74211/1‑149 (MQ=255)
         tgCGGTTACGGGTTCTGCATCGGCAGGGATGATGCGGAATAACCAGCCTTGGGCCATCGGGTCTTCGTTGACCAGTTCTGGGCTGTCTTCGAGCTGGCCGTTGAGGTCGACCACTTCACCCGTCAGTGGCATGTACACGCCGCTggcgg                                                                                                                                         <  4:212270/148‑1 (MQ=255)
             gTTACGGCTTCTGCATCGGCAGGGATGAAGCGGAAGAACCAGCCTTCGCCCATCGGGTCTTCGTTGACCAGTTCAGGGCTGTCTTCGAGCTGGCCGTTGACCTCGACCACTTCACCCGTCAGCGGCGTGTACATGCCACTgac                                                                                                                                           <  2:259983/143‑3 (MQ=255)
             gTTACGGCTTCTGCATCGGCAGGGATGAAGCGGAAGAACCAGCCTTCGCCCATCGGGTCTTCGTTGACCAGTTCAGGGCTGTCTTCGAGCTGGCCGTTGACCTCGACCACTTCACCCGTCAGCGGCGTGTACATGCCACTgac                                                                                                                                           >  1:259983/1‑141 (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                                                                                           >  6:234782/1‑144 (MQ=255)
                                                            gCCCATCGGGTCTTCGTTGACCAGTTCAGGGCTGTCTTCGAGCTGGCCGTTGACCTCGACCACTTCACCCGTCAGTGGCATGTACACGCCGCTGGCGGCCTTGACCGATTCGACGGTGGACGCTTCGTTGCCTTTTTCATActg                                                                                           <  5:234782/144‑1 (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                                                                 >  1:112388/1‑149 (MQ=255)
                                                                                 cAGTTCAGGGCTGTCTTCGAGCTGGCCGTTGACCTCGACCACTTCACCCGTCAGTGGCATGTACACGCCGCTGGCGGCCTTGACCGATTCGACGGTGGACGCTTCGTTGCCTTTTTCATACTGCTGCAACTCTGGCAGTTGCACGaaaa                                                                 >  4:134540/1‑149 (MQ=255)
                                                                                 cAGTTCAGGGCTGTCTTCGAGCTGGCCGTTGACCTCGACCACTTCACCCGTCAGTGGCATGTACACGCCGCTGGCGGCCTTGACCGATTCGACGGTGGACGCTTCGTTGCCTTTTTCATACTGCTGCAACTCTGGCAGTTGCACGaaaa                                                                 >  8:157564/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                                                                                                                             >  8:218019/1‑81 (MQ=255)
                                                                                         ggCTGTCTTCGAGCTGGCCGTTGACCTCGACCACTTCACCCGTCAGTGGCATGTACACGCCGCTGGCGGCCTTGACCGAtt                                                                                                                             <  7:218019/81‑1 (MQ=255)
                                                                                         ggCTGTCTTCGAGCTGGCCGTTGACCTCGACCACTTCACCCGTCAGTGGCATGTACACGCCGCTGGCGGCCTTGACCGATTCGACGGTGGACGCTTCGTTGCCTTTTTCATACTGCTGCAACTCTGGCAGTTGCACGAAAACCACGTCg                                                         >  7:341010/1‑149 (MQ=255)
                                                                                             gTCTTCGAGCTGGCCGTTGACCTCGACCACTTCACCCGTCAGTGGCATGTACACGCCGCTGGCGGCCTTGACCGATTCGACGGTGGACGCTTCGTTGCCTTTTTCATACTGCTGCAACTCTGGCAGTTGCACGAAAACCACGTCGcc                                                       >  1:47413/1‑147 (MQ=255)
                                                                                              tCTTCGAGCTGGCCGTTGACCTCGACCACTTCACCCGTCAGTGGCATGTACACGCCGCTGGCGGCCTTGACCGATTCGACGGTGGACGCTTCGTTGCCTTTTTCATACTGCTGCAACTCTGGCAGTTGCACGAAAACCACGTCGCCAAg                                                    <  7:157564/149‑1 (MQ=255)
                                                                                              tCTTCGAGCTGGCCGTTGACCTCGACCACTTCACCCGTCAGTGGCATGTACACGCCGCTGGCGGCCTTGACCGATTCGACGGTGGACGCTTCGTTGCCTTTTTCATACTGCTGCAACTCTGGCAGTTGCACGAAAACCACGTCGCCAAg                                                    >  1:174006/1‑149 (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                                          <  2:174006/148‑1 (MQ=255)
                                                                                                         gCCGTTGACCTCGACCACTTCACCCGTCAGTGGCATGTACACGCCGCTGGCGGCCTTGACCGATTCGACGGTGGACGCTTCGTTGCCTTTTTCATACTGCTGCAACTCTGGCAGTTGCACGAAAACCACGTCGCCAAGGGCATTCTgg                                          <  7:48237/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                                                <  1:223906/124‑1 (MQ=255)
                                                                                                                           ttCACCCGTCAGTGGCATGTACACGCCGCTGGCGGCCTTGACCGATTCGACGGTGGACGCTTCGTTGCCTTTTTCATACTGCTGCAACTCTGGCAGTTGCACGAAAACCACGTCGCCAAGGGCa                                                >  2:223906/1‑124 (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)
                                                                                                                                          cATGTACACGCCGCTGGCGGCCTTGACCGATTCGACGGTGGACGCTTCGTTGCCTTTTTCATACTGCTGCAACTCTGGCAGTTGCACGAAAACCACGTCGCCAAGGGCATTCTgg                                          <  7:52315/115‑1 (MQ=255)
                                                                                                                                          cATGTACACGCCGCTGGCGGCCTTGACCGATTCGACGGTGGACGCTTCGTTGCCTTTTTCATACTGCTGCAACTCTGGCAGTTGCACGAAAACCACGTCGCCAAGGGCATTCTgg                                          >  8:52315/1‑115 (MQ=255)
                                                                                                                                                cacGCCGCTGGCGGCCTTGACCGATTCGACGGTGGACGCTTCGTTGCCTTTTTCATACTGCTGCAACTCTGGCAGTTGCACGAAAACCACGTCGCCAAGGGCATTCTGGGCGTAGGCGGTGATGCCCACGGTGACGCTGCCATCGGCtt  <  3:134540/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GTCGAGCAGAGCGGTTACGGCTTCTGCATCGGCAGGGATGAAGCGGAAGAACCAGCCTTCGCCCATCGGGTCTTCGTTGACCAGTTCAGGGCTGTCTTCGAGCTGGCCGTTGACCTCGACCACTTCACCCGTCAGTGGCATGTACACGCCGCTGGCGGCCTTGACCGATTCGACGGTGGACGCTTCGTTGCCTTTTTCATACTGCTGCAACTCTGGCAGTTGCACGAAAACCACGTCGCCAAGGGCATTCTGGGCGTAGGCGGTGATGCCCACGGTGACGCTGCCATCGGCTT  >  NC_002947/1129607‑1129899

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