Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 371,124 T→G 9.5% F202C (TTC→TGC)  dgcA → dimethylglycine dehydrogenase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947371,1240TG 9.5% 184.7 / 10.2 63F202C (TTC→TGC) dgcAdimethylglycine dehydrogenase subunit
Reads supporting (aligned to +/- strand):  ref base T (28/29);  new base G (4/2);  total (32/31)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.20e-01

GCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGAAAGCCGTACGCGCCGAAGTGGGTGACGACTTCTGCGTGGGCATGCGTATCTGCGGTGACGAGTTCCACCCTGATGGCCTCAGCCACGAGGACATGAAGCAGATCGCTGCCTACTACGACGCCAC  >  NC_002947/370978‑371265
                                                                                                                                                  |                                                                                                                                             
gcAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCg                                                                                                                                             >  5:205548/1‑149 (MQ=255)
  aGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCTTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGc                                                                                                                                           >  8:345938/1‑149 (MQ=255)
       gcgcGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGa                                                                                                                                      <  7:306886/149‑1 (MQ=255)
        cgcgTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTg                                                                                                                                       >  8:23977/1‑147 (MQ=255)
             gttgcCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTAtt                                                                                                                                <  4:347328/146‑1 (MQ=255)
                 ccAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGa                                                                                                                              >  8:31199/1‑147 (MQ=255)
                  cAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGAAAg                                                                                                                           >  5:183590/1‑149 (MQ=255)
                        gggCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGAATTTCATGCGGCCTTCaaa                                                                                                                                    >  5:137250/1‑133 (MQ=255)
                        gggCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGAATTTCATGCGGCCTTCaaa                                                                                                                                    <  6:137250/134‑2 (MQ=255)
                          gcggcCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGAAAGCCGTAcgc                                                                                                                   >  3:240521/1‑149 (MQ=255)
                               cTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGAAAGCCGTACGCGCCg                                                                                                               >  2:207409/1‑148 (MQ=255)
                                  gACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGAAAGCCGTACGCGCCGAAGt                                                                                                           >  5:120838/1‑149 (MQ=255)
                                         tGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGAAAGCCGTACGCGCCGAAGTGGGTgac                                                                                                    <  4:240521/149‑1 (MQ=255)
                                         gggAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGAAAGCCGTACGCGCCGAAGTGGGTgac                                                                                                    <  5:157100/148‑1 (MQ=255)
                                           taaTGGTCGCCAGTGCACCAGCACATTATCTACCTGTGCTGTAGCCAACGGTTTCAGAAGCGTACACATATTTGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGAAAGCCGTACGCGCCGAAGTGGGTgacga                                                                                                  <  1:305407/148‑1 (MQ=255)
                                           gAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACTAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGaaa                                                                                                                            <  3:114707/123‑1 (MQ=255)
                                           gAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGaaa                                                                                                                            >  4:114707/1‑123 (MQ=255)
                                                 tCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGAAAGCCGTACGCGCCGAAGTGGGTGACGACTTCt                                                                                             <  6:161128/148‑1 (MQ=255)
                                                 tCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGAAAGCCGTACGCGCCGAAGTGGGTGACGACTTCTg                                                                                            >  4:78304/1‑149 (MQ=255)
                                                    gCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGAAAGCCGTACGCGCCGAAGTGGGTGACGACTTCTGCGt                                                                                         <  1:207409/149‑1 (MQ=255)
                                                      aGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGAAAGCCGTACGCGCCGAAGTGGGTGACGACTTCTGCGTgg                                                                                       <  7:345938/149‑1 (MQ=255)
                                                        tGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGAAAGCCGTACGCGCCGAAGTGGGTGACGACTTCTGCGTggg                                                                                      >  3:82791/1‑148 (MQ=255)
                                                            ccAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGAAAGCCGTACGCGCCGAAGTGGGTGACGACTTCTGCGTGGGCATGc                                                                                 <  7:86411/149‑1 (MQ=255)
                                                               gCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGAAAGCCGTACGCGCCGAAGTGGGTGACGACTTCTGCGTGGGCATGCGTa                                                                              <  6:183590/149‑1 (MQ=255)
                                                                cacaTGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGAAAGCCGTa                                                                                                                      >  8:13414/1‑108 (MQ=255)
                                                                cacaTGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGAAAGCCGTa                                                                                                                      <  7:13414/108‑1 (MQ=255)
                                                                   aTGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGAAAGCCGTACGCGCCGAAGTGGGTGACGACTTCTGCGTGGGCATGCGTATCTg                                                                          >  5:277681/1‑149 (MQ=255)
                                                                      aTCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGAAAGCCGTACGCGCCGAAGTGGGTGACGACTTCTGCGTGGGCATGCGTATCTGCgg                                                                       <  7:23977/149‑1 (MQ=255)
                                                                       tCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGAAAGCCGTACGCGCCGAAGTGGGTGACGACTTCTGCGTGGGCAt                                                                                   >  5:130863/1‑136 (MQ=255)
                                                                       tCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGAAAGCCGTACGCGCCGAAGTGGGTGACGACTTCTGCGTGGGCAt                                                                                   <  6:130863/136‑1 (MQ=255)
                                                                       tCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGAAAGCCGTACGCGCCGAAGTGGGTGACGACTTCTGCGTGGGCATGCGTATCTGCGGt                                                                      >  2:236689/1‑149 (MQ=255)
                                                                         gACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGAAAGCCGTACGCGCCGAAGTGGGTGACGACTTCTGCGTGGGCATGCGTATCTGCGGTGa                                                                    >  8:258569/1‑149 (MQ=255)
                                                                          aCCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGAAAGCCGTACGCGCCGAAGTGGGTGACGACTTCTGCGTGGGCATGCGTATCTGCGGTGAc                                                                   >  4:153714/1‑149 (MQ=255)
                                                                          aCCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGAAAGCCGTACGCGCCGAAGTGGGTGACGACTTCTGCGTGGGCATGCGTATCTGCGGTGAc                                                                   <  7:31199/149‑1 (MQ=255)
                                                                          aCCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGAAACCCGTACGCGCCGAAGTGGGTGACGACTTCTGCGTGGGCATGCGTATCTGCGGTGAc                                                                   <  1:288640/149‑1 (MQ=255)
                                                                             aGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGAAAGCCGTACGCGCCGAAGTGGGTGACGACTTCTGCGTGGGCATGCGTATCTGCGGTGACGAg                                                                >  6:320202/1‑149 (MQ=255)
                                                                              gTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGAAAGCCGTACGCGCCGAAGTGGGTGACGACTTCTGCGTGGGCATGCGTATCTGCGGTGACGAGt                                                               <  3:167569/149‑1 (MQ=255)
                                                                              gTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGAAAGCCGTACGCGCCGAAGTGGGTGACGACTTCTGCGTGGGCATGCGTATCTGCGGTGACGAGt                                                               >  6:6247/1‑149 (MQ=255)
                                                                                 cTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGAAAGCCGTACGCGCCGAAGTGGGTGACGACTTCTGCGTGGGCATGCGTATCTGCGGTGACGAGTTc                                                             >  8:157493/1‑148 (MQ=255)
                                                                                  tGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGAAAGCCGTACGCGCCGAAGTGGGTGACGACTTCTGCGTGGGCATGCGTATCTGCGGTGACGAGTTCCa                                                           >  5:89170/1‑149 (MQ=255)
                                                                                    gaccccTCGGGTCACCAAGCGTACCTACTCAGGGGGCGCCACATTTGCAGGCGGCATGAAATTCGGCCTTTACGTTTTGAAAGCCGTACGCGCCGTCGTGGGTGACGACTTCTGGGTGGGCATGCGTATCTGCGGTGACGAGTTcccc                                                          <  6:277681/145‑3 (MQ=255)
                                                                                     aGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGAAAGCCGTACGCGCCGAAGTGGGTGACGACTTCTGCGTGGGCATGCGTATCTGCGGTGACGAGTTCCAccc                                                        <  6:120838/149‑1 (MQ=255)
                                                                                         cACGGGTCAACAATCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTTAAAGCCGTACGAGCCGAAGTGGGTGACGACTTCTGCGTGGGCATGCGTATCTGCGGTGACGAGTTCCCACCTGAt                                                    >  3:209513/1‑149 (MQ=255)
                                                                                          aCGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGAAAGCCGTACGCGCCGAAGTGGGTGACGACTTCTGCGTGGGCATGCGTATCTGCGGTGACGAGTTCCACCCTGATg                                                   <  3:78304/149‑1 (MQ=255)
                                                                                            gggTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGAAAGCCGTACGCGCCGAAGTGGGTGACGACTTCTGCGTGGGCATGCGTATCTGCGGTGACGAGTTCCACCCTGATGGc                                                 <  4:82791/149‑1 (MQ=255)
                                                                                                 aacaaGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGAAAGCCGTACGCGCCGAAGTGGGTGACGACTTCTGCGTGGGCATGCGTATCTGCGGTGACGAGTTCCACCCTGATGGCCTCa                                             <  7:258569/148‑1 (MQ=255)
                                                                                                 aacaaGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGAAAGCCGTACGCGCCGAAGTGGGTGACGACTTCTGCGTGGGCATGCGTATCTGCGGTGACGAGTTCCACCCTGATGGCCTCa                                             <  6:205548/148‑1 (MQ=255)
                                                                                                 aacaaGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGAAAGCCGTACGCGCCGAAGTGGGTGACGACTTCTGCGTAGGCATGCGTATCTGCGGTGACGAGTTCCACCCTGATGGCCTCAg                                            >  8:334123/1‑149 (MQ=255)
                                                                                                       cGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGAAAGCCGTACGCGCCGAAGTGGGTGACGACTTCTGCGTGGGCATGCGTATCTGCGGTGACGAGTTCCACCCTGATGGCCTCAGCCACGa                                      <  3:153714/149‑1 (MQ=255)
                                                                                                          aCCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGAAAGCCGTACGCGCCGAAGTGGGTGACGACTTCTGCGTGGGCATGCGTATCTGCGGTGACGAGTTCCa                                                           <  4:301098/125‑1 (MQ=255)
                                                                                                          aCCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGAAAGCCGTACGCGCCGAAGTGGGTGACGACTTCTGCGTGGGCATGCGTATCTGCGGTGACGAGTTCCa                                                           >  3:301098/1‑125 (MQ=255)
                                                                                                            cgacgaATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGAAAGCCGTACGCGCCGAAGTGGGTGACGACTTCTGCGTGGGCATGCGTATCTGCGGTGACGAGTTCCACCCTGATGGCCTCAGCCACGAGGACa                                 <  4:17062/149‑1 (MQ=255)
                                                                                                            agaagaATGGGGCGGAAACTTTGAAGGCCTCATTAAATTCGGCCTTGAAGTATTGAAAGCAGTACGCGCCGAATTGTGTGACGACTTCTGCGTGGGAATGAGTATCTGCGGTGACGAGTTCCACCATGATTGCCTCAGCCACGAGGACa                                 >  2:65494/5‑149 (MQ=255)
                                                                                                                     gggCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGAAAGCCGTACGCGCCGAAGTGGGTGACGACTTCTGCGTGGGCATGCGTATCTGCGGTGACGAGTTCCACCCTGATGGCCTCAGCCACGAGGACATGAAGCAGa                        <  6:89170/149‑1 (MQ=255)
                                                                                                                           cACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGAAAGCCGTACGCGCCGAAGTGGGTGACGACTTCTGCGTGGGCATGCGTATCTGCGGTGACGAGTTCCACCCTGATGGCCTCAGCCACGAGGACATGAAGCAGATCGCTg                  >  7:289839/1‑149 (MQ=255)
                                                                                                                             tacttcaaGGCCGAATTTCATGCGGCCTTGAAGTATTGAAAGCCGTACGCGCCGAAGTGGGTGACGACTTCTGCGTGGGCATGCGTATCTGCGGTGACGAGTTCCACCCTGATGGCCTCAGCCACGAGGACATGAAGCAGATCGCTGcc                >  1:79357/7‑149 (MQ=255)
                                                                                                                             ccTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGAAAGCCGTACGCGCCGAAGTGGGTGACGACTTCTGCGTGGGCATGCGTATCTGCGGTGACGAGTTCCACCCTGATGGCCTCAGCCACGAGGACATGAAGCAGATCGCTGcc                >  5:353219/1‑149 (MQ=255)
                                                                                                                                ttcaaGGCCGAATTTCATGCGGCCTTGAAGTATTGAAAGCCGTACGCGCCGAAGTGGGTGACGACTTCTGCGTGGGCATGCGTATCTGCGGTGACGAGTTCCACCCTGAt                                                    <  7:221326/107‑1 (MQ=255)
                                                                                                                                ttcaaGGCCGAATTTCATGCGGCCTTGAAGTATTGAAAGCCGTACGCGCCGAAGTGGGTGACGACTTCTGCGTGGGCATGCGTATCTGCGGTGACGAGTTCCACCCTGAt                                                    >  8:221326/4‑110 (MQ=255)
                                                                                                                                   aaGGCCGAATTTCATGCGGCCTTGAAGTATTGAAAGCCGTACGCGCCGAAGTGGGTGACGACTTCTGCGTGGGCATGCGTATCTGCGGTGACGAGTTCCACCCTGATGGCCTCAGCCACGAGGACATGAAGCAGATCGCTGCCtactac          >  7:307098/1‑149 (MQ=255)
                                                                                                                                           aTGAAATTCGGCCTTGAAGTATTGAAAGCCGTACGCGCCGAAGTGGGTGACGACTTCTGCGTGGGCATGCGTATCTGCGGTGACGAGTTCCACCCTGATGGCCTCAGCCACGAGGACATGAAGCAGATCGCTGCCTACTACGACGCCAc  <  8:289839/149‑1 (MQ=255)
                                                                                                                                             gAAATTCGGCCTTGAAGTATTGAAAGCCGTACGCGCCGAAGTGGGTGACGACTTCTGCGTGGGCATGCGTATCTGCGGTGACGAGTTCCACCCTGATGGCCTCAGCCACGAGGACATGAAGCAGATCGCTGCctacta           >  3:162936/1‑138 (MQ=255)
                                                                                                                                             gAAATTCGGCCTTGAAGTATTGAAAGCCGTACGCGCCGAAGTGGGTGACGACTTCTGCGTGGGCATGCGTATCTGCGGTGACGAGTTCCACCCTGATGGCCTCAGCCACGAGGACATGAAGCAGATCGCTGCctacta           <  4:162936/138‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                             
GCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTTGAAGTATTGAAAGCCGTACGCGCCGAAGTGGGTGACGACTTCTGCGTGGGCATGCGTATCTGCGGTGACGAGTTCCACCCTGATGGCCTCAGCCACGAGGACATGAAGCAGATCGCTGCCTACTACGACGCCAC  >  NC_002947/370978‑371265

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