Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,750,090 C→A 6.3% G134C (GGT→TGT)  PP_2406 ← quinate/shikimate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,750,0900CA 6.3% 200.7 / 3.8 64G134C (GGT→TGT) PP_2406quinate/shikimate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (30/30);  new base A (2/2);  total (32/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.72e-01

CCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCATCTGCACCACCTGGCGCCGTGCAACGTCGGGCAGGCCACGGCGCAAGCCTTCGGCAAAGCCCAGGCAATCGGTGTTGTGGCCGACGCGCTTGCCGTCCTTGAGCACCACGGTGTTTACCGCGCCGATACCACGGGCCTCGTCGGAC  >  NC_002947/2749942‑2750238
                                                                                                                                                    |                                                                                                                                                    
ccAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCAcc                                                                                                                                                      >  6:71173/1‑149 (MQ=255)
 cAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCAcc                                                                                                                                                      >  2:251041/1‑148 (MQ=255)
 cAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCAcc                                                                                                                                                      >  7:272904/1‑148 (MQ=255)
    gCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCATc                                                                                                                                                  >  6:4387/1‑149 (MQ=255)
               tCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCATCTGCACCACCTg                                                                                                                                       >  5:199083/1‑149 (MQ=255)
                    gcCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCATCTGCACCACCTGGCGcc                                                                                                                                  >  5:272640/1‑149 (MQ=255)
                          aTGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCATCTGCACCACCTGGCGCCGTGCaa                                                                                                                            >  5:215560/1‑149 (MQ=255)
                            gggTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACGCAGCAGGGCATGGGCCACGGCCGAGCCCGCGCCACCGGCACCCATCTGCACCCCCTGGCGCCGTGCAACg                                                                                                                          >  4:106240/1‑149 (MQ=255)
                               tgtTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCATCTGCACCACCTGGCGCCGTGCAACGTCg                                                                                                                       >  4:30567/1‑149 (MQ=255)
                                 ttCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCATCTGCACCACCTGGCGCCGTGCAACGTCggg                                                                                                                     >  4:260278/1‑149 (MQ=255)
                                    aGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCATCTGCACCACCTGGCGCCGTGCAACGTCGGGCa                                                                                                                   <  7:167096/148‑1 (MQ=255)
                                    aGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCATCTGCACCACCTGGCGCCGTGCAACGTCGGGCa                                                                                                                   <  2:180502/148‑1 (MQ=255)
                                    aGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTGCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCATCTGCACCACCTGGCGCCGTGCAACGTCGGGCa                                                                                                                   <  1:58795/148‑1 (MQ=255)
                                        tGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCATCTGCACCACCTGGCGCCGTGCAACGTCGGGCAGGCCa                                                                                                              >  7:183568/1‑149 (MQ=255)
                                             accaACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCATCTGCACCACCTGGCGCCGTGCAACGTCGGGCAGGCCACGgcg                                                                                                         >  2:319445/1‑149 (MQ=255)
                                             accaACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCATTCACCCAGCAGGGCATGGGCAACGGCCGAGCCAGCGCCACCGGCACCCATCTGCACCACCTGGCGCCGTGCAACGTCGGGCAGGCCACGgcg                                                                                                         >  7:77166/1‑149 (MQ=255)
                                                aaCGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCATCTGCACCACCTGGCGCCGTGCAACGTCGGGCAGGCCACGGCGCaa                                                                                                      <  5:71173/149‑1 (MQ=255)
                                                aaCGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCATCTGCACCACCTGGCGCCGTGCAACGTCGGGCAGGCCACGGCGCaa                                                                                                      <  6:21053/149‑1 (MQ=255)
                                                     cTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCATCTGCACCACCTGGCGCCGTGCAACGTCGGGCAGGCCACGGCGCAAGCCtt                                                                                                 >  5:26674/1‑149 (MQ=255)
                                                        cgcACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCATCTGCACCACCTGGCGCCGTGCAACGTCGGGCAGGCCACGGCGCAAGCCTTCg                                                                                               <  3:260278/148‑1 (MQ=255)
                                                           aCGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCATCTGCACCACCTGGCGCCGTGCAACGTCGGGCAGGCCACGGCGCAAGCCTTCGGCa                                                                                            >  6:41598/1‑148 (MQ=255)
                                                            cGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCATCTGCACCACCTGGCGCCGTGCAACGTCGGGCAGGCCACGGCGCAAGCCTTCGGCaaa                                                                                          <  8:272904/149‑1 (MQ=255)
                                                             gggTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCATCTGCACCACCTGGCGCCGTGCAACGTCGGGCAGGCCACGGCGCAAGCCTTCGGCAAAg                                                                                         <  5:207636/149‑1 (MQ=255)
                                                               gTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCATCTGCACCACCTGGCGCCGTGCAACGTCGGGCAGGCCACGGCGCAAGCCTTCGGCAAAGcc                                                                                       >  4:129455/1‑149 (MQ=255)
                                                                 cGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCATCTGCACCTCCTGGCGCCGTGCAACGTCGGGCAGGCCACGGCGCAAGCCTTCGGCAAAGCCCa                                                                                     <  1:251041/149‑1 (MQ=255)
                                                                  gCATCCACTTCGAACAGCTCGAGCCGTTCTACCCCTTGTCCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCTTCTGCACCACCTGGCGCCGTGCAACGTCGGGCAGGCCACGGCGCAAGCCTTCGGCAAAGCCCAg                                                                                    <  8:77166/149‑1 (MQ=255)
                                                                  gCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCATCTGCACCACCTGGCGCCGTGCAACGTCGGGCAGGCCACGGCGCAAGCCTTCGGCAAAGCCCAg                                                                                    <  5:321928/149‑1 (MQ=255)
                                                                   cATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCATCTGCACCACCTGGCGCCGTGCAACGTCGGGCAGGCCACGGCGCAAGCCTTCGGCAAAGCCCAgg                                                                                   <  7:24847/149‑1 (MQ=255)
                                                                    aTCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACACATCTGCACCACCTGGCGCCGTGCAACGTCGGGCAGGCCACGGCGCAAGCCTTCGGCAAAGCCCAg                                                                                    <  1:100099/147‑1 (MQ=255)
                                                                    aTCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACACATCTGCACCACCTGGCGCCGTGCAACGTCGGGCAGGCCACGGCGCAAGCCTTCGGCAAAGCCCAg                                                                                    >  2:100099/1‑147 (MQ=255)
                                                                    aTCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACACATCTGCACCACCTGGCGCCGTGCAACGTCGGGCAGGCCACGGCGCAAGCCTTCGGCAAAGCCCAg                                                                                    >  2:104014/1‑147 (MQ=255)
                                                                    aTCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCCCCACCGGCACACATCTGCACCACCTGGCGCCGTGCAACGTCGGGCAGGCCACGGCGCAAGCCTTCGGCAAAGCCCAg                                                                                    <  1:104014/147‑1 (MQ=255)
                                                                    aTCCACTTCGAACAGCACCAGCCCTTCTACCACTTCACACAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCATCTGCACCACCTGGCGCCGTGCAACGTCGGGCAGGCCACGGCGCAAGCCTTCGGCAAAGCCCAgg                                                                                   >  1:314374/1‑148 (MQ=255)
                                                                     tacACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCATCTGCAACACCTGGCGCCGTGCAACGTCCGGCAGGCCACGGCGCAAGCCTTCGGCAAATCCCACGc                                                                                  >  1:181629/3‑148 (MQ=255)
                                                                       cACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCATCTGCACCACCTGGCGCCGTGCAACGTCGGGCAGGCCACGGCGCAAGCCTTCGGCAAAGCCCAGGCAAt                                                                               <  6:26674/149‑1 (MQ=255)
                                                                       cACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCATCTGCACCACCTGGCGCCGTGCAACGTCGGGCAGGCCACGGCGCAAGCCTTCGGCAAAGCCCAGGCAAt                                                                               <  6:199083/149‑1 (MQ=255)
                                                                          ttCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCATCTGCACCACCTGGCGCCGTGCAACGTCGGGCAGGCCACGGCGCAAGCCTTCGGCAAAGCCCAGGCAATCgg                                                                            >  4:135921/1‑149 (MQ=255)
                                                                              aaCAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCATCTGCACCACCTGGCGCCGTGCAACGTCGGGCAGGCCACGGCGCAAGCCTTCGGCAAAGCCCAGGCAATCGGtgtt                                                                        <  1:319445/149‑1 (MQ=255)
                                                                               aCAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCATCTGCACCACCTGGCGCCGTGCAACGTCTGGCAGGCCACGGCGCAAGCCTTCGGCAAACCCCAAGCAATCGGtgttg                                                                       >  8:203566/1‑149 (MQ=255)
                                                                                 aGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCATCTGCACCACCTGGCGCCGTGCAACGTCGGGCAGGCCACGGCGCAAGCCTTCGGCAAAGCCCAGGCAATCGGTGTtgtg                                                                     <  1:318629/149‑1 (MQ=255)
                                                                                    aCCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCATCTGCACCACCTGGCGCCGTGCAACGTCGGGCAGGCCACGGCGCAAGCCTTCGGCAAAGCCCAGGCAATCGGTGTTGTGGcc                                                                  >  8:30705/1‑149 (MQ=255)
                                                                                      cAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCATCTGCACCACCTGGCGCCGTGCAACGTCGGGCAGGCCACGGCGCAAGCCTTCGGCAAAGCCCAGGCAATCGGTGTTGTGGCCGa                                                                <  3:129455/149‑1 (MQ=255)
                                                                                        gCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCATCTGCACCACCTGGCGCCGTGCAACGTCGGGCAGGCCACGGCGCAAGCCTTCGGCAAAGCCCAGGCAATCGGTGTTGTGGCCGAcg                                                              >  5:216150/1‑149 (MQ=255)
                                                                                         ccGGTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCATCTGCCCCACCTGGCGCCGTGCAACGTCGGGCAGGCCACGGCGCAAGCCTTCGGCAAAGCCCAGGCAATCGGTGTTGTGGCCGAcgc                                                             <  5:4387/149‑1 (MQ=255)
                                                                                           gTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCATCTGCACCACCTGGCGCCGTGCAACGTCGGGCAGGCCACGGCGCAAGCCTTCGGCAAAGCCCAGGCAATCGGTGTTGTGGCCGACgcgc                                                           <  3:135921/149‑1 (MQ=255)
                                                                                            ttCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCATCTGCACCACCTGGCGCCGTGCAACGTCGGGCAGGCCACGGCGCAAGCCTTCGGCAAAGCCCAGGCAATCGGTGTTGTGGCCGACGCGCt                                                          >  3:192557/1‑149 (MQ=255)
                                                                                              cTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCATCTGCACCACCTGGCGCCGTGCAACGTCGGGCAGGCCACGGCGCAAGCCTTCGGCAAAGCCCAGGCAATCGGTGTTGTGGCCGACGCGCTTg                                                        >  6:22387/1‑149 (MQ=255)
                                                                                                    cTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCATCTGCACCACCTGGCGCCGTGCAACGTCGGGCAGGCCACGGCGCAAGCCTTCGGCAAAGCCCAGGCAATCGGTGTTGTGGCCGACGCGCTTGCCGTcc                                                  >  5:26527/1‑149 (MQ=255)
                                                                                                    cTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCATCTGCACCACCTGGCGCCGTGCAACGTCGGGCAGGCCACGGCGCAAGCCTTCGGCAAAGCCCAGGCAATCGGTGTTGTGGCCGACGCGCTTGCCGTcc                                                  >  3:246101/1‑149 (MQ=255)
                                                                                                    cTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCATCTGCACCACCTGGCGCCGTGCAACGTCGGGCAGGCCACGGCGCAAGCCTTCGGCAAAGCCCAGGCAATCGGTGTTGTGGCCGACGCGCTTGCCGTcc                                                  >  3:234477/1‑149 (MQ=255)
                                                                                                       ctcccAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCATCTGCACGACCTGGCGCCGTGCAACGTCGGGCAGGGCACGGCGGAAGCCTTCGGCAAAGCGCAGGCAATCGGTGTTGTGGCGGACGCGCTTGCCGTCCt                                                 <  7:203566/145‑1 (MQ=255)
                                                                                                        aCCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCATCTGcaccac                                                                                                                                          >  5:266716/1‑57 (MQ=255)
                                                                                                        aCCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCATCTGcaccac                                                                                                                                          <  6:266716/57‑1 (MQ=255)
                                                                                                        aCCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCATCTGCACCACCTGGCGCCGTGCAACGTCGGGCAGGCCACGGCGCAAGCCTTCGGCAAAGCCCAGGCAATCGGTGTTGTGGCCGACGCGCTTGCCGTCCTTGa                                              <  6:272640/149‑1 (MQ=255)
                                                                                                        aCCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCATCTGCACCACCTGGCGCCGTGCAACGTCGGGCAGGCCACGGCGCAAGCCTTCGGCAAAGCCCAGGCAATCGGTGTTGTGGCCGACGCGCTTGCCGTCCTTGa                                              <  3:30567/149‑1 (MQ=255)
                                                                                                             gcagGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCATCTGCCCCACCTGGCGCCGTGCAACGTCGGGCAGGCCACGGCGCAAGCCTTCGGCAAAGCCCAGGCAATCGGTGTTGTGGCCGACGCGCTTGCCGTCCTTGAGcacc                                         >  1:154753/1‑149 (MQ=255)
                                                                                                               agGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCATCTGCACCACCTGGCGCCGTGCAACGTCGGGCAGGCCACGGCGCAAGCCTTCGGCAAAGCCCAGGCAATCGGTGTTGTGGCCGACGCGCTTGCCGTCCTTGAGcaccac                                       >  4:256120/1‑149 (MQ=255)
                                                                                                                   cATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCATCTGCACCACCTGGCGCCGTGCAACGTCGGGCAGGCCACGGCGCAAGCCTTCGGCAAAGCCCAGGCAATCGGTGTTGTGGCCGACGCGCTTGCCGTCCTTGAGCACCACGgtg                                   <  6:26527/149‑1 (MQ=255)
                                                                                                                     tGGGCCACGGCCGAGCCAGCGCCACCGGCACCCATCTGCACCACCTGGCGCCGTGCAACGTCGGGCAGGCCACGGCGCAAGCCTTCGGCAAAGCCCAGGCAATCGGTGTTGTGGCCGACGCGCTTGCCGTCCTTGAGCACCACGGTGtt                                 >  3:290593/1‑149 (MQ=255)
                                                                                                                        gCCACGGCCGAGCCAGCGCCACCGGCACCCATCTGCACCACCTGGCGCCGTGCAACGTCGGGCAGGCCACGGCGCAAGCCTTCGGCAAAGCCCAGGCAATCGGTGTTGTGGCCGACGCGCTTGCCGTCCTTGAGCACCACGGTGTTTAc                              <  2:240373/149‑1 (MQ=255)
                                                                                                                         ccACGGCCGAGCCAGCGCCACCGGCACCCATCTGCACCACCTGGCGCCGTGCAACGTCGGGCAGGCCACGGCGCAAGCCTTCGGCAAAGCCCAGGCAATCGGTGTTGTGGCCGACGCGCTTGCCGTCCTTGAGCACCACGGTGTTTAcc                             <  4:234477/149‑1 (MQ=255)
                                                                                                                          cACGGCCGAGCCAGCGCCACCGGCACCCATCTGCACCACCTGGCGCCGTGCAACGTCGGGCAGGCCACGGCGCAAGCCTTCGGCAAAGCCCAGGCAATCGGTGTTGTGGCCGACGCGCTTGCCGTCCTTGAGCACCACGGTGTTTACcg                            >  6:100862/1‑149 (MQ=255)
                                                                                                                                cGAGCCAGCGCCACCGGCACCCATCTGCACCACCTGGCGCCGTGCAACGTCGGGCAGGCCACGGCGCAAGCCTTCGGCAAAGCCCAGGCAATCGGTGTTGTGGCCGACGCGCTTGCCGTCCTTGAGCACCACGGTGTTTACCGCGCCGa                      <  6:216150/149‑1 (MQ=255)
                                                                                                                                      aGCGCCACCGGCACCCATCTGCACCACCTGGCGCCGTGCAACGTCGGGCAGGCCACGGCGCAAGCCTTCGGCAAAGCCCAGGCAATCGGTGTTGTGGCCGACGCGCTTGCCGTCCTTGAGCACCACGGTGTTTACCGCGCCGATACCAc                >  4:28422/1‑149 (MQ=255)
                                                                                                                                        cgcCACCGGCACCCATCTGCACCTCCTGGCGCCGTGCAAAGTCGGGCAGGACACGGCGCAAGCCTTCGGCAAAGCCCAGGCAATCGGTGTTGTTGCCGACGCGCTTGCCGTCCTTGAGCACCACGGTGTTTACCGCGCCGATACCACgg              <  2:181629/149‑1 (MQ=255)
                                                                                                                                        cgcCACCGGCACCCATCTGCACCACCTGGCGCCGTGCAACGTCGGGCAGGCCACGGCGCAAGCCTTCGGCAAAGCCCAGGCAATCGGTGTTGTGGCCGACGCGCTTGCCGTCCTTGAGCACCACGGTGTTTACCGCGCCGATACCACgg              <  8:183568/149‑1 (MQ=255)
                                                                                                                                                 cACCCATCTGCACCACCTGGCGCCGTGCAACGTCGGGCAGGCCACGGCGCAAGCCTTCGGCAAAGCCCAGGCAATCGGTGTTGTGGCCGACGCGCTTGCCGTCCTTGAGCACCACGGTGTTTACCGCGCCGATACCACGGGCCtcgtcg     <  5:41598/149‑1 (MQ=255)
                                                                                                                                                    ccATCTGCACCACCTGGCGCCGTGCAACGTCGGGCAGGCCACGGCGCAAGCCTTCGGCAAAGCCCAGGCAATCGGTGTTGTGGCCGACGCGCTTGCCGTCCTTGAGCACCACGGTGTTTACCGCGCCGATACCACGGGCCTCGTCGGAc  <  3:28422/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCATCTGCACCACCTGGCGCCGTGCAACGTCGGGCAGGCCACGGCGCAAGCCTTCGGCAAAGCCCAGGCAATCGGTGTTGTGGCCGACGCGCTTGCCGTCCTTGAGCACCACGGTGTTTACCGCGCCGATACCACGGGCCTCGTCGGAC  >  NC_002947/2749942‑2750238

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: ‑