Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 353,133 C→T 8.3% intergenic (+97/+21) hisF → / ← cbcV imidazole glycerol phosphate synthase subunit HisF/choline/betaine/carnitine ABC transporter ATP binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947353,1330CT 8.3% 204.4 / 6.9 72intergenic (+97/+21)hisF/cbcVimidazole glycerol phosphate synthase subunit HisF/choline/betaine/carnitine ABC transporter ATP binding protein
Reads supporting (aligned to +/- strand):  ref base C (38/28);  new base T (4/2);  total (42/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 = 9.90e-01

AAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCGGTTTGATAACGAATCTGCAGTGCCTCACGCAT  >  NC_002947/352998‑353282
                                                                                                                                       |                                                                                                                                                     
aaGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCACTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCTGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCa                                                                                                                                         >  2:2746608/1‑150 (MQ=255)
                   gcgGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTc                                                                                                                      >  1:2911489/1‑150 (MQ=255)
                                              aCCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCg                                                                                           <  2:2911489/150‑1 (MQ=255)
                                                                        GGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGcc                                                                                                                 <  1:368301‑M1/102‑1 (MQ=255)
                                                                        GGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGcc                                                                                                                 >  2:368301‑M1/1‑102 (MQ=255)
                                                                          GCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGt                                                               <  2:621696‑M1/150‑1 (MQ=255)
                                                                            acgcccgctcccacagagatctcagcaaacctgtgggagcgggcgtgcccgcgaagggcCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCg                                                             >  2:2712317‑M2/60‑150 (MQ=255)
                                                                            ACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCg                                                             <  1:3016944‑M1/150‑1 (MQ=255)
                                                                            ACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCg                                                             <  1:1470726‑M1/150‑1 (MQ=255)
                                                                              GCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGagcag                                                                                                           <  2:1720754‑M1/102‑1 (MQ=255)
                                                                              GCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGagcag                                                                                                           >  1:1720754‑M1/1‑102 (MQ=255)
                                                                               CCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTa                                                          <  1:2045078‑M1/150‑1 (MQ=255)
                                                                                ccgctcccacagagatctcagcaaacctgtgggagcgggcgtgcccgcgaagggcCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTAt                                                         <  1:2712317‑M2/95‑1 (MQ=255)
                                                                                CCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTaaa                                                                                                <  1:1153639‑M1/111‑1 (MQ=255)
                                                                                CCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTaaa                                                                                                >  2:1153639‑M1/1‑111 (MQ=255)
                                                                                CCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCCCGTGGTAAAGCTCGGTATCCCCGAGTATCCCTACCACCCTGTCGTTAt                                                         <  1:2728605‑M1/150‑1 (MQ=255)
                                                                                 cgctcccacagagatctcagcaaacctgtgggagcgggcgtgcccgcgaagggcTGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCg                                                                                           >  2:1937271‑M2/55‑115 (MQ=255)
                                                                                 cgctcccacagagatctcagcaaacctgtgggagcgggcgtgcccgcgaagggcTGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCg                                                                                           <  1:1937271‑M2/61‑1 (MQ=255)
                                                                                  gctcccacagagatctcagcaaacctgtgggagcgggcgtgcccgcgaagggcCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATcc                                                       >  2:2309356‑M2/54‑150 (MQ=255)
                                                                                   ctcccacagagatctcagcaaacctgtgggagcgggcgtgcccgcgaagggcTGCAAAGCAGCCCCGATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCt                                                      >  1:2220896‑M2/53‑150 (MQ=255)
                                                                                    tcccacagagatctcagcaaacctgtgggagcgggcgtgcccgcgaagggcTGCGAAGCgg                                                                                                                                              <  2:2874279‑M2/10‑1 (MQ=255)
                                                                                    tcccacagagatctcagcaaacctgtgggagcgggcgtgcccgcgaagggcTGCGAAGCgg                                                                                                                                              >  1:2874279‑M2/52‑61 (MQ=255)
                                                                                    tcccacagagatctcagcaaacctgtgggagcgggcgtgcccgcgaagggcCGCGAAGCGGc                                                                                                                                             >  1:1559353‑M2/52‑62 (MQ=255)
                                                                                    tcccacagagatctcagcaaacctgtgggagcgggcgtgcccgcgaagggcCGCGAAGCGGc                                                                                                                                             <  2:1559353‑M2/11‑1 (MQ=255)
                                                                                    tcccacagagatctcagcaaacctgtgggagcgggcgtgcccgcgaagggcCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTg                                                     >  2:1767001‑M2/52‑150 (MQ=255)
                                                                                      ccacagagatctcagcaaacctgtgggagcgggcgtgcccgcgaagggcTGCAAAGCAGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTgca                                                   >  2:2859314‑M2/50‑150 (MQ=255)
                                                                                      ccacagagatctcagcaaacctgtgggagcgggcgtgcccgcgaagggcCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGcgcg                                                                                                       >  2:3043344‑M2/50‑98 (MQ=255)
                                                                                      ccacagagatctcagcaaacctgtgggagcgggcgtgcccgcgaagggcCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGcgcg                                                                                                       <  1:3043344‑M2/49‑1 (MQ=255)
                                                                                      ccacagagatctcagcaaacctgtgggagcgggcgtgcccgcgaagggcCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCg                                                                                           <  2:2978190‑M2/61‑1 (MQ=255)
                                                                                      ccacagagatctcagcaaacctgtgggagcgggcgtgcccgcgaagggcCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCg                                                                                           >  1:2978190‑M2/50‑110 (MQ=255)
                                                                                      CCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTg                                                                                                                          >  1:1960146‑M1/1‑79 (MQ=255)
                                                                                      CCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTg                                                                                                                          <  2:1960146‑M1/79‑1 (MQ=255)
                                                                                       cacagagatctcagcaaacctgtgggagcgggcgtgcccgcgaagggcCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTg                                                                                                     >  2:1648413‑M2/49‑99 (MQ=255)
                                                                                       cacagagatctcagcaaacctgtgggagcgggcgtgcccgcgaagggcCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTg                                                                                                     <  1:1648413‑M2/51‑1 (MQ=255)
                                                                                       CACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTgcag                                                  <  1:1834656‑M1/150‑1 (MQ=255)
                                                                                        acagagatctcagcaaacctgtgggagcgggagtgcccgcgaagggcCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTgcagc                                                 >  2:1894695‑M2/48‑150 (MQ=255)
                                                                                        ACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATcc                                                       >  1:1820181‑M1/1‑144 (MQ=255)
                                                                                        ACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATcc                                                       <  2:1820181‑M1/144‑1 (MQ=255)
                                                                                            gTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTa                                                          <  1:1223775/137‑1 (MQ=255)
                                                                                            gTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTa                                                          >  2:1223775/1‑137 (MQ=255)
                                                                                             tttGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGa                                                                                                               <  2:1367227/83‑1 (MQ=255)
                                                                                             tttGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGa                                                                                                               >  1:1367227/1‑83 (MQ=255)
                                                                                                 cTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCttg                                        <  1:621696/150‑1 (MQ=255)
                                                                                                  tGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCg                                                                                           <  2:1292291/98‑1 (MQ=255)
                                                                                                  tGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCg                                                                                           >  1:1292291/1‑98 (MQ=255)
                                                                                                  tGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCttgt                                       >  1:492318/1‑150 (MQ=255)
                                                                                                  tGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCttgt                                       >  2:192040/1‑150 (MQ=255)
                                                                                                   gagaTCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCtt                                         <  2:2315598/147‑1 (MQ=255)
                                                                                                   gagaTCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCtt                                         >  1:2315598/1‑147 (MQ=255)
                                                                                                   gagaTCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTACAGCACCAGCttgtt                                      >  2:699485/1‑150 (MQ=255)
                                                                                                   gagaTCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACAGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCttgtt                                      >  1:2283274/1‑150 (MQ=255)
                                                                                                     gaTCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGc                                    >  1:496605/1‑150 (MQ=255)
                                                                                                        ctctGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCgg                                 >  1:2075343/1‑150 (MQ=255)
                                                                                                         tctGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCGGt                                >  1:543951/1‑150 (MQ=255)
                                                                                                         tctGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCGGt                                >  1:2880040/1‑150 (MQ=255)
                                                                                                               ggagcgggcgtgcccgcgaagggcCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCGGTTTGATa                          >  1:1276030‑M2/25‑150 (MQ=255)
                                                                                                               ggagcgggcgtgcccgcgaagggcCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCGGTTTGATa                          >  1:604718‑M2/25‑150 (MQ=255)
                                                                                                               ggagcgggcgtgcccgcgaagggcCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCGGTTTGATa                          >  1:604530‑M2/25‑150 (MQ=255)
                                                                                                               ggagcgggcgtgcccgcgaagggcCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCGGTTTGATa                          >  2:2392814‑M2/25‑150 (MQ=255)
                                                                                                                 agcgggcgtgcccgcgaagggcCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCGGTTTGATAAc                        <  2:3016944‑M2/128‑1 (MQ=255)
                                                                                                                   cgggcgtgcccgcgaagggcCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCGGTTTGATAACGa                      >  2:1632287‑M2/21‑150 (MQ=255)
                                                                                                                   cgggcgtgcccgcgaagggcCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCGGTTTGATAACGa                      <  2:1834656‑M2/130‑1 (MQ=255)
                                                                                                                      gcgtgcccgcgaagggcCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTg                                                                                                     <  2:1345833‑M2/51‑1 (MQ=255)
                                                                                                                      gcgtgcccgcgaagggcCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTg                                                                                                     >  1:1345833‑M2/18‑68 (MQ=255)
                                                                                                                          gcccgcgaagggcCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCGGTTTGATAACGAATCTGCa               <  2:1470726‑M2/137‑1 (MQ=255)
                                                                                                                          gcccgcgaagggcCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCGGTTTGATAACGAATCTGCa               <  1:2309356‑M2/137‑1 (MQ=255)
                                                                                                                            ccgcgaagggcCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGGATCCCCACCACCCTGTc                                                              <  2:676919‑M2/90‑1 (MQ=255)
                                                                                                                            ccgcgaagggcCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCACGAGGAGCCCCACCACCCTGTc                                                              >  1:676919‑M2/12‑101 (MQ=255)
                                                                                                                             cgcgaagggcCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCGGTTTGATAACGAATCTGCAGTg            >  1:2155237‑M2/11‑150 (MQ=255)
                                                                                                                              gcgaagggcCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCGGTTTGATAACGAATCTGCAGTGc           >  2:1327602‑M2/10‑150 (MQ=255)
                                                                                                                               cgAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAACAGCGCGTGGTAAAGCTCGGTATCCCGGGGATTGGCCACCACCCTGTCGTTATCCTGAAGAATCACCTATTCCGTGTTTTGATTACAAATCTGCAGTGcc          >  1:533635/1‑150 (MQ=255)
                                                                                                                                     gcCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCGGTTTGATAACGAATCTGCAGTGCCTCACGc    >  1:1713401‑M2/3‑150 (MQ=255)
                                                                                                                                       cgcgAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGAGCTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCGGTTTGATAACGAATCTGCAGTGCCTCACGCAt  >  2:1200466/1‑150 (MQ=255)
                                                                                                                                       |                                                                                                                                                     
AAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCGGTTTGATAACGAATCTGCAGTGCCTCACGCAT  >  NC_002947/352998‑353282

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