Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,845,927 C→T 6.8% G257D (GGC→GAC)  PP_3395 ← LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,845,9270CT 6.8% 175.1 / 2.4 59G257D (GGC→GAC) PP_3395LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (29/26);  new base T (2/2);  total (31/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.19e-01

AGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCGTCAGTCAGCGCAATGACCCGCAGGTCGAACAACGCCGCATAGCGGGACACGGCCTCTTGCGGGAACACCGCCAGGCCAAGGCCAAACTGCACCAGACGGTATGCCGCGTCGAACGACGACACGTACGAACGAAAACGCAGCGCCTTGCCTGACTGTTCGGCGCTTTTACGCATGAAAGCGTTCAATGTGGCAAAGGCCGATAGCCCCAGGTGGTCGTAACCCAG  >  NC_002947/3845780‑3846075
                                                                                                                                                   |                                                                                                                                                     
aGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCGTCAGTCAGCGCAATGACCCGCAGGTCGAACAACGCCGCATAGCGGGACACGGCCTCTTGCGGGAACACCGCCaggc                                                                                                                                                       >  6:365569/1‑148 (MQ=255)
   cgcCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCGTCAGTCAGCGCAATGACCCGCAGGTCGAACAACGCCGCATAGCGGGACACGGCCTCTTGCGGGAACACCGCCaggccaag                                                                                                                                                   <  3:354475/149‑1 (MQ=255)
   cgcCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCGTCAGTCAGCGCAATGACCCGCAGGTCGAACAACGCCGCATAGCGGGACACGGCCTCTTGCGGGAACACCGCCaggccaag                                                                                                                                                   <  4:112527/149‑1 (MQ=255)
     ccgagccgaCAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCGTCAGTCAGCGCAATGACCCGCAGGTCGAACACCGCCGCATAGCGGGACACGGCCTCTTGCGGGAACACCGCCaggccaaggc                                                                                                                                                 >  4:3799/1‑149 (MQ=255)
      cgagccgaCAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCGTCAGTCAGCGCAATGACCCGCAGGTCGAACAACGCCGCATAGCGGGACACGGCCTCTTGCGGGAACACCGCCaggccaagg                                                                                                                                                  <  4:87211/147‑1 (MQ=255)
       gagccgaCAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCGTCAGTCAGCGCAATGACCCGCAGGTCGAACAACGCCGCATAGCGGGACTCGGCCTCTTGCGGGAACACCGCCaggccaaggcca                                                                                                                                               <  8:199093/149‑1 (MQ=255)
       gagccgaCAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCGTCAGTCAGCGCAATGACCCGCAGGTCGAACAACGCCGCATAGCGGGACACGGCCTCTTGCGGGAACACCGCCaggccaaggcc                                                                                                                                                >  4:255923/1‑148 (MQ=255)
           cgaCAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCGTCAGTCAGCGCAATGACCCGCAGGTCGAACAACGCCGCATAGCGGGACACGGCCTCTTGCGGGAACACCGCCAGGCCAAGGCCAAACt                                                                                                                                           >  2:409021/1‑149 (MQ=255)
             aCAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCGTCAGTCAGCGCAATGACCCGCAGGTCGAACAACGCCGCATAGCGGGACACGGCCTCTTGCGGGAACACCGCCAGGCCAAGGCCAAACTg                                                                                                                                          >  1:122201/1‑148 (MQ=255)
             aCAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCGTCAGTCAGCGCAATGACCCGCAGGTCGAACAACGCCGCATAGCGGGACACGGCCTCTTGCGGGAACACCGCCAGGCCAAGGCCAAACTGc                                                                                                                                         <  1:79365/149‑1 (MQ=255)
                  gagaGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCGTCAGTCAGCGCAATGACCCGCAGGTCGAACAACGCCGCATAGCGGGACACGGCCTCTTGCGGGAACACCGCCAGGTCAAGGCCAAACTGCAc                                                                                                                                       >  1:164419/1‑146 (MQ=255)
                  gagaGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCGTCAGTCAGAGCAATGACCCGCAGGTCGAACAACGCCGCATAGCGGGACACGGCCTCTTGCGGGAACACCGCCAGGTCAAGGCCAAACTGCAc                                                                                                                                       >  1:170502/1‑146 (MQ=255)
                          gtgCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCGTCAGTCAGCGCAATGACCCGCAGGTCGAACAACGCCGCATAGCGGGACACGGCCTCTTGCGGGAACACCGCCAGGCCAAGGCCAAACTGCACCAGACGGTATg                                                                                                                            >  2:354822/1‑149 (MQ=255)
                                gTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCGTCAGTCAGCGCAATGACCCGCAGGTCGAACAACGCCGCATAGCGGGACACGGCCTCTTGCGGGAACACCGCCAGGCCAAGGCCAAACTGCACCAGACGGTATGCCGCGt                                                                                                                      <  2:204645/149‑1 (MQ=255)
                                 tCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCGTCAGTCAGCGCAATGACCCGCAGGTCGAACAACGCCGCATAGCGGGACACGGCCTCTTGCGGGAACACCGCCAGGCCAAGGCCAAACTGCACCAGACGGTATGCCGCGTc                                                                                                                     >  6:365942/1‑149 (MQ=255)
                                    cgcATGCAGATGACGAATTCACCAAGCGCCCAGTCGTCAGTCAGCGCAATGACCCGCAGGTCGAACAACGCCGCATAGCGGGACACGGCCTCTTGCGGGAACACCGCCAGGCCAAGGCCAAACTGCACCAGACGGTATGCCGCGTCGaa                                                                                                                  <  3:255923/149‑1 (MQ=255)
                                    cgcATGCAGATGACGAATTCACCAAGCGCCCAGTCGTCAGTCAGCGCAATGACCCGCAGGTCGAACAACGCCGCATAGCGGGACACGGCCTCTTGCGGGAACACCGCCAGGCCAAGGCCAAACTGCACCAGACGGTATGCCGCGTCGaa                                                                                                                  <  4:36251/149‑1 (MQ=255)
                                     gcATGCAGATGACGAATTCACCAAGCGCCCAGTCGTCAGTCAGCGCAATGACCCGCAGGTCGAACAACGCCGCATAGCGGGACACGGCCTCTTGCGGGAACACCGCCAGGCCAAGGCCAAACTGCACCAGACGGTATGCCGCGTCGAac                                                                                                                 <  1:354822/149‑1 (MQ=255)
                                        tGCAGATGACGAATTCACCAAGCGCCCAGTCGTCAGTCAGCGCAATGACCCGCAGGTCGAACAACGCCGCATAGCGGGACACGGCCTCTTGCGGGAACACCGCCAGGCCAAGGCCAAACTGCACCAGACGGTATGCCGCGTCGAacgac                                                                                                              <  5:365569/149‑1 (MQ=255)
                                                 cGAATTCACCAAGCGCCCAGTCGTCAGTCAGCGCAATGACCCGCAGGTCGAACAACGCCGCATAGCGGGACACGGCCTCTTGCGGGAACACCGCCGGGCCAAGGCCAAACTGCACCAGACGGTATGCCGCGTCGAACGACGACACGTac                                                                                                    <  2:359881/149‑1 (MQ=255)
                                                  gAATTCACCAAGCGCCCAGTCGTCAGTCAGCGCAATGACCCGCAGGTCGAACAACGCCGCATAGCGGGACACGGCCTCTTGCGGGAACACCGCCAGGCCAAGGCCAAACTGCACCAGACGGTATGCCGCGTCGAACGACGACACGTac                                                                                                    >  1:158296/1‑148 (MQ=255)
                                                  gAATTCACCAAGCGCCCAGTCGTCAGTCAGCGCAATGACCCGCAGGTCGAACAACGCCGCATAGCGGGACACGGCCTCTTGCGGGAACACCGCCAGGCCAAGGCCAAACTGCACCAGACGGTATGCCGCGTCGAACGACGACACGTac                                                                                                    >  2:170963/1‑148 (MQ=255)
                                                      tCACCAAGCGCCCAGTCGTCAGTCAGCGCAATGACCCGCAGGTCGAACAACGCCGCATAGCGGGACACGGCCTCTTGCGGGAACACCGCCAGGCCAAGGCCAAACTGCACCAGACGGTATGCCGCGTCGAACGACGACACGTacgaac                                                                                                >  2:268419/1‑148 (MQ=255)
                                                          cAAGCGCCCAGTCGTCAGTCAGCGCAATGACCCGCAGGTCGAACAACGCCGCATAGCGGGACACGGCCTCTTGCGGGAACACCGCCAGGCCAAGGCCAAACTGCACCAGACGGTATGCCGCGTCGAACGACGACACGTACGAACGaaaa                                                                                           <  1:268419/149‑1 (MQ=255)
                                                              cgcCCAGTCGTCAGTCAGCGCAATGACCCGCAGGTCGAACAACGCCGCATAGCGGGACACGGCCTCTTGCGGGAACACCGCCAGGCCAAGGCCAAACTTCACCAGACGGTATGCCGCGTCGAACGACGACACGTACGAACGAAAACGCa                                                                                       <  8:266152/149‑1 (MQ=255)
                                                                cccAGTCGTCAGTCAGCGCAATGACCCGCAGGTCGAACAACGCCGCATAGCGGGACACGGCCTCTTGCGGGAACACCGCCAGGCCAAGGCCAAACTGCACCAGACGGTATGCCGCGTCGAACgacgac                                                                                                           >  1:87942/1‑128 (MQ=255)
                                                                cccAGTCGTCAGTCAGCGCAATGACCCGCAGGTCGAACAACGCCGCATAGCGGGACACGGCCTCTTGCGGGAACACCGCCAGGCCAAGGCCAAACTGCACCAGACGGTATGCCGCGTCGAACgacgac                                                                                                           <  2:87942/128‑1 (MQ=255)
                                                                     tcgtcAGTCAGCGCAATGACCCGCAGGTCGAACAACGCCGCATAGCGGGACACGGCCTCTTGCGGGAACACCGCCCGGCCAAGGCCAAACTGCACCAGACGGTATGCCGCGTCGAACGACGACACGTACGAACGAAAACGCAGCGCCtt                                                                                >  8:54204/1‑149 (MQ=255)
                                                                       gtcagtcaGCGCAATGACCCGCAGGTCGAACAACGCCGCATAGCGGGACACGGCCTCTTGCGGGAACACCGCCAGGCCAAGGCCAAACTGCACCAGACGGTATGCCGCGTCGAACGACGACACGTacgaa                                                                                                 <  1:117165/130‑1 (MQ=255)
                                                                       gtcagtcaGCGCAATGACCCGCAGGTCGAACAACGCCGCATAGCGGGACACGGCCTCTTGCGGGAACACCGCCAGGCCAAGGCCAAACTGCACCAGACGGTATGCCGCGTCGAACGACGACACGTacgaa                                                                                                 >  2:117165/1‑130 (MQ=255)
                                                                       gtcagtcaGCGCAATGACCCGCAGGTCGAACAACGCCGCATAGCGGGACACGGCCTCTTGCGGGAACACCGCCAGGCCAAGGCCAAACTGCACCAGACGGTATGCCGCGTCGAACGACGACACGTACGAACGAAAACGCAGCGCCTTGc                                                                              <  1:409021/149‑1 (MQ=255)
                                                                       gtcagtcaGCGCAATGACCCGCAGGTCGAACAACGCCGCATAGCGGGACACGGCCTCTTGCGGGAACACCGCCAGGCCAAGGCCAAACTGCACCAGACGGTATGCCGCGTCGAACGACGACACGTACGAACGAAAACGCAGCGCCTTGc                                                                              >  7:330529/1‑149 (MQ=255)
                                                                       gtcagtcaGCGCAATGACCCGCAGGTCGAACAACGCCGCATAGCGGGACACGGCCTCTTGCGGGAACACCGCCAGGCCAAGGCCAAACTGCACCAGACGGTATGCCGCGTCGAACGACGACACGTACGAACGAAAACGCAGCGCCTTGc                                                                              >  3:391279/1‑149 (MQ=255)
                                                                        tcagtcagCGCAATGACCCGCAGGTCGAACAACGCCGCATAGCGGGACACGGCCTCTTGCGGGAACACCGCCAGGTCAAGGCCAAACTGCACCAGACGGTATGCCGCGTCGAACGACGACACGTACGAACGAAAACGCAGCGCCTTGcc                                                                             <  2:164419/149‑1 (MQ=255)
                                                                        tcagtcagCGCAATGACCCGCAGGTCGAACAACGCCGCATAGCGGGACACGGCCTCTTGCGGGAACACCGCCAGGTCAAGGCCAAACTGCACCAGACGGTATGCCGCGTCGAACGACGACACGTACGAACGAAAACGCAGCGCCTTGcc                                                                             <  2:170502/149‑1 (MQ=255)
                                                                           gtcagCGCAATGACCCGCAGGTCGAACAACGCCGCATAGCGGGACACGGCCTCTTGCGGGAACACCGCCAGGCCAAGGCCAAACTGAACCAGACGGAATGCCGCGTCGAACGACGACACGTACGACAGAAAACGCAGCGCCTTGCCTGa                                                                          >  6:225500/1‑149 (MQ=255)
                                                                            tcagCGCAATGACCCGCAGGTCGAACAACGCCGCATAGCGGGACACGGCCTCTTGCGGGAACACCGCCAGGCCAAGGCCAAACTGCACCAGACGGTATGCCGCGTCGAACGACGACACGTACGAACGAAAACGCAGCGCCTTGCCTGAc                                                                         >  4:392546/1‑149 (MQ=255)
                                                                              agCGCAATGACCCGCAGGTCGAACAACGCCGCATAGCGGGACACGGCCTCTTGCGGGAACACCGCCAGGCCAAGGCCAAACTGCACCAGACGGTATGCCTCGTCGAACGACGACACGTACGAACGAAAACGCAgcgc                                                                                   <  3:275347/137‑1 (MQ=255)
                                                                              agCGCAATGACCCGCAGGTCGAACAACGCCGCATAGCGGGACACGGCCTCTTGCGGGAACACCGCCAGGCCAAGGCCAAACTGCACCAGACGGTATGCCGCGTCGAACGACGACACGTACGAACGAAAACGCAgcgc                                                                                   >  4:275347/1‑137 (MQ=255)
                                                                                cgcAATGACCCGCAGGTCGAACAACGCCGCATAGCGGGACACGGCCTCTTGCCGGAACACCGACAGGCCAAGGCCAAACTGCACCAGACGGTATCCCGCGACGAACGACGACACGTAAGAACGAAAACGAAGCGCCTTCCCTGACTGt                                                                      >  7:12406/1‑148 (MQ=255)
                                                                                         ccGCAGGTCGAACAACGCCGCATAGCGGGACACGGCCTCTTGCGGGAACACCGCCAGGCCAAGGCCAAACTGCACCAGACGGTATGCCGCGTCGAACGACGACACGTACGAACGAAAACGCAGCGCCTTGCCTGACTGTTCGGCGCttt                                                            >  7:139065/1‑149 (MQ=255)
                                                                                          cGCAGGTCGAACAACGCCGCATAGCGGGACACGGCCTCTTGCGGGAACACCGCCAGGCCAAGGCCAAACTGCACCAGACGCTATGCCGCGTCGAACCACGACAAGTAAGAACGAAAACGCAGCGCCTTGCATGACTGTTCGGCGCttt                                                            >  7:322564/1‑148 (MQ=255)
                                                                                           gCAGGTCGAACAACGCCGCATAGCGGGACACGGCCTCTTGCGGGAACACCGCCAGGCCAAGGCCAAACTGCACCAGACGGTATGCCGCGTCGAACGACGACACGTACGAACGAAAACGCAGCGCCTTGCCTGACTGTTCGGCGCtttt                                                           >  4:23549/1‑148 (MQ=255)
                                                                                            cAGGTCGAACAACGCCGCATAGCGGGACACGGCCTCTTGCGGGAACACCGCCAGGCCAAGGCCAAACTGCACCAGACGGTATGCCGCGTCGAACGACGACACGTACGAACGAAAACGCAGCGCCTTGCCTGACTGTTCGGCGCTTTTAc                                                         <  8:330529/149‑1 (MQ=255)
                                                                                               gTCGAACAACGCCGCATAGTGGGACACGGCCTCTTGCGGGAACACCGCCAGGCCAAGGCCAAACTGCACCAGACGGTATGCCGCGTCGAACGACGACACGTACGAACGAAAACGCAGCGCCTTGCCTGACTGTTCGACGCTTTTACGc                                                       >  3:365818/1‑148 (MQ=255)
                                                                                                  gAACAACGCCGCATAGCGGGACACGGCCTCTTGCGGGAACACCGCCAGGCCAAGGCCAAACTGCACCAGACGGTATGCCGCGTCGAACGACGACACGTACGAACGAAAACGCAGCGCCTTGCCTGACTGTTCGGCGCTTTTACGCATg                                                    >  5:88743/1‑148 (MQ=255)
                                                                                                   aacaacGCCGCATAGCGGGACACGGCCTCTTGCGGGAACACCGCCAGGCCAAGGCCAAACTGCACCAGACGGTATGCCGCGTCGAACGACGACACGTACGAACGAAAACGCAGCGCCTTGCCTGACTGTTCGGCGCTTTTACGCATGaa                                                  <  6:88743/149‑1 (MQ=255)
                                                                                                   aacaacGCCGCATAGCGGGACACGGCCTCTTGCGGGAACACCGCCAGGCCAAGGCCAAACTGCACCAGACGGTATGCCGCGTCGAACGACGACACGTACGAACGAAAACGCAGCGCCTTGCCTGACTGTTCGGCGCTTTTACGCATGa                                                   >  1:7318/1‑148 (MQ=255)
                                                                                                       acGCCGCATAGCGGGACACGGCCTCTTGCGGGAACACCGCCAGGCCAAGGCCAAACTGCACCAGACGGTATGCCGCGTCGAACGACGACACGTACGAACGAAAACGCAGCGCCTTGCCTGACTGTTCGGCGCTTTTACGCATGAAAGCg                                              <  8:367980/149‑1 (MQ=255)
                                                                                                        cgtcgcATAGGGGGACACGTCCTCTTGCGGTAACACCGCCAGGCCATGTGCTAACTGGACCAGACGGTATGCCGCGTCGAACGACGTCACGTACGAACTAAAACGCAGCGCGTTGCCTTACTGTTCGGCGCTTTTACTCATGAAAGCGt                                             <  8:12406/146‑1 (MQ=255)
                                                                                                        cgccgcATAGCGGGACACGGCCTCTTGCCGGAACACCGCCAGGCCAAGGCCAAACTGCACCAGACGGTATGCCGCGTCGAACGACGACACGTACGAACGAAAACGCAGCGCCTTGCCTGACTGTTCGGCGCTTTTACGCATGAAAGCGt                                             >  1:230747/1‑149 (MQ=255)
                                                                                                           cgcATAGCGGGACACGGCCTCTTGCGGGAACACCGCCAGGCCAAGGCCAAACTGCACCAGACGGTATGCCGCGTCGAACGACGACACGTACGAACGAAAACGCAGCGCCTTGCCTGACTGTTCGGCGCTTTTACGCATGAAAGCGTTc                                           >  3:15163/1‑148 (MQ=255)
                                                                                                                 gCGGGACACGGCCTCTTGCGGGAACACCGCCAGGCCAAGGCCAAACTGCACCAGACGGTATGCCGCGTCGAACGACGACACGTACGAACGAAAACGCAGCGCCTTGCCTGACTGTTCGGCGCTTTTACGCATGAAAGCGTTCAAtgtg                                     <  3:23549/148‑1 (MQ=255)
                                                                                                                 gCGGGACACGGCCTCTTGCGGGAACACCGCCAGGCCAAGGCCAAACTGCACCAGACGGTATGCCGCGTCGAACGACGACACGTACGAACGAAAACGCAGCGCCTTGCCTGACTGTTCGGCGCTTTTACGCATGAAAGCGTTCAAtgtg                                     <  6:389865/148‑1 (MQ=255)
                                                                                                                   gggACACGGCCTCTTGCGGGAACACCGCCAGGCCAAGGCCAAACTGCACCAGACGGTATGCCGCGTCGAACGACGACACGTACGAACGAAAACGCAGCGCCTTGCCTGACTGTTCGGCGCTTTTACGCATGAAAGCGTTCAAtg                                       <  1:170690/144‑1 (MQ=255)
                                                                                                                   gggACACGGCCTCTTGCGGGAACACCGCCAGGCCAAGGCCAAACTGCACCAGACGGTATGCCGCGTCGAACGACGACACGTACGAACGAAAACGCAGCGCCTTGCCTGACTGTTCGGCGCTTTTACGCATGAAAGCGTTCAAtg                                       >  2:170690/1‑144 (MQ=255)
                                                                                                                   gggACACGGCCTCTTGCGGGAACACCGCCAGGCCAAGGCCAAACTGCACCAGACGGTATGCCGCGTCGAACGACGACACGTACGAACGAAAACGCAGCGCCTTGCCTGACTGTTCGGCGCTTTTACGCATGAAAGCGTTCAATGTGGc                                   >  8:346281/1‑148 (MQ=255)
                                                                                                                           gCCTCTTGCGGGAACACCGCCAGGCCAAGGCCAAACTGCACCAGACGGTATGCCGCGTCGAACGACGACACGTACGAACGAAAACGCAGCGCCTTGCCTGACTGTTCGGCGCTTTTACGCATGAAAGCGTTCAATGTGGCAAAGGCCGa                          >  4:104850/1‑149 (MQ=255)
                                                                                                                             ctctTGCGGGAACACCGCCAGGCCAAGGCCAAACTGCACCAGACGGTATGCCGCGTCGAACGACGACACGTACGAACGAAAACGCAGCGCCTTGCCTGACTGTTCGGCGCTTTTACGCATGAAAGCGTTCAATGTGGCAAAGGCCGATa                        <  3:3799/149‑1 (MQ=255)
                                                                                                                              gctTGCGGGAACACCGCCAGGCCAAGGCCAAACTGCCCCAGACGGTATGCCGCGTCGAACGACGACACGTACGAACGAAACCGCAGCGCCTTGCCTGACTGTTCGGCGCTTTTACGCATGAAAGCGTTCAATGTGGCAAAGTCCGATAg                       <  2:128602/148‑1 (MQ=255)
                                                                                                                                   cGGGAACACCGCCAGGCCAAGGCCAAACTGCACCAGACGGTATGCCGCGTCGAACGACGACACGTACGAACGAAAACGCAGCGCCTTGCCTGACTGTTCGGCGCTTTTACGCATGAAAGCGTTCAATGTGGCAAAGGCCGATAGcccc                   <  8:139065/148‑1 (MQ=255)
                                                                                                                                    gggAACACCGCCAGGCCAAGGCCAAACTGCACCAGACGGTATGCCGCGTCGAACGACGACACGTACGAACGAAAACGCAGCGCCTTGCCTGACTGTTCGGCGCTTTTACGCATGAAAGCGTTCAATGTGGCAAAGGCCGATAGCCCCAg                 <  1:170963/149‑1 (MQ=255)
                                                                                                                                    gggAACACCGCCAGGCCAAGGCCAAACTGCACCAGACGGTATGCCGCGTCGAACGACGACACGTACGAACGAAAACGCAGCGCCTTGCCTGACTGTTCGGCGCTTTTACGCATGAAAGCGTTCAATGTGGCAAAGGCCGATAGCCCCAg                 <  2:158296/149‑1 (MQ=255)
                                                                                                                                          acCGCCAGGCCAAGGCCAAACTGCACCAGACGGTATGCCGCGTCGAACGACGACACGTACGAACGAAAACGCAGCGCCTTGCCTGACTGTTCGGCGCTTTTACGCATGAAAGCGTTCAATGTGGCAAAGGCCGATAGCCCCAGGTGGTc           <  2:7318/149‑1 (MQ=255)
                                                                                                                                            cGCCAGGCCAAGGCCAAACTGCACCAGACGGTATGCCGCGTCGAACGACGCCACGTACGAACTAAAACGCAGCGCCTTGCCTGACTCTTCGGCGCTTTTACGCATGAAAGCGTTCAATGAGGCAAAGGCCGATAGCCCCAGGTGGTCg          >  3:295874/1‑148 (MQ=255)
                                                                                                                                            cGCCAGGCCAAGGCCAAACTGCACCAGACGGTATGCCGCGTCGAACGACGACACGTACGAACGAAAACGCAGCGCCTTGCCTGACTGTTCGGCGCTTTTACGCATGAAAGCGTTCAATGTGGCAAAGGCCGATAGCCCCAGGTGGTCg          >  8:14103/1‑148 (MQ=255)
                                                                                                                                                   ccaaggccaaACTGCACCAGACGGTATGCCGCGTCGAACGACGACAGCGTACGAACGAAAACGCAGCGCCTTGCCTGACTGTTCGGCGCTTTTACGCATGAAAGCGTTCAATGTGGCAAAGGCCGATAGCCCCAGGTGGTCGTAAcc     >  1:89657/1‑147 (MQ=255)
                                                                                                                                                   ccaaggccaaACTGCACCAGACGGTATGCCGCGTCGAACGACGACACGTACGAACGAAAACGCAGCGCCTTGCCTGACTGTTCGGCGCTTTTACGCATGAAAGCGTTCAATGTGGCAAAGGCCGATAGCCCCAGGTGGTCGTAACCCAg  >  1:172670/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
AGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCGTCAGTCAGCGCAATGACCCGCAGGTCGAACAACGCCGCATAGCGGGACACGGCCTCTTGCGGGAACACCGCCAGGCCAAGGCCAAACTGCACCAGACGGTATGCCGCGTCGAACGACGACACGTACGAACGAAAACGCAGCGCCTTGCCTGACTGTTCGGCGCTTTTACGCATGAAAGCGTTCAATGTGGCAAAGGCCGATAGCCCCAGGTGGTCGTAACCCAG  >  NC_002947/3845780‑3846075

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: