Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,047,209 A→G 7.6% intergenic (+22/‑241) PP_3568 → / → quiA hypothetical protein/quinate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,047,2090AG 7.6% 187.8 / 3.8 68intergenic (+22/‑241)PP_3568/quiAhypothetical protein/quinate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (29/33);  new base G (3/2);  total (33/35)
Fisher's exact test for biased strand distribution p-value = 6.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGCCGGGCGATTGTCTGGACATCCCCCCGCACGACCGCCACCGCGTGGAATGGACCGCGCCCGGCACGGCCACCATCTGGCTCGCGGTGTTCTATAGCAGCAGCACCCCATGGCCTGCCTGAAGGGGCCATGCCTGCCCGTGCACAAGGCGCAGTTGTCCCGTCGCCGTGCAAGGTTGGGCTTACGACCAAAGCATGAATACAAATCGAAGACGCGTCATAAACGCCTCGTTCAATGCGTTTTTCAGCACCTTATCTATACTTTTGGCCACATCTTCGTAAGTTCAGCCCT  >  NC_002947/4047065‑4047356
                                                                                                                                                |                                                                                                                                                   
gcgcCGGGCGATTGTCTGGACATCCCCCCGCACGACCGCCACCGCGTGGAATGGACCGCGCCCGGCACGGCCACCATCTGGCTCGCGGTGTTCTATAGCAGCAGCACCCCATGGCCTCCCTGAAGGGGCCATGCCTGCCCGTGCACaa                                                                                                                                                  <  6:257767/148‑1 (MQ=255)
gcgcCGGGCGATTGTCTGGACATACCACCGCACGACCCCCACCGCGTGGAATGGACCGCGCCCGGCACGGCCACCATCTGGCTCGCGGTGTTCTATAGCAGCAGCACCCCATGGCCAGCCTGAACGGGCCATGCCTCCCCGTGCACAAg                                                                                                                                                 >  4:43098/1‑149 (MQ=255)
  gcCGGGCGATTGTCTGGACATCCCCCCGCACGACCGCCACCGCGTGGAATGGACCGCGCCCGGCACGGCCACCATCTGGCTCGCGGTGTTCTATAGCAGCAGCACCCCATGGCCTGCCTGAAGGGGCCATGCCTGCCCGTGCACAAGgc                                                                                                                                               <  5:312575/149‑1 (MQ=255)
     gggCGATTGTCTGGACATCCCCCCGCACGACCGCCACCGCGTGGAATGGACCGCGCCCGGCACGGCCACCATCTGGCTCGCGGTGTTCTATAGCAGCAGCACCCCATGGCCTGCCTGAAGGGGCCATGCCTGCCCGTGCACAAGGCGCa                                                                                                                                            <  5:39332/149‑1 (MQ=255)
     gggCGATTGTCTGGACATCCCCCCGCACGACCGCCACCGCGTGGAATGGACCGCGCCCGGCACGGCCACCATCTGGCTCGCGGTGTTCTATAGCAGCAGCACCCCATGGCCTGCCTGAAGGGGCCATGCCTGCCCGTGCACAAGGCGCa                                                                                                                                            <  6:43498/149‑1 (MQ=255)
      ggCGATTGTCTGGACATCCCCCCGCACGACCGCCACCGCGTGGAATGGACCGCGCCCGGCACGGCCACCATCTGGCTCGCGGTGTTCTATAGCAGCAGCACCCCATGGCCTGCCTGAAGGGGCCATGCCTGCCCGTGCACAAGGCGCAg                                                                                                                                           <  1:214413/149‑1 (MQ=255)
                     aTCCCCCCGCACGACCGCCACCGCGTGGAATGGACCGCGCCCGGCACGGCCACCATCTGGCTCGCGGTGTTCTATAGCAGCAGCACCCCATGGCCTGCCTGAAGGGGCCATGCCTGCCCGTGCACAAGGCGCAGTTGTCCCGTCGCCGt                                                                                                                            >  3:184355/1‑149 (MQ=255)
                     aTCCCCCCGCACGACCGCCACCGCGTGGAATGGACCGCGCCCGGCACGGCCACCATCTGGCTCGCGGTGTTCTATAGCAGCAGCACCCCATGGCCTGCCTGAAGGGGCCATGCCTGCCCGTGCACAAGGCGCAGTTGTCCCGTCGCCGt                                                                                                                            >  1:94133/1‑149 (MQ=255)
                       ccccccGCACGACCGCCACCGCGTGGAATGGACCGCGCCCGGCACGGCCACCATCTGGCTCGCGGTGTTCTATAGCAGCAGCACCCCATGGCCTGCCTGAAGGGGCCATGCCTGCCCGTGCACAAGGCGCAGTTGTCCCGTCGCCGTg                                                                                                                           >  2:68172/1‑148 (MQ=255)
                        cccccGCACGACCGCCACCGCGTGGAATGGACCGCGCCCGGCACGGCCACCATCTGGCTCGCGGTGTTCTATAGCAGCAGCACCCCATGGCCTGCCTGAAGGGGCCATGCCTGCCCGTGCACAAGGCGCAGTTGTCCCGTCGCCGTGc                                                                                                                          >  4:175669/1‑148 (MQ=255)
                         ccccGCACGACCGCCACCGCGTGGAATGGACCGCGCCCGGCACGGCCACCATCTGGCTCGCGGTGTTCTATAGCAGCAGCACCCCATGGCCTGCCTGAAGGGGCCATGCCTGCCCGTGCGCAAGGCGCAGTTGTCCCGTCGCCGTGc                                                                                                                          >  3:143883/1‑147 (MQ=255)
                         ccccGCACGACCGCCACCGCGTGGAATGGACCGCGCCCGGCACGGCCACCATCTGGCTCGCGGTGTTCTATAGCAGCAGCACCCCATGGCCTGCCTGAAGGGGCCATGCCTGCCCGTGCGCAAGGCGCAGTTGTCCCGTCGCCGTGc                                                                                                                          >  3:138917/1‑147 (MQ=255)
                         ccccGCACGACCGCCACCGCGTGGAATGGACCGCGCCCGGCACGGCCACCATCTGGCTCGCGGTGTTCTATAGCAGCAGCACCCCATGGCCTGCCTGAAGGGGCCATGCCTGCCCGTGCACAAGGCGCAGTTGTCCCGTCGCCGTGCaa                                                                                                                        <  1:122954/149‑1 (MQ=255)
                         ccccGCACGACCGCCACCGCGTGGAATGGACCGCGCCCGGCACGGCCACCATCTGGCTCGCGGTGTTCTATAGCAGCAGCACCCCATGGCCTGCCTGAAGGGGCCATGCCTGCCCGTGCACAAGGCGCAGTTGTCCCGTCGCCGTGCaa                                                                                                                        <  2:318476/149‑1 (MQ=255)
                          cccGCACGACCGCCACCGCGTGGAATGGACCGCGCCCGGCACGGCCACCATCTGGCTCGCGGTGTTCTATAGCAGCAGCACCCCATGGCCTGCCTGAAGGGGCCATGCCTGCCCGTGCACAAGGCGCAGTTGTCCCGTCGCCGTGCAAg                                                                                                                       >  6:44880/1‑149 (MQ=255)
                              cACGACCGCCACCGCGTGGAATGGACCGCGCCCGGCACGGCCACCATCTGGCTCGCGGTGTTCTATAGCAGCAGCACCCCATGGCCTGCCTGAAGGGGCCATGCCTGCCCGTGCACAAGGCGCAGTTGTCCCGTCGCCGTGCAAGGTTg                                                                                                                   <  7:70071/149‑1 (MQ=255)
                               aCGACCGCCACCGCGTGGAATGGACCGCGCCCGGCACGGCCACCATCTGGCTCGCGGTGTTCTATAGCAGCAGCACCCCATGGCCTGCCTGAAGGGGCCATGCCTGCCCGTGCACAAGGCGCAGTTGTCCCGTCGCCGTGCAAGGTTgg                                                                                                                  >  1:66106/1‑149 (MQ=255)
                                cGACCGCCACCGCGTGGAATGGACCGCGCCCGGCACGGCCACCATCTGGCTCGCGGTGTTCTATAGCAGCAGCACCCCATGGCCTGCCTGAAGGGGCCATGCCTGCCCGTGCACAAGGCGCAGTTGTCCCGTCGCCGTGCAAGGTTggg                                                                                                                 <  3:214069/149‑1 (MQ=255)
                                     gCCACCGCGTGGAATGGACCGCGCCCGGCACGGCCACCATCTGGCTCGCGGTGTTCTATAGCAGCAGCACCCCATGGCCTGCCTGAAGGGGCCATGCCTGCCCTTGCACAAGGCGCAGTTGTCCCGTCGCCGTGCAAGGTTGGTCTTAc                                                                                                            <  2:250732/149‑1 (MQ=255)
                                        aCCGCGTGGAATGGACCGCGCCCGGCACGGCCACCATCTGGCTCGCGGTGTTCTATAGCAGCAGCACCCCATGGCCTGCCTGAAGGGGCCATGCCTGCCCGTGCACAAGGCGCAGTTGTCCCGTCGCCGTGCAAGGTTGGGCTTACGAc                                                                                                         <  3:175669/149‑1 (MQ=255)
                                         ccGCGTGGAATGGACCGCGCCCGGCACGGCCACCATCTGGCTCGCGGTGTTCTATAGCAGCAGCACCCCATGGCCTGCCTGAAGGGGCCATGCCTGCCCGTGCACAAGGCGCAGTTGTCCCGTCGCCGTGCAAGGTTGGGCTTACGAc                                                                                                         >  6:183159/1‑148 (MQ=255)
                                                gAATGGACCGCGCCCGGCACGGCCACCATCTGGCTCGCGGTGTTCTATAGCAGCAGCACCCCATGGCCTGCCTGAAGGGGCCATGCCTGCCCGTGCACAAGGCGCAGTTCTCCCGTCGCCGTGCAAGGTTGGGCTTACGACCAAAGAAt                                                                                                 >  4:189204/1‑149 (MQ=255)
                                                    ggACCGCGCCCGGCACGGCCACCATCTGGCTCGCGGCGTTCTATAGCAGCAGCACCCCATGGCCTGCCTGAAGGGGCCATGCCTGCCCGTGCACAAGGCGCAGTTGTCCCGGCGCCGTGCAAGGTTGGGCTTACGACCAAAGCATGAAt                                                                                             <  7:218654/149‑1 (MQ=255)
                                                     gACCGCGCCCGGCACGGCCACCATCTGGCTCGCGGTGTTCTATAGCAGCAGCACCCCATGGCCTGCCTGAAGGGGCCATGCCTGCCCGTGCACAAGGCGCAGTTGTCCCGTCGCCGTGCAAGGTTGGGCTTACGACCAAAGCATGAAt                                                                                             <  4:184355/148‑1 (MQ=255)
                                                      aCCGCGCCCGGCACGGCCACCATCTGGCTCGCGGTGTTCTATAGCAGCAGCACCCCATGGCCTGCCTGAAGGGGCCATGCCTGCCCGTGCGCAAGGCGCAGTTGTCCCGTCGCCGTGCAAGGTTGGGCTTACGACCAAAGCATGAATa                                                                                            >  3:75361/1‑148 (MQ=255)
                                                       ccGCGCCCGGCACGGCCACCATCTGGCTCGCGGTGTTCTATAGCAGCAGCACCCCATGGCCTGCCTGAAGGGGCCATGCCTGCCCGTGCACAAGGCGCAGTTGTCCCGTCGCCGTGCAAGGTTGGGCTTACGACCAAAGCATGAATACa                                                                                          <  7:270604/149‑1 (MQ=255)
                                                        cgcgCCCGGCACGGCCACCATCTGGCTCGCGGTGTTCTATAGCAGCAGCACCCCATGGCCTGCCTGAAGGGGCCATGCCTGCCCGTGCACAAGGCGCAGTTGTCCCGTCGCCGTGCAAGGTTGGGCTTACGACCAAAGCATGAATACaa                                                                                         <  4:87932/149‑1 (MQ=255)
                                                              cGGCACGGCCACCATCTGGCTCGCGGTGTTCTATAGCAGCAGCACCCCATGGCCTGCCTGAAGGGGCCATGCCTGCCCGTGCACAAGGCGCAGTTGTCCCGTCGCCGTGCAAGGTTGGGCTTACGACCAAAGCATGAATACAAATCg                                                                                     >  5:98091/1‑147 (MQ=255)
                                                              cGGCACGGCCACCATCTGGCTCGCGGTGTTCTATAGCAGCAGCACCCCATGGCCTGCCTGAAGGGGCCATGCCTGCCCGTGCACAAGGCGCAGTTGTCCCGTCGCCGTGCAAGGTTGGGCTTACGACCAAAGCATGAATACAAATCGaa                                                                                   <  6:13025/149‑1 (MQ=255)
                                                                 cACGGCCACCATCTGGCTCGCGGTGTTCTATAGCAGCAGCACCCCATGGCCTGCCTGAAGGGGCCATGCCTGCCCGTGCACAAGGCGCAGTTGTCCCGTCGCCGTGCAAGGTTggg                                                                                                                 >  4:62035/1‑116 (MQ=255)
                                                                 cACGGCCACCATCTGGCTCGCGGTGTTCTATAGCAGCAGCACCCCATGGCCTGCCTGAAGGGGCCATGCCTGCCCGTGCACAAGGCGCAGTTGTCCCGTCGCCGTGCAAGGTTggg                                                                                                                 <  3:62035/116‑1 (MQ=255)
                                                                      ccaccaTCTGGCTCGCGGTGTTCTATAGCAGCAGCACCCCATGGCCTGCCTGAAGGGGCCATGCCTGCCCGTGCACAAGGCGCAGTTGTCCCGTCGCCGTGCAAGGTTGGGCTTACGACCAAAGCATGAATACAAATCGAAGACGCGt                                                                            >  3:6360/1‑148 (MQ=255)
                                                                       cacccTCTGGCTCGCGGTGTTCTATAGCAGCAGCCCCCCATGGCCTGCCTGAAGGGGCCATGCCTGCCCGTGCACAAGGCGCAGTTGTCCCGTCGCCGTGCAAGGTTGGGCTTACGACCAAAGCATGAATACAAATCGAAGACGCGTCa                                                                          >  1:189153/1‑149 (MQ=255)
                                                                        accaTCTGGCTCGCGGTGTTCTATAGCAGCAGCACCCCATGGCCTGCCTGAAGGGGCCATGCCTGCCCGTGCACAAGGCGCAGTTGTCCCGTCGCCGTGCAAGGTTGGGCTTACGACCAAAGCATGAATACAAATCGAAGACGCGTCAt                                                                         <  5:183159/149‑1 (MQ=255)
                                                                         ccaTCTGGCTCGCGGTGTTCTATAGCAGCAGCACCCCATGGCCTGCCTGAAGGGGCCATGCCTGCCCGTGCACAAGGCGCAGTTGTCCCGTCGCCGTGCAAGGTTGGGCTTACGACCAAAGCATGAATACAAATCGaa                                                                                   >  7:297677/1‑138 (MQ=255)
                                                                         ccaTCTGGCTCGCGGTGTTCTATAGCAGCAGCACCCCATGGCCTGCCTGAAGGGGCCATGCCTGCCCGTGCACAAGGCGCAGTTGTCCCGTCGCCGTGCAAGGTTGGGCTTACGACCAAAGCATGAATACAAATCGaa                                                                                   <  8:297677/138‑1 (MQ=255)
                                                                            tCTGGCTCGCGGTGTTCTATAGCAGCAGCACCCCATGGCCTGCCTGAAGGGGCCATGCCTGCCCGTGCACAAGGCGCAGTTGTCCCGTCGCCGTGCAAGGTTGGGCTTACGACCAAAGCATGAAt                                                                                             >  8:30853/1‑125 (MQ=255)
                                                                            tCTGGCTCGCGGTGTTCTATAGCAGCAGCACCCCATGGCCTGCCTGAAGGGGCCATGCCTGCCCGTGCACAAGGCGCAGTTGTCCCGTCGCCGTGCAAGGTTGGGCTTACGACCAAAGCATGAAt                                                                                             <  7:30853/125‑1 (MQ=255)
                                                                                gCTCGCGGTGTTCTATAGCAGCAGCACCCCATGGCCTGCCTGAAGGGGCCATGCCTGCCCGTGCACAAGGCGCAGTTGTCCCGTCGCCGTGCAAGGTTGGGCTTACGACCAAAGCATGAATACAAATCGAAGACGCGTCATAAACGc                                                                   <  5:44880/147‑1 (MQ=255)
                                                                                    gcgGTGTTCTATAGCAGCAGCACCCCATGGCCTGCCTGAAGGGGCCATGCCTGCCCGTGCACAAGGCGCAGTTGTCCCGTCGCCGTGCAAGGTTGGGCTTACGACCAAAGCATGAATACAAATCGAAGACGCGTCATAAACGCCTCGtt                                                             <  3:58215/149‑1 (MQ=255)
                                                                                    gcgGTGTTCTATAGCAGCAGCACCCCATGGCCTGCCTGAAGGGGCCATGCCTGCCCGTGCACAAGGCGCAGTTGTCCCGTCGCCGTGCAAGGTTGGGCTTACGACCAAAGCATGAATACAAATCGAAGACGCGTCATAAACGCCTCGtt                                                             <  1:200900/149‑1 (MQ=255)
                                                                                            cTATAGCAGCAGCACCCCATGGCCTGCCTGAAGGGGCCATGCCTGCCCGTGCACAAGGCGCAGTTGTCCCGTCGCCGTGCAAGGTTGGGCTTACGACCAAAGCATGAATACAAATCGAAGACGCGTCATAAACGCCTCGTTCAATGCGt                                                     >  7:168430/1‑149 (MQ=255)
                                                                                             tataGCAGCAGCACCCCATGGCCTGCCTGAAGGGGCCATGCCTGCCCGTGCACAAGGCGCAGTTGTCCCGTCGCCGTGCAAGGTTGGGCTTACGACCAAAGCATGAATACAAATCGAAGACGCGTCATAAACGCCTCGTTCAATGCGtt                                                    >  8:51737/1‑149 (MQ=255)
                                                                                              ataGCAGCAGCACCCCATGGCCTGCCTGAAGGGGCCATGCCTGCCCGTGCTCAAGGCGCAGTTGTGCCGTCGCCGTGCAAGGTTGGGGTTACGGCCAAAGCAGGATTTCAAATCGATGACGCGTCATAAACGCGTCGTTGATTGCGttt                                                   >  2:185235/1‑149 (MQ=255)
                                                                                              ataGCAGCAGCACCCCATGGCCTGCCTGAAGGGGCCATGCCTGCCCGTGCACAAGGCGCAGTTGTCCCGTCGCCGTGCAAGGTTGGGCTTACGACCAAAGCATGAATACAAATCGAAGACGCGTCAt                                                                         >  1:146324/1‑127 (MQ=255)
                                                                                              ataGCAGCAGCACCCCATGGCCTGCCTGAAGGGGCCATGCCTGCCCGTGCACAAGGCGCAGTTGTCCCGTCGCCGTGCAAGGTTGGGCTTACGACCAAAGCATGAATACAAATCGAAGACGCGTCAt                                                                         <  2:146324/127‑1 (MQ=255)
                                                                                               taGCAGCAGCACCCCATGGCCTGCCTGAAGGGGCCATGCCTGCCCGTGCACAAGGCGCAGTTGTCCGGTCGGCGTGCAAGGTTGGGCTTACGACCAAAGCATGAATACAAATCGAAGACGCGTCATAAACGCCTCGTTCAATGCGttt                                                   <  3:43098/148‑1 (MQ=255)
                                                                                               taGCAGCAGCACCCCATGGCCTGCCTGAAGGGGCCATGCCTGCCCGTGCACAAGGCGCAGTTGTCCCGTCGCCGTGCAAGGTTGGGCTTACGACCAAAGCAGGAATACAAATCGAAGACGCGTCATAAACGCCTCGTTCAATGCGtttt                                                  >  3:285311/1‑149 (MQ=255)
                                                                                                 gcagcagcaCCCCATGGCCTGCCTGAAGGGGCCATGCCTGCCCGTGCACAAGGCGCAGTTGTCCCGTCGCCGTGCAAGGTTGGGCTTACGACCAAAGCATGAATACAAATCGAAGACGCGTCATAAACGCCTCGTTCAATGCGttttt                                                 <  1:68172/148‑1 (MQ=255)
                                                                                                 gcagcagcaCCCCATGGCCTGCCTGAAGGGGCCATGCCTGCCCGTGCACAAGGCGCAGTTGTCCCGTCGCCGTGCAAGGTTGGGCTTACGACCAAAGCATGAATACAAATCGAAGACGCGTCATAAACGCCTCGTTCAATGCGttttt                                                 >  3:121249/1‑148 (MQ=255)
                                                                                                  cagcagcaCCCCATGGCCTGCCTGAAGGGGCCATGCCTGCCCGTGCACAAGGCGCAGTTGTCCCGTCGCCGTGCAAGGTTGGGCTTACGACCAAAGCATGAATACAAATCGAAGACGCGTCATAAACGCCTCGTTCAATGCGTTTTTCa                                               <  3:189204/149‑1 (MQ=255)
                                                                                                     cagcaCCCCATGGCCTGCCTGAAGGGGCCATGCCTGCCCGTGCACAAGGCGCAGTTGTCCCGTCGCCGTGCAAGGTTGGGCTTACGACCAAAGCATGAATACAAATCGAAGACGCGTCATAAACGCCTCGTTCAATGCGTTTTTCAGCa                                            <  4:6360/149‑1 (MQ=255)
                                                                                                       gcaCCCCATGGCCTGCCTGAAGGGGCCATGCCTGCCCGTGCACAAGGCGCAGTTGTCCCGTCGCCGTGCAAGGTTGGGCTTACTACCAAAGCATGAATACAAATCGAAGACGCGTCATAAACGCCTCGTTCAATGCGTTTTTCAGCAcc                                          <  2:66106/149‑1 (MQ=255)
                                                                                                         aCCCCATGGCCTGCCTGAAGGGGCCATGCCTGCCCGTGCACAAGGCGCAGTTGTCCCGTCGCCGTGCAAGGTTGGGCTTACGACCAAAGCATGAATACAAATCGAAGACGCGTCATAAACGCCTCGTTCAATGCGTTTTTCAGCACCt                                         >  1:204035/1‑148 (MQ=255)
                                                                                                         aCCCCATGGCCTGCCTGAAGGGGCCATGCCTGCCCGTGCACAAGGCGCAGTTGTCCCGTCGCCGTGCAAGGTTGGGCTTACGACCAAAGCATGAATACAAATCGAAGACGCGTCATAAACGCCTCGTTCAATGCGTTTTTCAGCACCt                                         >  8:141273/1‑148 (MQ=255)
                                                                                                         aCCCCATGGCCTGCCTGAAGGGGCCATGCCTGCCCGTGCACAAGGCGCAGTTGTCCCGTCGCCGTGCAAGGTTGGGCTTACGACCAAAGCATGAATACAAATCGAAGACGCGTCATAAACGCCTCGTTCAATGCGGTTTTCAGCACCt                                         >  8:117251/1‑148 (MQ=255)
                                                                                                              aTGGCCTGCCTGAAGGGGCCATGCCTGCCCGTGCACAAGGCGCAGTTGTCCCGTCGCCGTGCAAGGTTGGGCTTACGACCAAAGCATGAATACAAATCGAAGACGCGTCATAAACGCCTCGTTCAATGCGTTTTTCAGCACCTTATCta                                   <  2:94133/149‑1 (MQ=255)
                                                                                                                  cctgcctgAAGGGGCCATGCCTGCCCGTGCACAAGGCGCAGTTGTCCCGTCGCCGTGCAAGGTTGGGCTTACGACCAAAGCATGAATACAAATCGAAGACGCGTCATAAACGCCTCGTTCAATGCGTTTTTCAGCACCTTATCTATACt                               <  2:189153/149‑1 (MQ=255)
                                                                                                                   ctgcctgAAGGGGCCATGCCTGCCCGTGCACAAGGCGCAGTTGTCCCGTCGCCGTGCAAGGTTGGGCTTACGACCAAAGCATGAATACAAATCGAAGACGCGTCATAAACGCCTCGTTCAATGCGTTTTTCAGCACCTTATCTATACtt                              >  4:241316/1‑149 (MQ=255)
                                                                                                                        tgAAGGGGCCATGCCTGCCCGTGCACAAGGCGCAGTTGTCCCGTCGCCGTGCAAGGTTGGGCTTACGACCAAAGCATGAATACAAATCGAAGACGCGTCATAAACGCCTCGTTCAATGCGTTTTTCAGCACCTTATCTATACTTTTGGc                         >  5:134285/1‑149 (MQ=255)
                                                                                                                            ggggCCATGCCTGCCCGTGCACAAGGCGCAGTTGTCCCGTCGCCGTGCAAGGTTGGGCTTACGACCAAAGCATGAATACAAATCGAAGACGCGTCATAAACGCCTCGTTCAATGCGTTTTTCAGCACCTTATCTATACTTTTGGCcac                      >  6:3374/1‑148 (MQ=255)
                                                                                                                              ggCCATGCCTGCCCGTGCACAAGGCGCAGTTGTCCCGTCGCCGTGCAAGGTTGGGCTTACGACCAAAGCATGAATACAAATCGAAGACGCGTCATAAACGCCTCGTTCAAt                                                         <  4:7738/111‑1 (MQ=255)
                                                                                                                              ggCCATGCCTGCCCGTGCACAAGGCGCAGTTGTCCCGTCGCCGTGCAAGGTTGGGCTTACGACCAAAGCATGAATACAAATCGAAGACGCGTCATAAACGCCTCGTTCAAt                                                         >  3:7738/1‑111 (MQ=255)
                                                                                                                                 cATGCCTGCCCGTGCGCAAGGCGCAGTTGTCCCGTCGCCGTGCAAGGTTGGGCTTACGACCAAAGCATGAATACAAATCGAAGACGCGTCATAAACGCCTCGTTCAATGCGTTTTTCAGCACCTTATCTATACTTTTGGCCACATCTTc                <  4:138917/149‑1 (MQ=255)
                                                                                                                                 cATGCCTGCCCGTGCGCAAGGCGCAGTTGTCCCGTCGCCGTGCAAGGTTGGGCTTACGACCAAAGCATGAATACAAATCGAAGACGCGTCATAAACGCCTCGTTCAATGCGTTTTTCAGCACCTTATCTATACTTTTGGCCACATCTTc                <  4:143883/149‑1 (MQ=255)
                                                                                                                                 cATGCCTGCCCGTGCACAAGGCGCAGTTGTCCCGTCGCCGTGCAAGGTTGGGCTTACGACCAAAGCATGAATACAAATCGAAGACGCGTCATAAACGCCTCGTTCAATGCGTTTTTCAGCACCTTAt                                      <  2:40703/127‑1 (MQ=255)
                                                                                                                                 cATGCCTGCCCGTGCACAAGGCGCAGTTGTCCCGTCGCCGTGCAAGGTTGGGCTTACGACCAAAGCATGAATACAAATCGAAGACGCGTCATAAACGCCTCGTTCAATGCGTTTTTCAGCACCTTAt                                      >  1:40703/1‑127 (MQ=255)
                                                                                                                                          ccGTGCACAAGGCGCAGTTGTCCCGTCGCCGTGCAAGGTTGGGCTTACGACCAAAGCATGAATACAAATCGAAGACGCGTCATAAACGCCTCGTTCAATGCGTTTTTCAGCACCTTATCTATACTTTTGGCCACATCTTCGTAAGtt         >  2:52619/1‑147 (MQ=255)
                                                                                                                                               cacaAGGCGCAGTTGTCCCGTCGCCGTGCAAGGTTGGGCTTACGACCAAAGCATGAATACAAATCGAAGACGCGTCATAAACGCCTCGTTCAATGCGTTTTTCAGCACCTTATCTATACTTTTGGCCACATCTTCGTAAGTTCAGCCCt  >  7:165625/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
GCGCCGGGCGATTGTCTGGACATCCCCCCGCACGACCGCCACCGCGTGGAATGGACCGCGCCCGGCACGGCCACCATCTGGCTCGCGGTGTTCTATAGCAGCAGCACCCCATGGCCTGCCTGAAGGGGCCATGCCTGCCCGTGCACAAGGCGCAGTTGTCCCGTCGCCGTGCAAGGTTGGGCTTACGACCAAAGCATGAATACAAATCGAAGACGCGTCATAAACGCCTCGTTCAATGCGTTTTTCAGCACCTTATCTATACTTTTGGCCACATCTTCGTAAGTTCAGCCCT  >  NC_002947/4047065‑4047356

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: