Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,127,417 C→T 6.6% G712S (GGT→AGT)  gcvP‑I ← glycine dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,127,4170CT 6.6% 174.5 / 2.6 61G712S (GGT→AGT) gcvP‑Iglycine dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (27/29);  new base T (2/2);  total (30/31)
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.42e-01

ATCCATGTGATGGGCAGGATGCTGGCGCTACCGAACGGCGCAGCGCATACCGCGCCCTGGGTGTTTTCCAGCTGTGCATGGCCTGGCAGGAACGGCGCCAGGTGCGACTTGACGCCAATCGGGCCAACACCCGGGCCACCACCAC‑CGTGCGGGATGCAGAAGGTCTTGTGCAGGTTCAGGTGGGAAACGTCGCCGCCGAACTTGCCTGGCGCGCACAGTCCGACCATGGCGTTCATGTTGGCGCCATCGATGTACACCTGGCCGCCGTTGTCGTGGATGATCGCG  >  NC_002947/1127278‑1127562
                                                                                                                                           |                                                                                                                                                  
aTCCATGTGATGGGCAGGATGCTGGCGCTACCGAACGGCGCAGCGCATACCGCGCCCTGGGTGTTTTCCAGCTGTGCATGGCCTGGCAGGAACGGCGCCAGGTGCGACTTGACGCCAATCGGGCCAACACCCGGGCCACCACCAC‑CGTg                                                                                                                                          >  5:274401/1‑149 (MQ=255)
 tCCATGTGATGGGCAGGATGCTGGCGCTACCGAACGGCGCAGCGCATACCGCGCCCTGGGTGTTTTCCAGCTGTGCATGGCCTGGCAGGAACGGCGCCAGGTGCGACTTGACGCCAATCGGGCCAACACCCGGGCCACCACCAC‑CGTGc                                                                                                                                         <  2:327930/149‑1 (MQ=255)
    aTGTGATGGGCAGGATGCTGGCGCTACCGAACGGCGCAGCGCATACCGCGCCCTGGGTGTTTTCCAGCTGTGCATGGCCTGGCAGGAACGGCGCCAGGTGCGACTTGACGCCAATCGGGCCAACACCCGGGCCACCACCAC‑CGTGCggg                                                                                                                                      <  1:276388/149‑1 (MQ=255)
             gCAGGATGCTGGCGCTACCGAACGGCGCAGCGCATACCGCGCCCTGGGTGTTTTCCAGCTGTGCATGGCCTGGCAGGAACGGCGCCAGGTGCGACTTGACGCCAATCGGGCCAACACCCGGGCCACCACCAC‑CGTGCGGGATGCAGAAg                                                                                                                             >  8:100797/1‑149 (MQ=255)
             gCAGGATGCTGGCGCTACCGAACGGCGCAGCGCATACCGCGCCCTGGGTGTTTTCCAGCTGTGCATGGCCTGGCAGGAACGGCGCCAGGTGCGACTTGACGCCAATCGGGCCAACACCCGGGCCACCACCAC‑CGTGCGGGATGCAGAAg                                                                                                                             <  5:201350/149‑1 (MQ=255)
                ggATGCTGGCGCTACCGAACGGCGCAGCGCATACCGCGCCCTGGGTGTTTTCCAGCTGTGCATGGCCTGGCAGGAACGGCGCCAGGTGCGACTTGACGCCAATCGGGCCAACACCCGGGCCACCACCAC‑CGTGCGGGATGCAGAAGGTc                                                                                                                          <  1:155967/149‑1 (MQ=255)
                  aTGCTGGCGCTACAGAACGGCGCAGCGCATACCGCGCCCTGGGTGTTTTCCAGCTGTGCATGGCCTGGCAGGAACGGCGCCAGGTGCGACTTGACGCCAATCGGGCCAACACCCGGGCCACCACCAC‑CGTGCGGGATGCAGAAGGTCtt                                                                                                                        <  2:304708/149‑1 (MQ=255)
                     cTGGCGCTACCGAACGGCGCAGCGCATACCGCGCCCTGGGTGTTTTCCAGCTGTGCATGGCCTGGCAGGAACGGCGCCAGGTGCGACTTGACGCCAATCGGGCCAACACCCGGGCCACCACCAC‑CGTGCGGGATGCAGAAGGTCTtgt                                                                                                                      <  6:274401/148‑1 (MQ=255)
                     cTGGCGCTACCGAACGGCGCAGCGCATACCGCGCCCTGGGTGTTTTCCAGCTGTGCATGGCCTGGCAGGAACGGCGCCAGGTGCGACTTGACGCCAATCGGGCCAACACCCGGGCCACCACCAC‑CGTGCGGGATGCAGAAGGTCTtgt                                                                                                                      <  7:100797/148‑1 (MQ=255)
                        gcgcTACCGAACGGCGCAGCGCATACCGCGCCCTGGGTGTTTTCCAGCTGTGCATGGCCTGGCAGGAACGGCGCCAGGTGCGACTTGACGCCAATCGGGCCAACACCCGGGCCACCACCAC‑CGTGCGGGATGCAGAAGGTCTTGTGCa                                                                                                                   <  1:179836/148‑1 (MQ=255)
                         cgcTACCGAACGGCGCAGCGCATACCGCGCCCTGGGTGTTTTCCAGCTGTGCATGGCCTGGCAGGAACGGCGCCAGGTGCGACTTGACGCCAATCGGGCCAACACCCGGGCCACCACCAC‑CGTGCGGGATGCAGAAGGTCTTGTGCAgg                                                                                                                 >  5:253630/1‑149 (MQ=255)
                              ccGAACGGCGCAGCGCATACCGCGCCCTGGGTGTTTTCCAGCTGTGCATGGCCTGGCAGGAACGGCGCCAGGTGCGACTTGACGCCAATCGGGCCAACACCCGGGCCACCACCAC‑CGTGCGGGATGCAGAAGGTCTTGTGCAGGTTCAg                                                                                                            >  8:220911/1‑149 (MQ=255)
                              ccGAACGGCGCAGCGCATACCGCGCCCTGGGTGTTTCCCAGCTGTGCATGGCCTGGCAGGAACGGCGCCAGGTGCGACTTGACGCCAATCGGGCCAACACCCGGGCCACTACCAC‑CGTGCGGGATGCAGAAGGTCTTGTGCAGGTTCAg                                                                                                            >  6:250198/1‑149 (MQ=255)
                              ccGAACGGCGCAGCGCATACCGCGCCCTGGGTGTTTCCCAGCTGTGCATGGCCTGGCAGGAACGGCGCCAGGTGCGACTTGACGCCAATCGGGCCAACACCCGGGCCACTACCAC‑CGTGCGGGATGCAGAAGGTCTTGTGCAGGTTCAg                                                                                                            >  6:244834/1‑149 (MQ=255)
                                 aaCGGCGCAGCGCATACCGCGCCCTGGGTGTTTTCCAGCTGTGCATGGCCTGGCAGGAACGGCGCCAGGTGCGACTTGACGCCAATCGGGCCAACACCCGGGCCACCACCACACGTGCGGGATGCAGAAGGTCTTGTGCAGGTTCAGGt                                                                                                          <  1:249606/149‑1 (MQ=255)
                                     gcgcagcgcaTACCGCGCCCTGGGTGTTTTCCAGCTGTGCATGGCCTGGCAGGAACGGCGCCAGGTGCGACTTGACGCCAATCGGGCCAACACCCGGGCCACCACCAC‑CGTGCGGGATGCAGAAGGTCTTGTGCAGGTTCAGGTGGGa                                                                                                      >  2:347006/1‑148 (MQ=255)
                                       gcagcgcaTACCGCGCCCTGGGTGTTTTCCAGCTGTGCATGGCCTGGCAGGAACGGCGCCAGGTGCGACTTGACGCCAATCGGGCCAACACCCGGGCCACCACCAC‑CGTGCGGGATGCAGAAGGTCTTGTGCAGGTTCAGGTGGGAAAc                                                                                                   <  6:214592/149‑1 (MQ=255)
                                         agcgcaTACCGCGCCCTGGGTGTTTTCCAGCTGTGCATGGCCTGGCAGGAACGGCGCCAGGTGCGACTTGACGCCAATCGGGCCAACACCCGGGCCACCACCAC‑CGTGCGGGATGCAGAAGGTCTTGTGCAGGTTCAGGTGGGAAACGt                                                                                                 >  3:154781/1‑149 (MQ=255)
                                            gcaTACCGCGCCCTGGGTGTTTTCCAGCTGTGCATGGCCTGGCAGGAACGGCGCCAGGTGCGACTTGACGCCAATCGGGCCAACACCCGGGCCACCACCAC‑CGTGCGGGATg                                                                                                                                   <  6:206459/112‑1 (MQ=255)
                                            gcaTACCGCGCCCTGGGTGTTTTCCAGCTGTGCATGGCCTGGCAGGAACGGCGCCAGGTGCGACTTGACGCCAATCGGGCCAACACCCGGGCCACCACCAC‑CGTGCGGGATg                                                                                                                                   >  5:206459/1‑112 (MQ=255)
                                            gcaTACCGCGCCCTGGGTGTTTTCCAGCTGTGCATGGCCTGGCAGGAACGGCGCCAGGTGCGACTTGACGCCAATCGGGCCAACACCCGGGCCACCACCAC‑CGTGCGGGATGCAGAAGGTCTTGTGCAGGTTCAGGTGGGAAACGTcgc                                                                                              >  1:78526/1‑149 (MQ=255)
                                              aTACCGCGCCCTGGGTGTTTTCCAGCTGTGCATGGCCTGGCAGGAACGGCGCCAGGTGCGACTTGACGCCAATCGGGCCAACACCCGGGCCACCACCAC‑CGTGCGGGATGCAGAAGGTCTTGTGCAGGTTCAGGTGGGAAACGTcgccg                                                                                            >  8:82160/1‑149 (MQ=255)
                                               tACCGCGCCCTGGGTGTTTTCCAGCTGTGCATGGCCTGGCAGGAACGGCGCCAGGTGCGACTTGACGCCAATCGGGCCAACACCCGGGCCACCACCAC‑CGTGCGGGATGCAGAAGGTCTTGTGCAGGTTCAGGTGGGAAACGTcgccgc                                                                                           >  8:186886/1‑149 (MQ=255)
                                                   gcgcCCTGGGTGTTTTCCAGCTGTGCATGGCCTGGCAGGAACGGCGCCAGGTGCGACTTGACGCCAATCGGGCCAACACCCGGGCCACCACCAC‑CGTGCGGGATGCAGAAGGTCTTGTGCAGGTTCAGGTGGGAAACGTCGCCGCCGaa                                                                                       >  7:205596/1‑149 (MQ=255)
                                                          gggTGTTTTCCAGCTGTGCATGGCCTGGCAGGAACGGCGCCAGGTGCGACTTGACGCCAATCGGGCCAACACCCGGGCCACCACCAC‑CGTGCGGGATGCAGAAGGTCTTGTGCAGGTTCAGGTGGGAAACGTCGCCGCCGAACTTGCCt                                                                                >  4:120229/1‑149 (MQ=255)
                                                           ggTGTTTTCCAGCTGTGCATGGCCTGGCAGGAACGGCGCCAGGTGCGACTTGACGCCAATCGGGCCAACACCCGGGCCACCACCAC‑CGTGCGGGATGCAGAAGGTCTTGTGCAGGTTCAGGTGGGAAACGTCGCCGCCGAACTTGCCTg                                                                               >  7:39754/1‑149 (MQ=255)
                                                            gtgtTTTCCAGCTGTGCATGGCCTGGCAGGAACGGCGCCAGGTGCGACTTGACGCCAATCGGGCCAACACCCGGGCCACCACCAC‑CGTGCGGGATGCAGAAGGTCtt                                                                                                                        <  7:81165/107‑1 (MQ=255)
                                                            gtgtTTTCCAGCTGTGCATGGCCTGGCAGGAACGGCGCCAGGTGCGACTTGACGCCAATCGGGCCAACACCCGGGCCACCACCAC‑CGTGCGGGATGCAGAAGGTCtt                                                                                                                        >  8:81165/1‑107 (MQ=255)
                                                             tgtTTTCCAGCTGTGCATGGCCTGGCAGGAACGGCGCCAGGTGCGACTTGACGCCAATCGGGCCAACACCCGGGCCACCACCAC‑CGTGCGGGATGCAGAAGGTCTTGTGCAGGTTCAGGTGGGa                                                                                                      <  6:357783/124‑1 (MQ=255)
                                                             tgtTTTCCAGCTGTGCATGGCCTGGCAGGAACGGCGCCAGGTGCGACTTGACGCCAATCGGGCCAACACCCGGGCCACCACCAC‑CGTGCGGGATGCAGAAGGTCTTGTGCAGGTTCAGGTGGGAAACGTCGCCGCCGAACTTGCCTGgc                                                                             >  5:132267/1‑149 (MQ=255)
                                                               ttttCCAGCTGTGCATGGCCTGGCAGGAACGGCGCCAGGTGCGACTTGTCGCCAATCGGGCCAACACCCGGGCCACCACCAC‑CGTGCGGGATGCAGAAGGTCTTGTGCAGGTTCAGGTGGGAAACGTCGCCGCCGAAGTTGCCTGgcgc                                                                           >  2:145487/1‑149 (MQ=255)
                                                                tttCCAGCTGTGCATGGCCTGGCAGGAACGGCGCCAGGTGCGACTTGACGCCAATCGGGCCAACACCCGGGCCACCACCAC‑CGTGCGGGATGCAGAAGGTCTTGTGCAGGTTCAGGTGGGAAACGTCGCCGCCGAACTTGCCTGgcgcg                                                                          >  6:70544/1‑149 (MQ=255)
                                                                  tCCAGCTGTGCATGGCCTGGCAGGAACGGCGCCAGGTGCGACTTGACGCCAATCGGGCCAACACCCGGGCCACCACCAC‑CGTGCGGGATGCAGAAGGTCTTGTGCAGGTTCAGGTGGGAAACGTCGCCGCCGAACTTGCCTGgcgcgc                                                                         <  3:120229/148‑1 (MQ=255)
                                                                       cTGTGCATGGCCTGGCAGGAACGGCGCCAGGTGCGACTTGACGCCAATCGGGCCAACACCCGGGCCACCACCAC‑CGTGCGGGATGCAGAAGGTCTTGTGCAGGTTCAGGTGGGAAACGTCGCCGCCGAACTTGCCTGGCGCGCACAGTc                                                                   <  8:253034/149‑1 (MQ=255)
                                                                           gCATGGCCTGGCAGGAACGGCGCCAGGTGCGACTTGACGCCAATCGGGCCAACACCCGGGCCACAACCAC‑CGTGCGGGATGCAGAAGGTCTTGTGCAGGTTCAGGTGGGAAACGTCGCAGCCGAACTTGCCTGGCGCGCACAGTCCGAc                                                               >  8:6168/1‑149 (MQ=255)
                                                                            cATGGCCTGGCAGGAACGGCGCCAGGTGCGACTTGACGCCAATCGGGCCAACACCCGGGCCACCACCAC‑CGTGCGGGATGCAGAAGGTCTTGTGCAGGTTCAGGTGGGAAACGTCGCCGCCGAACTTGCCTGGCGCGCACAGTCCGAc                                                               >  2:216084/1‑148 (MQ=255)
                                                                              tGGCCTGGCAGGAACGGCGCCAGGTGCGACTTGACGCCAATCGGGCCAACACCCGGGCCACCACCAC‑CGTGCGGGATGCAGAAGGTCTTGTGCAGGTTCAGGTGGGAAACGTCGCCGCCGAACTTGCCTGGCGCGCACAGTCCGACCAt                                                            >  6:11257/1‑149 (MQ=255)
                                                                                gCCTGGCAGGAACGGCGCCAGGTGCGACTTGACGCCAATCGGGCCAACACCCGGGCCACCACCAC‑CGTGCGGGATGCAGAAGGTCTTGTGCAGGTTCAGGTGGGAAACGTCGCCGCCGAACTTGCCTGGCGCGCACAGTCCGACCATg                                                           <  6:132267/148‑1 (MQ=255)
                                                                                       aGGAACGGCGCCAGGTGCGACTTGACGCCAATCGGGCCAACACCCGGGCCACCACCAC‑CGTGCGGGATGCAGAAGGTCTTGTGCAGGTTCAGGTGGGAAACGTCGCCGCCGAACTTGCCTGGCGCGCACAGTCCGACCATGGCGTTCAt                                                   <  8:205596/149‑1 (MQ=255)
                                                                                                 ccAGGTGCGACTTGACGCCAATCGGGCCAACACCCGGGCCACCACCAC‑CGTGCGGGATGCAGAAGGTCTTGTGCAGGTTCAGGTGGGAAACGTCGCCGCCGAACTTGCCTGGCGCGCACAGTCCGACCATGGCGTTCATGTTGGCGCCa                                         >  4:181041/1‑149 (MQ=255)
                                                                                                  cAGGTGCGACTTGACGCCAATCGGGCCAACACCCGGGCCACCACCAC‑CGTGCGGGATGCAGAAGGTCTTGTGCAGGTTCAGGTGGGAAACGTCGCCGCCGAACTTGCCTGGCGCGCACAGTCCGACCATGGCGTTCATGTTGGCGCCAt                                        <  1:216084/149‑1 (MQ=255)
                                                                                                      tGCGACTTGACGCCAATCGGGCCAACACCCGGGCCACCACCAC‑CGTGCGGGATGCAGAAGGTCTTGTGCAGGTTCAGGTGGGAAACGTCGCCGCCGAACTTGCCTGGCGCGCACAGTCCGACCATGGCGTTCATGTTGGCGCCATCGAt                                    <  4:154781/149‑1 (MQ=255)
                                                                                                      tGCGACTTGACGCCAATCGGGCCAACACCCGGGCCACCACCAC‑CGTGCGGGATGCAGAAGGTCTTGTGCAGGTTCAGGTGGGAAACGTCGCCGCCGAACTTGCCTGGCGCGCACAGTCCGACCATGGCGTTCATGTTGGCGCCATCGAt                                    <  5:11257/149‑1 (MQ=255)
                                                                                                         gACTTGACGCCAATCGGGCCAACACCCGGGCCACCACCAC‑CGTGCGCGATGCAGAAGGTCTTGTGCAGGTTCAGGTGGGAAACGTCGCCGCCGAACTTGCCTGGCGCGCACAGTCCGACCATGGCGTTCATGTTGGCGCCATCGATGTa                                 >  3:217085/1‑149 (MQ=255)
                                                                                                          aCTTGACGCCAATCGGGCCAACACCCGGGCCACCACCAC‑CGTGCGGGATGCAGAAGGTCTTGTGCAGGTTCAGGTGGGAAACGTCGCCGCCGAACTTGCCTGGCGCGCACAGTCCGACCATGGCGTTCATGTTGGCGCCATCGATGTac                                <  1:347006/149‑1 (MQ=255)
                                                                                                          aCTGGACGCCAATCCGGCCAACACCCGGGCCACCACCAC‑CTCGCGGTATGCAGAAGGTCTTGTGCAGGTGCAGGCGGGAAACGTCGCCGCCGAACTTGCCTGGCGCGCACAGTCCGACCATGGCGTTCATGTGGGCGCCATCGATGTac                                <  6:253630/149‑1 (MQ=255)
                                                                                                              gACGCCAATCGGGCCAACACCCGGGCCACCACCAC‑CGTGCGGGATGCAGAAGGTCTTGTGCAGGTTCAGGTGGGAAACGTcgccg                                                                                            <  3:29883/85‑1 (MQ=255)
                                                                                                              gACGCCAATCGGGCCAACACCCGGGCCACCACCAC‑CGTGCGGGATGCAGAAGGTCTTGTGCAGGTTCAGGTGGGAAACGTcgccg                                                                                            >  4:29883/1‑85 (MQ=255)
                                                                                                              gACGCCAATCGGGCCAACACCCGGGCCACCACCAC‑CGTGCGGGATGCAGAAGGTCTTGTGCAGGTTCAGGTGGGAAACGTCGCCGCCGAACTTGCCTGGCGCGCACAGTCCGACCATGGCGTTCATGTTGGCGCCATCGATGTACACCt                            >  2:32076/1‑149 (MQ=255)
                                                                                                                    gggagccccAACACCCGGGCCACCACCAC‑CGTGCGGGATGCAGAAGGTCTTGTGCAGGTTCAGGTGGGAAACGTCGCCGCCGAACTTGCCTGGCGCGCACAGTCCGACCATGGCGTTCATGTTGGCGCCATCGATGTACACCTGgcc                        <  1:130239/140‑1 (MQ=255)
                                                                                                                     aTCGGGCCAACACCCGGGCCACCACCAC‑CGTGCGGGATGCAGAAGGTCTTGTGCAGGTTCAGGTGGGAAACGTCGCCGCCGAACTTGCCTGGCGCGCACAGTCCGACCATGGCGTTCATGTTGGCGCCATCGATGTACACCTGgccgcc                     >  4:328415/1‑149 (MQ=255)
                                                                                                                      tCGGGCCAACACCCGGGCCACCACCAC‑CGTGCGGGATGCAGAAGGTCTTGTGCAGGTTCAGGTGGGAAACGTCGCCGCCGAACTTGCCTGGCGCGCACAGTCCGACCATGGCGTTCATGTTGGCGCCATCGATGTACACCTGGccgccg                    <  1:332737/149‑1 (MQ=255)
                                                                                                                        gggCCAACACCCGGGCCACCACCAC‑CGTGCGGGATGCAGAAGGTCTTGTGCAGGTTCAGGTGGGAAACGTCGCCGCCGAACTTGCCTGGCGCGCACAGTCCGACCATGGCGTTCATGTTGGCGCCATCGATGTACACCTGGCCGCCGtt                  <  7:82160/149‑1 (MQ=255)
                                                                                                                          gCCAACACCCGGGCCACCACCAC‑CGTGCGGGATGCAGAAGGTCTTGTGCAGGTTCAGGTGGGAAACGTCGCCGCCGAACTTGCCTGGCGCGCACAGTCCGACCATGGCGTTCATGTTGGCGCCATCGATGTACACCTGGCCGCCGTTGt                <  2:78526/149‑1 (MQ=255)
                                                                                                                            cAACACCCGGGCCACCACCAC‑CGTGCGGGATGCAGAAGGTCTTGTGCAGGTTCAGGTGGGAAACGTCGCCGCCGAACTTGCCTGGCGCGCACAGTCCGACCATGGCGTTCATGTTGGCGCCATCGATGTACACCTGgcc                        >  2:130239/1‑139 (MQ=255)
                                                                                                                            cAACACCCGGGCCACCACCAC‑CGTGCGGGATGCAGAAGGTCTTGTGCAGGTTCAGGTGGGAAACGTCGCCGCCGAACTTGCCTGGCGCGCACAGTCCGACCATGGCGTTCATGTTGGCGCCATCGATGTACACCTGGCCGCCGTTGTCg              >  3:227407/1‑149 (MQ=255)
                                                                                                                             aaCACCCGGGCCACTACCAC‑CGTGCGGGATGCAGAAGGTCTTGTGCAGGTTCAGGTGGGAAACGTCGCCGCCGAACTTGCCTGGCGCGCACAGTCCGACCATGGCGTTCATGTTGGCGCCATCGATGTACACCTGGCCGCCGTTGTCGt             <  5:244834/149‑1 (MQ=255)
                                                                                                                             aaCACCCGGGCCACTACCAC‑CGTGCGGGATGCAGAAGGTCTTGTGCAGGTTCAGGTGGGAAACGTCGCCGCCGAACTTGCCTGGCGCGCACAGTCCGACCATGGCGTTCATGTTGGCGCCATCGATGTACACCTGGCCGCCGTTGTCGt             <  5:250198/149‑1 (MQ=255)
                                                                                                                             aaCACCCGGGCCACCACCAC‑CGTGCGGGATGCAGAAGGTCTTGAGCAGGTTCAGCTGGGAAACGTCGCCGCCGAACTTGCCTGGCGCGCACAGTCCGACCATGGCGTTCATGTTGGCGCCATCGATGTACACCTGGCCGCCGTTGTCGt             >  5:146325/1‑149 (MQ=255)
                                                                                                                              acacCCGGGCCACCACCAC‑CGTGCGGGATGCAGAAGGTCTTGTGCAGGTTCAGGTGGGAAACGTCGCCGCCGAACTTGCCTGGCGCGCACAGTCCGACCATGGCGTTCATGTTGGCGCCATCGATGTACACCTGGCCGCCGTTGTCGTg            <  7:220911/149‑1 (MQ=255)
                                                                                                                                acCCGGGCCACCACCAC‑CGTGCGGGATGCAGAAGGTCTTGTGCAGGTTCAGGTGGGAAACGTCGCCGCCGAACTTGCCTGGCGCGCACAGTCCGACCATGGCGTTCATGTTGGCGCCATCGATGTACACCTGGCCGCCGTTGTCGTGga          <  7:7575/149‑1 (MQ=255)
                                                                                                                                           caccac‑cGTGCGGGATGCAGAAGGTCTTGTGCAGGTTCAGGTGGGAAACGTCGCCGCCGAACTTGCCTGGCGCGCACAGTCCGACCATGGCGTTCATGTTGGCGCCATCGATGTACACCTGGCCGCCGTTGTCGTGGATGATcgcg  >  3:67348/1‑146 (MQ=255)
                                                                                                                                           |                                                                                                                                                  
ATCCATGTGATGGGCAGGATGCTGGCGCTACCGAACGGCGCAGCGCATACCGCGCCCTGGGTGTTTTCCAGCTGTGCATGGCCTGGCAGGAACGGCGCCAGGTGCGACTTGACGCCAATCGGGCCAACACCCGGGCCACCACCAC‑CGTGCGGGATGCAGAAGGTCTTGTGCAGGTTCAGGTGGGAAACGTCGCCGCCGAACTTGCCTGGCGCGCACAGTCCGACCATGGCGTTCATGTTGGCGCCATCGATGTACACCTGGCCGCCGTTGTCGTGGATGATCGCG  >  NC_002947/1127278‑1127562

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