Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_pae_pgi 3,560,455 +G intergenic (‑2/+1) glpR ← / ← glpR pseudogene, DNA‑binding transcriptional repressor;regulator; Energy metabolism, carbon: Anaerobic respiration; repressor of the glp operon/pseudogene, DNA‑binding transcriptional repressor;regulator; Energy metabolism, carbon: Anaerobic respiration; repressor of the glp operon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_pgi3,560,4551.G96.5% 296.1 / 4.5 87intergenic (‑2/+1)glpR/glpRpseudogene, DNA‑binding transcriptional repressor;regulator; Energy metabolism, carbon: Anaerobic respiration; repressor of the glp operon/pseudogene, DNA‑binding transcriptional repressor;regulator; Energy metabolism, carbon: Anaerobic respiration; repressor of the glp operon
Reads supporting (aligned to +/- strand):  ref base . (1/2);  new base G (36/48);  total (37/50)
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.87e-01

TGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑C‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGT  >  NC_000913_3_pae_pgi/3560306‑3560594
                                                                                                                                                        |                                                                                                                                           
tgCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑c                                                                                                                                              >  1:532021/1‑150 (MQ=255)
tgCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑c                                                                                                                                              >  1:41354/1‑150 (MQ=255)
   cGACATCGGTCAACAGCGTCGAGCCAATGGGGCTTTGCTCCGCCCCTTTGCGGTCGATGCGCTCTTGTTGTTCGGCCTGGGGGGCCTTGCGATCGCGCCACGGCTTGTTCACCGAACTGGAAGGCAGCGCCGCACAGCCATGAT‑GG‑cgca                                                                                                                                           >  2:534182/1‑147 (MQ=255)
    gATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑cgcag                                                                                                                                          <  2:386172/150‑5 (MQ=255)
     atatCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑cgcagg                                                                                                                                         >  2:71890/1‑145 (MQ=255)
       atCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGAt                                                                                                                                       >  4:201319/1‑150 (MQ=255)
       atCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGAt                                                                                                                                       >  1:430639/1‑150 (MQ=255)
         cGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCa                                                                                                                                     <  1:71903/150‑1 (MQ=255)
          gATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAg                                                                                                                                    >  2:400648/1‑150 (MQ=255)
              aaCAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGttt                                                                                                                                <  2:244305/150‑1 (MQ=255)
              aaCAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGttt                                                                                                                                <  1:437795/150‑1 (MQ=255)
               aCAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGtttt                                                                                                                               <  2:121147/150‑1 (MQ=255)
                 cccctCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATCCGCTCTTTTTCTTCCGTCTGCGTGGCCTTGCGATCGTGCCACCGCGTGTTAACCGAACTGTAACGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGc                                                                                                                             <  1:308229/146‑1 (MQ=255)
                   cGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTc                                                                                                                           <  2:353930/150‑1 (MQ=255)
                        aGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GGTG‑CAGGATCAGGTTTTGCTCCGCCa                                                                                                                      >  2:14666/1‑150 (MQ=255)
                         gCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAg                                                                                                                     >  2:426825/1‑150 (MQ=255)
                            aTTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑cgcag                                                                                                                                          <  3:261798/126‑5 (MQ=255)
                            aTTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑cgcag                                                                                                                                          >  4:261798/1‑122 (MQ=255)
                             ttGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCa                                                                                                                 >  2:620998/1‑150 (MQ=255)
                                gggATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTg                                                                                                              >  1:393064/1‑150 (MQ=255)
                                  gATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCa                                                                                                                      <  2:96134/140‑1 (MQ=255)
                                  gATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCa                                                                                                                      >  1:96134/1‑140 (MQ=255)
                                  gATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAg                                                                                                            <  3:388679/150‑1 (MQ=255)
                                   aTTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCGCCGCCAGCTCATTGAgg                                                                                                           <  1:534182/150‑1 (MQ=255)
                                    tttGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGt                                                                                                          >  2:419875/1‑150 (MQ=255)
                                     ttGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTc                                                                                                         >  1:694370/1‑150 (MQ=255)
                                     ttGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTc                                                                                                         >  2:250696/1‑150 (MQ=255)
                                       gCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTcgc                                                                                                       >  1:940492/1‑150 (MQ=255)
                                       gCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTcgc                                                                                                       >  1:940487/1‑150 (MQ=255)
                                       gCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTcgc                                                                                                       >  3:13496/1‑150 (MQ=255)
                                           cGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTcgcg                                                                                                      <  1:777679/147‑1 (MQ=255)
                                           cGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTcgcg                                                                                                      >  2:777679/1‑147 (MQ=255)
                                            gCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGa                                                                                                                                        <  1:139819/112‑1 (MQ=255)
                                            gCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGa                                                                                                                                        >  2:139819/1‑112 (MQ=255)
                                            gCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGa                                                                                                  <  2:655867/150‑1 (MQ=255)
                                                cTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGGGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGTTGGG‑G‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAg                                                                                              <  2:143813/150‑1 (MQ=255)
                                                  ttGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCGCGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTc                                                                                            <  1:17831/150‑1 (MQ=255)
                                                     cGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGc                                                                                         <  1:678990/150‑1 (MQ=255)
                                                       ggCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCgg                                                                                       <  1:250696/150‑1 (MQ=255)
                                                        gCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCggg                                                                                      <  4:204631/150‑1 (MQ=255)
                                                         cGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGc                                                                                     <  2:940492/150‑1 (MQ=255)
                                                         cGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGc                                                                                     <  2:940487/150‑1 (MQ=255)
                                                         cGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGc                                                                                     <  1:426825/150‑1 (MQ=255)
                                                             gcgcTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGAc                                                                                 <  2:694370/150‑1 (MQ=255)
                                                               gcTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑cg                                                                                                                                             <  1:874656/88‑2 (MQ=255)
                                                               gcTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑cg                                                                                                                                             >  2:874656/1‑87 (MQ=255)
                                                                 tctTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCgcgg                                                                                                     >  1:813456/1‑126 (MQ=255)
                                                                 tctTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCgcgg                                                                                                     <  2:813456/126‑1 (MQ=255)
                                                                   tttttCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGaa                                                                           >  1:921250/1‑150 (MQ=255)
                                                                   tttttCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGaa                                                                           >  1:920956/1‑150 (MQ=255)
                                                                    ttttCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGaaa                                                                          <  2:502767/150‑1 (MQ=255)
                                                                        cttcGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGaa                                                                           >  1:634685/1‑145 (MQ=255)
                                                                        cttcGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGaa                                                                           <  2:634685/145‑1 (MQ=255)
                                                                           cGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGc                                                                                                                        <  4:444520/97‑1 (MQ=255)
                                                                           cGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGc                                                                                                                        >  3:444520/1‑97 (MQ=255)
                                                                           cGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTa                                                                   <  2:393064/150‑1 (MQ=255)
                                                                           cGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTa                                                                   <  4:520167/150‑1 (MQ=255)
                                                                             gTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTAcc                                                                 <  2:41354/150‑1 (MQ=255)
                                                                              tCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCa                                                                <  1:400648/150‑1 (MQ=255)
                                                                              tCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCa                                                                <  1:931582/150‑1 (MQ=255)
                                                                                tGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGc                                                              <  3:201319/150‑1 (MQ=255)
                                                                                tGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGc                                                              <  2:430639/150‑1 (MQ=255)
                                                                                    tGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCtt                                                          >  2:54197/1‑150 (MQ=255)
                                                                                    tGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCtt                                                          <  1:620998/150‑1 (MQ=255)
                                                                                     ggCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCa                                                                >  3:428215/1‑143 (MQ=255)
                                                                                     ggCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCa                                                                <  4:428215/143‑1 (MQ=255)
                                                                                      gCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCg                                                        >  3:148002/1‑150 (MQ=255)
                                                                                        cTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGAt                                                                                                                                       <  4:134046/69‑1 (MQ=255)
                                                                                        cTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGAt                                                                                                                                       >  3:134046/1‑69 (MQ=255)
                                                                                         ttGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTa                                                     <  1:843362/150‑1 (MQ=255)
                                                                                          tGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTAc                                                    >  2:344608/1‑150 (MQ=255)
                                                                                            cGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTg                                                  <  1:419875/150‑1 (MQ=255)
                                                                                             gATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGc                                                                                                                    <  2:321550/83‑1 (MQ=255)
                                                                                             gATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGc                                                                                                                    >  1:321550/1‑83 (MQ=255)
                                                                                             gATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGa                                                 <  4:148002/150‑1 (MQ=255)
                                                                                             gATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGa                                                 <  2:532021/150‑1 (MQ=255)
                                                                                                  tGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGAACCGCCATGAT‑GG‑CGCAGCATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATaa                                            >  4:596989/1‑150 (MQ=255)
                                                                                                   gacacaGCGTGTTAAACGAACTGGAAGGCAGAGCCGCACAGCAATGAT‑GG‑CGAAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTAGCGGCGAATAGTCTGCGGGCTGAAGGAGAAATGCTCTACTAGCTCTTAGGTACTGCCATaaa                                           >  2:1004314/3‑149 (MQ=255)
                                                                                                      aCGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGc                                                                                                                    >  3:221013/1‑74 (MQ=255)
                                                                                                      aCGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGc                                                                                                                    <  4:221013/74‑1 (MQ=255)
                                                                                                       cGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCctg                                       >  2:42028/1‑150 (MQ=255)
                                                                                                            tgtTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGtt                                  <  2:755502/150‑1 (MQ=255)
                                                                                                                 aCCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAAcc                             >  2:935802/1‑150 (MQ=255)
                                                                                                                  ccGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCa                            <  1:71890/150‑1 (MQ=255)
                                                                                                                   cGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAg                           >  2:379178/1‑150 (MQ=255)
                                                                                                                     aaCTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCg                                                        >  2:266134/1‑119 (MQ=255)
                                                                                                                     aaCTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCg                                                        <  1:266134/119‑1 (MQ=255)
                                                                                                                          gatagCAGCGCCGCACCGCCATGGA‑GC‑CGCAGTGTCGGATTTTGCTCAGCCAGCTCATTGGGGTCGCCGCGCCTAGTCTTATAGCTGACGGCGAACTGCTCTACCAGCTCTTCGGTACTTACCTAACCCTGCTATTTAACCCGTTCGata                    <  4:10024/146‑1 (MQ=255)
                                                                                                                          gAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGata                    >  1:416416/1‑150 (MQ=255)
                                                                                                                           aaGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGataa                   <  1:344608/150‑1 (MQ=255)
                                                                                                                             ggCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGat                     <  2:374633/146‑1 (MQ=255)
                                                                                                                             ggCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGat                     >  1:374633/1‑146 (MQ=255)
                                                                                                                                  cgccgcACCGCCATGAT‑GGTG‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGtt            <  1:14666/150‑1 (MQ=255)
                                                                                                                                  cgccgcACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGtt            <  1:54197/150‑1 (MQ=255)
                                                                                                                                      gcACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGa        <  3:333248/150‑1 (MQ=255)
                                                                                                                                            ccATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTtgt  >  2:352440/1‑150 (MQ=255)
                                                                                                                                                        |                                                                                                                                           
TGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑C‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGT  >  NC_000913_3_pae_pgi/3560306‑3560594

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: ‑