Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,171,481 G→A 5.9% pseudogene (247/299 nt) gatR ← pseudogene, repressor for gat operon; interrupted by IS3; split galactitol utilization operon repressor, fragment 2; split galactitol utilization operon repressor, interrupted

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,171,4810GA 5.9% 204.7 / 3.2 68pseudogene (247/299 nt)gatRpseudogene, repressor for gat operon; interrupted by IS3; split galactitol utilization operon repressor, fragment 2; split galactitol utilization operon repressor, interrupted
Reads supporting (aligned to +/- strand):  ref base G (36/27);  new base A (2/2);  total (39/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.06e-01

GCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTATTACCAGACATTATTTTCGCCGCACCGCCATGAAAGCGCGTAACCACGCCTTTTTGTTCGAGAAAGCGCAAATCGGCACGGATTGTCGCT  >  NC_000913/2171335‑2171629
                                                                                                                                                  |                                                                                                                                                    
gcggcggcACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGgcg                                                                                                                                                    >  3:416165/1‑149 (MQ=255)
gcggcggcACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGgcg                                                                                                                                                    >  3:705475/1‑149 (MQ=255)
      gcACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTgg                                                                                                                                              >  3:488251/1‑149 (MQ=255)
        acaGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCg                                                                                                                                            >  1:206951/1‑149 (MQ=255)
            ttGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCt                                                                                                                                        >  1:35932/1‑149 (MQ=255)
            ttGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCt                                                                                                                                        >  1:507980/1‑149 (MQ=255)
                        ccTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCtt                                                                                                                            >  1:369950/1‑149 (MQ=255)
                         cTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCttt                                                                                                                           >  1:123574/1‑149 (MQ=255)
                             gTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAAc                                                                                                                       <  4:180203/149‑1 (MQ=255)
                              tCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAAcc                                                                                                                      >  4:783702/1‑149 (MQ=255)
                                 ggAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGa                                                                                                                   <  2:474797/149‑1 (MQ=255)
                                 ggAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGa                                                                                                                   >  3:277750/1‑149 (MQ=255)
                                        tttGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTg                                                                                                            >  4:561034/1‑149 (MQ=255)
                                        tttGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGACGCGCTGGCGAGCTGAAATCGCTCttt                                                                                                                           <  3:292793/134‑1 (MQ=255)
                                        tttGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGACGCGCTGGCGAGCTGAAATCGCTCttt                                                                                                                           >  2:169897/1‑134 (MQ=255)
                                        tttGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGACGCGCTGGCGAGCTGAAATCGCTCttt                                                                                                                           >  4:292793/1‑134 (MQ=255)
                                        tttGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGACGCGCTGGCGAGCTGAAATCGCTCttt                                                                                                                           <  1:169897/134‑1 (MQ=255)
                                         ttGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTgg                                                                                                           <  3:561034/149‑1 (MQ=255)
                                         ttGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTgg                                                                                                           >  3:22355/1‑149 (MQ=255)
                                         ttGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTgg                                                                                                           <  2:35932/149‑1 (MQ=255)
                                           gAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCt                                                                                                             >  1:703348/1‑145 (MQ=255)
                                           gAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCt                                                                                                             <  2:703348/145‑1 (MQ=255)
                                            aaaaCAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGtc                                                                                                        >  4:104603/1‑149 (MQ=255)
                                              aaCAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGtctc                                                                                                      >  4:460714/1‑149 (MQ=255)
                                                 aaGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGt                                                                                                   >  3:138674/1‑149 (MQ=255)
                                                   gAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGttt                                                                                                 >  4:73582/1‑149 (MQ=255)
                                                    aaCCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGa                                                                                                                   >  2:343272/1‑130 (MQ=255)
                                                    aaCCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGa                                                                                                                   <  1:343272/130‑1 (MQ=255)
                                                       cTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGTGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACt                                                                                             >  2:416313/1‑149 (MQ=255)
                                                       cTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACt                                                                                             <  3:783702/149‑1 (MQ=255)
                                                          gCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTAtt                                                                                          >  2:644752/1‑149 (MQ=255)
                                                            ttAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTATTAc                                                                                        <  4:416924/149‑1 (MQ=255)
                                                             tAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTATTAcc                                                                                       <  4:277750/149‑1 (MQ=255)
                                                               ggCCTGTGTCCATAGTACGTGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTATTACCAg                                                                                     >  2:163933/1‑149 (MQ=255)
                                                                gCCTTTGTCCCTATTACGTGGGTAGGATCAGTCATCCCTTGGTGGATCATTTTGACTGCCGCCTGCGAAATTCTGTGTTTTGGCGCGCTGGCGCGCTGAAATCGCTCTTTAAGGCCGACTTCCTGGGTCTCGGTTTCACTATTACCAGa                                                                                    >  2:698286/1‑149 (MQ=255)
                                                                  cTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTATTAc                                                                                        <  3:560462/143‑1 (MQ=255)
                                                                  cTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTATTAc                                                                                        >  4:560462/1‑143 (MQ=255)
                                                                   tgtgTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTATTACCAGACat                                                                                 >  1:197237/1‑149 (MQ=255)
                                                                         cATATTACGTGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTATTACCAGACATTAtttt                                                                           >  1:6613/1‑149 (MQ=255)
                                                                         cATATTACGTGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTATTACCAGACATTAtttt                                                                           >  4:561081/1‑149 (MQ=255)
                                                                           tatTACGTGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGTCGCGCTGGCGAGCTGAAATCGCTCTTTCAAccc                                                                                                                     >  3:360270/1‑105 (MQ=255)
                                                                           tatTACGTGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAAccc                                                                                                                     <  4:360270/105‑1 (MQ=255)
                                                                           tatTACGTGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTATTACCAGACATTATTTTcg                                                                         >  4:333824/1‑149 (MQ=255)
                                                                                         gATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTATTACCAGACATTATTTTCGCCGCACCGCCATGa                                                           >  3:632927/1‑149 (MQ=255)
                                                                                               tCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTATTACCAGACATTATTTTCGCCGCACCGCCATGAAAgcgc                                                     <  2:69676/149‑1 (MQ=255)
                                                                                                 aTCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTATTACCAGACATTATTTTCGCCGCACCGCCATGAAAGCGCGt                                                   <  1:166223/149‑1 (MQ=255)
                                                                                                     cTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTATTACCAGACATTATTTTCGCCGCACCGCCATGAAAGCGCGTAAcc                                               >  2:263470/1‑149 (MQ=255)
                                                                                                      ttCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTATTACCAGACATTATTTTCGCCGCACCGCCATGAAAGCGCGTAACCa                                              >  1:200229/1‑149 (MQ=255)
                                                                                                      ttCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTATTACCAGACATTATTTTCGCCGCACCGCCATGAAAGCGCGTAACCa                                              >  4:290880/1‑149 (MQ=255)
                                                                                                       tCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTATTACCAGCCATTATTCTCGCCGCACCGCCATGAAAGCGCGTAACCAc                                             <  2:123574/149‑1 (MQ=255)
                                                                                                         gTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTATTACCAGACATTATTTTCGCCGCACCGCCATGAAAGCGCGTAACCACGc                                           >  2:71204/1‑149 (MQ=255)
                                                                                                         gTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTATTACCAGACATTATTTTCGCCGCACCGCCATGAAAGCGCGTAACCACGc                                           <  2:206951/149‑1 (MQ=255)
                                                                                                            gATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTATTACCAGACATTATTTTCGCCGCACCGCCATGAAAGCGCGTAACCACGCCtt                                        >  3:700399/1‑149 (MQ=255)
                                                                                                             aTCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTATTACCAGACATTATTTTCGCCGCACCGCCATGAAAGCGCGTAACCACGCCttt                                       >  4:774617/1‑149 (MQ=255)
                                                                                                             aTCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTATTACCAGACATTATTTTCGCCGCACCGCCATGAAAGCGCGTAACCACGCCttt                                       >  2:522797/1‑149 (MQ=255)
                                                                                                              tCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTATTACCAGACATTATTTTCGCCGCACCGCCATGAAAGCGCGTAACCACGCCtttt                                      >  1:639056/1‑149 (MQ=255)
                                                                                                               cATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTATTACCAGACATTATTTTCGCCGCACCGCCATGAAAGCGCGTAACCACGCCttttt                                     <  4:22355/149‑1 (MQ=255)
                                                                                                                 ttttGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTATTACCAGACATTATTTTcgcc                                                                       <  2:779090/113‑1 (MQ=255)
                                                                                                                 ttttGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTATTACCAGACATTATTTTcgcc                                                                       >  1:779090/1‑113 (MQ=255)
                                                                                                                   ttGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTATTACCAGACATTATTTTCGCCGCACCGCCATGAAAGCGCGTAACCACGCCTTTTTGTTc                                 <  4:416165/149‑1 (MQ=255)
                                                                                                                   ttGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTATTACCAGACATTATTTTCGCCGCACCGCCATGAAAGCGCGTAACCACGCCTTTTTGTTc                                 <  4:705475/149‑1 (MQ=255)
                                                                                                                     gACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTATTACCAGACATTATTTTCGCCGCACCGCCATGAAAGCGCGTAACCACGCCTTTTTGTTCga                               <  4:488251/149‑1 (MQ=255)
                                                                                                                     gACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTATTACCAGACATTATTTTCGCCGCACCGCCATGAAAGCGCGTAACCACGCCTTTTTGTTCga                               <  2:369950/149‑1 (MQ=255)
                                                                                                                          ccgccTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTATTACCAGACATTATTTTCGCCGCACCGCCATGAAAGCGCGTAACCACGCCTTTTTGTTCGAGAAAg                          <  3:104603/149‑1 (MQ=255)
                                                                                                                                   cTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTATTACCAGACATTATTTTCGCCGCACCGCCATGAAAGCGCGTAACCACGCCTTTTTGTTCGAGAAAGCGCAAATCg                 >  2:653891/1‑149 (MQ=255)
                                                                                                                                    tATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTATTACCAGACATTATTTTCGCCGCACCGCCATGAAAGCGCGTAACCACGCCTTTTTGTTCGAGAAAGCGCAAATCgg                <  4:632927/149‑1 (MQ=255)
                                                                                                                                    tATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTATTACCAGACATTATTTTCGCCGCACCGCCATGAAAGCGCGTAACCACGCCTTTTTGTTCGAGAAAGCGCAAATCgg                <  4:773152/149‑1 (MQ=255)
                                                                                                                                           ttttttGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTATTACCAGACATTATTTTCGCCGCACCGCCATGAAAGCGCGTAACCACGCCTTTTTGTTCGAGAAAGCGCAAATCGGCACGGAt         <  2:200229/149‑1 (MQ=255)
                                                                                                                                             ttttGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTATTACCAGACATTATTTTCGCCGCACCGCCATGAAAGCGCGTAACCACGCCTTTTTGTTCGAGAAAGCGCAAATCGGCACGGATTg       <  1:644752/149‑1 (MQ=255)
                                                                                                                                               ttGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTATTACCAGACATTATTTTCGCCGCACCGCCATGAAAGCGCGTAACCACGCCTTTTTGTTCGAGAAAGCGCAAATCGGCACGGATTGTc     <  1:71204/149‑1 (MQ=255)
                                                                                                                                                  gcgcgcTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTATTACCAGACATTATTTTCGCCGCACCGCCATGAAAGCGCGTAACCACGCCTTTTTGTTCGAGAAAGCGCAAATCGGCACGGATTGTCGCt  >  3:397371/1‑149 (MQ=255)
                                                                                                                                                  gcgcgcTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTATTACCAGACATTATTTTCGCCGCACCGCCATGAAAGCGCGTAACCACGCCTTTTTGTTCGAGAAAGCGCAAATCGGCACGGATTGTCGCt  >  1:437373/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTATTACCAGACATTATTTTCGCCGCACCGCCATGAAAGCGCGTAACCACGCCTTTTTGTTCGAGAAAGCGCAAATCGGCACGGATTGTCGCT  >  NC_000913/2171335‑2171629

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: