Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,010,976 C→A 8.5% R13S (CGC→AGC)  bcsG → endoglucanase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,010,9760CA 8.5% 215.1 / 6.7 69R13S (CGC→AGC) bcsGendoglucanase
Reads supporting (aligned to +/- strand):  ref base C (33/30);  new base A (2/4);  total (35/34)
Fisher's exact test for biased strand distribution p-value = 4.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTGCTAGTGCGCCTGTGTTCGGCATTACGTGGCGTGTTGCGCCGTCGCCTTCCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGCCGGGCAAGTGGCTGCGTATT  >  NC_002947/3010828‑3011122
                                                                                                                                                    |                                                                                                                                                   
ttgCTAGTGCGCCTGTGTTCGGCATTACGTGGCGTGTTGCGCCGTCGCCTTCCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGccc                                                                                                                                                     >  3:353472/1‑149 (MQ=255)
  gCTAGTGCGCCTGTGTTCGGCATTACGTGGCGTGTTGCGCCGTCGCCTTCCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGc                                                                                                                                                   >  7:90742/1‑149 (MQ=255)
  gCTAGTGCGCCTGTGTTCGGCATTACGTGGCGTGTTGCGCCGTCGCCTTCCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGc                                                                                                                                                   >  2:48559/1‑149 (MQ=255)
          gcCTGTGTTCGGCATTACGTGGCGTGTTGCGCCGTCGCCTTCCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTgg                                                                                                                                           >  4:28970/1‑149 (MQ=255)
             tgtgTTCGGCATTACGTGGCGTGTTGCGCCGTCGCCTTCCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCAGTCCAGGCCAGCGGGccagtacc                                                                                                                                        >  3:40247/1‑143 (MQ=255)
                 ttCGGCATTACGTGGCGTGTTGCGCCGTCGCCTTCCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTgg                                                                                                                                      >  3:180023/1‑147 (MQ=255)
                 ttCGGCATTACGTGGCGTGTTGCGCCGTCGCCTTCCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCAGTCCAGGCCAGCGGGCCAGTACCGGGa                                                                                                                                    <  6:211970/149‑1 (MQ=38)
                 ttCGGCATTACGTGGCGTGTTGCGCCGTCGCCTTCCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCAGTCCAGGCCAGCGGGCCAGTACCGGGa                                                                                                                                    <  4:40247/149‑1 (MQ=38)
                          aCGTGGCGTGTTGCGCCGTCGCCTTCCCCCCCGTTAGCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAAc                                                                                                                           >  1:246489/1‑149 (MQ=255)
                             tGGCGTGTTGCGCCGTCGCCTTCCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACctc                                                                                                                        >  8:168752/1‑149 (MQ=255)
                                    ttGCGCCGTCGCCTTCCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTAttttt                                                                                                                 <  1:48559/149‑1 (MQ=255)
                                        gcCGTCGCCTTCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGc                                                                                                             >  6:193640/1‑148 (MQ=255)
                                           gTCGCCTTCCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCaaa                                                                                                          <  3:28970/149‑1 (MQ=255)
                                              gCCTTCCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGtt                                                                                                       >  2:245216/1‑149 (MQ=255)
                                              gCCTTCCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGtt                                                                                                       <  4:353472/149‑1 (MQ=255)
                                              gCCTTCCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGtt                                                                                                       >  7:326527/1‑149 (MQ=255)
                                                 ttCCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTct                                                                                                    >  4:333638/1‑149 (MQ=255)
                                                    ccccccGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTctgc                                                                                                  <  1:245216/148‑1 (MQ=255)
                                                     cccccGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTctgctg                                                                                                <  7:49640/149‑1 (MQ=255)
                                                      ccccGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGc                                                                                                                                               <  6:104519/101‑1 (MQ=255)
                                                      ccccGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGc                                                                                                                                               >  5:104519/1‑101 (MQ=255)
                                                      ccccGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTggg                                                                                                                                     >  8:287398/1‑111 (MQ=255)
                                                      ccccGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTggg                                                                                                                                     <  7:287398/111‑1 (MQ=255)
                                                      ccccGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTgg                                                                                               <  3:333638/149‑1 (MQ=255)
                                                      ccccGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTgg                                                                                               >  5:102246/1‑149 (MQ=255)
                                                       cccGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGt                                                                                              >  1:125511/1‑149 (MQ=255)
                                                               cTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTgg                                                                                      <  4:180023/149‑1 (MQ=255)
                                                               cTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTgg                                                                                      >  1:257436/1‑149 (MQ=255)
                                                                tGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTgg                                                                                                                                      >  3:287204/1‑100 (MQ=255)
                                                                tGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTgg                                                                                                                                      <  4:287204/100‑1 (MQ=255)
                                                                tGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTggg                                                                                     >  8:187053/1‑149 (MQ=255)
                                                                tGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTggg                                                                                     >  8:192244/1‑149 (MQ=255)
                                                                 gAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTggg                                                                                     >  7:29113/1‑148 (MQ=255)
                                                                 gAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTggg                                                                                     >  3:131365/1‑148 (MQ=255)
                                                                     ccGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGc                                                                               <  7:187053/149‑1 (MQ=255)
                                                                     ccGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGc                                                                               >  4:185998/1‑149 (MQ=255)
                                                                     ccGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGc                                                                               <  7:192244/149‑1 (MQ=255)
                                                                       gCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTc                                                                             <  7:168752/149‑1 (MQ=255)
                                                                            tcGTACGCCTCCGCTCTTCACGCTATGAACCCCATAATTTCTAATACCATCTTCTGCCTCCCTGTACTGGCTCTCTGGCCTGTACTTGTATGGTTGAACCTCAATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCTACt                                                                         <  2:332215/147‑1 (MQ=255)
                                                                                 cgccgccgcGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAgg                                                                   <  3:185998/149‑1 (MQ=255)
                                                                                 cgccgccgcGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAg                                                                    <  6:102246/148‑1 (MQ=255)
                                                                                       cgcgCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGc                                                             <  2:125511/149‑1 (MQ=255)
                                                                                               aCGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCtg                                                     <  2:257436/149‑1 (MQ=255)
                                                                                                          cccAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTg                                                                                                                                       >  3:183617/1‑57 (MQ=255)
                                                                                                          cccAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTg                                                                                                                                       <  4:183617/57‑1 (MQ=255)
                                                                                                                   aCTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGc                                  <  1:263684/149‑1 (MQ=255)
                                                                                                                    cTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCtt                                <  2:246489/149‑1 (MQ=255)
                                                                                                                     tAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCt                                    <  8:225815/144‑1 (MQ=255)
                                                                                                                     tAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCt                                    >  7:225815/1‑144 (MQ=255)
                                                                                                                      aaTACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGc                                          >  1:7250/1‑137 (MQ=255)
                                                                                                                      aaTACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGc                                          <  2:7250/137‑1 (MQ=255)
                                                                                                                       aTACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGt                              <  8:29113/148‑1 (MQ=255)
                                                                                                                           caccaTCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGc                                  >  7:45706/1‑140 (MQ=255)
                                                                                                                           caccaTCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGc                                  <  8:45706/140‑1 (MQ=255)
                                                                                                                           caccaTCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCtt                         >  8:327755/1‑149 (MQ=255)
                                                                                                                            accaTCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTc                                                                                                     >  6:271292/1‑73 (MQ=255)
                                                                                                                            accaTCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTc                                                                                                     <  5:271292/73‑1 (MQ=255)
                                                                                                                            accaTCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGttt                                      >  4:319768/1‑135 (MQ=255)
                                                                                                                            accaTCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCttt                        >  3:24325/1‑149 (MQ=255)
                                                                                                                               aTCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAATGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTGGTTCGGTGCGTTTCTGCTTGTCCCTTTGcc                     >  4:231939/1‑149 (MQ=255)
                                                                                                                                tCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGCCg                    <  8:326527/149‑1 (MQ=255)
                                                                                                                                 cACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGCCgg                   >  2:109273/1‑149 (MQ=255)
                                                                                                                                  aCCCTCCCAGTCCAGGCCAGCGGGCCAGTACCGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGc                                  <  6:17893/133‑1 (MQ=37)
                                                                                                                                    ccTCCCAGTCCAGGCCAGCGGGCCAGTACCGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCt                                            <  1:249269/121‑1 (MQ=37)
                                                                                                                                    ccTCCCAGTCCAGGCCAGCGGGCCAGTACCGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCt                                            >  2:249269/1‑121 (MQ=37)
                                                                                                                                              cTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGCCGGGCAAGTGGCTGc       <  4:24325/148‑1 (MQ=255)
                                                                                                                                                 gCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTg                                   <  8:325656/117‑1 (MQ=255)
                                                                                                                                                 gCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTg                                   >  7:325656/1‑117 (MQ=255)
                                                                                                                                                 gCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGCCgg                   >  5:363030/1‑133 (MQ=255)
                                                                                                                                                 gCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGCCgg                   <  6:363030/133‑1 (MQ=255)
                                                                                                                                                  cccGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGCCGGGCAAGTGGCTGCGTAtt  >  6:366824/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TTGCTAGTGCGCCTGTGTTCGGCATTACGTGGCGTGTTGCGCCGTCGCCTTCCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGCCGGGCAAGTGGCTGCGTATT  >  NC_002947/3010828‑3011122

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: