Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,717,845 T→G 10.8% N250H (AAC→CAC)  paaG ← 2‑(1,2‑epoxy‑1,2‑dihydrophenyl)acetyl‑CoA isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,717,8450TG10.8% 173.1 / 13.9 65N250H (AAC→CAC) paaG2‑(1,2‑epoxy‑1,2‑dihydrophenyl)acetyl‑CoA isomerase
Reads supporting (aligned to +/- strand):  ref base T (28/30);  new base G (4/3);  total (32/33)
Fisher's exact test for biased strand distribution p-value = 7.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.00e-01

GACGGGGTGGCCGGCCTGGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTC‑CTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCTGCGCTTGATCAGTGTCAGGCCGTAGGTGGGCTGGGCGGCGAGGT  >  NC_002947/3717700‑3717991
                                                                                                                                                 |                                                                                                                                                   
gACGGGGTGGCCGGCCTGGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCAc                                                                                                                                                  >  5:245146/1‑149 (MQ=255)
 aCGGGGTGGCCGGCCTGGGCCGGGACCTGGGCGTTACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACGTGCACATTGCTTGCGAGTGCTGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTGATGAAGGCGTTCACg                                                                                                                                                 <  7:71101/149‑1 (MQ=255)
  cGGGGTGGCCGGCCGGGGGCGCGACCTGGGCGATACCGGCGCCCCTGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGc                                                                                                                                                <  2:333596/149‑1 (MQ=255)
      ggggccggccTGGGCCGCGACCTGGGCGATCCCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTGAACGCCTTc                                                                                                                                            <  7:344736/147‑1 (MQ=255)
         gccggccTGGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTGAACGCCTTcat                                                                                                                                          >  1:82324/1‑146 (MQ=255)
         gccggccTGGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTGAACGCCTTcat                                                                                                                                          <  2:82324/148‑3 (MQ=255)
          gcggccCAGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCATCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCgg                                                                                                                                        >  7:158196/2‑149 (MQ=255)
          ccggccTGGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCgg                                                                                                                                        <  5:328352/149‑1 (MQ=255)
          ccggccTGGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCgg                                                                                                                                        <  2:47497/149‑1 (MQ=255)
           ctgccTGGGCCGCGACCTGGGCGATACCGTCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTGAACGCCTTcatgaa                                                                                                                                       <  5:39982/147‑6 (MQ=255)
              ccTGGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTcgc                                                                                                                                          >  1:265027/1‑143 (MQ=255)
              ccTGGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTcgc                                                                                                                                          <  2:265027/143‑1 (MQ=255)
               cTGGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTc                                                                                                                                   >  8:14841/1‑149 (MQ=255)
                   gCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTC‑CTCg                                                                                                                              <  5:29082/149‑1 (MQ=255)
                     cgcgACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTC‑CTCGCt                                                                                                                            >  5:277879/1‑149 (MQ=255)
                        gACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTC‑CTCGCTgcg                                                                                                                         >  6:283805/1‑149 (MQ=255)
                           cTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTC‑CTCGCTGCGccc                                                                                                                      <  5:231/149‑1 (MQ=255)
                              ggCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTC‑CTCGCTGCGCCCTGc                                                                                                                   <  3:142490/149‑1 (MQ=255)
                              ggCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTC‑CTCGCTGCGCCCTGc                                                                                                                   <  6:258347/149‑1 (MQ=255)
                                     ccGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTC‑CTCGCTGCGCCCTGCCAGgcgc                                                                                                            <  6:345475/149‑1 (MQ=255)
                                       ggCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGTCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTC‑CTCGCTGCGCCCTGCCAGGCGCTg                                                                                                          <  8:12495/149‑1 (MQ=255)
                                        gcgcCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTC‑CTCGCTGCGCCCTGCCAGGCGCTGc                                                                                                         <  2:188528/149‑1 (MQ=255)
                                            ccATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTC‑CTCGCTGCGCCCTGCCAGGCGCTGCAGGt                                                                                                     >  2:136293/1‑149 (MQ=255)
                                                     cGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTC‑CTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCgctcca                                                                                            <  1:147542/149‑1 (MQ=255)
                                                      ggCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTC‑CTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCgctccag                                                                                           >  4:171553/1‑149 (MQ=255)
                                                       gCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTC‑CTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCgctccagc                                                                                          <  6:277879/149‑1 (MQ=255)
                                                         aCCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGCTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTC‑CTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCgctccagctc                                                                                        <  4:26551/149‑1 (MQ=255)
                                                              tCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTCACTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCtccagctccagc                                                                                    <  1:284616/149‑1 (MQ=255)
                                                                aCCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTC‑CTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTgc                                                                                 >  5:109913/1‑149 (MQ=255)
                                                                              aTTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTC‑CTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCAtt                                                                   <  3:171553/149‑1 (MQ=255)
                                                                               ttgcttgcGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTC‑CTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTg                                                                  <  6:245146/149‑1 (MQ=255)
                                                                                 gcttgcGAGTGAGGCCATTTTCAGCGCCCCTTGAAGGGTGGCGTGCGTTTGTTGATGAAGGCGTTCACGCCTTCGCGGTAGTC‑CTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGt                                                                 <  6:364360/148‑1 (MQ=255)
                                                                                  cttgcGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTC‑CTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTc                                                                >  8:340913/1‑148 (MQ=255)
                                                                                  cttgcGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTC‑CTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACATTTGTc                                                                >  3:153325/1‑148 (MQ=255)
                                                                                     gcGAGGGCGGCCATTGTCCGCGCCCCTTGAAGGCTGGCGTGCGTTGTGTCATGAAGGCGTTCACGCCTTCGCGGTAGTC‑CTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCg                                                               <  8:158196/146‑1 (MQ=255)
                                                                                       gAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTC‑CTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACt                                                          >  8:164677/1‑149 (MQ=255)
                                                                                              gCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTC‑CTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCAtt                                                    >  3:263531/1‑148 (MQ=255)
                                                                                                cATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTC‑CTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAg                                                 >  6:39219/1‑149 (MQ=255)
                                                                                                 aTTTTCAGCGCCCCTTGAAGGCTGGCGTGGGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTC‑CTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAgg                                                >  2:311261/1‑149 (MQ=255)
                                                                                                 aTTTTCAGCGCCCCTTGAAGGCT‑GCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTC‑CTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGc                                               >  8:106115/1‑149 (MQ=255)
                                                                                                 aTTTTCAGCGCCCCTTGAAGGCT‑GCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTC‑CTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGc                                               >  8:100564/1‑149 (MQ=255)
                                                                                                          gcCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTC‑CTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCTGCGCtt                                        >  1:148843/1‑148 (MQ=255)
                                                                                                           cccgTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTC‑CTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCTGCGCTTGa                                      <  7:164677/149‑1 (MQ=255)
                                                                                                                 gAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTC‑CTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCTGCGCTTGATCAgt                                 <  7:340913/148‑1 (MQ=255)
                                                                                                                  aaGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCTGTAGTC‑CTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCTGCGCtt                                        <  5:250444/140‑1 (MQ=255)
                                                                                                                  aaGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTC‑CTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCTGCGCtt                                        >  6:250444/1‑140 (MQ=255)
                                                                                                                       tGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTC‑CTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCTGCGCTTGATCAGTGTCAGGc                          >  6:191251/1‑149 (MQ=255)
                                                                                                                         gcgtgcgtTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTC‑CTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCTGCGCTTGATCa                                   <  2:164516/138‑1 (MQ=255)
                                                                                                                         gcgtgcgtTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTC‑CTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCTGCGCTTGATCa                                   >  1:164516/1‑138 (MQ=255)
                                                                                                                            tgggtTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTA‑CTCGCTGCGCCCTGCCAGGGGCTGCAGGTCGCGCTCGATCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCTGCGCTTGATCAGTGTCTGTCCGTAg                     >  7:256567/4‑149 (MQ=255)
                                                                                                                            tgcgtTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTC‑CTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCTGCGCTTGATCAgt                                 >  1:154301/1‑137 (MQ=255)
                                                                                                                               gtTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTC‑CTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCTGCGCTTGATCAGTGTCAGGCCGTAGGtg                  >  5:302199/1‑149 (MQ=255)
                                                                                                                                    ttCATGAAGGCGTTCACGCCTTCGCGGTAGTC‑CTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCaaaa                                                                        >  3:24285/1‑90 (MQ=255)
                                                                                                                                    ttCATGAAGGCGTTCACGCCTTCGCGGTAGTC‑CTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCaaaa                                                                        <  4:24285/90‑1 (MQ=255)
                                                                                                                                       gcgAAGGCGTGAACGCCTTCGCGGTAGTC‑CTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCTGCGCTTGATCAGTGTCAGGCCGTAGGtgggctgggc          >  7:123264/3‑149 (MQ=255)
                                                                                                                                       aTGAAGGCGTTCACGCCTTCGCGGTAGTC‑CTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCTGCGCTTGATCAGTGTCAGGCCGTAGGtgggctggg           <  7:100564/148‑1 (MQ=255)
                                                                                                                                       aTGAAGGCGTTCACGCCTTCGCGGTAGTC‑CTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCTGCGCTTGATCAGTGTCAGGCCGTAGGtgggctggg           <  7:106115/148‑1 (MQ=255)
                                                                                                                                       aTGAAGGCGTTCACGCCTTCGCGGTAGTC‑CTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCTGCGCTTGATCAGTGTCAGGCCGTAGGtgggctggg           <  5:39219/148‑1 (MQ=255)
                                                                                                                                        cgAAGGCGTGAACGCCTTCGCGGTAGTC‑CTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCTGCGCTTGATCAGTGTCAGGCCGTAGGTGGGCTGggcg         >  2:109905/2‑149 (MQ=255)
                                                                                                                                         gAAGGCGTTCACGCCTTCGCGGTAGTC‑CTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCTGCGCTTGATCAGTGTCAGGCCGTAGGTGGGCTGggcgg        >  2:100357/1‑149 (MQ=255)
                                                                                                                                         gAAGGCGTGAACGCCTTCGCGGTAGTC‑CTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCTGCGCTTGATCAGTGTCAGGCCGTAGGTGGGCTGggcgg        >  3:237975/1‑149 (MQ=255)
                                                                                                                                          aaGGCGTTCACGCCTTCGCGGTAGTC‑CTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCTGCGCTTGATCAGTGTCAGGCCGTAGGTGGGCTGggcggc       <  5:185284/149‑1 (MQ=255)
                                                                                                                                           aGGCGTTCACGCCTTCGCGGTAGTC‑CTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCTGCGCTTGATCAGTGTCAGGCCGTAGGTGGGCTGggcggc       >  6:353138/1‑148 (MQ=255)
                                                                                                                                             gCGTTCACGCCTTCGCGGTAGTC‑CTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCTGCGCTTGATCAGTGTCAGGCCGTAGGTGGGCTGGGCGGCGAg    >  2:178923/1‑149 (MQ=255)
                                                                                                                                               gTTCACGCCTTCGCGGTAGTC‑CTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCTGCGCTTGATCAGTGTCAGGCCGTAGGTGGGCTGGGCGGCGAGGt  <  5:283805/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
GACGGGGTGGCCGGCCTGGGCCGCGACCTGGGCGATACCGGCGCCCATGGCGCCGGCACCGATCACCGCCACCTGCACATTGCTTGCGAGTGCGGCCATTTTCAGCGCCCCTTGAAGGCTGGCGTGCGTTTGTTCATGAAGGCGTTCACGCCTTCGCGGTAGTC‑CTCGCTGCGCCCTGCCAGGCGCTGCAGGTCGCGCTCCAGCTCCAGCTGCGCCTCAAAACCATTGTCGAAACTGGCATTGAGGCTGCGCTTGATCAGTGTCAGGCCGTAGGTGGGCTGGGCGGCGAGGT  >  NC_002947/3717700‑3717991

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑