Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,149,172 G→A 5.6% A605V (GCG→GTG)  yegI ← protein kinase‑related putative non‑specific DNA‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,149,1720GA 5.6% 225.5 / 2.3 72A605V (GCG→GTG) yegIprotein kinase‑related putative non‑specific DNA‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (37/30);  new base A (2/2);  total (40/32)
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.21e-01

AATTGAAAAATGGCTCAGCAGCGGCTTAAATCTGCCTGCGCCTGAGCCAGTTTTTCTGCTGCCTGACGTAACGGTTCCATCAACGGCATGGTCCGTGCTGGAGCATGAAGGCGGAATCGCTGCAACTCTGTCGCGCCGACTGTCAACGCCGATTCCAGCCGATGTCGTTTGGCGGTCATTTTCGCCATTACCGCTTGTCTGTCTGCTGGCGTTATCGCTTCATTCGGCCTGAACACAAAACGGCGTTCACAGCTCGCTTTCCAGTCGATAACCGCCTGGGTCAGATGGTCGCCA  >  NC_000913/2149025‑2149318
                                                                                                                                                   |                                                                                                                                                  
aaTTGAAAAATGGCTCAGCAGCGGCTTAAATCTGCCTGCGCCTGAGCCAGTTTTTCTGCTGCCTGACGTAACGGTTCCATCAACGGCATGGTCCGTGCTGGAGCATGAAGGCGGAATCGCTGCAACTCTGTCGCGCCGACTGTCAACGc                                                                                                                                                   >  6:267099/1‑149 (MQ=255)
            gCTCAGCAGCGGCTTAAATCTGCCTGCGCCTGAGCCAGTTTTTCTGCTGCCTGACGTAACGGTTCCATCAACGGCATGGTCCGTGCTGGAGCATGAAGGCGGAATCGCTGCAACTCTGTCGCGCCGACTGTCAACGCCGATTCCAGcc                                                                                                                                        <  1:14452/148‑1 (MQ=255)
              tCAGCAGCGGCTTAAATCTGCCTGCGCCTGAGCCAGTTTTTCTGCTGCCTGACGTAACGGTTCCATCAACGGCATGGTCCGTGCTGGAGCATGAAGGCGGAATCGCTGCAACTCTGTCGCGCCGACTGTCAACGCCGATTCCAGCCGAt                                                                                                                                     >  1:395902/1‑149 (MQ=255)
              tCAGCAGCGGCTTAAATCTGCCTGCGCATGAGCCAGTTTTTCTGCTGCCTGACGTAACGGTTCCATCAACTGCATGGTCCGTGCTGGAGCATGAAGGCGGAATCGCTGCAACTCTGTCGCGCCGACTGTCAACGCCGATTCCAGCCGAt                                                                                                                                     >  1:401898/1‑149 (MQ=255)
                   agcGGCTTAAATCTGCCTGCGCCTGAGCCAGTTTTTCTGCTGCCTGACGTAACGGTTCCATCAACGGCATGGTCCGTGCTGGAGCATGAAGGCGGAATCGCTGCAACTCTGTCGCGCCGACTGTCAACGCCGATTCCAGCCGATGTCGt                                                                                                                                >  8:319343/1‑149 (MQ=255)
                       gCTTAAATCTGCCTGCGCCTGAGCCAGTTTTTCTGCTGCCTGACGTAACGGTTCCATCAACGGCATGGTCCGTGCTGGAGCATGAAGGCGGAATCGCTGCAACTCTGTCGCGCCGACTGTCAACGCCGATTCCAGCCGATGTCGTTTgg                                                                                                                            <  5:49948/149‑1 (MQ=255)
                        cTTAAATCTGCCTGCGCCTGAGCCAGTTTTTCTGCTGCCTGACGTAACGGTTCCATCAACGGCATGGTCCGTGCTGGAGCATGAAGGCGGAATCGCTGCAACTCTGTCGCGCCGACTGTCAACGCCGATTCCAGCCGATGTCGTTTGGc                                                                                                                           <  7:164884/149‑1 (MQ=255)
                         ttAAATCTGCCTGCGCCTGAGCCAGTTTTTCTGCTGCCTGACGTAACGGTTCCATCAAGGGCATGGTCCGTGCGGGAGCATGAAGGCGGAATCGCTGCAACTCTGTCGCGCCGCCTGTCAACGCCGATACCAGCCGATGTCGTTTGCCg                                                                                                                          >  4:14459/1‑149 (MQ=255)
                           aaaTCTGCCTGCGCCTGAGCCAGTTTTTCTGCTGCCTGACGTAACGGTTCCATCAACGGCATGGTCCGTGCTGGAGCATGAAGGCGGAATCGCTGCAACTCTGTCGCGCCGACTGTCAACGc                                                                                                                                                   >  8:133595/1‑122 (MQ=255)
                           aaaTCTGCCTGCGCCTGAGCCAGTTTTTCTGCTGCCTGACGTAACGGTTCCATCAACGGCATGGTCCGTGCTGGAGCATGAAGGCGGAATCGCTGCAACTCTGTCGCGCCGACTGTCAACGc                                                                                                                                                   <  7:133595/122‑1 (MQ=255)
                             aTCTGCCTGCGCCTGAGCCAGTTTTTCTGCTGCCTGACGTAACGGTTCCATCAACGGCATGGTCCGTGCTGGAGCATGAAGGCGGAATCGCTGCAACTCTGTCGCGCCGACTGTCAACGCCGATTCCAGCCGATGTCGTTTGGCGGTCa                                                                                                                      >  3:163914/1‑149 (MQ=255)
                                tgcctgcGCCTGAGCCAGTTTTTCTGCTGCCTGACGTAACGGTTCCATCAACGGCATGGTCCGTGCTGGAGCATGAAGGCGGAATCGCTGCAACTCTGTCGCGCCGACTGTCAACGCCGATTCCAGCCGATGTCGTTTGGCGGTCAttt                                                                                                                   >  1:206974/1‑149 (MQ=255)
                                  cctgcGCCTGAGCCAGTTTTTCTGCTGCCTGACGTAACGGTTCCATCAACGGCATGGTCCGTGCTGGAGCATGAAGGCGGAATCGCTGCAACTCTGTCGCGCCGACTGTCAACGCCGATTCCAGCCGATGTCGTTTGGCGGTCAttt                                                                                                                   <  1:202301/147‑1 (MQ=255)
                                      cgcCTGAGCCAGTTTTTCTGCTGCCTGACGTAACGGTTCCATCAACGGCATGGTCCGTGCTGGAGCATGAAGGCGGAATCGCTGCAACTCTGTCGCGCCGACTGTCAACGCCGATTCCAGCCGATGTCGTTTGGCGGTCATTTTCGCCa                                                                                                             <  3:375506/149‑1 (MQ=255)
                                        ccTGAGCCAGTTTTTCTGCTGCCTGACGTAACGGTTCCATCAACGGCATGGTCCGTGCTGGAGCATGAAGGCGGAATCGCTGCAACTCTGTCGCGCCGACTGTCAACGCCGATTCCAGCCGATGTCGTTTGGCGGTCATTTTCGCCAt                                                                                                            >  7:159695/1‑148 (MQ=255)
                                               cAGTTTTTCTGCTGCCTGACGTAACGGTTCCATCAACGGCATGGTCCGTGCTGGAGCATGAAGGCGGAATCGCTGCAACTCTGTCGCGCCGACTGTCAACGCCGATTCCAGCCGATGTCGTTTGGCGGTCATTTTCGCCATTACCGCtt                                                                                                    >  7:160697/1‑149 (MQ=255)
                                               cAGTTTTTCTGCTGCCTGACGTAACGGTTCCATCAACGGCATGGTCCGTGCTGGAGCATGAAGGCGGAATCGCTGCAACTCTGTCGCGCCGACTGTCAACGCCGATTCCAGCCGATGTCGTTTGGCGGTCATTTTCGCCATTACCGCtt                                                                                                    >  5:302823/1‑149 (MQ=255)
                                               cAGTTTTTCTGCTGCCTGACGTAACGGTTCCATCAACGGCATGGTCCGTGCTGGAGCATGAAGGCGGAATCGCTGCAACTCTGTCGCGCCGACTGTCAACGCCGATTCCAGCCGATGTCGTTTGGCGGTCATTTTCGCCATTACCGCtt                                                                                                    <  1:220670/149‑1 (MQ=255)
                                                       ctgctgCCTGACGTAACGGTTCCATCAACGGCATGGTCCGTGCTGGAGCATGAAGGCGGAATCGCTGCAACTCTGTCGCGCCGACTGTCAACGCCGATTCCAGCCGATGTCGTTTGGCGGTCATTTTCGCCATTACCGCTtgtctgtc                                                                                             <  7:319343/148‑1 (MQ=255)
                                                         gctgcCTGACGTAACGGTTCCATCAACGGCATGGTCCGTGCTGGAGCATGAAGGCGGAATCGCTGCAACTCTGTCGCGCCGACTGTCAACGCCGATTCCAGCCGATGTCGTTTGGCGGTCATTTTCGCCATTACCGCTtgtctgaccc                                                                                           >  1:233825/1‑144 (MQ=255)
                                                          ctgcCTGACGTAACGGTTCCATCAACGGCATGGTCCGTGCTGGAGCATGAAGGCGGAATCGCTGCAACTCTGTCGCGCCGACTGTCAACGCCGATTCCAGCCGATGTCGTTTGGCGGTCATTTTCGCCATTACCGCTTGtctgtctg                                                                                           <  1:196045/147‑1 (MQ=255)
                                                               tGACGTAACGGTTCCATCAACGGCATGGTCCGTGCTGGAGCATGAAGGCGGAATCGCTGCAACTCTGTCGCGCCGACTGTCAACGCCGATTCCAGCCGATGTCGTTTGGCGGTCATTTTCGCCATTACCGCTtgt                                                                                                  >  7:198115/1‑135 (MQ=255)
                                                               tGACGTAACGGTTCCATCAACGGCATGGTCCGTGCTGGAGCATGAAGGCGGAATCGCTGCAACTCTGTCGCGCCGACTGTCAACGCCGATTCCAGCCGATGTCGTTTGGCGGTCATTTTCGCCATTACCGCTtgt                                                                                                  <  8:198115/135‑1 (MQ=255)
                                                               tGACGTAACGGTTCCATCAACGGCATGGTCCGTGCTGGAGCATGAAGGCGGAATCGCTGCAACTCTGTCGCGCCGACTGTCAACACCGATTCCAGCCGATGTCGTTTGGCGGTCATTTTCGCCATTAc                                                                                                         <  6:39602/128‑1 (MQ=255)
                                                               tGACGTAACGGTTCCATCAACGGCATGGTCCGTGCTGGAGCATGAAGGCGGAATCGCTGCAACTCTGTCGCGCCGACTGTCAACACCGATTCCAGCCGATGTCGTTTGGCGGTCATTTTCGCCATTAc                                                                                                         >  5:39602/1‑128 (MQ=255)
                                                                gACGTAACGGTTCCATCAACGGCATGGTCCGTGCTGGAGCATGAAGGCGGAATCGCTGCAACTCTGTCGCGCCGACTGTCAACGCCGATTCCAGCCGATGTCGTTTGGCGGTCATTTTCGCCATTACCGCTTGTCTGTCTGCTGGCGtt                                                                                   >  7:9352/1‑149 (MQ=255)
                                                                gACGTAACGGTTCCATCAACGGCATGGTCCGTGCTGGAGCATGAAGGCGGAATCGCTGCAACTCTGTCGCGCCGACTGTCAACCCCGATTCCAGCCGATGTCGTTTGGCCGTCATTTTCGCCATAACCGCTTGTCTGTCTGCAGGCGt                                                                                    >  6:370979/1‑148 (MQ=255)
                                                                   gTAACGGTTCCATCAACGGCATGGTCCGTGCTGGAGCATGAAGGCGGAATCGCTGCAACTCTGTCGCGCCGACTGTCAACGCCGATTCCAGCCGATGTCGTTTGGCGGTCATTTTCGCCATTACCGCTTGTCTGTCTGCTGGCGTTATc                                                                                >  4:171109/1‑149 (MQ=255)
                                                                   gTAACGGTTCCATCAACGGCATGGTCCGTGCTGGAGCATGAAGGCGGAATCGCTGCAACTCTGTCGCGCCGACTGTCAACGCCGATTCCAGCCGATGTCGTTTGGCGGTCATTTTCGCCATTACCGCTTGTCTGTCTGCTGGCGTTATc                                                                                >  7:117841/1‑149 (MQ=255)
                                                                    tAACGGTTCCATCAACGGCATGGTCCGTGCTGGAGCATGAAGGCGGAATCGCTGCAACTCTGTCGCGCCGACTGTCAACGCCGATTCCAGCCGATGTCGTTTGGCGGTCATTTTCGCCATTACCGCTTGTCTGTCTGCTGGCGTTATCg                                                                               <  8:46252/149‑1 (MQ=255)
                                                                    tAACGGTTCCATCAACGGCATGGTCCGTGCTGGAGCATGAAGGCGGAATCGCTGCAACTCTGTAGCGCCGACTGTCAACGCCGATTCCAGCCGATGTCGTTTGGCGGTCATTTTCGCCATTACCGCTAGTCTGTCTGCTGGCGTTATCg                                                                               >  7:313072/1‑149 (MQ=255)
                                                                       cGGTTCCATCAACGGCATGGTCCGTGCTGGAGCATGAAGGCGGAATCGCTGCAACTCTGTCGCGCCGACTGTCAACGCCGATTCCAGCCGATGTCGTTTGGCGGTCATTTTCGCCATTACCGCTTGTCTGTCTGCTGGTGTTATCGCtt                                                                            <  8:160697/149‑1 (MQ=255)
                                                                        ggTTCCATCAACGGCATGGTCCGTGCTGGAGCATGAAGGCGGAACCGCTGCAACTCTGTCGCGCCGACTGTCAACGCCGATTCCAGCCGATGTCGTTTGGCGGTCATTTTCGCCATTACCGCTTGTCTGTCTGCTGGCGTTa                                                                                  >  2:401781/1‑142 (MQ=255)
                                                                        ggTTCCATCAACGGCATGGTCCGTGCTGGAGCATGAAGGCGGAACCGCTGCAACTCTGTCGCGCCGACTGTCAACGCCGATTCCAGCCGATGTCGTTTGGCGGTCATTTTCGCCATTACCGCTTGTCTGTCTGCTGGCGTTa                                                                                  <  1:401781/142‑1 (MQ=255)
                                                                         gTTCCATCAACGGCATGGTCCGTGCTGGAGCATGAAGGCGGAATCGCTGCAACTCTGTCGCGCCGACTGTCAACGCCGATTCCAGCCGATGTCGTTTGGCGGTCATTTTCGCCATTACCGCTTGTCTGTCTGCTGGCGTTATCGCttaa                                                                          <  8:9352/149‑3 (MQ=255)
                                                                             cATCAACGGCATGGTCCGTGCTGGAGCATGAAGGCGGAATTGCTGCAACTCTGTCGCGCCGACTGTCAACGCCGATTCCAGCCGATGTCGTTTGGCGGTCATTTTCGCCATTACCGCTTGTCTGTCTGCTGGCGTTATCGCTTCATTCg                                                                      >  4:239402/1‑149 (MQ=255)
                                                                             cATCAACGGCATGGTCCGTGCTGGAGCATGAAGGCGGAATCGCTGCAACTCTGTCGCGCCGACTGTCAACGCCGATTCCCGCCGATGTCGTTTGGCGGTCATTTTCGCCATTACCGCTTGTCTGTCTGCTGGCGTTATCGCTTCATTCg                                                                      >  6:348598/1‑149 (MQ=255)
                                                                             cATCAACGGCATGGTCCGTGCTGGAGCATGAAGGCGGAATCGCTGCAACTCTGTCGCGCCGACTGTCAACGCCGATTCCAGCCGATGTCGTTTGGCGGTCATTTTCGCCATTACCGCTTGTCTGTCTGCTGGCGTTATCGCTTCATTCg                                                                      >  6:280057/1‑149 (MQ=255)
                                                                              aTCAACGGCATGGTCCGTGCTGGAGCATGAAGGCGGAATCGCTGCAACTCTGTCGCGCCGACTGTCAACGCCGATTCCAGCCGATGTCGTTTGGCGGTCATTTTCGCCATTACCGCTTGTCTGTCTGCTGGCGTTATCGCTTCATTCgg                                                                     <  5:267099/149‑1 (MQ=255)
                                                                              aTCAACGGCATGGTCCGTGCTGGAGCATGAAGGCGGAATCGCTGCAACTCTGTCGCGCCGACTGTCAACGCCGATTCCAGCCGATGTCGTTTGGCGGTCATTTTCGCCATTACCGCTTGTCTGTCTGCTGGCGTTATCGCTTCATTCgg                                                                     >  7:313318/1‑149 (MQ=255)
                                                                                 aaCGGCATGGTCCGTGCTGGAGCATGAAGGCGGAATCGCTGCAACTCTGTCGCGCCGACTGTCAACGCCGATTCCAGCCGATGTCGTTTGGCGGTCATTTTCGCCATTACCGCTTGTCTGTCTGCTGGCGTTATCGCTTCATTCGGCCt                                                                  <  8:66210/149‑1 (MQ=255)
                                                                                 aaCGGCATGGTCCGTGCTGGAGCATGAAGGCGGAATCGCTGCAACTCTGTCGCGCCGACTGTCAACACCGATTc                                                                                                                                             >  3:344983/1‑74 (MQ=255)
                                                                                 aaCGGCATGGTCCGTGCTGGAGCATGAAGGCGGAATCGCTGCAACTCTGTCGCGCCGACTGTCAACACCGATTc                                                                                                                                             <  4:344983/74‑1 (MQ=255)
                                                                                    ggCATGGTCCGTGCTGTAGCATGAAGGCGGAATCGCTGCAACTCTGTCGCGCCGACTGTCAACGCCGATTCCAGCCGATGTCGTTTGGCGGTCATTTTCGCCATTACCGCTTGGCTGTCTGCTGGCGTTATCGCTTCATTCGGCCTGaa                                                               <  6:27631/149‑1 (MQ=255)
                                                                                      cATGGTCCGTGCTGGAGCATGAAGGCGGAATCGCTGCAACTCTGTCGCGCCGACTGTCAACGCCGATTCCAGCCGATGTCGTTTGGCGGTCATTTTCGCCATTACCGCTTGTCTGTCTGCTGGCGTTATCGCTTCATTCGGCCTGAaca                                                             <  2:401898/149‑1 (MQ=255)
                                                                                      cATGGTCCGTGCTGGAGCATGAAGGCGGAATCGCTGCAACTCTGTCGCGCCGACTGTCAACGCCGATTCCAGCCGATGTCGTTTGGCGGTCATTTTCGCCATTACCGCTTGTCTGTCTGCTGGCGTTATCGCTTCATTCGGCCTGAaca                                                             <  2:395902/149‑1 (MQ=255)
                                                                                           tCCGTGCTGGAGCATGAAGGCGGAATCGCTGCAACTCTGTCGCGCCGACTGTCAACGCCGATTCCAGCCGATGTCGTTTGGCGGTCATTTTCGCCATTACCGCTTGTCTGTCTGCTGGCGTTATCGCTTCATTCGGCCTGAACACaaaa                                                        <  2:390911/149‑1 (MQ=255)
                                                                                               tGCTGGAGCATGAAGGCGGAATCGCTGCAACTCTGTCGCGCCGACTGTCAACGCCGATTCCAGCCGATGTCGTTTGGCGGTCATTTTCGCCATTACCGCTTGTCTGTCTGCTGGCGTTATCGCTTCATTCGGCCTGAACACAAAACGGc                                                    <  8:117841/149‑1 (MQ=255)
                                                                                                   ggAGCATGAAGGCGGAATCGCTGCAACTCTGTCTCGCCGACTGTCAACGCCGATTCCAGCCGATGTCGTTTGGCGGTCATTTTCGCCATTACCGCTTGTCTGTCTGCTGGCGTTATCGCTTCATTCGGCCTGAACACAAAACGGCGtt                                                 >  1:252186/1‑148 (MQ=255)
                                                                                                   ggAGCATGAAGGCGGAATCGCTGCAACTCTGTCGCGCCGACTGTCAACGCCGATTCCAGCCGATGTCGTTTGGCGGTCATTTTCGCCATTACCGCTTGTCTGTCTGCTGGCGTTATCGCTTCATTCGGCCTGAACACAAAACGGCGtt                                                 >  4:280744/1‑148 (MQ=255)
                                                                                                    gAGCATGAAGGCGGAATCGCTGCAACTCTGTCGCGCCGACTGTCAACGCCGATTCCAGCCGATGTCGTTTGGCGGTCATTTTCGCCATTACCGCTTGTCTGTCTGCTGGCGTTATCGCTTCATTCGGCCTGAACACAAAACGGCGTTca                                               >  8:117627/1‑149 (MQ=255)
                                                                                                      gCATGAAGGCGGAATTGCTGCAACTCTGTCGCGCCGACTGTCAACGCCGATTCCAGCCGATGTCGTTTGGCGGTCATTTTCGCCATTACCGCTTGTCTGTCTGCTGGCGTTATCGCTTCATTCGGCCTGAACACAAAACGGCGTTcaca                                             <  3:239402/149‑1 (MQ=255)
                                                                                                             ggCGGAATCGCTGCAACTCTGTCGCGCCGACTGTCAACGCCGATTCCAGCCGATGTCGTTTGGCGGTCATTTTCGCCATTACCGCTTGTCTGTCTGATGGCGTTATCGCTTCATTCGGCCTGAACACAAAACGGCGTTCAAAGCTCGCt                                      >  6:32654/1‑149 (MQ=255)
                                                                                                              gCGGAATCGCTGCAACTCTGTCGCGCCGACTGTCAACGCCGATTCCAGCCGATGTCGTTTGGCGGTCATTTTCGCCATTACCGCTTGTCTGTCTGCTGGCGTTATCGCTTCATTCGGCCTGAACACAAAACGGCGTTCACAGCTCGCtt                                     <  4:163914/149‑1 (MQ=255)
                                                                                                                 gAATCGCTGCAACTCTGTCGCGCCGACTGTCAACGCCGATTCCAGCCGATGTCGTTTGGCGGTCATTTTCGCCATTACCGCTTGTCTGTCTGCTGGCGTTATCGCTTCATTCGGCCTGAACACAAAACGGCGTTCACAGCTCGCTTTc                                   >  7:65818/1‑148 (MQ=255)
                                                                                                                  aaTCGCTGCAACTCTGTCGCGCCGACTGTCAACGCCGATTCCAGCCGATGTCGTTTGGCGGTCATTTTCGCCATTACCGCTTGTCTGTCTGCTGGc                                                                                      <  8:354983/96‑1 (MQ=255)
                                                                                                                  aaTCGCTGCAACTCTGTCGCGCCGACTGTCAACGCCGATTCCAGCCGATGTCGTTTGGCGGTCATTTTCGCCATTACCGCTTGTCTGTCTGCTGGc                                                                                      >  7:354983/1‑96 (MQ=255)
                                                                                                                  aaTCGCTGCAACTCTGTCGCGCCGACTGTCAACGCCGATTCCAGCCGATGTCGTTTGGCGGTCATTTTCGCCATTACCGCTTGTCTGTCTGCTGGCGTTATCGCTTCATTCGGCCTGAACACAAAACGGCGTTCACAGCTCGCTTTCCa                                 <  4:238197/149‑1 (MQ=255)
                                                                                                                    tCGCTGCAACTATGTCGCTCCGACTGTCAACGCCGATTCCAGCCGATGTCGTTTGGCGGTCATTTTCGCCATTACCGCTTGTCTGTCTGCTGGCGTTATCGCTTCATTCGGCCTGAACACAAAACGGCGTTCACAGCTCGCTTTCCAGt                               <  8:313072/149‑1 (MQ=255)
                                                                                                                      gCTGCAACTCTGTCGCGCCGACTGTCAACGCCGATTCCAGCCGATGTCGTTTGGCGGTCATTTTCGCCATTACCGCTTGTCTGTCTGCTGGCGTTATCGCTTCATTCGGCCTGAACACAAAACGGCGTTCACAGCTCGCTTTCCAGTc                              <  8:159695/148‑1 (MQ=255)
                                                                                                                          cAACTCTGTCGCGCCGACTGTCAACGCCGATTCCAGCCGATGTCGTTTGGCGGTCATTTTCGCCATTACCGCTTGTCTGTCTGCTGGCGTTATCGCTTCATTCGGCCTGAACACAAAACGGCGTTCACAGCTCGCTTTCCAGTCGATaa                         >  4:128133/1‑149 (MQ=255)
                                                                                                                            aCTCTGTCGCGCCGAATGTCAACGCCGATTCCAGCCGATGTCGTTTGGCGGTCATTTTCGCCATTACCGCTTGTCTGTCTGCTGGCGTTATCGCTTCATTCGTCCTGAACACAAAACGGCGTTCACAGCTCGCTTTCCAGTCGATAAcc                       <  5:370979/149‑1 (MQ=255)
                                                                                                                                  tCGCGCCGACTGTCAACGCCGATTCCAGCCGATGTCGTTTGGCGGTCATTTTCGCCATTACCGCTTGTCTGTCTGCTGGCGTTATCGCTTCATTCGGCCTGAACACAAAACGGCGTTCACAGCTCGCTTTCCAGTCGATAACCGCCTgg                 <  3:171109/149‑1 (MQ=255)
                                                                                                                                        cGACTGTCAACGCCGATTCCAGCCGATGTCGTTTGGCGGTCATTTTCGCCATTACCGCTTGTCTGTCTGCTGGCGTTATCGCTTCATTCGGCCTGAACACAAAACGGCGTTCACAGCTCGCTTTCCAGTCGATAACCGCCTGGGt               >  3:399895/1‑145 (MQ=255)
                                                                                                                                            tGTCAACGCCGATTCCAGCCGATGTCGTTTGGCGGTCATTTTCGCCATTACCGCTTGTCTGTCTGCTGGCGTTATCGCTTCATTCGGCCTGAACACAAAACGGCGTTCACAGCTCGCTTTCCAGTCGATAACCGCCTGGGTCAGATGGt       >  3:160665/1‑149 (MQ=255)
                                                                                                                                              tCAACGCCGATTCCAGCCGATGTCGTTTGGCGGTCATTTTCGCCATTACCGCTTGTCTGTCTGCTGGCGTTATCGCTTCATTCGGCCTGAACACAAAACGGCGTTCACAGCTCGCTTTCCAGTCGATAACCGCCTGGGTCAGATGGTc      <  2:206974/148‑1 (MQ=255)
                                                                                                                                              tCAACGCCGATTCCAGCCGATGTCGTTTGGCGGTCATTTTCGCCATTACCGCTTGTCTGTCTGCTGGCGTTATCGCTTCATTCGGCCTGAACACAAAACGGCGTTCACAGCTCGCTTTCCAGTCGATAACCGCCTGGGTCAGATGGTCg     >  2:317393/1‑149 (MQ=255)
                                                                                                                                              tCAACGCCGATTCCAGCCGATGTCGTTTGGCGGTCATTTTCGCCATTACCGCTTGTCTGTCTGCTGGCGTTATCGCTTCATTCGGCCTGAACACAAAACGGCGTTCACAGCTCGCTTTCCAGTCGATAACCGCCTGGGTCAGAAGGTCg     >  4:280422/1‑149 (MQ=255)
                                                                                                                                                aaCGCCGATTCCAGCCGATGTCGTTTGGCGGTCATTTTCGCCATTACCGCTTGTCTGTCTGCTGGCGTTATCGCTTCATTCGGCCTGAACACAAAACGGCGGTCACAGCTCTCTTTCCAGTCGATAACCGCCTGGGTCAGATGGTCGc    >  6:280063/1‑148 (MQ=255)
                                                                                                                                                 aCGCCGATTCCAGCCGATGTCGTTTGGCGGTCATTTTCGCCATTACCGCTTGTCTGTCTGCTGGCGTTATCGCTTCATTCGGCCTGAACACAAAACGGCGTTCACAGCTCGCTTTCCAGTCGATAACCGCCTGGGTCAGATGGTCGcc   >  8:302050/1‑148 (MQ=255)
                                                                                                                                                 aCGCCGATTCCAGCCGATGTCGTTTGGCGGTCATTTTCGCCATTACCGCTTGTCTGTCTGCTGGCGTTATCGCTTCATTCGGCCTGAACACAAAACGGCGTTCACAGCTCGCTTTCCAGTCGATAACCGCCTGGGTCAGATGGTCGcc   >  4:143134/1‑148 (MQ=255)
                                                                                                                                                  cGCCGATTCCAGCCGATGTCGTTTGGCGGTCATTTTCGCCATTACCGCTTGTCTGTCTGCTGGCGTTATCGCTTCATTCGGCCTGAACACAAAACGGCGTTCACAGCTCGCTTTCCAGTCGATAACCGCCTGGGTCAGATGGTCGCCa  >  8:344284/1‑148 (MQ=255)
                                                                                                                                                   gCCGATTCCAGCCGATGTCGTTTGGCGGTCATTTTCGCCATTACCGCTTGTCTGTctgctg                                                                                        <  6:94264/61‑1 (MQ=255)
                                                                                                                                                   gCCGATTCCAGCCGATGTCGTTTGGCGGTCATTTTCGCCATTACCGCTTGTCTGTctgctg                                                                                        >  5:94264/1‑61 (MQ=255)
                                                                                                                                                   gCCGATTCCAGCCGATGTCGTTTGGCGGTCATTTTCGCCATTACCGCTTGTCTGTCTGCTGGCGTTATCGCTTCATTCGGCCTGAACACAAAACg                                                      >  6:396500/1‑95 (MQ=255)
                                                                                                                                                   gCCGATTCCAGCCGATGTCGTTTGGCGGTCATTTTCGCCATTACCGCTTGTCTGTCTGCTGGCGTTATCGCTTCATTCGGCCTGAACACAAAACg                                                      <  5:396500/95‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
AATTGAAAAATGGCTCAGCAGCGGCTTAAATCTGCCTGCGCCTGAGCCAGTTTTTCTGCTGCCTGACGTAACGGTTCCATCAACGGCATGGTCCGTGCTGGAGCATGAAGGCGGAATCGCTGCAACTCTGTCGCGCCGACTGTCAACGCCGATTCCAGCCGATGTCGTTTGGCGGTCATTTTCGCCATTACCGCTTGTCTGTCTGCTGGCGTTATCGCTTCATTCGGCCTGAACACAAAACGGCGTTCACAGCTCGCTTTCCAGTCGATAACCGCCTGGGTCAGATGGTCGCCA  >  NC_000913/2149025‑2149318

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