Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,509,703 T→C 100% N120S (AAC→AGC)  serS ← serine‑‑tRNA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,509,7030TC95.4% 203.9 / ‑0.7 65N120S (AAC→AGC) serSserine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base T (2/1);  new base C (27/35);  total (29/36)
Fisher's exact test for biased strand distribution p-value = 5.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.51e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CCGGACCACGCAGCACAGCAAAACGGGCACCGGACAGCTTGGCAGCCGCTTCGAAGTCCAGGCCACCGCTGATTTCACCGAGGGCCACATGGTCCTTGATTTCGAAGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGTTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGCGTCCGGCAGGTTGGGAATGGTCAGCAGGATGCCATCCAGCTCGGCCTGAATGGCGTCCAGCTCAGCCTTGCCGGCAGCCAGTTCGTTGGCCATGCGCTCGACGTCA  >  NC_002947/4509557‑4509843
                                                                                                                                                  |                                                                                                                                            
ccGGACCACGCAGCACAGCAAAACGGGCACCGGACAGCTTGGCAGCCGCTTCGAAGTCCAGGCCACCGCTGATTTCACCGAGGGCCACATGGTCCTTGATTTCGAAGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTg                                                                                                                                            <  6:21188/149‑1 (MQ=255)
ccGGACCACGCAGCACAGCAAAACGGGCACCGGACAGCTTGGCAGCCGCTTCGAAGTCCAGGCCACCGCTGATTTCACCGAGGGCCACATGGTCCTTGATTTCGAAGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTg                                                                                                                                            <  1:225556/149‑1 (MQ=255)
 cGGACCACGCAGCACAGCAAAACGGGCACCGGACAGCTTGGCAGCCGCTTCGAAGTCCAGGCCACCGCTGATTTCACCGAGGGCCACATGGTCCTTGATTTCGAAGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTgt                                                                                                                                           <  1:156447/149‑1 (MQ=255)
       aCGCAGCACAGCAAAACGGGCACCGGACAGCTTGGCAGCCGCTTCGAAGTCCAGGCCACCGCTGATTTCACCGAGGGCCACATGGTCCTTGATTTCGAAGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCtt                                                                                                                                     <  5:29438/149‑1 (MQ=255)
       aCGCAGCACAGCAAAACGGGCACCGGACAGCTTGGCAGCCGCTTCGAAGTCCAGGCCACCGCTGATTTCACCGAGGGCCACATGGTCCTTGATTTCGAAGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCt                                                                                                                                      >  5:649/1‑148 (MQ=255)
       aCGCAGCACAGCAAAACGGGCACCGGACAGCTTGGCAGCCGCTTCGAAGTCCAGGCCACCGCTGATTTCACCGAGGGCCACATGGTCCTTGATTTCGAAGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTGGCGTCtt                                                                                                                                     <  7:95135/149‑1 (MQ=255)
               cagcaAAACGGGCACCGGACAGCTTGGCAGCCGCTTCGAAGTCCAGGCCACCGCTGATTTCACCGAGGGCCACATGGTCCTTGATTTCGAAGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTCACTGgcg                                                                                                                             <  4:56514/149‑1 (MQ=255)
               cagcaAAACGGGCACCGGACAGCTTGGCAGCCGCTTCGAAGTCCAGGCCACCGCTGATTTCACCGAGGGCCACATGGTCCTTGATTTCGAAGGCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTCACTGgcg                                                                                                                             <  1:42635/149‑1 (MQ=255)
                  caAAACGGGCACCGGACAGCTTGGCAGCCGCTTCGAAGTCCAGGCCACCGCTGATTTCACCGAGGGCCACATGGTCCTTGATTTCGAAGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTgtcgtc                                                                                                                                       >  5:62634/1‑136 (MQ=255)
                  caAAACGGGCACCGGACAGCTTGGCAGCCGCTTCGAAGTCCAGGCCACCGCTGATTTCACCGAGGGCCACATGGTCCTTGATTTCGAAGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTgtcgtc                                                                                                                                       <  6:62634/136‑1 (MQ=255)
                    aaaCGGGCACCGGACAGCTTGGCAGCCGCTTCGAAGTCCAGGCCACCGCTGATTTCACCGAGGGCCACATGGTCCTTGATTTCGAAGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTCACTGGCGCCGAc                                                                                                                        <  1:132289/149‑1 (MQ=255)
                    aaaCGGGCACCGGACAGCTTGGCAGCCGCTTCGAAGTCCAGGCCACCGCTGATTTCACCGAGGGCCACATGGTCCTTGATTTCGAAGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTCACTGGCGCCGAc                                                                                                                        >  6:42681/1‑149 (MQ=255)
                    aaaCGGGCACCGGACAGCTTGGCAGCCGCTTCGAAGTCCAGGCCACCGCTGATTTCACCGAGGGCCACATGGTCCTTGATTTCGAAGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTCACTGGCGCCGAc                                                                                                                        <  8:205070/149‑1 (MQ=255)
                      aCGGGCACCGGACAGCTTGGCAGCCGCTTCGAAGTCCAGGCCACCGCTGATTTCACCGAGGGCCACATGGTCCTTGATTTCGAAGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTCACTGGCGCCGACCg                                                                                                                      <  6:273875/149‑1 (MQ=255)
                          gCACCGGACAGCTTGGCAGCCGCTACGAAGTCCAGGCCACCGCTGATTTCACCGAGGGCCACATGGTCCTTGATTTCGAAGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTCACTGGCGCCGACCGGCAc                                                                                                                  >  6:106384/1‑149 (MQ=255)
                               cgaactggctgccggcaaggctgagctggacgccattcaggccgagctggatggcatcctgctgaccattcccaacctgccggacgccaGCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTCACTGGCGCCGACCGGCACGCTgg                                                                                                             >  3:32605‑M1/90‑149 (MQ=255)
                                 aCAGCTTGGCAGCCGCTTCGAAGTCCAGGCCACCGCTGATTTCACCGAGGGCCACATGGTCCTTGATTTCGAAGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGCg                                                                                                           >  7:11566/1‑149 (MQ=255)
                                    gCTTGGCAGCCGCTTCGAAGTCCAGGCCACCGCTGATTTCACCGAGGGCCACATGGTCCTTGATTTCGAAGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGCGtgc                                                                                                        >  1:105333/1‑147 (MQ=255)
                                     cTTGGCAGCCGCTTCGAAGTCCAGGCCACCGCTGATTTCACCGAGGGCCACATGGTCCTTGATTTCGAAGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGCGTCCg                                                                                                       >  3:181921/1‑149 (MQ=255)
                                     cTTGGCAGCCGCTTCGAAGTCCAGGCCACCGCTGATTTCACCGAGGGCCACATGGTCCTTGATTTCGAAGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGCGTCCg                                                                                                       >  4:130203/1‑149 (MQ=255)
                                       ctgccggcaaggctgagctggacgccattcaggccgagctggatggcatcctgctgaccattcccaacctgccggacgccaGCGTGCCCCAACGGCGCACTTCGACGTTGTCGTCTTCACTGGCGCCGACCg                                                                                                                      >  7:231351‑M1/82‑132 (MQ=255)
                                       ctgccggcaaggctgagctggacgccattcaggccgagctggatggcatcctgctgaccattcccaacctgccggacgccaGCGTGCCCCAACGGCGCACTTCGACGTTGTCGTCTTCACTGGCGCCGACCg                                                                                                                      <  8:231351‑M1/51‑1 (MQ=255)
                                         gCAGCCGCTTCGAAGTCCAGGCCACCGCTGATTTCACCGAGGGCCACATGGTCCTTGATTTCGAAGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGCGTc                                                                                                         >  6:138868/1‑143 (MQ=255)
                                         gCAGCCGCTTCGAAGTCCAGGCCACCGCTGATTTCACCGAGGGCCACATGGTCCTTGATTTCGAAGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGCGTc                                                                                                         <  5:138868/143‑1 (MQ=255)
                                          cAGCCGCTTCGAAGTCCAGGCCACCGCTGATTTCACCGAGGGCCACATGGTCCTTGATTTCGAAGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTg                                                                                                                                            >  8:130327/1‑107 (MQ=255)
                                          cAGCCGCTTCGAAGTCCAGGCCACCGCTGATTTCACCGAGGGCCACATGGTCCTTGATTTCGAAGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTg                                                                                                                                            <  7:130327/107‑1 (MQ=255)
                                              cGCTTCGAAGTCCAGGCCACCGCTGATTTCACCGAGGGCCACATGGTCCTTGATTTCGAAGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGCGTCCGGCAGGTTgg                                                                                              <  8:252603/149‑1 (MQ=255)
                                                 ttCGAAGTCCAGGCCACCGCTGATTTCACCGAGGGCCACATGGTCCTTGATTTCGAAGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGAAGCTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGCGTCCGGCAGGTTGGGaa                                                                                           >  4:8929/1‑149 (MQ=255)
                                                          cAGGCCACCGCTGATTTCACCGAGGGCCACATGGTCCTTGATTTCGAAGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGCGTCCGGCAGGTTGGGAATGGTcagc                                                                                   <  5:42681/148‑1 (MQ=255)
                                                          cAGGCCACCGCTGATTTCACCGAGGGCCACATGGTCCTTGATTTCGAAGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGCGTCCGGCAGGTTGGGAATGGTcag                                                                                    >  2:83830/1‑147 (MQ=255)
                                                             gCCACCGCTGATTTCACCGAGGGCCACATGGTCCTTGATTTCGAAGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGCGTCCGGCAGGTTGGGAATGGTCAGCAGGa                                                                               <  3:130203/149‑1 (MQ=255)
                                                                   gCTGATTTCACCGAGGGCCACATGGTCCTTGATTTCGAAGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGCGTCCGGCAGGTTGGGAATGGTCAGCAGGATGCCAt                                                                         >  7:268438/1‑149 (MQ=255)
                                                                        tttCACCGAGGGCCACATGGTCCTTGATTTCGAAGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTCACTGGCGCCGCCCGGCACGCTGGCGTCCGGCAGGTTGGGAATGGTCAGCAGGATGCCATCCAGc                                                                    >  3:268623/1‑149 (MQ=255)
                                                                                aGGGCCACATGGTCCTTGATTTCGAAGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGCGTCCGGCAGGTTGGGAATGGTCAGCAGGATg                                                                             >  2:264373/1‑132 (MQ=255)
                                                                                aGGGCCACATGGTCCTTGATTTCGAAGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGCGTCCGGCAGGTTGGGAATGGTCAGCAGGATg                                                                             <  1:264373/132‑1 (MQ=255)
                                                                                  ggCCACATGGTCCTTGATTTCGAAGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGCGTCCGGCAGGTTGGGAATGGTCAGCAGGATGCCATCCAGCTCGGCCTGaa                                                          <  3:8929/149‑1 (MQ=255)
                                                                                   gCCACATGGTCCTTGATTTCGAAGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGCGTCCGGCAGGTTGGGAATGGTCAGCAGGATGCCATCCAGCTCGGCCTGAAt                                                         <  8:11566/149‑1 (MQ=255)
                                                                                    ccACATGGTCCTTGATTTCGAAGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGCGTCCGGCAGGTTgg                                                                                              >  7:142914/1‑111 (MQ=255)
                                                                                    ccACATGGTCCTTGATTTCGAAGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGCGTCCGGCAGGTTgg                                                                                              <  8:142914/111‑1 (MQ=255)
                                                                                    ccACATGGTCCTTGATTTCGAAGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGCGTCCGGCAGGTTGGGAATGGTCAGCAGGATGCCATCCAGCTCGGCCTGAATg                                                        <  1:319927/149‑1 (MQ=255)
                                                                                      agaTGTTCCTTGATTTCGAAGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACTCTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGCGTCCGGCAGGTTGTGAATGGTCAGCAGGATGCCATCCAGCTCGGCCTGAATGGc                                                      <  4:177382/147‑1 (MQ=255)
                                                                                        aTGGTCCTTGATTTCGAAGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGCGTCCGGCAGGTTGGGAATGGTCAGCAGGATGCCATCCAGCTCGGCCTGAATGGCGt                                                    >  5:235193/1‑149 (MQ=255)
                                                                                         tGGTCCTTGATTTCGAAGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGc                                                                                                            >  1:234191/1‑92 (MQ=255)
                                                                                         tGGTCCTTGATTTCGAAGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGc                                                                                                            <  2:234191/92‑1 (MQ=255)
                                                                                                tGATTTCGAAGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGCGTCCGGCAGGTTGGGAATGGTCAGCAGGATGCCATCCAGCTCGGCCTGAATGGCg                                                     >  5:59982/1‑140 (MQ=255)
                                                                                                tGATTTCGAAGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGCGTCCGGCAGGTTGGGAATGGTCAGCAGGATGCCATCCAGCTCGGCCTGAATGGCg                                                     <  6:59982/140‑1 (MQ=255)
                                                                                                    gccGAAGGCGTTCGACTTCGGCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGCGTCCGGCAGGTTGGGAATGGTCAGCAGGATGCCATCCAGCTCGGCCTGAATGGCGTCCAGCTCAGCCt                                        >  3:173464/3‑149 (MQ=255)
                                                                                                        aaGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGCGTCCGGCAGGTTGGGAATGGTCAGCAGGATGCCATCCAGCTCGGCCTGAATGGCGTCCAGCTCAGCCTTGcc                                    >  6:121527/1‑149 (MQ=255)
                                                                                                           tCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGCGTCCGGCAGGTTGGGAATGGTCAGCAGGATGCCATCCAGCTCGGCCTGAATGGCGTCCAGCTCAGCCTTGc                                     <  4:314327/145‑1 (MQ=255)
                                                                                                           tCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGCGTCCGGCAGGTTGGGAATGGTCAGCAGGATGCCATCCAGCTCGGCCTGAATGGCGTCCAGCTCAGCCTTGc                                     >  3:314327/1‑145 (MQ=255)
                                                                                                                cGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGCGTCCGGCAGGTTGGGAATGGTCAGCAGGATGCCATCCAGCTCGGCCTGAATGGCGTCCAGCTCAGCCTTGCCGGCAGCCa                            <  6:235193/149‑1 (MQ=255)
                                                                                                                cGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGCGTCCGGCAGGTTGGGAATGGTCAGCAGGATGCCATCCAGCTCGGCCTGAATGGCGTCCAGCTCAGCCTTGCCGGCAGCCa                            <  3:152879/149‑1 (MQ=255)
                                                                                                                      cGGCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGCGTCCGGCAGGTTGGGAATGGTCAGCAGGATGCCATCCAGCTCGGCCTGAATGGCGTCCAGCTCAGCCTTGCCGGCAGCCAGTTCg                       <  4:181921/148‑1 (MQ=255)
                                                                                                                             CCCCAACGGCGCACTTCGACGCTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGCGTCCGGCAGGTTGGGAATGGTCAGCAGGATGCCATCCAGCTCGGCCTGAATGGCGTCCAGCTCAGCCTTGCCGGCAGCCAGTTCGTTGGCCAt               <  4:227924‑M1/149‑1 (MQ=255)
                                                                                                                             CCCCAACGGCGCACTTCGACGCTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGCGTCCGGCAGGTTGGGAATGGTCAGCAGGATGCCATCCAGCTCGGCCTGAATGGCGTCCAGCTCAGCCTTGCCGGCAGCCAGTTCGTTGGCCAt               <  7:259192‑M1/149‑1 (MQ=255)
                                                                                                                                 aaCGGCGCACTTCGACGCTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGCGTCCGGCAGGTTGGGAATGGTCAGCAGGATGCCATCCAGCTCGGCCTGAATGGCGTCCAGCTCAGc                                          >  8:9785/1‑118 (MQ=255)
                                                                                                                                 aaCGGCGCACTTCGACGCTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGCGTCCGGCAGGTTGGGAATGGTCAGCAGGATGCCATCCAGCTCGGCCTGAATGGCGTCCAGCTCAGc                                          <  7:9785/118‑1 (MQ=255)
                                                                                                                                     gcgcACTTCGACGCTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGCGTCCGGCAGGTTGGGAATGGTCAGCAGGATGCCATCCAGCTCGGCCTGAATGGCGTCCAGCTCAGCCTTGCCGGCAGCCAGTTCGTTGGCCATGCGCTCGa       <  8:268438/149‑1 (MQ=255)
                                                                                                                                      cgcACTTCGACGCTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGCGTCCGGCAGGTTGGGAATGGTCAGCAGGATGCCATCCAGCTCGGCCTGAATGGCGTCCAGCTCAGCCTTGCCGGCAGCCAGTTCGTTGGCCATGCGCTCGAc      >  8:224303/1‑149 (MQ=255)
                                                                                                                                       gcACTTCGACGCTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGCGTCCGGCAGGTTGGGAATGGTCAGCAGGATGCCATCCAGCTCGGCCTGAATGGCGTCCAGCTCAGCCTTGCCGGCAGCCAGTTCGTTGGCCATGCGCTCGACg     >  8:109263/1‑149 (MQ=255)
                                                                                                                                          cTTCGACGTTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGCGTCCGGCAGGTTGGGAATGGTCAGCAGGATGCCATCCAGCTCGGCCTGAATGGCGTCCAGCTCAGCCTTGCCGGCAGCCAGTTCGTTGGCCATGCGCTCGACGt    >  3:18/1‑147 (MQ=255)
                                                                                                                                          cTTCGACGCTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGCGTCCGGCAGGTTGGGAATGGTCAGCAGGATGCCATCCAGCTCGGCCTGAATGGCGTCCAGCTCAGCCTTGCCGGCAGCCAGTTCGTTGGCCATGCGCTCGACGTCa  <  3:227924/149‑1 (MQ=255)
                                                                                                                                          cTTCGACGCTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGCGTCCGGCAGGTTGGGAATGGTCAGCAGGATGCCATCCAGCTCGGCCTGAATGGCGTCCAGCTCAGCCTTGCCGGCAGCCAGTTCGTTGGCCATGCGCTCGACGTCa  <  4:173464/149‑1 (MQ=255)
                                                                                                                                              gACGCTGTCTTCTTCACTGGCGCCGACCGGCACGCTGGCGTCCGGCAGGTTGGGAATGGTCAGCAGGATGCCATCCAGCTCGGCCTGAATGGCGTCCAGCTCAGc                                          <  1:300264/105‑1 (MQ=255)
                                                                                                                                              gACGCTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGCGTCCGGCAGGTTGGGAATGGTCAGCAGGATGCCATCCAGCTCGGCCTGAATGGCGTCCAGCTCAGc                                          >  2:300264/1‑105 (MQ=255)
                                                                                                                                                  |                                                                                                                                            
CCGGACCACGCAGCACAGCAAAACGGGCACCGGACAGCTTGGCAGCCGCTTCGAAGTCCAGGCCACCGCTGATTTCACCGAGGGCCACATGGTCCTTGATTTCGAAGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGTTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGCGTCCGGCAGGTTGGGAATGGTCAGCAGGATGCCATCCAGCTCGGCCTGAATGGCGTCCAGCTCAGCCTTGCCGGCAGCCAGTTCGTTGGCCATGCGCTCGACGTCA  >  NC_002947/4509557‑4509843

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: