Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,331,048 G→C 6.3% F318L (TTC→TTG) ‡ nspC ← carboxynorspermidine/carboxyspermidine decarboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,331,0480GC 6.3% 206.3 / 4.1 63F318L (TTC→TTG) ‡nspCcarboxynorspermidine/carboxyspermidine decarboxylase
Reads supporting (aligned to +/- strand):  ref base G (26/33);  new base C (2/2);  total (28/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.23e-01

TCTTACGACAGGCTGGAAACGTAGTCTTCGAAACCGAACTCGCGGACGACTTCGACGCTGCCGTCGAGCTGCTTGACCACGATGGAAGGCATCTTCAGCCCGTTGAACCAGTTTTTCTTGACCATGGTGTAACCCGCCGTGTCGATGAACGACA‑GGCGATCGCCAATGGCCAGCGGACGGTCGAATTGGTACTCGCCGAAGATGTCACCGGCCAGGCATGACTTGCCGCAGACCATGTAGGTATGCTCGCCGTCATTGGGGGCCATCTTGGCGTTGAGGCGGTAGATCAGCAG  >  NC_002947/3330902‑3331194
                                                                                                                                                  |                                                                                                                                                   
tCTTACGACAGGCTGGAAACGTAGTCTTCGAAACCGAACTCGCGGACGACTTCGACGCTGCCGTCGAGCTGCTTGACCACGATGGAAGGCATCTTCAGCCCGTTGAACCAGTTTTTCTTGACCATGGTGTAACCCGCCGTGTCGATGaa                                                                                                                                                   >  8:411494/1‑149 (MQ=255)
    aCGACAGGCTGGAAACGTAGTCTTCGAAACCGAACTCGCGGACGACTTCGACGCTGCCGTCGAGCTGCTTGACCACGATGGAAGGCATCTTCAGCCCGTTGAACCAGTTTTTCTTGACCATGGTGTAACCCGCCGTGTCGATGAAc                                                                                                                                                  <  3:141043/146‑1 (MQ=255)
     cGACAGGCTGGAAACGTAGTCTTCGAAACCGAACTCGCGGACGACTTCGACGCTGCCGTCGAGCTGCTTGACCACGATGGAAGGCATCTTCAGCCCGTTGAACCAGTTTTTCTTGACCATGGTGTAACCCGCCGTGTCGATGAACGACa                                                                                                                                              <  4:211027/149‑1 (MQ=255)
     cGACAGGCTGGAAACGTAGTCTTCGAAACCGAACTCGCGGACGACTTCGACGCTGCCGTCGAGCTGCTTGACCACGATGGAAGGCATCTTCAGCCCGTTGAACCAGTTTTTCTTGACCATGGTGTAACCCGCCGTGTCGATGAACGACa                                                                                                                                              <  6:120809/149‑1 (MQ=255)
        cAGGCTGGAAACGTAGTCTTCGAAACCGAACTCGCGGACGACTTCGACGCTGCCGTCGAGCTGCTTGACCACGATGGAAGGCATCTTCAGCCCGTTGAACCAGTTTTTCTTGACCATGGTGTAACCCGCCGTGTCGTTCATCGACacgg                                                                                                                                           <  2:282790/149‑4 (MQ=255)
         aGGCTGGAAACGTAGTCTTCGAAACCGAACTCGCGGACGACTTCGACGCTGCCGTCGAGCTGCTTGACCACGATGGAAGGCATCTTCAGCCCGTTGAACCAGTTTTTCTTGACCATGGTGTAACCCGCCGTGTCGATGAACGACA‑GGc                                                                                                                                          >  8:52632/1‑148 (MQ=255)
             tGGAAACGTAGTCTTCGAAACCGAACTCGCGGACGACTTCGACGCTGCCGTCGAGCTGCTTGACCACGATGGAAGGCATCTTCAGCCCGTTGAACCAGTTTTTCTTGACCATGGTGTAACCCGCCGTGTCGATGAACGACA‑GGCGATCg                                                                                                                                     >  3:132701/1‑149 (MQ=255)
               gAAACGTAGTCTTCGAAACCGAACTTGCGGACGACTTCGACGCTGCCGTCGAGCTGCTTGACCACGATGGAAGGCATCTTCAGCCCGTTGAACCAGTTTTTCTTGACCATGGTGTAACCCGCCGTGTCGATGAACGACA‑GGCGATCGc                                                                                                                                    <  7:14280/148‑1 (MQ=255)
               gAAACGTAGTCTTCGAAACCGAACTCGCGGACGACTTCGACGCTGCCGTCGAGCTGCTTGACCACGATGGAAGGCATCTTCAGCCCGTTGAACCAGTTTTTCTTGACCATGGTGTAACCCGCCGTGTCGATGAACGACA‑GGCGATCGc                                                                                                                                    <  7:52632/148‑1 (MQ=255)
                  aCGTAGTCTTCGAAACCGAACTCGCGGACGACTTCGACGCTGCCGTCGAGCTGCTTGACCACGATGGAAGGCATCTTCAGCCCGTTTAACCAGTTTTTCTTGACCATGGTGTAAACCGCCGTGTCGATGAACGACA‑GGCGATCGCCAAt                                                                                                                                >  1:392809/1‑149 (MQ=255)
                   cGTAGTCTTCGAAACCGAACTCGCGGACGACTTCGACGCTGCCGTCGAGCTGCTTGACCACGATGGAAGGCATCTTCAGCCCGTTGAACCAGTTTTTCTTGACCATGGTGTAACCCGCCGTGTCGATGAACGACA‑GGCGATCGCCAATg                                                                                                                               <  4:372919/149‑1 (MQ=255)
                        tCTTCGAACCCGAACTCGCGGACGACTTCGACGCTGCCGTCGAGCTGCTTGACCACGATGGAAGGCATCTTCAGCCCGTTGAACCAGTTTTTCTTGACCATGGTGTAACCCGCCGTGTCGATGAACGACA‑GGCGATCGCCAATGGCCAg                                                                                                                          <  4:235227/149‑1 (MQ=255)
                          ttCGAAACCGAACTCGCGGACGACTTCGACGCTGCCGTCGAGCTGCTTGACCACGATGGAAGGCATCTTCAGCCCGTTGAACCAGTTTTTCTTGACCATGGTGTAACCCGCCGTGTCGATGAACGACA‑GGCGATCGCCAATGGCCAGCg                                                                                                                        >  5:240101/1‑149 (MQ=255)
                                     aCTCGCGGACGACTTCGACGCTGCCGTCGAGCTGCTTGACCACGATGGAAGGCATCTTCAGCCCGTTGAACCAGTTTTTCTTGACCATGGTGTAACCCGCCGTGTCGATGAACGACA‑GGCGATCGCCAATGGCCAGCGGACGGTCGAAt                                                                                                             <  2:407555/149‑1 (MQ=255)
                                         gcgGACGACTTCGACGCTGCCGTCGAGCTGCTTGACCACGATGGAAGGCATCTTCAGCCCGTTGAACCAGTTTTTCTTGACCATGGTGTAACCCGCCGTGTCGATGaa                                                                                                                                                   <  4:57966/108‑1 (MQ=255)
                                         gcgGACGACTTCGACGCTGCCGTCGAGCTGCTTGACCACGATGGAAGGCATCTTCAGCCCGTTGAACCAGTTTTTCTTGACCATGGTGTAACCCGCCGTGTCGATGaa                                                                                                                                                   >  3:57966/1‑108 (MQ=255)
                                         gcgGACGACTTCGACGCTGCCGTCGAGCTGCTTGACCACGATGGAAGGCATCTTCAGCCCGTTGAACCAGTTTTTCTTGACCATGGTGTAACCCGCCGTGTCGATGAACGACA‑GGCGATCGCCAATGGCCAGCGGACGGTCGAATTg                                                                                                           <  1:167389/147‑1 (MQ=255)
                                         gcgGACGACTTCGACGCTGCCGTCGAGCTGCTTGACCACGATGGAAGGCATCTTCAGCCCGTTGAACCAGTTTTTCTTGACCATGGTGTAACCCGCCGTGTCGATGAACGACA‑GGCGATCGCCAATGGCCAGCGGACGGTCGAATTGGt                                                                                                         <  5:249344/149‑1 (MQ=255)
                                                 cTTCGACGCTGCCGTCGAGCTGCTTGACCACGATGGAAGGCATCTTCAGCCCGTTGAACCAGTTTTTCTTGACCATGGTGTAACCCGCCGTGTCGATGAACGACA‑GGCGATCGCCAATGGCCAGCGGACGGTCGAATTGGTACTCGcc                                                                                                  >  2:232153/1‑148 (MQ=255)
                                                  ttCGACGCTGCCGTCGAGCTGCTTGACCACGATGGAAGGCATCTTCAGCCCGTTGAACCAGTTTTTCTTGACCATGGTGTAACCCGCCGTGTCGATGAACGACA‑GGCGATc                                                                                                                                      >  3:338790/1‑111 (MQ=255)
                                                  ttCGACGCTGCCGTCGAGCTGCTTGACCACGATGGAAGGCATCTTCAGCCCGTTGAACCAGTTTTTCTTGACCATGGTGTAACCCGCCGTGTCGATGAACGACA‑GGCGATc                                                                                                                                      <  4:338790/111‑1 (MQ=255)
                                                  ttCGACGCTGCCGTCGAGCTGCTTGACCACGATGGAAGGCATCTTCAGCCCGTTGAACCAGTTTTTCTTGACCATGGTGTAACCCGCCGTGTCGATGAACGACA‑GGCGATCGCCAATGGCCAGCGGACGGTCGAATTg                                                                                                           <  7:126789/138‑1 (MQ=255)
                                                  ttCGACGCTGCCGTCGAGCTGCTTGACCACGATGGAAGGCATCTTCAGCCCGTTGAACCAGTTTTTCTTGACCATGGTGTAACCCGCCGTGTCGATGAACGACA‑GGCGATCGCCAATGGCCAGCGGACGGTCGAATTg                                                                                                           >  8:126789/1‑138 (MQ=255)
                                                   tCGACGCTGCCGTCGAGCTGCTTGACCACGATGGAAGGCATCTTCAGCCCGTTGAACCAGTTTTTCTTGACCATGGTGTAACCCGCCGTGTCGATGAACGACA‑GGAGATCGCCAATGGCCAGCGGACGGTCGAATTGGTACTCGCCg                                                                                                 >  1:240691/1‑147 (MQ=255)
                                                    cGACGCTGCCGTCGAGCTGCTTGACCACGATGGAAGGCATCTTCAGCCCGTTGAACCAGTTTTTCTTGACCATGGTGTAACCCGCCGTGTCGATGAACGACA‑GGCGATCGCCAATGGCCAGCGGACGGTCGAATTGGTACTCGCCGa                                                                                                <  1:221772/147‑1 (MQ=255)
                                                     gACGCTGCCGTCGAGCTGCTTGACCACGATGGAAGGCATCTTCAGCCCGTTGAACCAGTTTTTCTTGACCATGGTGTAACCCGCCGTGTCGATGAACGACA‑GGCGATCGCCAATGGCCAGCGGACGGTCGAATTGGTACTCGCCGAAg                                                                                              <  7:411494/148‑1 (MQ=255)
                                                      aCGCTGCCGTCGAGCTGCTTGACCACGATGGAAGGCATCTTCAGCCCGTTGAACCAGTTTTTCTTGACCATGGTGTAACCCGCCGTGTCGATGAACGACA‑GGCGATCGCCAATGGCCAGCGGACGGTCGAATTGGTACTCGCCGAAGAt                                                                                            <  6:240101/149‑1 (MQ=255)
                                                           gCCGTCGAGCTGCTTGACCACGATGGAAGGCATCTTCAGCCCGTTGAACCAGTTTTTCTTGACCATGGTGTAACCCGCCGTGTCGATGAACGACA‑GGCGATCGCCAATGGCCAGCGGACGGTCGAATTGGTACTCGCCGAAGATGTCAc                                                                                       >  8:109826/1‑149 (MQ=255)
                                                              gTCGAGCTGCTTGACCACGATGGAAGGCATCTTCAGCCCGTTGAACCAGTTTTTCTTGACCATGGTGTAACCCGCCGTGTCGATGAACGACA‑GGCGATCGCCAATGGCCAGCGGACGGTCGAATTGGTACTCGCCGAAGATGTCACCgg                                                                                    <  3:52968/149‑1 (MQ=255)
                                                               tCGAGCTGCTTGACCACGATGGAAGGCATCTTCAGCCCGTTGAACCAGTTTTTCTTGACCATGGTGTAACCCGCCGTGTCGATGAACGACA‑GGCGATCGCCAATGGCCAGCGGACGGTCGAATTGGTACTCGCCGAAGATGTCACCGGc                                                                                   <  5:386157/149‑1 (MQ=255)
                                                                  aGCTGCTTGACCACGATGGAAGGCATCTTCAGCCCGTTGAACCAGTTTTTCTTGACCATGGTGTAACCCGCCGTGTCGATGAACGACA‑GGCGATCGCCAATGGCCAGCGGACGGTCGAATTGGTACTCGCCGAAGATGTCACCGGCCAg                                                                                <  7:20564/149‑1 (MQ=255)
                                                                  aGCTGCTTGACCACGATGGAAGGCATCTTCAGCCCGTTGAACCAGTTTTTCTTGACCATGGTGTAACCCGCCGTGTCGATGAACGACA‑GGCGATCGCCAATGGCCAGCGGACGGTCGAATTGGTACTCGCCGAAGATGTCACCGGCCAg                                                                                <  5:372352/149‑1 (MQ=255)
                                                                        ttGACCACGATGGAAGGCATCTTCAGCCCGTTGAACCAGTTTTTCTTGACCATGGTGTAACCCGCCGTGTCGATGAACGACA‑GGCGATCGCCAATGGCCAGCGGACGGTCGAATTGGTACTCGCCGAAGATGTCACCGGCCAGGCATGa                                                                          >  7:96329/1‑149 (MQ=255)
                                                                             cccgctcgAAGGCCCCTTCACCCCCTTCACCCAGTTTTTCTTGCCCCTGGTGTAACCCGCCGTGTCGTTGCACGCCC‑GGCGATCGCCAATTGCCAGCGGACGGTCTAATTGGTACTCGCCGAAGATGTCACCGGCCAGGCATGACTTGc                                                                     <  2:392809/142‑1 (MQ=255)
                                                                              aCGATGGAAGGCATCTTCAGCCCGTTGAACCAGTTTTTCTTGACCATGGGGTAACCCGCCGTGTCGATGAACGACA‑GGCGATCGCCAATGGCCAGCGGACGGTCGAATTGGTACTCGCCGAAGATGTCACCGGCCAGGCATGACTTGcc                                                                    <  8:131147/149‑1 (MQ=255)
                                                                                  tGGAAGGCATCTTCAGCCCGTTGAACCAGTTTTTCTTGACCATGGTGTAACCCGCCGTGTCGATGAACGACA‑GGCGATCGCCAATGGCCAGCGGACGGTCGAATTGGTACTCGCCGAAGATGTCACCGGCCAGGCATGACTTGCCg                                                                   >  5:386097/1‑146 (MQ=255)
                                                                                  tGGAAGGCATCTTCAGCCCGTTGAACCAGTTTTTCTTGACCATGGTGTAACCCGCCGTGTCGATGAACGACA‑GGCGATCGCCAATGGCCAGCGGACGGTCGAATTGGTACTCGCCGAAGATGTCACCGGCCAGGCATGACTTGCCGCAg                                                                <  7:271440/149‑1 (MQ=255)
                                                                                    gAAGGCATCTTCAGCCCGTTGAACCAGTTTTTCTTGACCATGGTGTAACCCGCCGTGTCGATGAACGACA‑GGCGATCGCCAATGGCCAGCGGACGGTCGAATTGGTACTCGCCGAAGATGTCACCGGCCAGGCATGACTTGCCGCAGAc                                                              <  3:206450/149‑1 (MQ=255)
                                                                                          aTCTTCAGCCCGTTGAACCAGTTTTTCTTGACCATGGTGTAACCCGCCGTGTCGATGAACGACA‑GGCGATCGCCAATGGCCAGCGGACGGTCGAATTGGTACTCGCCGAAGATGTCACCGGCCAGGCATGACTTGCCGCAGACCATGTa                                                        <  5:55836/149‑1 (MQ=255)
                                                                                               cAGCCCGTTGAACCAGTTTTTCTTGACCATGGTGTAACCCGCCGTGTCGATGAACGACA‑GGCGATCGCCAATGGCCAGCGGACGGTCGAATTg                                                                                                           <  3:353733/93‑1 (MQ=255)
                                                                                               cAGCCCGTTGAACCAGTTTTTCTTGACCATGGTGTAACCCGCCGTGTCGATGAACGACA‑GGCGATCGCCAATGGCCAGCGGACGGTCGAATTg                                                                                                           >  4:353733/1‑93 (MQ=255)
                                                                                               cAGCCCGTTGAACCAGTTTTTCTTGACCATGGTGTAACCCGCCGTGTCGATGAACGACA‑GGCGATCGCCAATGGCCAGCGGACGGTCGAATTGGTACTCGCCGAAGATGTCACCGGCCAGGCATGACTTGCCGCAGACCATGTAgg                                                      <  2:61658/146‑1 (MQ=255)
                                                                                                 gCCCGTTGAACCAGTTTTTCTTGACCATGGTGTAACCCGCCGTGTCGATGAACGACA‑GGCGATCGCCAATGGCCAGCGGACGGTCGAATTGGTACTCGCCGAAGATGTCACCGGCCAGGCATGACTTGCCGCAGACCATGTAGGTATGc                                                 >  8:269733/1‑149 (MQ=255)
                                                                                                 gCCCGTTGAACCAGTTTTTCTTGACCATGGTGTAACCCGCCGTGTCGATGAACGACA‑GGCGATCGCCAATGGCCAGCGGACGGTCGAATTGGTACTCGCCGAAGATGTCACCGGCCAGGCATGACTTGCCGCAGACCATGTAGGTATGc                                                 >  2:272661/1‑149 (MQ=255)
                                                                                                     gTTGAACCAGTTTTTCTTGACCATGGTGTAACCCGCCGTGTCGATGAACGACA‑GGCGATCGCCAATGGCCAGCGGACGGTCGAATTGGTACTCGCCGAAGATGTCACCGGCCAGGCATGACTTGCCGCAGACCATGTAGGTATGCTCGc                                             >  3:248620/1‑149 (MQ=255)
                                                                                                         aaCCAGTTTTTCTTGACCATGGTGTAACCCGCCGTGTCGATGAACGACA‑GGCGATCGCCAATGGCCAGCGGACGGTCGAATTGGTACTCGCCGAAGATGTCACCGGCCAGGCATGACTTGCCGCAGACCATGTAGGTATGCTCGCCGt                                          >  3:382696/1‑148 (MQ=255)
                                                                                                         aaCCAGTTTTTCTTGACCATGGTGTAACCCGCCGTGTCGATGAACGACA‑GGCGATCGCCAATGGCCAGCGGACGGTCGAATTGGTACTCGCCGAAGATGTCACCGGCCAGGCATGACTTGCCGCAGACCATGTAGGTATGCTCGCCGt                                          >  3:389043/1‑148 (MQ=255)
                                                                                                            cAGTTTTTCTTGACCATGGTGTAACCCGCCGTGTCGATGAACGACA‑GGCGATCGCCAATGGCCAGCGGACGGTCGAATTGGTACTCGCCGAAGATGTCACCGGCCAGGCATGACTTGCCGCAGACCATGTAGGTATGCTCGCCGTCAtt                                      >  6:49180/1‑149 (MQ=255)
                                                                                                                 tttCTTGACCATGGTGTAACCCGCCGTGTCGATGAACGACA‑GGCGATCGCCAATGGCCAGCGGACGGTCGAATTGGTACTCGCCGAAGATGTCACCGGCCAGGCATGACTTGCCGCAGACCATGTAGGTATGCTCGCCGTCATTggggg                                 >  6:20927/1‑149 (MQ=255)
                                                                                                                 tttCTTGACCATGGTGTAACCCGCCGTGTCGATGAACGACA‑GGCGATCGCCAATGGCCAGCGGACGGTCGAATTGGTACTCGCCGAAGATGTCACCGGCCAGGCATGACTTGCCGCAGACCATGTAGGTATGCTCGCCGTCATTgggg                                  <  1:232153/148‑1 (MQ=255)
                                                                                                                    cTTGACCATGGTGTAACCCGCCGTGTCGATGAACGACA‑GGCGATCGCCAATGGCCAGCGGACGGTCGAATTGGTACTCGCCGAAGATGTCACCGGCCAGGCATGACTTGCCGCAGACCATGTAGGTATGCTCGCCGTCATTGGGGGcc                               <  1:204841/148‑1 (MQ=255)
                                                                                                                         ccATGGTGTAACCCGCCGTGTCGATGAACGACA‑GGCGATCGCCAATGGCCAGCGGACGGTCGAATTGGTACTCGCCGAAGATGTCACCGGCCAGGCATGACTTGCCGCAGACCATGTAGGTATGCTCGCCGTCATTGGGGGCCATCTTg                         >  3:165956/1‑149 (MQ=255)
                                                                                                                               tgtAACCCGCCGTGGCGATGAACGACA‑GGCGATCGCCAATGGCCAGCGGACGGTCGAATTGGTACTCGCCGAAGATGTCACCGGCCTGGCATGACTTGCCGCAGACCATGTAGGTATGCTCGCCGTCATTGGGGGCCATCTTGGCGGTg                   >  7:368150/1‑149 (MQ=255)
                                                                                                                                   aCCCGCCGTGTCGATGAACGACA‑GGCGATCGCCAATGGCGAGCGGACGGTCGAATTGGTACTCGCCGAAGATGTCACCTGCCAGGCATGACTTGCCGCAGACCATTTAGGTATGCTCGCCGTCATTGGGGTCCATCTTGGCGTTGAgg                <  4:382696/148‑1 (MQ=255)
                                                                                                                                   aCCCGCCGTGTCGATGAACGACA‑GGCGATCGCCAATGGCCAGCGGACGGTCGAATTGGTACTCGCCGAAGATGTCACCGGCCAGGCATGACTTGCCGCAGACCATGTAGGTATGCTCGCCGTCATTGGGGGCCATCTTGGCGTTGAgg                <  4:389043/148‑1 (MQ=255)
                                                                                                                                    cccGCCGTGTCGATGAACGACA‑GGCGATCGCCAATGGCCAGCGGACGGTCGAATTGGTACTCGCCGAAGATGTCACCGGCCAGGCATGACTTGCCGCAGACCATGTAGGTATGCTCGCCGTCATTGGGGGCCATCTTGGCGTTGAGGCg              >  8:49801/1‑149 (MQ=255)
                                                                                                                                     ccgccgTGTCGATGAACGACA‑GGCGATCGCCAATGGCCAGCGGACGGTCGAATTGGTACTCGCCGAAGATGTCACCGGCCAGGCATGACTTGCCGCAGACCATGTAGGTATGCTCGCCGTCATTGGGGGCCATCTTGGCGTTGAGGCgg             <  2:360447/149‑1 (MQ=255)
                                                                                                                                      tcgcctgtCGTTCATCGACACGGCGATCGCCAATGGCCAGCGGACGGTCGAATTGGTACTCGCCGAAGATGTCACCGGCCAGGCATGACTTGCCGCAGACCATGTAGGTATGCTCGCCGTCATTGGGGGCCATCTTGGCGTTGAGGCgg             >  4:190257/6‑149 (MQ=255)
                                                                                                                                         cctgtCGTTCATCGACACGGCGATCGCCAATGGCCAGCGGACGGTCGAATTGGTACTCGCCGAAGATGTCACCGGCCAGGCATGACTTGCCGCAGACCATGTAGGTATGCTCGCCGTCATTGGGGGCCATCTTGGCGTTGAGGCGGt            >  4:108906/3‑147 (MQ=255)
                                                                                                                                         cctgtCGTTCATCGACACGGCGATCGCCAATGGCCAGCGGACGGTCGAATTGGTACTCGCCGAAGATGTCACCGGCCAGGCATGACTTGCCGCAGACCATGTAGGTATGCTCGCCGTCATTGGGGGCCATCTTGGCGTTGAGGCGGt            <  3:108906/145‑1 (MQ=255)
                                                                                                                                               gATGAACGACA‑GGCGATCGCCAATGGCCAGCGGACGGTCGAATTGGTACTCGCCGAAGATGTCACCGGCCAGGCATGACTTGCCGCAGACCATGTAGGTATGCTCGCCGTCATTGGGGGCCATCTTGGCGTTGAGGCGGTAGATcagc    <  8:96329/148‑1 (MQ=255)
                                                                                                                                               gATGAACGACA‑GGCGATCGCCAATGGCCAGCGGACGGTCGAATTGGTACTCGCCGAAGATGTCACCGGCCAGGCATGACTTGCCGCAGACCATGTAGGTATGCTCGCCGTCATTGGGGGCCATCTTGGCGTTGAGGCGGTAGATcag     >  7:31833/1‑147 (MQ=255)
                                                                                                                                                aTGAACGACA‑GGCGATCGCCAATGGCCAGCGGACGGTCGAATTGGTACTCGCCGAAGATGTCACCGGCCAGGCATGACTTGCCGCAGACCATGTAGGTATGCTCGCCGTCATTGGGGGCCATCTTGGCGTTGAGGCGGTAGATcagcag  >  8:72938/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
TCTTACGACAGGCTGGAAACGTAGTCTTCGAAACCGAACTCGCGGACGACTTCGACGCTGCCGTCGAGCTGCTTGACCACGATGGAAGGCATCTTCAGCCCGTTGAACCAGTTTTTCTTGACCATGGTGTAACCCGCCGTGTCGATGAACGACA‑GGCGATCGCCAATGGCCAGCGGACGGTCGAATTGGTACTCGCCGAAGATGTCACCGGCCAGGCATGACTTGCCGCAGACCATGTAGGTATGCTCGCCGTCATTGGGGGCCATCTTGGCGTTGAGGCGGTAGATCAGCAG  >  NC_002947/3330902‑3331194

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