Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,527,439 C→T 6.4% R318C (CGT→TGT)  PP_4860 → mannose‑6‑phosphate isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,527,4390CT 6.4% 184.2 / 2.3 62R318C (CGT→TGT) PP_4860mannose‑6‑phosphate isomerase
Reads supporting (aligned to +/- strand):  ref base C (27/31);  new base T (2/2);  total (29/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.64e-01

CGTCAATAGACCGGGCTTTCGCCTTTGCCGAGCAATTCGGGGTGAAGGATGCGGCCGTACTGGCCATGCTGAATGTGGATGGCCGGGTGCTTGATGCCTCCCAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTCAGCGCGGTTCCTGCGTGACGCCGGCTGGTATGAGTGCCGAGGGGGTGACGGCGCGGTGAGCAGGGATGACATGCCATCGACCACGCCTTATCACTTG  >  NC_002947/5527293‑5527582
                                                                                                                                                  |                                                                                                                                               
cGTCAATAGACCGGGCTTTCGCCTTTGCCGAGCAATTCGGGGTGAAGGATGCGGCCGTACTGGCCATGCTGAATGTGGATGGCCGGGTGCTTGATGCCTCCCAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGt                                                                                                                                               <  2:332902/149‑1 (MQ=255)
         aCCGGGCTTTCGCCTTTGCCGAGCAATTCGGGGTGAAGGATGCGGCCGTACTGGCCATGCTGAATGTGGATGGCCGGGTGCTTGATGCCTCCCAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGt                                                                                                                                      <  4:152512/149‑1 (MQ=255)
           cGGGCTTTCGCCTTTGCCGAGCAATTCGGGGTGAAGGATGCGGCCGTACTGGCCATGCTGAATGTGGATGGCCGGGTGCTTGATGCCTCCCAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTg                                                                                                                                     >  6:336731/1‑148 (MQ=255)
           cGGGCTTTCGCCTTTGCCGAGCAATTCGGGGTGAAGGATGCGGCCGTACTGGCCATGCTGAATGTGGATGGCCGGGTGCTTGATGCCTCCCAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTg                                                                                                                                     >  7:326111/1‑148 (MQ=255)
                tttCGCCTTTGCCGAGCAATTCGGGGTGAAGGATGCGGCCGTACTGGCCATGCTGAATGTGGATGGCCGGGTGCTTGATGCCTCCCAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGTGTGAAGGCAGTGAAGCa                                                                                                                                >  6:382432/1‑148 (MQ=255)
                tttCGCCTTTGCCGAGCAATTCGGGGTGAAGGATGCGGCCGTACTGGCCATGCTGAATGTGGATGGCCGGGTGCTTGATGCCTCCCAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGTGTGAAGGCAGTGAAGCa                                                                                                                                >  6:377302/1‑148 (MQ=255)
                   cGCCTTTGCCGAGCAATTCGGGGGGAAGGATGCGGCCGTACTTGCCCTGCGGAATGGGGCTGGCCGGGGGCTTGCTGCCTCCCAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGGGAAGGCAGTGAAGCAAAg                                                                                                                             <  1:265664/148‑1 (MQ=255)
                            cGAGCAATTCGGGGTGAAGGATGCGGCCGTACTGGCCATGCTGAATGTGGATGGCCGGGTGCTTGATGCCTCCCAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGa                                                                                                                                    <  5:241333/132‑1 (MQ=255)
                            cGAGCAATTCGGGGTGAAGGATGCGGCCGTACTGGCCATGCTGAATGTGGATGGCCGGGTGCTTGATGCCTCCCAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGa                                                                                                                                    >  6:241333/1‑132 (MQ=255)
                            cGAGCAATTCGGGGTGAAGGATGCGGCCGTACTGGCCATGCTGAATGTGGATGGCCGGGTGCTTGATGCCTCCCAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCgggg                                                                                                                    >  5:286917/1‑148 (MQ=255)
                            cGAGCAATTCGGGGTGAAGGATGCGGCCGTACTGGCCATGCTGAATGTGGATGGCCGGGTGCTTGATGCCTCCCAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGc                                                                                                                   <  2:22574/149‑1 (MQ=255)
                                cAATTCGGGGTGAAGGATGCGGCCGTACTGGCCATGCTGAATGTGGATGGCCGGGTGCTTGATGCCTCCCAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGtt                                                                                                               <  6:286917/149‑1 (MQ=255)
                                 aaTTCGGGGTGAAGGATGCGGCCGTACTGGCATGCTGAATGTGGATGGCCGGGTGCTTGATGCCTCCCAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTgg                                                                                                             >  3:184890/1‑149 (MQ=255)
                                      ggggTGAAGGATGCGGCCGTACTGGCCATGCTGAATGTGGATGGCCGGGTGCTTGCTGCCTCCCAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGgcgc                                                                                                         <  8:242796/149‑1 (MQ=255)
                                      ggggTGAAGGATGCGGCCGTACTGGCCATGCTGAATGTGGATGGCCGGGTGCTTGATGCCTCCCAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGgcgc                                                                                                         <  6:78158/149‑1 (MQ=255)
                                          tGAAGGATGCGGCCGTACTGGCCATGCTGAATGTGGATGGCCGGGTGCTTGATGCCTCCCAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCtt                                                                                                     >  6:114028/1‑149 (MQ=255)
                                                aTGCGGCCGTACTGGCCATGCTGAATGTGGATGGCCGGGTGCTTGATGCCTCCCAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGt                                                                                                                                      <  7:80460/110‑1 (MQ=255)
                                                aTGCGGCCGTACTGGCCATGCTGAATGTGGATGGCCGGGTGCTTGATGCCTCCCAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGt                                                                                                                                      >  8:80460/1‑110 (MQ=255)
                                                           cTGGCCATGCTGAATGTGGATGGCCGGGTGCTCGAGGCCTCCCAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTCAGCGCGGTTCCTGc                                                                                     <  8:248574/148‑1 (MQ=255)
                                                               ccATGCTGAATGTGGATGGCCGGGTGCTTGATGCCTCCCAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTCAGCGCGGTTCCTGCGTGAc                                                                                >  4:87954/1‑149 (MQ=255)
                                                               ccATGCTGAATGTGGATGGCCGGGTGCTTGATGCCTCACAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCACTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTCAGCGCGGTTCCTGCGTGAc                                                                                >  8:141242/1‑149 (MQ=255)
                                                                  tGCTGAATGTGGATGGCCGGGTGCTTGATGCCTCCCAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTCAGCGCGGTTCCTGCGTGACGcc                                                                             >  4:297109/1‑149 (MQ=255)
                                                                    cTGAATGTGGATGGCCGGGTGCTTGATGCCTCCCAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTCAGCGCGGTTCCTGCGTGACGCCgg                                                                           >  7:31507/1‑149 (MQ=255)
                                                                     tGAATGTGGATGGCCGGGTGCTTGATGCCTCCCAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTCAGCGCGGTTCCTGCGTGACGCCGGc                                                                          >  4:84212/1‑149 (MQ=255)
                                                                       aaTGTGGATGGCCGGGTGCTTGATGCCTCCCAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGTGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTCAGCGCGGTTCCTGCGTGACGCCGGCTg                                                                        <  5:382432/149‑1 (MQ=255)
                                                                       aaTGTGGATGGCCGGGTGCTTGATGCCTCCCAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGTGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTCAGCGCGGTTCCTGCGTGACGCCGGCTg                                                                        <  5:377302/149‑1 (MQ=255)
                                                                       aaTGTGGATGGCCGGGTGCTTGATGCCTCCCAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTCAGCGCGGTTCCTGCGTGACGCCGGCTg                                                                        <  5:336731/149‑1 (MQ=255)
                                                                         tgtgGATGGCCGGGTGCTTGATGCCTCCCAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTCAGCGCGGTTCCTGCGTGACGCCGGCTGGt                                                                      <  5:65637/149‑1 (MQ=255)
                                                                           tgGATGGCCGGGTTCTTGATGCCTCCCAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAATCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTCAGCGCGGTTCCTGCGTGACGCCGGCTGGTAt                                                                    <  4:37083/149‑1 (MQ=255)
                                                                           tgGATGGCCGGGTGCTTGATGCCTCCCAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTCAGCGCGGTTCCTGCGTGACGCCGGCTGGTAt                                                                    >  2:296960/1‑149 (MQ=255)
                                                                                     ggTGCTTGATGCCTCCCAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTCAGCGCGGTTCCTGCGTGACGCCGGCTGGTATGAGTGCCGAg                                                          <  5:114028/149‑1 (MQ=255)
                                                                                      gTGCTTGATGCCTCCCAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTCAGCGCGGTTCCTGCGTGACGCCGGCTGGTATGAGTGCCGAgg                                                         <  4:184890/149‑1 (MQ=255)
                                                                                         cTTGATGCCTCCCAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTCAGCGCGGTTCCTGCGTGACGCCGGCTGGTATGAGTTCCGAgggg                                                       <  7:141242/148‑1 (MQ=255)
                                                                                          ttGATGCCTCCCAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTCAGCGCGGTTCCTGCGTGACGCCGGCTGGTATGAGTGCCGAGGGGGt                                                     >  7:118910/1‑149 (MQ=255)
                                                                                                  tCCCAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTCAGCGCGGTTCCTGCGTGACGCCGGCTGGTATGAGTGCCGAGGGGGTGACGgcgc                                             >  6:14123/1‑149 (MQ=255)
                                                                                                  tCCCAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTCAGCGCGGTTCCTGCGTGACGCCGGCTGGTATGAGTGCCGAGGGGGTGACGgcgc                                             >  6:19640/1‑149 (MQ=255)
                                                                                                         gTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTCAGCGCGGTTCCTGCGTGACGCCGGCTGGTATGAGTGCCGAGGGGGTGACGGCGCGGTGAg                                       >  6:227402/1‑148 (MQ=255)
                                                                                                            cgagcgcgcAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTCAGCGCGGTTCCTGCGTGACGCCGGCTGGTATGAGTGCCGAGGGGGTGACGGCGCGGTGAGCAgg                                   <  1:135361/144‑1 (MQ=255)
                                                                                                              tGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTCAGCGCGGTTCCTGCGTGACGCCGGCTGGTATGAGTGCCGAGGGGGTGACGGCGCGGTGAGCAgg                                   >  2:135361/1‑147 (MQ=255)
                                                                                                               gggCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTCAGCGCGGTTCCTGCGTGACGCCGGCTGGTATGAGTGCCGAGGGGGTGACGGCGCGGTGAGCAGGGAt                                >  4:216535/1‑149 (MQ=255)
                                                                                                                ggCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTCAGCGCGGTTCCTGCGTGACGCCGGCTGGTATGAGTGCCGAGGGGGTGACGGCGCGGTGAGCAGGGAt                                <  5:55779/148‑1 (MQ=255)
                                                                                                                ggCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTCAGCGCGGTTCCTGCGTGACGCCGGCTGGTATGAGTGCCGAGGGGGTGACGGCGCGGTGAGCAGGGAt                                <  3:84212/148‑1 (MQ=255)
                                                                                                                ggCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTCAGCGCGGTTCCTGCGTGACGCCGGCTGGTATGAGTGCCGAGGGGGTGACGGCGCGGTGAGCAGGGAt                                <  8:31507/148‑1 (MQ=255)
                                                                                                                    cAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTCAGCGCGGTTCCTGCGTGACGCCGGCTGGTATGAGTGCCGAGGGGGTGACGGCGCGGTGAGCAGGGATGACAt                           <  5:227402/149‑1 (MQ=255)
                                                                                                                      ggCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCGTTCAGCGCGGTTCCTGCGTGACGCCGGCTGGTATGAGTGCCGAGGGGGTGACGGCGCGGTGAGCAGGGATGACATGc                         <  8:326111/149‑1 (MQ=255)
                                                                                                                          gAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTCAGCGCGGTTCCTGCGTGACGCCGGCTGGTATGAGTGCCGAGGGGGTGACGGCGCGGTGAGCAGGGATGACATGCCATc                     <  1:296960/149‑1 (MQ=255)
                                                                                                                               ccTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTCAGCGCGGTTCCTGCGTGACGCCGGCTGGTATGAGTGCCGAGGGGGTGACGGCGCGGTGAGCa                                     <  1:15711/128‑1 (MQ=255)
                                                                                                                               ccTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTCAGCGCGGTTCCTGCGTGACGCCGGCTGGTATGAGTGCCGAGGGGGTGACGGCGCGGTGAGCa                                     >  2:15711/1‑128 (MQ=255)
                                                                                                                               ccTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTCAGCGCGGTTCCTGCGTGACGCCGGCTGGTATGAGTGCCGAGGGGGTGACGGCGCGGTGAGCAGGGATGACATGCCATCGACCa                <  5:19640/149‑1 (MQ=255)
                                                                                                                               ccTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTCAGCGCGGTTCCTGCGTGACGCCGGCTGGTATGAGTGCCGAGGGGGTGACGGCGCGGTGAGCAGGGATGACATGCCATCGACCa                <  5:14123/149‑1 (MQ=255)
                                                                                                                                cTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTCAGCGCGGTTCCTGCGTGACGCCGGCTGGTATGAGTGCCGAGGGGGTGACGGCGCGGTGAGCAGGGATGACATGCCATCGACCAc               >  3:129429/1‑149 (MQ=255)
                                                                                                                                  gCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTCAGCGCGGTTCCTGCGTGACGCCGGCTGGTATGAGTGCCGAGGGGGTGACGGcgcg                                            <  7:88098/118‑1 (MQ=255)
                                                                                                                                  gCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTCAGCGCGGTTCCTGCGTGACGCCGGCTGGTATGAGTGCCGAGGGGGTGACGGcgcg                                            >  8:88098/1‑118 (MQ=255)
                                                                                                                                        gcTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTCAGCGCGGTTCCTGCGTGACGCCGGCTGGTATGAGTGCCGAGGGGGTGACGgc                                               >  7:68285/1‑109 (MQ=255)
                                                                                                                                        gcTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTCAGCGCGGTTCCTGCGTGACGCCGGCTGGTATGAGTGCCGAGGGGGTGACGgc                                               <  8:68285/109‑1 (MQ=255)
                                                                                                                                          ttACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTAAGCGCGGTTCCTGCGTGACGCCGGCTGGTATGAGTGCCGAGGGGGTGACGGCGCGGTGAGCAGGGATGACATGCCATCGACCACGCCTTATCAc     >  4:116583/1‑149 (MQ=255)
                                                                                                                                          ttACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCATTTCAGCGCGGTTCCTGCGTGACGCCGGCTGGTATGAGTGCCGAGGGGGTGACGGCGCGGTGAGCAGGGATGACATGCCATCGACCACGCCTTATCAc     >  6:144545/1‑149 (MQ=255)
                                                                                                                                             cATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTCAGCGCGGTTCCTGCGTGACGCCGGCTGGTATGAGTGCCGAGGGGGTGACGGCGCGGTGAGCAGGGATGACATGCCATCGACCACGCCTTATCACTTg  <  3:216535/149‑1 (MQ=255)
                                                                                                                                              aTTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTCAGCGCGGTTCCTGCGTGACGCCGGCTGGTATGAGTGCCGAGGGGGTGACGGCGCGGTGAGCAGGGATGAc                             <  4:345620/121‑1 (MQ=255)
                                                                                                                                              aTTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTCAGCGCGGTTCCTGCGTGACGCCGGCTGGTATGAGTGCCGAGGGGGTGACGGCGCGGTGAGCAGGGATGAc                             <  4:340094/121‑1 (MQ=255)
                                                                                                                                              aTTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTCAGCGCGGTTCCTGCGTGACGCCGGCTGGTATGAGTGCCGAGGGGGTGACGGCGCGGTGAGCAGGGATGAc                             >  3:345620/1‑121 (MQ=255)
                                                                                                                                              aTTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTCAGCGCGGTTCCTGCGTGACGCCGGCTGGTATGAGTGCCGAGGGGGTGACGGCGCGGTGAGCAGGGATGAc                             >  3:340094/1‑121 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
CGTCAATAGACCGGGCTTTCGCCTTTGCCGAGCAATTCGGGGTGAAGGATGCGGCCGTACTGGCCATGCTGAATGTGGATGGCCGGGTGCTTGATGCCTCCCAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTCAGCGCGGTTCCTGCGTGACGCCGGCTGGTATGAGTGCCGAGGGGGTGACGGCGCGGTGAGCAGGGATGACATGCCATCGACCACGCCTTATCACTTG  >  NC_002947/5527293‑5527582

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: