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.4% 266.2 / 2.0 85intergenic (‑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 . (2/1);  new base G (37/44);  total (39/46)
Fisher's exact test for biased strand distribution p-value = 5.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

GCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑C‑CA‑GGAT‑CAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTT  >  NC_000913_3_pae_pgi/3560307‑3560597
                                                                                                                                                      |                                                                                                                                                
gCCGATATCGATAAACCGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑c                                                                                                                                                   <  1:862654/149‑1 (MQ=255)
   gATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑cgca                                                                                                                                                >  1:100226/1‑146 (MQ=255)
     tatCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑cgcagg                                                                                                                                              >  1:520378/1‑144 (MQ=255)
         gATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGAGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCA‑GGAT‑Ca                                                                                                                                        >  1:375358/1‑149 (MQ=255)
          aTAAACAGCGTCGAGCCATTGGGGATTTGCTCCGTCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTG‑CACGGCGTGGTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCA‑GGAT‑CAgg                                                                                                                                      <  2:1088831/149‑1 (MQ=255)
          aTAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAAACGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGAA‑GGAT‑CAg                                                                                                                                       >  1:592695/1‑149 (MQ=255)
           tAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCA‑GGAT‑CAgg                                                                                                                                      <  1:1354831/149‑1 (MQ=255)
              aCAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCA‑GGAT‑CAgg                                                                                                                                      >  1:1508837/1‑146 (MQ=255)
              aCAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGCGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCA‑GGAT‑CAgg                                                                                                                                      <  2:1508837/146‑1 (MQ=255)
              aCAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTACAAACGCCGGTGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCA‑GGAT‑CAGGttt                                                                                                                                   <  1:1577251/149‑1 (MQ=255)
               cAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCA‑GGAT‑CAGGtttt                                                                                                                                  >  2:906961/1‑149 (MQ=255)
                  cgttcaGCCATTTGGGAGTTGCTCCGCCACTGTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCA‑GGAT‑CAGGTTTTGCt                                                                                                                               <  2:805826/144‑1 (MQ=255)
                   gTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCCGCGCCGCACCGCCATGATGG‑CGCA‑GGAT‑ACGGTTTTGCTc                                                                                                                              >  1:248989/1‑149 (MQ=255)
                    tcgtgtcATTGGGGCTTTTCTCCGCCGCTGTGCGGGCGATGCGCTCCTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTTGAAGGCAGCGCCGACCCGCCATGATGG‑CGCA‑GGAT‑CAGGTTTTGCTcc                                                                                                                             <  2:592695/143‑1 (MQ=255)
                    tCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCA‑GGAT‑CAGGTTTTGCTcc                                                                                                                             <  1:1047972/149‑1 (MQ=255)
                     cgcgCCATTGGGGATTTGCTGCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCA‑GGAT‑CAGGTTTTGCTCCg                                                                                                                            <  2:131749/146‑1 (MQ=255)
                      gAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCA‑GGAT‑CAGGTTTTGCTCCGc                                                                                                                           >  1:976301/1‑149 (MQ=255)
                        gCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCA‑GGAT‑CAGGTTTTGCTCCGCCa                                                                                                                         <  2:976301/149‑1 (MQ=255)
                         ccATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCA‑GGAT‑CAGGTTTTGCTCCGCCAg                                                                                                                        >  2:128690/1‑149 (MQ=255)
                          cATTGGGGATTTGCCCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTGCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCA‑GGAT‑CAGGTTTTGCTCCGCCAGc                                                                                                                       <  1:1021959/149‑1 (MQ=255)
                            ttGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCA‑GGAT‑CAGGTTTTGCTCCGCCAGCTc                                                                                                                     >  1:1408776/1‑149 (MQ=255)
                            ttGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCA‑GGAT‑CAGGTTTTGCTCCGCCAGCTc                                                                                                                     >  2:1566093/1‑149 (MQ=255)
                               gggATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGTG‑CA‑GGAT‑CAGGTTTTGCTCCGCCAGCTCAtt                                                                                                                  >  1:1367049/1‑149 (MQ=255)
                               gggATGTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCACACCGCCATGATGG‑CGCA‑GGAT‑CAGGTTTTGCTCCGCCAGCTCAtt                                                                                                                  <  1:1426714/149‑1 (MQ=255)
                                ggATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCA‑GGAT‑CAGGTTTTGCTCCGCCAGCTCATTg                                                                                                                 <  2:100226/149‑1 (MQ=255)
                                ggATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCA‑GGAT‑CAGGTTTTGCTCCGCCAGCTCATTg                                                                                                                 <  1:402376/149‑1 (MQ=255)
                                  aTTTGCTCCGCCACTTTGCGGGCGATGCGCTCAATTCATCAGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCA‑GGAT‑CAGGTTTTGCTCCGCCAGCTCATTGAg                                                                                                               <  1:1107153/149‑1 (MQ=255)
                                      gCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCA‑GGAT‑CAGGTTTTGCTCCGCCAGCTCATTGAGGTcg                                                                                                           >  2:1447153/1‑149 (MQ=255)
                                        tCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCA‑GGAT‑CAGGTTTTGCTCCGCCAGCTCATTGAGGTcgcg                                                                                                         >  1:1066259/1‑149 (MQ=255)
                                          cGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCCGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCA‑GGAT‑CAGGTTTTGCTCCGCCAGCTCATTGAGGTCgcggc                                                                                                       <  2:562248/149‑1 (MQ=255)
                                           gCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCA‑GGAT‑CAGGTTTTGCTCCGCCAGCTCATTGAGGTCgaggcg                                                                                                      >  2:1210914/1‑149 (MQ=255)
                                            ccACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCA‑GGAT‑CAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGa                                                                                                     >  1:485110/1‑149 (MQ=255)
                                            ccACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCA‑GGAT‑CAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGa                                                                                                     >  1:123106/1‑149 (MQ=255)
                                                       gCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCA‑GGAT‑CAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCgg                                                                                          <  1:146781/149‑1 (MQ=255)
                                                        cGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAGCCGAACTGGAAGGCAGCGCCGCACCGCCATGCTGG‑CGGA‑GGGT‑CAGGTTTTGCTCCGCCAGTTCATTGAGGTCGCGGCGAATAGTCTGCggg                                                                                         <  1:864509/149‑1 (MQ=255)
                                                            gcgcTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCA‑GGAT‑CAGGTTTTGCTCCGCCAGCTAATTGAGGGCGCGGCGAATAGTCTGCGGGCTgc                                                                                     >  2:1292672/1‑148 (MQ=255)
                                                              gcTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCA‑GGAT‑CAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACg                                                                                   <  1:1411844/149‑1 (MQ=255)
                                                              gcTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCA‑GGAT‑CAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACg                                                                                   >  1:1409231/1‑149 (MQ=255)
                                                                tctTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCA‑GGAT‑CAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGga                                                                                 >  2:810012/1‑149 (MQ=255)
                                                                  tttttCTTCGGTCTGGGTGGCCTTGCGATCGTGCCAAGGCGTGTTAACCGGAGTGGGAGGCAGCGACGCACCGCCATGCTGG‑CGCA‑AGAT‑CATGTCTCGCTCCGCCATCTCATTGAGGTCACGCCGAATAGCCTCCGGGCTGACgtgca                                                                               >  2:789351/1‑145 (MQ=255)
                                                                    tttCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCA‑GGAT‑CAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGCAGaaa                                                                             >  2:809807/1‑149 (MQ=255)
                                                                      tcctcGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGCTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CTCA‑GGAT‑CAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGACTAGTCTGCGGGCTGACGGAGAAATg                                                                           <  1:1226731/146‑1 (MQ=255)
                                                                       cttcGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGACGGCAGCGCCGCCCCGCCATGATGG‑CGCA‑GGAT‑CAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGc                                                                          <  1:1541013/149‑1 (MQ=255)
                                                                        ttcGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCA‑GGAT‑CAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGct                                                                         <  1:1335880/149‑1 (MQ=255)
                                                                         ccGGTCTGGGTGGCCCTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCA‑GGAT‑CAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGctc                                                                        <  1:488010/148‑1 (MQ=255)
                                                                           ggTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAGCTGGAAGGCAGCGCCGCACCGCCATTATGG‑CGCA‑GGAT‑CAGGTTTTGCTCCGCCAGCTCATTGAGGTAGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTa                                                                      >  2:1308625/1‑149 (MQ=255)
                                                                           ggTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCA‑GGAT‑CAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTa                                                                      <  2:191422/149‑1 (MQ=255)
                                                                            gTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCA‑GGAT‑CAGGTCACGGCTGGCCGTCGCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTAc                                                                     >  1:780976/1‑149 (MQ=255)
                                                                             tCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCA‑GGAT‑CAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTAcc                                                                    >  2:890754/1‑149 (MQ=255)
                                                                             tCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCA‑GGAT‑CAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTAcc                                                                    <  2:1066259/149‑1 (MQ=255)
                                                                             tCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCA‑GGAT‑CAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTAcc                                                                    >  1:1351332/1‑149 (MQ=255)
                                                                                      ccTTGCGATCGGGCCACGGCGGGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCA‑GGAT‑CAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCTGGTTTTCg                                                           <  2:1003240/149‑1 (MQ=255)
                                                                                       cTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCA‑GGAT‑CAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACTAGCTCTTCgg                                                          <  1:1566093/149‑1 (MQ=255)
                                                                                       cTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCA‑GGAT‑CAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCgg                                                          >  1:384776/1‑149 (MQ=255)
                                                                                       cTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGA‑A‑GG‑GAATGCAGGGTTTGCTCCGTCTGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGCCGGAGAAATGCCCTGCTAGctctcccg                                                          >  2:642/1‑145 (MQ=255)
                                                                                         tGCGATCGTGGCACGGCGGTTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCA‑GGAT‑CAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTa                                                        <  2:1284187/149‑1 (MQ=255)
                                                                                         tGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCA‑GGAT‑CAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTa                                                        <  2:520378/149‑1 (MQ=255)
                                                                                          gCGATCGTGCCACGGCGGTTTAACCGAACTGGAAGGCCGCGCCGCACCGCCATGAGGG‑CGCA‑GGAT‑CAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTAc                                                       <  2:730578/149‑1 (MQ=255)
                                                                                             aTCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCA‑GGAT‑CAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGa                                                    >  1:1037694/1‑149 (MQ=255)
                                                                                                     aCGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCA‑GGAT‑CAGGTTTTGCCCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCCCTTCGGTACTGACATAAccc                                            >  2:418483/1‑149 (MQ=255)
                                                                                                       ggCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCA‑GGAT‑CAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCctg                                          <  1:95587/149‑1 (MQ=255)
                                                                                                       ggCGTCTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCA‑GGAT‑CAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGCACTGACATAACCctg                                          <  1:1366531/149‑1 (MQ=255)
                                                                                                         cGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCA‑GGAT‑CAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCctgct                                        <  1:890754/149‑1 (MQ=255)
                                                                                                           tgtTAACCGAACTGGAAGGCAGCGCCGCACCGCCTTGCTGG‑CGCA‑GCGT‑CATTCTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGt                                      <  2:1541092/149‑1 (MQ=255)
                                                                                                           ggtTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCA‑GGAT‑CAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGt                                      <  1:1032199/148‑1 (MQ=255)
                                                                                                             ttAACCGAACTGGAAGGCAGCGCCGCACCCCCATGATGG‑CGCA‑GGAT‑CAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGttt                                    <  1:529852/149‑1 (MQ=255)
                                                                                                                aCCGAACTGGAAGGCAGCGCCGGACTTCCATGATGG‑G‑GATGAAT‑TTCGAAAATCGCTGTCAGCTCATTGAGTTCGCGTCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAAc                                 <  1:1176348/149‑1 (MQ=255)
                                                                                                                aCCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCA‑GGAT‑CAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAAc                                 >  2:329894/1‑149 (MQ=255)
                                                                                                                 ccGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCA‑GGAT‑CAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGct                                                                         >  2:842932/1‑108 (MQ=255)
                                                                                                                 ccGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCA‑GGAT‑CAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGct                                                                         <  1:842932/108‑1 (MQ=255)
                                                                                                                     aCTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCA‑GGAT‑CAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGtt                            <  1:193754/149‑1 (MQ=255)
                                                                                                                      cTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCA‑GGAT‑CAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTc                           <  1:1210914/149‑1 (MQ=255)
                                                                                                                       tGGAAGGCAGCGCCGCACCGCCATGCTGG‑CGCA‑GGAA‑CAGATTTAAATCCGCCAGCACATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGCTTAACCaccgcg                          >  1:232107/1‑144 (MQ=255)
                                                                                                                        ggAAGGCAGCGCCGCACCGCCATGATGG‑CGCA‑GGAT‑CAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGa                         <  2:1549249/149‑1 (MQ=255)
                                                                                                                        ggAAGGCAGCGCCGCACCGCCATGATGG‑CGCA‑GGAT‑CAGGTTTGGCTCCGCCAGCTCATTGAGGTCGCTGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGa                         <  2:375358/149‑1 (MQ=255)
                                                                                                                         gAAGGCAGCGCCGCACCGCCATGATGG‑CGCA‑GGAT‑CAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGAc                                                   >  1:769758/1‑122 (MQ=255)
                                                                                                                         gAAGGCAGCGCCGCACCGCCATGATGG‑CGCA‑GGAT‑CAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGAc                                                   <  2:769758/122‑1 (MQ=255)
                                                                                                                             gCAGCGCCGCACCGCCATGATGG‑CGCA‑GGAT‑CAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGataata                    <  1:128690/149‑1 (MQ=255)
                                                                                                                               aGCGCCGCACCGCCATGATGG‑CGCA‑GGAT‑AAGGTTTTGCTCCGCGAGCTCATTGCGGACGCAGCGAATAGGCTGCGGGCTGACGGATAAATGCTCTACCATCGCTTCGGTACTGACATAACCCTGTCGTTTAACCAGTTCGATAGTAcc                  >  1:1361010/1‑149 (MQ=255)
                                                                                                                               aGCGCCGCACCGCCATGATGG‑CGCA‑GGAT‑AAGGTTTTGCTCCGCAAGCTCATTGCGGACGCCGCGAATAGGCTGCGGGCTGACGGAGAAATGCTCTACCAGCTATTCGGTACTGACATAACCCTGGTGTTTAACCAGTTCGATAGTAcc                  >  1:1360991/1‑149 (MQ=255)
                                                                                                                                  gccgcACCGCCATGATGG‑CGCA‑GGAT‑CAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGtt               <  1:810012/149‑1 (MQ=255)
                                                                                                                                    cgcACCGCCATGATGG‑CGCA‑GGAT‑CAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCa                                                                   >  2:1597007/1‑95 (MQ=255)
                                                                                                                                    cgcACCGCCATGATGG‑CGCA‑GGAT‑CAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCa                                                                   <  1:1597007/95‑1 (MQ=255)
                                                                                                                                          gCCATGATGG‑CGGA‑GGAT‑CAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAAAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGtt       <  1:906961/149‑1 (MQ=255)
                                                                                                                                           ccATGATGG‑CGCA‑GGAT‑CAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTtg      >  1:71394/1‑149 (MQ=255)
                                                                                                                                             atgatgG‑CGCA‑GGAT‑CAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGCAATGCTCTACCAGCTCTTCGGTCCTGACATAACCCTGCTGTTTAACCAGTTAGATAATACCGTTGTGACGTtgtg    >  2:128674/1‑149 (MQ=255)
                                                                                                                                             atgatgG‑CGCA‑GGAT‑CAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACAAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCCTTCTGACGTtgtg    >  1:659234/1‑149 (MQ=255)
                                                                                                                                               gatgG‑CGCA‑GGAT‑CAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGtgtt  >  2:848380/1‑149 (MQ=255)
                                                                                                                                                      |                                                                                                                                                
GCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑C‑CA‑GGAT‑CAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTT  >  NC_000913_3_pae_pgi/3560307‑3560597

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