Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,591,760 C→G 5.0% P158P (CCC→CCG mccB → methylcrotonyl‑CoA carboxylase biotin‑containing subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,591,7600CG 5.0% 604.1 / 15.8 179P158P (CCC→CCGmccBmethylcrotonyl‑CoA carboxylase biotin‑containing subunit beta
Reads supporting (aligned to +/- strand):  ref base C (99/71);  new base G (4/5);  total (103/76)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 187 total aligned reads displayed.

GACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATCCCGCAGATCGCCGTGGTGATGGGCTCGTGCACCGCTGGCGGCGCGTACGTGCCGGCCATGGCCGACGAAGCGATCATGGTGCGCCAGC  >  NC_002947/4591611‑4591908
                                                                                                                                                     |                                                                                                                                                    
gacgCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTcccc                                                                                                                                                      <  1:3152132/150‑1 (MQ=255)
  cgCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGa                                                                                                                                                    <  1:213036/150‑1 (MQ=255)
   gCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGAc                                                                                                                                                   >  2:1036860/1‑150 (MQ=255)
        ggTGAAGGGCGGCACCTATTATCCGGTCACGGTGAAGAAGCACCTGCGTGCGCAGGAAATCGCGGGCGAGAACCGTCTGCCCTGCATTTATCTGGTCGATTCCGGCGGCGCGAACCTGCCCAATCAGGACGAGGTTTTTCCGGATCGCGa                                                                                                                                              <  1:3831603/150‑1 (MQ=255)
        ggTGAAGGGCGGCACCTATTATCCGGTCACGGTGAAGAAGCACCTGCGTGCGCAGGAAATCGCGGGCGAGAACCGTCTGCCCTGCATTTATCTGGTCGATTCCGGCGGCGCGAACCTGCCCAATCAGGACGAGGTTTTTCCGGATCGCGa                                                                                                                                              >  1:379051/1‑150 (MQ=255)
        ggTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGa                                                                                                                                              <  1:4051909/150‑1 (MQ=255)
               ggCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTc                                                                                                                                       <  1:3620683/150‑1 (MQ=255)
                 cGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCgg                                                                                                                                     <  2:2990054/150‑1 (MQ=255)
                 cGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCgg                                                                                                                                     <  2:40776/150‑1 (MQ=255)
                   gTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGgc                                                                                                                                   >  1:872426/1‑150 (MQ=255)
                   gTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGgc                                                                                                                                   <  2:2926794/150‑1 (MQ=255)
                   gTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGgc                                                                                                                                   >  1:22833/1‑150 (MQ=255)
                   gTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGgc                                                                                                                                   <  2:527189/150‑1 (MQ=255)
                           tacCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTttc                                                                                                                           >  1:2353004/1‑150 (MQ=255)
                           tacCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTttc                                                                                                                           >  1:1383363/1‑150 (MQ=255)
                             cccGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTttctt                                                                                                                         <  1:3489728/150‑1 (MQ=255)
                             cccGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTttctt                                                                                                                         <  1:917125/150‑1 (MQ=255)
                             cccGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTttctt                                                                                                                         <  2:3369002/150‑1 (MQ=255)
                             cccGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTttctt                                                                                                                         >  2:553543/1‑150 (MQ=255)
                              ccGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTttcttc                                                                                                                        >  1:2119786/1‑150 (MQ=255)
                                ggtcaCGGTGAAGAAGCACCTGCGTGCGCAGGAAATCGCGGGCGCGAACCGTCTGCCCTGCATTTAGCTGGTCGATTCCGGCGGCGCGAACCTGCCCACTCAGGACGAGGGTTTTCCGGATCGCGATCCGTTCGGGCGCATCTTCTCCaa                                                                                                                      >  2:3065701/5‑150 (MQ=255)
                                    aCGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAg                                                                                                                  <  1:1951649/150‑1 (MQ=255)
                                      ggTGAAGAAGCACCTGCGTGCGCAGGAAATCGCGGGCGAGAACCGTCTGCCCTGCATTTATCTGGTCGATTCCGGCGGCGCGAACCTGCCCAATCAGGACGAGGTTTTTCCGGATCGCGATCATTTCGGGCGCATCTTCTACAATCAGGc                                                                                                                >  2:43774/1‑150 (MQ=255)
                                      ggTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCCTTTCGGGCGCATTTTCTTCAACCAGGc                                                                                                                >  1:966313/1‑150 (MQ=255)
                                        tGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCa                                                                                                              >  1:3243435/1‑150 (MQ=255)
                                        tGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCa                                                                                                              <  2:4947051/150‑1 (MQ=255)
                                          aagaagCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAAc                                                                                                            <  2:4629442/150‑1 (MQ=255)
                                            gaagCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCCTTTCGGGCGCATTTTCTTCAACCAGGCCAACAt                                                                                                          >  1:4724939/1‑150 (MQ=255)
                                            gaagCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCaa                                                                                                             <  2:1908466/147‑1 (MQ=255)
                                            gaagCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACAt                                                                                                          <  2:3173867/150‑1 (MQ=255)
                                            gaagCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACGTGCCCCGCCAGGCCGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACag                                                                                                          >  1:520636/1‑149 (MQ=255)
                                              agCACCTGCGTTCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAATTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGa                                                                                                        <  1:1115123/150‑1 (MQ=255)
                                              agCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGa                                                                                                        <  2:1337453/150‑1 (MQ=255)
                                                 aCCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAgcg                                                                                                     >  2:3129730/1‑150 (MQ=255)
                                                   cTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAgcgcg                                                                                                   >  2:1103547/1‑150 (MQ=255)
                                                      cGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCaa                                                                                                <  2:4714743/150‑1 (MQ=255)
                                                      cGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCaa                                                                                                <  2:3371247/150‑1 (MQ=255)
                                                      cGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCaa                                                                                                <  2:2590505/150‑1 (MQ=255)
                                                      cGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCaa                                                                                                <  2:2353004/150‑1 (MQ=255)
                                                         gCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCACGGc                                                                                             >  1:2987624/1‑150 (MQ=255)
                                                         gCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGc                                                                                             >  1:624934/1‑150 (MQ=255)
                                                           acaAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCAt                                                                                           >  2:866479/1‑150 (MQ=255)
                                                             aaaCCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATcc                                                                                         >  1:4458141/1‑150 (MQ=255)
                                                             aaaCCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATcc                                                                                         >  1:1093935/1‑150 (MQ=255)
                                                             aaaCCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATcc                                                                                         >  1:3902976/1‑150 (MQ=255)
                                                               aCCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAATGCATCCCg                                                                                       >  1:3650765/1‑150 (MQ=255)
                                                               aCCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATCCCg                                                                                       >  2:3760499/1‑150 (MQ=255)
                                                                 cATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATCCCGCa                                                                                     >  2:4393607/1‑150 (MQ=255)
                                                                 cATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATCCCGCa                                                                                     >  2:3688030/1‑150 (MQ=255)
                                                                 cATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATCCCGCa                                                                                     >  2:3688016/1‑150 (MQ=255)
                                                                   tCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATCCCGCAGa                                                                                   >  1:990400/1‑150 (MQ=255)
                                                                     gCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATCCCGCAGATc                                                                                 <  2:988662/150‑1 (MQ=255)
                                                                       ccTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAGGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGCGCGCGCAAGGCAGCCCGCAGATCGc                                                                               >  2:1698767/1‑150 (MQ=255)
                                                                         tGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCCTTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATCCCGCAGATCGCCg                                                                             >  1:4626983/1‑150 (MQ=255)
                                                                         tGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATCCCGCAGATCGCCg                                                                             >  2:2471328/1‑150 (MQ=255)
                                                                          gCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATCCCGCAGATCGCCgt                                                                            >  1:2751610/1‑150 (MQ=255)
                                                                            aGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATCCCGCAGATCGCCgtgg                                                                          <  1:1237762/150‑1 (MQ=255)
                                                                             gAACCGCCTGCCATGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATCCCGCAGATCGCCgtggt                                                                         <  1:4235937/150‑1 (MQ=255)
                                                                               aCCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATCCCGCAGATCGCCGTGGTGa                                                                       <  1:1822354/150‑1 (MQ=255)
                                                                               aCCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATCCCGCAGATCGCCGTGGTGa                                                                       <  2:3951413/150‑1 (MQ=255)
                                                                                ccGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATCCCGCAGATCGCCGTGGTGAt                                                                      <  2:3509226/150‑1 (MQ=255)
                                                                                 cGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATCCCGCAGATCGCCGTGGTGATg                                                                     >  2:2908615/1‑150 (MQ=255)
                                                                                  gCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATCCCGCAGATCGCCGTGGTGATgg                                                                    >  1:526516/1‑150 (MQ=255)
                                                                                        cGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATCCCGCAGATCGCCGTGGTGATGGGCTCGt                                                              <  1:131835/150‑1 (MQ=255)
                                                                                          tGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATCCCGCAGATCGCCGTGGTGATGGGCTCGTGc                                                            <  1:1815365/150‑1 (MQ=255)
                                                                                           gCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATCCCGCAGATCGCCGTGGTGATGGGCTCGTGCa                                                           <  2:5091205/150‑1 (MQ=255)
                                                                                             atctatctGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATCCCGCAGATCGCCGTGGTGATGGGCTCGTGCAcc                                                         >  1:3084462/1‑150 (MQ=255)
                                                                                              tctatctGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATCCCGCAGATCGCCGTGGTGATGGGCTCGTGCACCg                                                        >  1:4396449/1‑150 (MQ=255)
                                                                                              tcgctctGGTGGACTCCGGCGGCGCCACCTGGCCCCGCCTGGACGACGTGGTCCCCGACCGCGAGCATTTCGGGCGCATGTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCAGCCCGCAGATCGCCGTGGTGAGGGGCTCGTGCACCg                                                        <  2:3371626/146‑1 (MQ=255)
                                                                                               ctatctGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAATATGAGCGCGCAAGGCATCCCGCAGATCGCCGTGGTGATGGGCTCGTGCACCGc                                                       >  2:1307220/1‑150 (MQ=255)
                                                                                                 atctGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATCCCGCAGATCGCCGTGGTGATGGGCTCGTGCACCGCTg                                                     >  1:4274727/1‑150 (MQ=255)
                                                                                                  tctGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATCCCGCAGATCGCCGTGGTGATGGGCTCGTGCACCGCTgg                                                    <  1:2780863/150‑1 (MQ=255)
                                                                                                  tctGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATCCCGCAGATCGCCGTGGTGATGGGCTCGTGCACCGCTgg                                                    >  2:3393813/1‑150 (MQ=255)
                                                                                                  tctGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATCCCGCAGATCGCCGTGGTGATGGGCTCGTGCACCGCTgg                                                    >  2:4985172/1‑150 (MQ=255)
                                                                                                  tctGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATCCCGCAGATCGCCGTGGTGATGGGCTCGTGCACCGCTgg                                                    <  2:4982087/150‑1 (MQ=255)
                                                                                                    tggtggACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCCACATGAGCGCGCAAGGCATCCCGCAGATCGCCGTGGTGATGGGCTCGTGCACCGCTggcg                                                  >  2:2002923/1‑150 (MQ=255)
                                                                                                    tggtggACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATCCCGCAGATCGCCGTGGTGATGGGCTCGTGCACCGCTggcg                                                  >  2:1435996/1‑150 (MQ=255)
                                                                                                    tggtggACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATCCCGCAGATCGCCGTGGTGATGGGCTCGGGCACCGCTggcg                                                  >  2:2507661/1‑150 (MQ=255)
                                                                                                         gACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATCCCGCAGATCGCCGTGGTGATGGGCTCGTGCACCGCTGGCGGcgcg                                             >  1:1421724/1‑150 (MQ=255)
                                                                                                          aCTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATCCCGCAGATCGCCGTGGTGATGGGCTCGTGCACCGCTGGCGGCGCGt                                            <  2:3210448/150‑1 (MQ=255)
                                                                                                          aCTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATCCCGCAGATCGCCGTGGTGATGGGCTCGTGCACCGCTGGCGGCGCGt                                            <  2:4281709/150‑1 (MQ=255)
                                                                                                           cTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATCCCGCAGATCGCCGTGGTGATGGGCTCGTGCACCGCTGGCGGCGCGTa                                           >  1:593751/1‑150 (MQ=255)
                                                                                                            tCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATCCCGCAGATCGCCGTGGTGATGGGCTCGTGCACCGCTGGCGGCGCGTAc                                          <  1:3615729/150‑1 (MQ=255)
                                                                                                                  ggcgCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATCCCGCAGATCGCCGTGGTGATGGGCTCGTGCACCGCTGGCGGCGCGTACGTGCCg                                    <  2:680551/150‑1 (MQ=255)
                                                                                                                        aaCTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATCCCGCAGATCGCCGTGGTGATGGGCTCGTGCACCGCTGGCGGCGCGTACGTGCCGGCCATg                              >  1:723310/1‑150 (MQ=255)
                                                                                                                        aaCTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATCCCGCAGATCGCCGTGGTGATGGGCTCGTGCACCGCTGGCGGCGCGTACGTGCCGGCCATg                              >  1:4195459/1‑150 (MQ=255)
                                                                                                                               cccGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTttc                                                                                                                           <  2:1684837/50‑1 (MQ=255)
                                                                                                                               cccGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTttc                                                                                                                           >  1:1684837/1‑50 (MQ=255)
                                                                                                                               cccGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATCCCGCAGATCGCCGTGGTGATGGGCTCGTGCACCGCTGGCGGCGCGTACGTGCCGGCCATGGCcgacg                       >  2:4848208/1‑150 (MQ=255)
                                                                                                                                ccGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATCCCGCAGATCGCCGTGGTGATGGGCTCGTGCACCGCTGGCGGCGCGTACGTGCCGGCCATGGCcgacga                      >  2:4246305/1‑150 (MQ=255)
                                                                                                                                ccGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATCCCGCAGATCGCCGTGGTGATGGGCTCGTGCACCGCTGGCGGCGCGTACGTGCCGGCCATGGCcgacga                      >  2:1624652/1‑150 (MQ=255)
                                                                                                                                      ggACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATCCCGCAGATCGCCGTGGTGATGGGCTCGTGCACCGCTGGCGGCGCGTACGTGCCGGCCATGGCCGACGAAGCGAt                >  2:30973/1‑150 (MQ=255)
                                                                                                                                      ggACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATCCCGCAGATCGCCGTGGTGATGGGCTCGTGCACCGCTGGCGGCGCGTACGTGCCGGCCATGGCCGACGAAGCGAt                >  1:3975408/1‑150 (MQ=255)
                                                                                                                                       gACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATCCCGCAGATCGCCGTGGTGATGGGCTCGTGCACCGCTGGCGGCGCGTACGTGCCGGCCATGGCCGACGAAGCGATc               >  2:3779292/1‑150 (MQ=255)
                                                                                                                                         cGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATCCCGCAGATCGCCGTGGTGATGGGCTCGTGCACCGCTGGCGGCGCGTACGTGCCGGCCATGGCCGACGAAGCGATCAt             >  2:4524311/1‑150 (MQ=255)
                                                                                                                                             gtgtTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATCCCGCAGATCGCCGTGGTGATGGGCTCGTGCACCGCTGGCGGCGCGTACGTGCCGGCCATGGCCGACGAAGCGATCATGGTg         >  2:3741490/1‑150 (MQ=255)
                                                                                                                                                  ccccGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATCCCGCAGATCGCCGTGGTGATGGGCTCGt                                                              <  2:3147674/92‑1 (MQ=255)
                                                                                                                                                  ccccGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATCCCGCAGATCGCCGTGGTGATGGGCTCGTGCACCGCTGGCGGCGCGTACGTGCCGGCCATGGCCGACGAAGCGATCATGGTGCGCca    >  2:3661822/1‑150 (MQ=255)
                                                                                                                                                   cccGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATCCCGCAGATCGCCGTGGTGATGGGCTCGTGCACCGCTGGCGGCGCGTACGTGCCGGCCATGGCCGACGAAGCGATCATGGTGCGCcag   >  1:480885/1‑150 (MQ=255)
                                                                                                                                                   cccGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATCCCGCAGATCGCCGTGGTGATGGGCTCGTGCACCGCTGGCGGCGCGTACGTGCCGGCCATGGCCGACGAAGCGATCATGGTGCGCcag   >  1:2786559/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATCCCGCAGATCGCCGTGGTGATGGGCTCGTGCACCGCTGGCGGCGCGTACGTGCCGGCCATGGCCGACGAAGCGATCATGGTGCGCCAGC  >  NC_002947/4591611‑4591908

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: