Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,006,596 G→C 7.6% A427G (GCC→GGC)  PP_3535 ← gamma‑glutamyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,006,5960GC 7.6% 204.3 / 6.2 66A427G (GCC→GGC) PP_3535gamma‑glutamyltransferase
Reads supporting (aligned to +/- strand):  ref base G (37/24);  new base C (3/2);  total (40/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.87e-01

GTTTGCCCGGCTGCACCCGATTGGCCACAGGCTTGCCGTGCTCGTGGGGGATGAACGAGAAGTCGGTCAGCTGATTGTTGAGCAGGAAGCCCTTGACCATCAGGTGCGAGCCAAATGC‑‑GGCCTCCACCGAGGTGGTCATGGCCAGGGCCTGGCCGTTGTCATCCACCGCACTGAGGTGAGAGGTGGAGATGCGCAGTGGCGAGCGGTCGGGCGCAAGCGCCAGGTCCGCTCCTTGCGGCTGGCCCGGACGGGCGCGCTTCATGCTGTGCTCGCCGATCAGCGCCGCGCGGCTGG  >  NC_002947/4006450‑4006743
                                                                                                                                                    |                                                                                                                                                    
gTTTGACCGGCTGCACCCGATTGGCCACAGGCTTGCCGTGCTCGTGGGGGATGAACGAGAAGTCGGTCAGCTGATTGTTGAGCAGGAAGCCCTTGACAATCAGGTGCGAGCCAAATGC‑‑GGCCTCCACCGAGGTGGTCATGGCCAGGGcc                                                                                                                                                    >  8:208840/1‑149 (MQ=255)
      ccGGCTGCACCCGATTGGCCACAGGCTTGCCGTGCTCGTGGGGGATGAACGAGAAGTCGGTCAGCTGATTGTTGAGCAGGAAGCCCTTGACCATCAGGTGCGAGCCAAATGC‑‑GGCCTCCACCGAGGTGGTCATGGCCAGGGCCTGGCCg                                                                                                                                              >  1:327437/1‑149 (MQ=255)
      ccGGCTGCACCCGATTGGCCACAGGCTTGCCGTGCTCGTGGGGGATGAACGAGAAGTCGGTCAGCTGATTGTTGAGCAGGAAGCCCTTGACCATCAGGTGCGAGCCAAATGC‑‑GGCCTCCACCGAGGTGGTCATGGCCAGGCCCTGGcca                                                                                                                                              >  2:373451/1‑148 (MQ=255)
        ggCTGCACCCGATTGGCCACAGGCTTGCCGTGCTCGTGGGGGATGAACGAGAAGTCGGTCAGCTGATTGTTGAGCAGGAAGCCCTTGACCATCAGGTGCGAGCCAAATGC‑‑GGCCTCCACCGAGGTGGTCATGGCCAGGCCCTGGccatg                                                                                                                                            <  1:373451/149‑4 (MQ=255)
              aCCCGATTGGCCACAGGCTTGCCGTGCTCGTGGGGGATGAACGAGAAGTCGGTCAGCTGATTGTTGAGCAGGAAGCCCTTGACCATCAGGTACGAGCCAAATGC‑‑GGCCTCCACCGAGGTGGTCATGGCCAGGGCCTGGCCGTTGTCATc                                                                                                                                      >  6:217556/1‑149 (MQ=255)
                 cGATTGGCCACAGGCTTGCCGTGCTCGTGGGGGATGAACGAGAAGTCGGTCAGCTGATTGTTGAGCAGGAAGCCCTTGACCATCAGGTGCGAGCCAAATGC‑‑GGCCTCCACCGAGGTGGTCATGGCCAGGGCCTGGCCGTTGTCATCCAc                                                                                                                                   >  2:55459/1‑149 (MQ=255)
                   aTTGGCCACAGGCTTGCCGTGCTCGTGGGGGATGAACGAGAAGTCGGTCAGCTGATTGTTGAGCAGGAAGCCCTTGACCATCAGGTGCGAGCCAAATGC‑‑GGCCTCCACCGAGGTGGTCATGGCCAGGGCCTGGCCGTTGt                                                                                                                                          >  5:7104/1‑140 (MQ=255)
                   aTTGGCCACAGGCTTGCCGTGCTCGTGGGGGATGAACGAGAAGTCGGTCAGCTGATTGTTGAGCAGGAAGCCCTTGACCATCAGGTGCGAGCCAAATGC‑‑GGCCTCCACCGAGGTGGTCATGGCCAGGGCCTGGCCGTTGt                                                                                                                                          <  6:7104/140‑1 (MQ=255)
                     tGGCCACAGGCTTGCCGTGCTCGTGGGGGATGAACGAGAAGTCGGTCAGCTGATTGTTGAGCAGGAAGCCCTTGACCATCAGGTGCGAGCCAAATGC‑‑GGCCTCCACCGAGGTGGTCATGGCCAGGGCCTGGCCGTTGTCATCCACCGCa                                                                                                                               >  5:256150/1‑149 (MQ=255)
                         cacaGGCTTGCCGTGCTCGTGGGGGATGAACGAGAAGTCGGTCAGCTGATTGTTGAGCAGGAAGCCCTTGACCATCAGGTGCGAGCCAAATGC‑‑GGCCTCCACCGAGGTGGTCATGGCCAGGGCCTGGCCGTTGTCATCCACCGCACTGa                                                                                                                           <  6:175670/149‑1 (MQ=255)
                          acaGGCTTGCCGTGCTCGTGGGGGATGAACGAGAAGTCGGTCAGCTGATTGTTGAGCAGGAAGCCCTTGACCATCAGGTGCGAGCCAAATGC‑‑GGCCTCCACCGAGGTGGTCATGGCCAGGGCCTGGCCGTTGTCATCCACCGCACTGAg                                                                                                                          <  6:220345/149‑1 (MQ=255)
                           caGGCTTGCCGTGCTCGTGGGGGATGAACGAGAAGTCGGTCAGCTGATTGTTGAGCAGGAAGCCCTTGACCATCAGGTACGAGCCAAATGC‑‑GGCCTCCACCGAGGTGGTCATGGCCAGGGCCTGGCCGTTGTCATCCACCGCACTGAgg                                                                                                                         <  5:217556/149‑1 (MQ=255)
                                  gCCGTGCTCGTGGGGGATGAACGAGAAGTCGGTCAGCTGATTGTTGAGCAGGAAGCCCTTGACCATCAGGTGCGAGCCAAATGC‑‑GGCCTCCACCGAGGTGGTCATGGCCAGGGCCTGGCCGTTGTCATCCACCGCACTGAGGTgaga                                                                                                                    <  5:384913/147‑1 (MQ=255)
                                  gCCGTGCTCGTGGGGGATGAACGAGAAGTCGGTCAGCTGATTGTTGAGCAGGAAGCCCTTGACCATCAGGTGCGAGCCAAATGC‑‑GGCCTCCACCGAGGTGGTCATGGCCAGGGCCTGGCCGTTGTCATCCACCGCACTGAGGTGAGAgg                                                                                                                  >  8:25300/1‑149 (MQ=255)
                                  gCCGTGCTCGTGGGGGATGAACGAGAAGTCGGTCAGCTGATTGTTGAGCAGGAAGCCCTTGACCATCAGGTGCGAGCCAAATGC‑‑GGCCTCCACCGAGGTGGTCATGGCCAGGGCCTGGCCGTTGTCATCCACCGCACTGAGGTGAGAgg                                                                                                                  >  4:87393/1‑149 (MQ=255)
                                   ccGTGCTCGTGGGGGATGAACGAGAAGTCGGTCAGCTGATTGTTGAGCAGGAAGCCCTTGACCATCAGGTGCGAGCCAAATGCTTGGCCTCCACCGAGGTGGTCATGGCCAGGGCCTGGCCGTTGTCATCCACCGCACTGAGGTGagag                                                                                                                   <  1:112584/149‑1 (MQ=255)
                                    cGTGCTCGTGGGGGATGAACGAGAAGTCGGTCAGCTGATTGTTGAGCAGGAAGCCCTTGACCATCAGGTGCGAGCCAAATGC‑‑GGCCTCCACCGAGGTGGTCATGGCCAGGGCCTGGCCGTTGTCATCCACCGCACTGAGGTGAGAGGTg                                                                                                                <  3:87393/149‑1 (MQ=255)
                                    cGTGCTCGTGGGGGATGAACGAGAAGTCGGTCAGCTGATTGTTGAGCAGGAAGCCCTTGACCATCAGGTGCGAGCCAAATGC‑‑GGCCTCCACCGAGGTGGTCATGGCCAGGCCCTGGccat                                                                                                                                             <  5:1538/120‑3 (MQ=255)
                                    cGTGCTCGTGGGGGATGAACGAGAAGTCGGTCAGCTGATTGTTGAGCAGGAAGCCCTTGACCATCAGGTGCGAGCCAAATGC‑‑GGCCTCCACCGAGGTGGTCATGGCCAGGCCCTGGccat                                                                                                                                             >  6:1538/1‑118 (MQ=255)
                                          cGTGGGGGATGAACGAGAAGTCGGTCAGCTGATTGTTGAGCAGGAAGCCCTTGACCATCAGGTGCGAGCCAAATGC‑‑GGCCTCCACCGAGGTGGTCATGGCCAGGGCCTGGCCGTTGTCATCCACCGCACTGAGGTGAGAGGTGGAGATg                                                                                                          >  4:53948/1‑149 (MQ=255)
                                             gggggATGAACGAGAAGTCGGTCAGCTGATTGTTGAGCAGGAAGCCCTTGACCATCAGGTGCGAGCCAAATGC‑‑GGCCTCCACCGAGGTGGTCATGGCCAGGGCCTGGCCGTTGTCATCCACCGCACTGAGGTGAGAGGTGGAGATgcg                                                                                                        >  7:24525/1‑148 (MQ=255)
                                              ggggATGAACGAGAAGTCGGTCAGCTGATTGTTGAGCAGGAAGCCCTTGACCATCAGGTGCGAGCCAAATGC‑‑GGCCTCCACCGAGGTGGTCATGGCCAGGGCCTGGCCGTTGTCATCCACCGCACTGAGGTGAGAGGTGGAGATGCGCa                                                                                                      <  1:55459/149‑1 (MQ=255)
                                                     aaCGAGAAGTCGGTCAGCTGATTGTTGAGCAGGAAGCCCTTGACCATCAGGTGCGAGCCAAATGC‑‑GGCCTCCACCGAGGTGGTCATGGCCAGGGCCTGGCCGTTGTCATCCACCGCACTGAGGTGAGAGGTGGAGATGCGCAGTGGc                                                                                                 >  1:171463/1‑147 (MQ=255)
                                                        gagaAGTCGGTCAGCTGATTGTTGAGCAGGAAGCCCTTGACCATCAGGTGCGAGCCAAATGC‑‑GGCCTCCACCGAGGTGGTCATGGCCAGGGCCTGGCCGTTGTCATCCACCGCACTGAGGTGAGAGGTGGAGATGCGCAGTGGCGAGCg                                                                                            >  7:89389/1‑149 (MQ=255)
                                                        gagaAGTCGGTCAGCTGATTGTTGAGCAGGAAGCCCTTGACCATCAGGTGCGAGCCAAATGC‑‑GGCCTCCACCGAGGTGGTCATGGCCAGGGCCTGGCCGTTGTCATCCACCGCACTGAGGTGAGAGGTGGAGATGCGCAGTGGCGAGCg                                                                                            <  2:15939/149‑1 (MQ=255)
                                                                ggTCAGCTGATTGTTGAGCAGGAAGCCCTTGACCATCAGGTGCGAGCCAAATGC‑‑GGCCTCCACCGAGGTGGTCATGGCCAGGGCCTGGCCGTTGTCATCCACCGCACTGAGGTGAGAGGTGGAGATGCGCAGTGGCGAGCGGTCGGgcg                                                                                    <  7:25300/149‑1 (MQ=255)
                                                                  tCAGCTGATTGTTGAGCAGGAAGCCCTTGACCATCAGGTGCGAGCCAAATGC‑‑GGCCTCCACCGAGGTGGTCATGGCCAGGGCCTGGCCGTTGTCATCCACCGCACTGAGGTGAGAGGTGGAGATGCGCAGTGGCGAGCGGTCGGgcgc                                                                                   >  2:230226/1‑148 (MQ=255)
                                                                    aGCTGATTGTTGAGCAGGAAGCCCTTGACCATCAGGTGCGAGCCAAATGC‑‑GGCCTCCACCGAGGTGGTCATGGCCAGGGCCTGGCCGTTGTCATCCACCGCACTGAGGTGAGAGGTGGAGATGCGCAGTGGCGAGCGGTCGGGCGCaa                                                                                 <  3:53948/148‑1 (MQ=255)
                                                                       tGATTGTTGAGCAGGAAGCCCTTGACCATCAGGTGCGAGCCAAATGC‑‑GGCCTCCACCGAGGTGGTCATGGCCAGGGCCTGGCCGTTGTCATCCACCGCACTGAGGTGAGAGGTGGAGATGCGCAGTGGCGAGCGGTCGGGCGCAAgcgc                                                                            >  3:39531/1‑149 (MQ=255)
                                                                       tGATTGTTGAGCAGGAAGCCCTTGACCATCAGGAGCCAGCCAAATGC‑‑GGCCTCCACCGAGGTGGTCATGGCCAGGGCCTGGCCGTTGTCATCCACCGCACTGAGGTGAGAGGTGGAGATGCGCAGTGGCGAGCGGTCGGGCGCAAgcgc                                                                            >  3:44453/1‑149 (MQ=255)
                                                                        gATTGTTGAGCAGGAAGCCCTTGACCATCAGGTGCGAGCCAAATGC‑‑GGCCTCCACCGAGGTGGTCATGGCCAGGGCCTGGCCGTTGTCATCCACCGCACTGAGGTGAGAGGTGGAGATGCGCAGTGGCGAGCGGTCGGGCGCAAgcgc                                                                            >  6:232513/1‑148 (MQ=255)
                                                                         aTTGTTGAGCAGGAAGCCCTTGACCATCAGGTGCGAGCCAAATGC‑‑GGCCTCCACCGAGGTGGTCATGGCCAGGGCCTGGCCGTTGTCATCCACCGCACTGAGGTGAGAGGTGGAGATGCGCAGTGGCGAGCGGTCGGGCGCAAgcgc                                                                            <  1:230226/147‑1 (MQ=255)
                                                                             ttgAGCAGGAAGCCCTTGACCATCAGGTGCGAGCCAAATGC‑‑GGCCTCCACCGAGGTGGTCATGGCCAGGGCCTGGCCGTTGTCATCCACCGCACTGAGGTGAGAGGTGGAGATGCGCAGTCGCGAGCGGTCGGGCGCAAGCGCCAGGTc                                                                      >  5:285406/1‑149 (MQ=255)
                                                                                aGCAGGAAGCCCTTGACCATCAGGTGCGAGCCAAATGC‑‑GGCCTCCACCGAGGTGGTCATGGCCAGGGCCTGGCCGTTGTCATCCACCGCACTGAGGTGAGAGGTGGAGATGCGCAGTGGCGAGCGGTCGGGCGCAAGCGCCAGGTCCGc                                                                   <  2:327437/149‑1 (MQ=255)
                                                                                  cAGGAAGCCCTTGACCATCAGGTGCGAGCCAAATGC‑‑GGCCTCCACCGAGGTGGTCATGGCCAGGGCCTGGCCGTTGTCATCCACCGCACTGAGGTGAGAGGTGGAGATGCGCAGTGGCGAGCGGTCGGGCGCAAGCGCCAGGTCCGCTc                                                                 >  7:178083/1‑149 (MQ=255)
                                                                                   aGGAAGCCCTTGACCATCAGGTGCGAGCCAAATGC‑‑GGCCTCCACCGAGGTGGTCATGGCCAGGGCCTGGCCGTTGTCATCCACCGCACTGAGGTGAGAGGTGGAGATGCGCAGTGGCGAGCGGTCGGGCGCAAGCGCCAGGTCCGCTcc                                                                >  2:254051/1‑149 (MQ=255)
                                                                                     gAAGCCCTTGACCATCAGGTGCGAGCCAAATGC‑‑GGCCTCCACCGAGGTGGTCATGGCCAGGGCCTGGCCGTTGTCATCCACCGCACTGAGGTGAGAGGTGGAGATGCGCAGTGGCGAGCGGTCGGGCGCAAGCGCCAGGTCCGCTCCtt                                                              >  4:125296/1‑149 (MQ=255)
                                                                                       aGCCCTTGACCATCAGGTGCGAGCCAAATGC‑‑GGCCTCCACCGAGGTGGTCATGGCCAGGGCCTGGCCGTTGTCATCCACCGCACTGAGGTGAGAGGTGGAGATGCGCAGTGGCGAGCGGTCGGGCGCAAGCGCCAGGTCCGCTCCTTGc                                                            <  5:109339/149‑1 (MQ=255)
                                                                                        gCCCTTGACCATCAGGTGCGAGCCAAATGC‑‑GGCCTCCACCGAGGTGGTCATGGCCAGGGCCTGGCCGTTGTCATCCACCGCACTGAGGTGAGAGGTGGAGATGCGCAGTGGCGAGCGGTCGGGCGCAAGCGCCAGGTCCGCTCCTTGCg                                                           >  2:228930/1‑149 (MQ=255)
                                                                                           cTTGACCATCAGGTGCGAGCCAAATGC‑‑GGCCTCCACCGAGGTGGTCATGGCCAGGGCCTGGCCGTTGTCATCCACCGCACTGAGGTGAGAGGTGGAGATGCGCAGTGGCGAGCGGTCGGGCGCAAGCGCCAGGTCCGCTCCTTGCGGCt                                                        >  6:78624/1‑149 (MQ=255)
                                                                                               acgaTCAGGTGCGAGCCATATGC‑‑GGCCTCCACCGTGGTGGTCATGGGGATGGCCTGTCCGTTGTCATCCACCGCACTGAGGTGAGAGGTGGAGATGCGCAGTGGCGAGCGGTCGGGCGCAAGCGCCAGGTCCGCTCCTTGCGGCTGGc                                                     <  6:285406/145‑1 (MQ=255)
                                                                                                 cATCAGGTGCTAGGAAAATGC‑‑GGCCTCCACCGAGGTGGTAATGGCCAGGGCCTGGCCGTTGTCATCCACCGCACTGAGGTGAGAGTTGGAGATGCGCAGTGGCGAGCGGTCGGGCGCAAGCGCCAGGTCCGCTCCTTGCGGCTGGCCCg                                                  <  4:44453/149‑1 (MQ=255)
                                                                                                 cATCAGGTGCGAGCCAAATGC‑‑GGCCTCCACCGAGGTGGTCATGGCCAGGGCCTGGCCGTTGTCATCCACCGCACTGAGGTGAGAGGTGGAGATGCGCAGTGGCGAGCGGTCGGGCGCAAGCGCCAGGTCCGCTCCTTGCGGCTGGCCCg                                                  <  4:39531/149‑1 (MQ=255)
                                                                                                     aGGTGCGAGCCAAATGC‑‑GGCCTCCACCGAGGTGGTCATGGCCAGGGCCTGGCCGTTGTCATCCACCGCACTGAGGTGAGAGGTGGAGATGCGCAGTGGCGAGCGGTCGGGCGCAAGCGCCAGGTCCGCTCCTTGCGGCTGGCCCGGAc                                               <  8:24525/148‑1 (MQ=255)
                                                                                                     aGGTGCGAGCCAAATGC‑‑GGCCTCCACCGAGGTGGTCATGGCCAGGGCCTGGCCGTTGTCATCCACCGCACTGAGGTGAGAGGTGGAGATGCGCAGTGGCGAGCGGTCGGGCGCAAGCGCCAGGTCCGCTCCTTGCGGCTGGCCCGGACg                                              >  8:348880/1‑149 (MQ=255)
                                                                                                             gCCAAATGC‑‑GGCCTCCACCGAGGTGGTCATGGCCAGGGCCTGGCCGTTGTCATCCACCGCACTGAGGTGAGAGGTGGAGATGCGCAGTGGCGAGCGGTCGGGCGCAAGCGCCCAGGTCCGCTCCTTGCGGCTGGCCCGGACGGgcgcgc                                       <  1:254051/149‑1 (MQ=255)
                                                                                                             gCCAAATGC‑‑GGCCTCCACCGAGGTGGTCATGGCCAGGGCCTGGCCGTTGACATCCACCGCACTGAGGTGAGAGGTGGAGATGCGCAGTGGCGAGCGGTCGGGCGCAAGCGCCAGGTCCGCTCCTTGCGGCTGGCCCGGACGGGCGCGCt                                      >  8:62355/1‑149 (MQ=255)
                                                                                                                 aaTGC‑‑GGCCTCCACCGAGGTGGTCATGGCCAGGGCCTGGCCGTTGTCATCCACCGCACTGAGGTGAGAGGTGGAGATGCGCAGTGGCGAGCGGTCGGGCGCAAGCGCCAGGTCCGCTCCTTGCGGCTGGCCCGGACGGGCGCGCTTCAt                                  >  6:383484/1‑149 (MQ=255)
                                                                                                                        ggCCTCCACCGAGGTGGTCATGGCCAGGGCCTGGCCGTTGTCATCCACCGCACTGAGGTGAGAGGTGGAGATGCGCAGTGGCGAGCGGTCGGGCGCAAGCGCCAGGTCCGCTCCTTGCGGCTGGCCCGGACGGGCGCGCTTCATGCtgt                             <  7:329059/149‑1 (MQ=255)
                                                                                                                             ccACCGAGGTGGTCATGGCCAGGGCCTGGCCGTTGTCATCCACCGCACTGAGGTGAGAGGTGGAGATGCGCAGTGGCGAGCGGTCGGGCGCAAGCGCCAGGTCCGCTCCTTGCGGCTGGCCCGGACGGGCGCGCTTCATGCTGTGCTCg                        >  2:97663/1‑149 (MQ=255)
                                                                                                                              cACCGAGGTGGTCATGGCCAGGGCCTGGCCGTTGTCATCCACCGCACTGAGGTGAGAGGTGGAGATGCGCAGTGGCGAGCGGTCGGGCGCAAGCGCCAGGTCCGCTCCTTGCGGCTGGCCCGGACGGGCGCGCtt                                     <  8:161581/135‑1 (MQ=255)
                                                                                                                              cACCGAGGTGGTCATGGCCAGGGCCTGGCCGTTGTCATCCACCGCACTGAGGTGAGAGGTGGAGATGCGCAGTGGCGAGCGGTCGGGCGCAAGCGCCAGGTCCGCTCCTTGCGGCTGGCCCGGACGGGCGCGCtt                                     >  7:161581/1‑135 (MQ=255)
                                                                                                                              cACCGAGGTGGTCATGGCCAGGGCCTGGCCGTTGTCATCCACCGCACTGAGGTGAGAGGTGGAGATGCGCAGTGGCGAGCGGTCGGGCGCAAGCGCCAGGTCCGCTCCTTGCGGCTGGCCCGGACGGGCGCGCTTCATGCTGTGCTCGc                       >  3:250131/1‑149 (MQ=255)
                                                                                                                              cACCGAGCTGGTCATGGCCAGGGCCTGGCCGTTGTCATCCACCGCAGTGAGGTGAGAGGTGGAGGTGCGCAGTGGCGAGCGGTCGGGCGCAAGCGCCAGGTCCGCTGCTGGCGGCTGGCCCGGACGGGCGGGCGTCATGCTGTGCTCGc                       >  4:101036/1‑149 (MQ=255)
                                                                                                                                ccGAGGTGGTCATGGCCAGGGCCTGGCCGTTGTCATCCACCGCACTGAGGTGAGAGGTGGAGATGCGCAGTGGCGAGCGGTCGGGCGCAAGCGCTAGGTCCGCTCCTTGCGGCTGGCCCGGACGGGCGCGCTTCATGCTGTGCTCGCCg                     <  1:97663/149‑1 (MQ=255)
                                                                                                                                   aGGTGGTCATGGCCAGGGCCTGGCCGTTGTCATCCACCGCACTGAGGTGAGAGGTGGAGATGCGCAGTGGCGAGCGGTCGGGCGCAAGCGCCAGGTCCGCTCCTTGCGGCTGGCCCGGACGGGCGCGCTTCATGCTGTGCTCGCCGAt                   >  4:89755/1‑148 (MQ=255)
                                                                                                                                           aTGGCCAGGGCCTGGCCGTTGTCATCCACCGCACTGAGGTGAGAGGTGGAGATGCGCAGTGGCGAGCGGTCGGGCGCAAGCGCCAGGTCCGCTCCTTGCGGCTGGCCCGGACGGGCGCGCTTCATGCTGTGCTCGCCGATCAGcgccgc          >  5:352175/1‑149 (MQ=255)
                                                                                                                                            cggCCAGGCCCTGGCCGTTGTCATCCACCGCACTGAGGTGAGAGGTGGAGATGCGCAGTGGCGAGCGGTCGGGCGCAAGCGCCAGGTCCGCTCCTTGCGGCTGGCCCGGACGGGCGCGCTTCATGCTGTGCTCGCCGATCAGCGCcgcg         >  4:303561/2‑149 (MQ=255)
                                                                                                                                               ccAGGGCCTGGCCGTTGTCATCCACCGCACTGAGTTGAGAGGTGGAGATGCGCAGTGGCGAGCGGTCGGGCGCAAGCGCCAGGTCCGCTCCTTGCGGCTGGCCCGGACGGGCGCGCTTCATGCTGTGCTCGCCGATCAGCGCCGCGCgg      >  1:50453/1‑149 (MQ=255)
                                                                                                                                               ccAGGGCCTGGCCGTTGTCATCCACCGCACTGAGGTGAGAGGTGGAGATGCGCAGTGGCGAGCGGTCGGGCGCAAGCGCCAGGTCCGCTCCTTGCGGCTGGCCCGGACGGGCGCGCTTCATGCTGTGCTCGCCGATCAGcgcc            >  2:170471/1‑143 (MQ=255)
                                                                                                                                               ccAGGGCCTGGCCGTTGTCATCCACCGCACTGAGGTGAGAGGTGGAGATGCGCAGTGGCGAGCGGTCGGGCGCAAGCGCCAGGTCCGCTCCTTGCGGCTGGCCCGGACGGGCGCGCTTCATGCTGTGCTCGCCGATCAGcgcc            <  1:170471/143‑1 (MQ=255)
                                                                                                                                               ccAGGGCCTGGCCGTTGTCATCCACCGCACTGAGGTGAGAGGTGGAGATGCGCAGTGGCGAGCGGTCGGGCGCAAGCGCCAGGTCCGCTCCTTGCGGCTGGCCCGGACGGGCGCGCTTCATGCTGTGCTCGCCGATCAGCGCCGCGCgg      >  5:307906/1‑149 (MQ=255)
                                                                                                                                                cAGGGCCTGGCCGTTGTCATCCACCGCACTGAGGTGAGAGGTGGAGATGCGCAGTGGCGAGCGGTCGGGCGCAAGCGCCAGGTCCGCTCCTTGCGGCTGGCCCGGACGGGCGCGCTTCATGCTGTGCTCGCCGATCAGCGCCGCGCGGc     >  3:352269/1‑149 (MQ=255)
                                                                                                                                                 aGGGCCTGGCCGTTGTCATCCACCGCACTGAGGTGAGAGGTGGAGATGCGCAGTGGCGAGCGGTCGGGCGCAAGCGCCAGGTCCGCTCCTTGCGGCTGGCCCGGACGGGCGCGCTTCa                                   <  3:208084/118‑1 (MQ=255)
                                                                                                                                                 aGGGCCTGGCCGTTGTCATCCACCGCACTGAGGTGAGAGGTGGAGATGCGCAGTGGCGAGCGGTCGGGCGCAAGCGCCAGGTCCGCTCCTTGCGGCTGGCCCGGACGGGCGCGCTTCa                                   >  4:208084/1‑118 (MQ=255)
                                                                                                                                                 aGGGCCTCCCAGTTGTCATCCACCGCACAGAGCTCAGAGGAGCAGATGCGAAGAGGCGAGCGGTCGGCCGCAAGCGCCAGGTCCGCTCCTTGCGCCTGGCCCGGACTGCCGCGCATCATACTGTGCTCGCCGATCAGCGCCGCGCGGCt    >  6:134546/1‑149 (MQ=255)
                                                                                                                                                   ggCCTGGCCGTTGTCATCCACCGCACTGAGGTGAGAGGTGGAGATGCGCAGTGGCGAGCGGTCGGGCGCAAGCGCCAGGTCCGCTCCTTGCGGCTGGCCCGGACGGGCGCGCTTCATGCTGTGCTCGCCGATCAGCGCCGCGCGGCTgg  <  3:125296/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GTTTGCCCGGCTGCACCCGATTGGCCACAGGCTTGCCGTGCTCGTGGGGGATGAACGAGAAGTCGGTCAGCTGATTGTTGAGCAGGAAGCCCTTGACCATCAGGTGCGAGCCAAATGC‑‑GGCCTCCACCGAGGTGGTCATGGCCAGGGCCTGGCCGTTGTCATCCACCGCACTGAGGTGAGAGGTGGAGATGCGCAGTGGCGAGCGGTCGGGCGCAAGCGCCAGGTCCGCTCCTTGCGGCTGGCCCGGACGGGCGCGCTTCATGCTGTGCTCGCCGATCAGCGCCGCGCGGCTGG  >  NC_002947/4006450‑4006743

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: