Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,835,320 C→T 5.9% Q18* (CAG→TAG)  PP_5114 → adenine methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,835,3200CT 5.9% 192.6 / 2.5 68Q18* (CAG→TAG) PP_5114adenine methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (36/28);  new base T (2/2);  total (38/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.15e-01

TCGTCACCGTTGGCCCGAAAGTGCCACAAAAACCGCTGCCAGCACCCACTGACAAGCCCTCCGAGCAACCGCTCGGGGTACCGGAGCACTAATGCCAAGATCCACCCCTCCCGCCCGCCCGCAAGCGGGCCAAAGCAAGGGCCAGGGCCACCTGCGCATCATCGCCGGCGAATGGCGCAGCCGCCGCCTGGCGGTGCCGGAAGGCGAAGGCCTGCGGCCAACGCCTGACCGCGTGCGGGAAACCGTATTCAACTGGCTGGCCCCCTACATCGAGGGCGCCCGCGTGCTGGA  >  NC_002947/5835178‑5835468
                                                                                                                                              |                                                                                                                                                    
tCGTCACCGTTGGCCCGAAAGTGCCACAAAAACCGCTGCCAGCACCCACTGACAAGCCCTCCGAGCAACCGCTCGGGGTACCGGAGCACTAATGCCAAGATCCACCCCTCCCGCCCGCCCGCAAGCGGGCCAAAGCAAgggccagggcca                                                                                                                                               <  2:702395/150‑1 (MQ=255)
  gTCACCGTTGGCCCGAAAGTGCCACAAAAACCGCTGCCAGCACCCACTGACAAGCCCTCCGAGCAACCGCTCGGGGTACCGGAGCACTAATGCCAAGATCCACCCCTCCCGCCCGCCCGCAAGCGGGCCAAAGCAAGGGCTAGGGCCAcc                                                                                                                                             >  1:246570/1‑150 (MQ=255)
                gAAAGTGCCACAAAAACCGCTGCCAGCACCCACTGACAAGCCCTCCGAGCAACCGCTCGGGGTACCGGAGCACTAATGCCAAGATCCACCCCTCCCGCCCGCCCGCAAGCGGGCCAAAGCAAGGGCCAGGGCCACCTGCGCATCATCGcc                                                                                                                               >  1:818272/1‑150 (MQ=255)
                         acaAAAACCGCTGCCAGCACCCACTGACAAGCCCTCCGAGCAACCGCTCGGGGTACCGGAGCACTAATGCCAAGATCCACCCCTCCCGCACGACCACCGGCACCGACAAGCCTGGCACACGGCCACCGGCGCATCATCGCCGGCGAATgg                                                                                                                      >  2:655168/1‑150 (MQ=255)
                           aaaaaCCGCTGCCAGCACCCACTGACAAGCCCTCCGAGCAACCGCTCGGGGTACCGGAGCACTAATGCCAAGATCCACCCCTCCCGCCCGCCCGCAAGCGGGCCGAAGCAAGGGCCAGGGCCACCTGCGCATCATCGCCGGCGAATGgcg                                                                                                                    >  1:579928/1‑150 (MQ=255)
                             aaaCCGCTGCCAGCACCCACTGACAAGCCCTCCGAGCAACCGCTCGGGGTACCGGAGCACTAATGCCAAGATCCACCCCTCCCGCCCGCCCGCAAGCGGGCCAAAGCAAGGGCCAGGGCCACCTGCGCATCATCGCCGGCGAATGGCGCa                                                                                                                  >  1:1015054/1‑150 (MQ=255)
                                 cGCTGCCAGCACCCACTGACAAGCCCTCCGAGCAACCGCTCGGGGTACCGGAGCACTAATGCCAAGATCCACCCCTCCCGCCCGCCCGCAAGCGGGCCAAAGCAAGGGCCAGGGCCACCTGCGCATCATCGCCGGCGAAt                                                                                                                        >  2:625552/1‑140 (MQ=255)
                                 cGCTGCCAGCACCCACTGACAAGCCCTCCGAGCAACCGCTCGGGGTACCGGAGCACTAATGCCAAGATCCACCCCTCCCGCCCGCCCGCAAGCGGGCCAAAGCAAGGGCCAGGGCCACCTGCGCATCATCGCCGGCGAAt                                                                                                                        <  1:625552/140‑1 (MQ=255)
                                  gCTGCCAGCACCCACTGACAAGCCCTCCGAGCAACCGCTCGGGGTACCGGAGCACTAATGCCAAGATCCACCCCTCCCGCCCGCCCGCAAGCGGGCCGAAGCAAGGGCCAGGGCCACCTGCGCATCATCGCCGGCGAATGGCGCAgccgc                                                                                                             <  2:579928/150‑1 (MQ=255)
                                      ccAGCACCCACTGACAAGCCCTCCGAGCAACCGCTCGGGGTACCGGAGCACTAATGCCAAGATCCACCCCTCCCGCCCGCCCGCAAGCGGGCCAAAGCAAGGGCCAGGGCCACCTGCGCATCATCGCCGGCGAATGGCGCAgccgccgcc                                                                                                         >  1:692111/1‑150 (MQ=255)
                                      ccAGCACCCACTGACAAGCCCTCCGAGCAACCGCTCGGGGTACCGGAGCACTAATGCCAAGATCCACCCCTCCCGCCCGCCCGCAAGCGGGCCAAAGCAAGGGCCAGGGCCACCTGCGCATCATCGCCGGCGAATGGCGCAgccgccgcc                                                                                                         >  2:150110/1‑150 (MQ=255)
                                        aGCACCCACTGACAAGCCCTCCGAGCAACCGCTCGGGGTACCGGAGCACTAATGCCAAGATCCACCCCTCCCGCCCGCCCGCAAGCGGGCCAAAGCAAGGGCCAGGGCCACCTGCGCATCATCGCCGGCGAATGGCGCAGCCGCCGCCTg                                                                                                       >  1:193959/1‑150 (MQ=255)
                                         gCACCCACTGACAAGCCCTCCGAGCAACCGCTCGGGGTACCGGAGCACTAATGCCAAGATCCACCCCTCCCGCCCGCCCGCAAGCGGGCCAAAGCAAGGGCCAGGGCCACCTGCGCATCATCGCCGGCGAATGGCGCAGCCGCCGCCTgg                                                                                                      >  1:239060/1‑150 (MQ=255)
                                             ccACTGACAAGCCCTCCGAGCAACCGCTCGGGGTACCGGAGCACTAATGCCAAGATCCACCCCTCCCGCCCGCCCGCAAGCGGGCCAAAGCAAGGGCCAGGGCCACCTGCGCATCATCGCCGGCGAATGGCGCAGCCGCCGCCTGGCGGt                                                                                                  >  1:824812/1‑150 (MQ=255)
                                             ccACTGACAAGCCCTCCGAGCAACCGCTCGGGGTACCGGAGCACTAATGCCAAGATCCACCCCTCCCGCCCGCCCGCAAGCGGGCCAAAGCAAGGGCCAGGGCCACCTGCGCATCATCGCCGGCGAATGGCGCAGCCGCCGCCTGGCGGt                                                                                                  >  1:953256/1‑150 (MQ=255)
                                              cACTGACAAGCCCTCCGAGCAACCGCTCGGGGTACCGGAGCACTAATGCCAAGATCCACCCCTCCCGCCCGCCCGCAAGCGGGCCAAAGCAAGGGCTAGGGCCACCTGCGCATCATCGCCGGCGAATGGCGCAGCCGCCGCCTGGCGGTg                                                                                                 <  2:246570/150‑1 (MQ=255)
                                              cACTGACAAGCCCTCCGAGCAACCGCTCGGGGTACCGGAGCACTAATGCCAAGATCCACCCCTCCCGCCCGCCCGCAAGCGGGCCAAAGCAAGGGCTAGGGCCACCTGCGCATCATCGCCGGCGAATGGCGCAGCCGCCGCCTGGCGGTg                                                                                                 >  2:599513/1‑150 (MQ=255)
                                               attgcgaaGCCCTCTGAGCAACCGCGTGGGGTACCGGAGCACTGATGCCAAGATCGACCGCTCCCGCCCGGCCGCAGGCGAGCCAAAGCAAGGGCCAGGGCCACCTGCGCATCATCGCCGGCTAATGACAAAGCCGCCGCCTGGCGGTGc                                                                                                >  2:152025/7‑150 (MQ=255)
                                               aCTGACAAGCCCTCCGAGCAACCGCTCGGGGTACCGGAGCACTAATGCCAAGATCCACCCCTCCCGCCCGCCCGCAAGCGGGCCAAAGCAAGGGCCAGGGCCACCTGCGCATCATCGCCGGCGAATGGCGCAGCCGCCGCCTGGCGGTGc                                                                                                >  2:99966/1‑150 (MQ=255)
                                                 tGACAAGCCCTCCGAGCAACCGCTCGGGGTACCGGAGCACTAATGCCAAGATCCACCCCTCCCGCCCGCCCGCAAGCGGGCCAAAGCAAGGGCCAGGGCCACCTGCGCATCATCGCCGGCGAATGGCGCAGCCGCCGCCTGGCGGTGCCg                                                                                              >  1:38914/1‑150 (MQ=255)
                                                  gACAAGCCCTCCGAGCAACCGCTCGGGGTACCGGAGCACTAATGCCAAGATCCACCCCTCCCGCCCGCCCGCAAGCGGGCCAAAGCAAGGGCCAGGGCCACCTGCGCATCATCGCCGGCGAATGGCGCAGCCGCCGCCTGGCGGTGCCgg                                                                                             <  2:309896/150‑1 (MQ=255)
                                                    cAAGCCCTCCGAGCAACCGCTCGGGGTACCGGAGCACTAATGCCAAGATCCACCCCTCCCGCCCGCCCGCAAGCGGGCCAAAGCAAGGGCCAGGGCCACCTGCGCATCATCGCCGGCGAATGGCGCAGCCGCCGCCTGGCGGTGCCGgaa                                                                                           <  2:765552/150‑1 (MQ=255)
                                                    cAAGCCCTCCGAGCAACCGCTCGGGGTACCGGAGCACTAATGCCAAGATCCACCCCTCCCGCCCGCCCGCAAGCGGGCCAAAGCAAGGGCCAGGGCCACCTGCGCATCATCGCCGGCGAATGGCGCAGCCGCCGCCTGGCGGTGCCGgaa                                                                                           <  2:818272/150‑1 (MQ=255)
                                                     aaGCCCTCCGAGCAACCGCTCGGGGTACCGGAGCACTAATGCCAAGATCCACCCCTCCCGCCCGCCCGCAAGCGGGCCAAAGCAAGGGCCAGGGCCACCTGCGCATCATCGCCGGCGAATGGCGCAGCCGCCGCCTGGCGGTGCCGgaag                                                                                          >  1:1126940/1‑150 (MQ=255)
                                                     aaGCCCTCCGAGCAACCGCTCGGGGTACCGGAGCACTAATGCCAAGATCCACCCCTCCCGCCCGCCCGCAAGCGGGCCAAAGCAAGGGCCAGGGCCACCTGCGCATCATCGCCGGCGAATGGCGCAGCCGCCGCCTGGCGGTGCCGgaag                                                                                          <  1:129619/150‑1 (MQ=255)
                                                      aGCCCTCCGAGCAACCGCTCGGGGTACCGGAGCACTAATGCCAAGATCCACCCCTCCCGCCCGCCCGCAAGCGGGCCAAAGCAAGGGCCAGGGCCACCTGCGCATCATCGCCGGCGAATGGCGCAGCCGCCGCCTGGCGGTGCCGgaagg                                                                                         <  2:1285980/150‑1 (MQ=255)
                                                       gCCCTCCGAGCAACCGCTCGGGGTACCGGAGCACTAATGCCAAGATCCACCCCTCCCGCCCGCCCGCAAGCGGGCCAAAGCAAGGGCTAGGGCCACCTGCGCATCATCGCCGGCGAATGGCGCAGCCGCCGCCTGGCGGTGCCGgaaggc                                                                                        <  1:599513/150‑1 (MQ=255)
                                                       gCCCTCCGAGCAACCGCTCGGGGTACCGGAGCACTAATGCCAAGATCCACCCCTCCCGCCCGCCCGCAAGCGGGCCAAAGCAAGGGCCAGGGCCACCTGCGCATCATCGCCGGCGAATGGCGCAGCCGCCGCCTGGCGGTGCCGgaaggc                                                                                        >  2:928144/1‑150 (MQ=255)
                                                        cccTCCGAGCAACCGCTCGGGGTACCGGAGCACTAATGCCAAGATCCACCCCTCCCGCCCGCCCGCAAGCGGGCCAAAGCAAGGGCCAGGGCCACCTGCGCATCATCGCCGGCGAATGGCGCAGCCGCCGCCTGGCGGTGCCGgaaggcg                                                                                       >  1:1148768/1‑150 (MQ=255)
                                                          cTCCGAGCAACCGCTCGGGGTACCGGAGCACTAATGCCAAGATCCACCCCTCCCGCCCGCCCGCAAGCGGGCCAAAGCAAGGGCCAGGGCCACCTGCGCATCATCGCCGGCGAATGGCGCAGCCGCCGCCTGGCGGTGCCGgaaggcgaa                                                                                     <  1:99966/150‑1 (MQ=255)
                                                           tCCGAGCAACCGCTCGGGGTACCGGAGCACTAATGCCAAGATCCACCCCTCCCGCCCGCCCGCAAGCGGGCCAAAGCAAGGGCCAGGGCCACCTGCGCATCATCGCCGGCGAATGGCGCAGCCGCCGCCTGGCGGTGCCGgaaggcgaag                                                                                    <  2:239060/150‑1 (MQ=255)
                                                           tCCGAGCAACCGCTCGGGGTACCGGAGCACTAATGCCAAGATCCACCCCTCCCGCCCGCCCGCAAGCGGGCCAAAGCAAGGGCCAGGGCCACCTGCGCATCATCGCCGGCGAATGGCGCAGCCGCCGCCTGGCGGTGCCGgaaggcgaag                                                                                    <  2:692111/150‑1 (MQ=255)
                                                               agcggttgCTCGGAGTACCGGAGCACTAATGCCAAGATCCACCCCTCCCGCCCGCCCGCAAGCGGGCCAAAGCAAGGGCCAGGGCCACCTGCGCATCATCGCCGGCGAATGGCGCAGCCGCCGCCTGGCGGTGCCGGAAGGCGAAGGCCt                                                                                <  1:957291/143‑1 (MQ=255)
                                                                 cAACCGCTCGGGGTACCGGAGCACTAATGCCAAGATCCACCCCTCCCGCCCGCCCGCAAGCGGGCCAAAGCAAGGGCCAGGGCCACCTGCGCATCATCGCCGGCGAATGGCGCAGCCGCCGCCTGGCGGTGCCGGACGGCGAAGGCCTGc                                                                              >  2:1248295/1‑150 (MQ=255)
                                                                   gttgCTCGGGGTACCGGAGCACTAATGCCAAGATCCACCCCTGCCGCCCGCCCGCAAGCGGGCCAAAGCAAGGGCCAGGGCCACCTGCGCATCATCGCCGGCGAATGGCGCAGCCGCCGCCTGGCGGTGCCGGAATGCGAAGGCCTGCgg                                                                            >  2:490688/4‑150 (MQ=255)
                                                                     cGCTCGGGGTACCGGAGCACTAATGCCAAGATCCACCCCTCCCGCCCGCCCGCAAGCGGGCCAAAGCAAGGGCCAGGGCCACCTGCGCATCATCGCCGGCGAATGGCGCAGCCGCCGCCTGGCGGTGCCGGAAGGCGAAGGCCTGCGGcc                                                                          <  1:711708/150‑1 (MQ=255)
                                                                       cTCGGGGTACCGGAGCACTAATGCCAAGATCCACCCCTCCCGCCCGCCCGCAAGCGGGCCAAAGCAAGGGCCAGGGCCACCTGCGCATCATCGCCGGCGAATGGCGCAGCCGCCGCCTGGCGGTGCCGGAAGGCGAAGGCCTGCGGCCaa                                                                        <  2:801541/150‑1 (MQ=255)
                                                                            ggTACCGGAGCACTAATGCCAAGATCCACCCCTCCCGCTACCGCGCAAGCGGGCCAAAGCAAGGGCCAGGGCCACCTGCGCATCATCGCCGGCGAATGGCGCAGCCGCCGCCTGGCGGTGCCGGAAGGCGAAGGCCTGCGGCCAACGCCt                                                                   <  1:655168/150‑1 (MQ=255)
                                                                            ggTACCGGAGCACTAATGCCAAGATCCACCCCTCCCGCCCGCCCGCAAGCGGGCCAAAGCAAGGGCCAGGGCCACCTGCGCATCATCGCCGGCGAATGGCGCAGCCGCCGCCTGGCGGTGCCGGAAGGCGAAGGCCTGCGGCCAACGCCt                                                                   <  1:916841/150‑1 (MQ=255)
                                                                            ggTACCGGAGCACTAATGCCAAGATCCACCCCTCCCGCCCGCCCGCAAGCGGGCCAAAGCAAGGGCCAGGGCCACCTGCGCATCATCGCCGGCGAATGGCGCAGCCGCCGCCTGGCGGTGCCGGAAGGCGAAGGCCTGCGGCCAACGCCt                                                                   <  1:807790/150‑1 (MQ=255)
                                                                              tACCGGAGCACTAATGCCAAGATCCACCCCTCCCGCCCGCCCGCAAGCGGGCCAAAGCAAGGGCCAGGGCCACCTGCGCATCATCACCGGCGAATGGCGCAGCCGCCGCCTGGCGGTGCCGGAAGGCGAAGGCCTGCGGCCAACGCCTGa                                                                 >  2:123141/1‑150 (MQ=255)
                                                                                  ggAGCACTAATGCCAAGATCCACCCCTCCCGCCCGCCCGCAAGCGGGCCAAAGCAAGGGCCAGGGCCACCTGCGCATCATCGCCGGCGAATGGCGCAGCCGCCGCCTGGCGGTGCCGGAAGGCGAAGGCCTGCGGCCAACGCCTGACcgc                                                             <  2:1148768/150‑1 (MQ=255)
                                                                                        cTAATGCCAAGATCCACCCCTCCCACCCGCCCGCAAGCGGGCCAAAGCAAGGGCCAGGGCCACCTGCGCATCATCGCCGGCGAATGGCGCAGCCGCCGCCTGGCGGTGCCGGAAGGCGAAGGCCTGCGGCCAACGCCTGa                                                                 >  2:757040/1‑140 (MQ=255)
                                                                                        cTAATGCCAAGATCCACCCCTCCCACCCGCCCGCAAGCGGGCCAAAGCAAGGGCCAGGGCCACCTGCGCATCATCGCCGGCGAATGGCGCAGCCGCCGCCTGGCGGTGCCGGAAGGCGAAGGCCTGCGGCCAACGCCTGa                                                                 <  1:757040/140‑1 (MQ=255)
                                                                                           aTGCCAAGATCCACCCCTCCCGCCCGCCCGCAAGCGGGCCAAAGCAAGGGCCAGGGCCACCTGCGCATCATCGCCGGCGAATGGCGCAgccgccgc                                                                                                          <  2:406910/96‑1 (MQ=255)
                                                                                           aTGCCAAGATCCACCCCTCCCGCCCGCCCGCAAGCGGGCCAAAGCAAGGGCCAGGGCCACCTGCGCATCATCGCCGGCGAATGGCGCAgccgccgc                                                                                                          >  1:406910/1‑96 (MQ=255)
                                                                                           aTGCCAAGATCCACCCCTCCCGCCCGCCCGCAAGCGGGCCAAAGCAAGGGCCAGGGCCACCTGCGCATCATCGCCGGCGAATGGCGCAGCCGCCGCCTGGCGGTGCCGGAAGGCGAAGGCCTGCGGCCAACGCCTGACCGCGTGCGGGaa                                                    >  2:234438/1‑150 (MQ=255)
                                                                                             gCCAAGATCCACCCCTCCCGCCCGCCCNCAAGCGGGCCAANGCAAGGNCCAGGGCCACCTGCGCATCATCGCCGGCGAATGGCGCAGCCGCCGCCTGGCGGTGCCGGAAGGCGAAGGCCTGCGGCCAACGCCTGACCGCGTGCGGGAAAc                                                  >  1:702057/1‑150 (MQ=255)
                                                                                                   aTCCACCCCTCCCGCCCGCCCGCAAGCGGGCCAAAGCAAGGGCCAGGGCCACCTGCGCATCATCGCCGGCGAATGGCGCAGCCGCCGCCTGGCGGTGCCGGAAGGCGAAGGCCTGCGGCCAACGCCTGACCGCGTGCGGGAAACCGTAtt                                            >  1:894147/1‑150 (MQ=255)
                                                                                                       aCCCCTCCCGCCCGCCCGCAAGCGGGCCAAAGCAAGGGCCAGGGCCACCTGCGCATCATCGCCGGCGAATGGCGCAGCCGCCGCCTGGCGGTGCCGGAAGGCGAAGGCCTGCGGCCAACGCCTGACCGCGTGCGGGAAACCGTATTCAAc                                        >  1:671230/1‑150 (MQ=255)
                                                                                                         cccTCCCGCCCGCCCGCAAGCGGGCCAAAGCAAGGGCCAGGGCCACCTGCGCATCATCGCCGGCGAATGGCGCAGCCGCCGCCTGGCGGTGCCGGAAGGCGAAGGCCTGCGGCCAACGCCTGACCGCGTGCGGGAAACCGTATTCAActg                                      >  1:1129715/1‑150 (MQ=255)
                                                                                                           cTCCCGCCCGCCCGCAAGCGGGCCAAAGCAAGGGCCAGGGCCACCTGCGCATCATCGCCGGCGAATGGCGCAGCCGCCGCCTGGCGGTGCCGGAAGGCGAAGGCCTGCGGCCAACGCCt                                                                   <  2:906233/119‑1 (MQ=255)
                                                                                                           cTCCCGCCCGCCCGCAAGCGGGCCAAAGCAAGGGCCAGGGCCACCTGCGCATCATCGCCGGCGAATGGCGCAGCCGCCGCCTGGCGGTGCCGGAAGGCGAAGGCCTGCGGCCAACGCCt                                                                   >  1:906233/1‑119 (MQ=255)
                                                                                                              gctttggcccgcAAGCGGGCCAAAGCAAGGGCCAGGGCCACCTGCGCATCATCGCCGGCGAATGGCACAGCCGCCGCCTGGCGGTGCCGGAAGGCGAAGGCCTGCGGCCAACGCCTGACCGCGTGCGGGAAACCGTATTCAActggct                                   >  2:264927/7‑148 (MQ=255)
                                                                                                              gctttggcccgcAAGCGGGCCAAAGCAAGGGCCAGGGCCACCTGCGCATCATCGCCGGCGAATGGCACAGCCGCCGCCTGGCGGTGCCGGAAGGCGAAGGCCTGCGGCCAACGCCTGACCGCGTGCGGGAAACCGTATTCAActggct                                   <  1:264927/142‑1 (MQ=255)
                                                                                                                tttggcccgcTTGCGGGCCAAAGCAAGGGCCAGGGCCACCTGCGCATCATCGCCGGCGAATGGCGCAGCCGCCGCCTGGCGGTGCCGGAAGGCGAAGGCCTGCGGCCAACGCCTGACCGCGTGCGGGAAACCGTATTCAACTGGCTGGcc                               >  1:877385/5‑150 (MQ=255)
                                                                                                                       cgcAAGCGTGCCAAAGCAAGGGCCAGGGCCACCTGCGCATCATCGCCGGCGAATGGCGCAGCCGCCGCCTGGCGGTGCCGGAAGGCGAAGGCCTGCGGCCAACGCCTGACCGCGTGCGGGAAACCGTATTCAACTGGCTGGCCCCCTACa                        >  1:893156/1‑150 (MQ=255)
                                                                                                                        gcAAGCGGGCCAAAGCAAGGGCCAGGGCCACCTGCGCATCATCGCCGGCGAATGGCGCAGCCGCCGCCTGGCGGTGCCGGAAGGCGAAGGCCTGCGGCCAACGCCTGCCCGCGTGCGGGCACCCCTACTCACCTGGCTCGCCCCCCAccc                       >  2:710521/1‑148 (MQ=255)
                                                                                                                        gcAAGCGGGCCAAAGCAAGGGCCAGGGCCACCTGCGCATCATCGCCGGCGAATGGCGCAGCCGCCGCCTGGCGGTGCCGGAAGGCGAAGGCCTGCGGCCAACGCCTGACCGCGTGCGGGAAACCGTATTCAACtggctggc                                >  1:1147372/1‑141 (MQ=255)
                                                                                                                        gcAAGCGGGCCAAAGCAAGGGCCAGGGCCACCTGCGCATCATCGCCGGCGAATGGCGCAGCCGCCGCCTGGCGGTGCCGGAAGGCGAAGGCCTGCGGCCAACGCCTGACCGCGTGCGGGAAACCGTATTCAACtggctggc                                <  2:1147372/141‑1 (MQ=255)
                                                                                                                            gtgggCCAAAGCAAGGGCCAGGGCCTCCTGCGCATCATCGCCGGCGAATGGCGCAGCCGCCGCCTGGCGGTGCCGGAAGGCGAAGGCCTGCGGCCAACGCCTGACCGCGTGCGGGAAACCGTATTCAACTGGCTGGCCCCCTACATCGAg                   <  2:1155863/148‑1 (MQ=255)
                                                                                                                            gCGGGCCAAAGCAAGGGCCAGGGCCACCTGCGCATCATCGCCGGCGAATGGCGCAGCCGCCGCCTGGCGGTGCCGGAAGGCGAAGGCCTGCGGCCAACGCCTGACCGCGTGCGGGAAACCGTATTCAACTGGCTGGCCCCCTACATCGAg                   <  2:1126940/150‑1 (MQ=255)
                                                                                                                             cGGGCCAAAGCAAGGGCCAGAGCCACATGCGAAATAATGCCGGCGAATGTCGCAGCCGCCGCCTGGCGGTAGCGGAAGGCGAAGGCCTGCGGCCAACGCCTGACCGCGTGCGGGAAACCGTATTCAACTGGCTGGCCCCCTACATCGAgg                  <  1:1248295/150‑1 (MQ=255)
                                                                                                                                   aaaGCAAGGGCCAGGGCCACCTGCGCGAAGTCGCCGGCGAATGGCGCAGCCGCCGCCTGGCGGTGCCGGAAGGCGAAGGCCTGCGGCCAACGCCTGACCGCGTGCGGGAAACCGTATTCAACTGGCTGGCCCCCTACATCGAGGGCGccc            >  1:674939/1‑150 (MQ=255)
                                                                                                                                     aGCAAGGGCCAGGGCCACCTGCGCATCATCGCCGGCGAATAATGCAGCAGTAGGCAAACGTTGCTACGAGATGAAGGCCTGCGGCCAACGCCTGACCGCTTGCGGGAAACCGTATTCAACTGGCTGGCCCCCTACATCGAGGGCGCCcgc          <  1:152025/150‑1 (MQ=255)
                                                                                                                                       cAAGGGCCAGGGCCACCTGCGCATCATCGCCGGCGAATGGCGCAGCCGCCGCCTGGCGGTGCCGGAAGGCGAAGGCCTGCGGCCAACGCCTGACCGCGTGCGGGAAACCGTATTCAACTGGCTGGCCCCCTACATCGAg                   <  2:724161/139‑1 (MQ=255)
                                                                                                                                       cAAGGGCCAGGGCCACCTGCGCATCATCGCCGGCGAATGGCGCAGCCGCCGCCTGGCGGTGCCGGAAGGCGAAGGCCTGCGGCCAACGCCTGACCGCGTGCGGGAAACCGTATTCAACTGGCTGGCCCCCTACATCGAg                   >  1:724161/1‑139 (MQ=255)
                                                                                                                                        aaGGGCCAGGTCCACCTGCGCATCATCGCCGGCGAATGGCGCAGCCGCCGCCTGGCGGTGCCGGAAGGCGAAGGCCTGCGGCCAACGCCTGACCGCGTGCGGGAAACCGTa                                              <  2:140693/111‑1 (MQ=255)
                                                                                                                                        aaGGGCCAGGGCCACCTGCGCATCATCGCCGGCGAATGGCGCAGCCGCCGCCTGGCGGTGCCGGAAGGCGAAGGCCTGCGGCCAACGCCTGACCGCGTGCGGGAAACCGTa                                              >  1:140693/1‑111 (MQ=255)
                                                                                                                                         agggccagggccACCTGCGCATCATCGCCGGCGAATGGCGCAGCCTCCGCCGGGCGGTGCCGGAAGGCGAAGGCCTGCGGCCAACGCCTGACCGCGTGCg                                                        <  2:314964/100‑1 (MQ=255)
                                                                                                                                         agggccagggccACCTGCGCATCATCGCCGGCGAATGGCGCAGCCGCCGCCTGGCGGTGCCGGAAGGCGAAGGCCTGCGGCCAACGCCTGACCGCGTGCg                                                        >  1:314964/1‑100 (MQ=255)
                                                                                                                                           ggccagggccaCCTGCGCATCATCGCCGGCGAATGGCGCAGCCGCCGCCTGGCGGTGCCGGAAGGCGAAGGCCTGCGGCCAACGCCTGACCGCGTGCGGGAAACCGTATTCAACTGGCTGGCCCCCTACATCGAGGGCGCCCGCGTGCTg    <  2:877385/150‑1 (MQ=255)
                                                                                                                                             ccagggccaCCTGCGCATCATCGCCGGCGAATGGCGCAGCCGCCGCCTGGCGGTGCCGGAAGGCGAAGGCCTGCGGCCAACGCCTGACCGCGTGCGGGAAACCGTATTCAACTGGCTGGCCCCCTACATCGAGGGCGCCCGCGTGCTGGa  <  2:1015054/150‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
TCGTCACCGTTGGCCCGAAAGTGCCACAAAAACCGCTGCCAGCACCCACTGACAAGCCCTCCGAGCAACCGCTCGGGGTACCGGAGCACTAATGCCAAGATCCACCCCTCCCGCCCGCCCGCAAGCGGGCCAAAGCAAGGGCCAGGGCCACCTGCGCATCATCGCCGGCGAATGGCGCAGCCGCCGCCTGGCGGTGCCGGAAGGCGAAGGCCTGCGGCCAACGCCTGACCGCGTGCGGGAAACCGTATTCAACTGGCTGGCCCCCTACATCGAGGGCGCCCGCGTGCTGGA  >  NC_002947/5835178‑5835468

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