Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,614,751 T→C 5.0% intergenic (‑48/‑53) sad ← / → yneJ succinate semialdehyde dehydrogenase, NAD(P)+‑dependent/putative DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,614,7510TC 5.0% 246.0 / 2.9 79intergenic (‑48/‑53)sad/yneJsuccinate semialdehyde dehydrogenase, NAD(P)+‑dependent/putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (42/33);  new base C (2/2);  total (44/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.96e-01

CGTTTTCGATATCGTCAGCGCCAGCCCACGGCAGCACAGAAAGTTGTTCACCCGTGGCAGGATTTATCGAAATTGCATGAGTTGCCGGAGTAATGGTCATCGGGGTATCTCCTTTATGAGTCATGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCACGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAAGCTGGCAGCATAACCCAGGCTGCAGCAAAAGTGCATCGCGTGCCGTCGAATCT  >  NC_000913/1614604‑1614898
                                                                                                                                                   |                                                                                                                                                   
cGTTTTCGATATCGTCAGCGCCAGCCCACGGCAGCACAGAAAGTTGTTCACCCGTGGCAGGATTTATCGAAATTGCATGAGTTGCCGGAGTAATGGTCATCGGGGTATCGCCTTTATTAGTCATGGTATGAATATACGCAGATTTAct                                                                                                                                                     >  2:221272/1‑148 (MQ=255)
cGTGTTCGATTTCGTCTGCGCCAGCCCACGGCAGCACAGAAAGTTGTTCACCCGTGGCAGGATTTATCGAAATTGCATGAGTTGCCGGAGTAAGGGTCATCGGGGTATCTCCTTTATGAGTCATGGTATGAAGATACGCAGATTTAct                                                                                                                                                     <  6:109149/148‑1 (MQ=255)
   tttCGATATCGTCAGCGCCAGCCCACGGTAGCACAGAAAGTTGTTCACCCGTGGCAGGATTTATCGAAATTGCATGAGTTGCCGGAGTAATGGTCATCGGGGTATCTCCTTTATGAGTCATGGTATGAAGATACGCAGATTTACTCtt                                                                                                                                                  >  4:173379/1‑148 (MQ=255)
         tatCGTCAGCGCCAGCCCACGGCAGCACAGAAAGTTGTTCACCCGTGGCAGGATTTATCGAAATTGCATGAGTTGCCGGAGTAATGGTCATCGGGGTATCTCCTTTATGAGTCATGGTATGAAGATACGCAGATTTACTCTTGCttt                                                                                                                                             >  8:320008/1‑147 (MQ=255)
         tatCGTCAGCGCCAGCCCACGGCAGCACAGAAAGTTGTTCACCCGTGGCAGGATTTATCGAAATTGCATGAGTTGCCGGAGTAATGGTCATCGGGGTATCTCCTTTATGAGTCATGGTATGAAGATAAGCAGATTTACTCTTGCttt                                                                                                                                             >  6:251271/1‑147 (MQ=255)
          atCGTCAGCGCCAGCCCACGGCAGCACAGAAAGTTGTTCACCCGTGGCAGGATTTATCGAAATTGCATGGGTTGCCGGAGTAATGGTCATCGGGGTATCTCCTTTATGAGTCATGGTATGAAGAGACGCAGATTGTCTCGTGCTTTa                                                                                                                                            >  6:145561/1‑147 (MQ=255)
           tCGTCAGCGCCAGCCCACGGCAGCACAGAAAGTTGTTCACCCGTGGCAGGATTTATCGAAATTGCATGAGTTGCCGGAGTAATGGTCATCGGGGTATCTCCTTTATGAGTCATGGTATGAAGATACGCAGATTTACTCTTGCTTTaaa                                                                                                                                          <  7:236532/148‑1 (MQ=255)
            cGTCAGCGCCAGCCCACGGCAGCACAGAAAGTTGTTCACCCGTGGCAGGATTTATCGAAATTGCATGAGTTGCCGGAGTAATGGTCATCGGGGTATCTGCTTTATGAGTCATGGTATGAAGATACGCAGATTTACTCTTGCTTTaaaa                                                                                                                                         >  5:356734/1‑148 (MQ=255)
              tCAGCGCCAGCCCACGGCAGCACAGAAAGTTGTTCACCCGTGGCAGGATTTATCGAAATTGCATGAGTTGCCGGAGTAATGGTCATCGGGGTATCTCCTTTATGAGTCATGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAAt                                                                                                                                        >  2:207468/1‑147 (MQ=255)
              tCAGCGCCAGCCCACGGCAGCACAGAAAGTTGTTCACCCGGGGCAGGATTTATCGAAATTGCATGAGTTGCCGGAGTAATGGTCATCGGGGTATCTCCTTTATGAGTCATGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATg                                                                                                                                       <  6:302282/148‑1 (MQ=255)
                    ccAGCCCACGTCAGCACAGAATTTTGTTCACCCGTGGCAGGATTTATCGAAATTGCATGAGTTGCCGGAGTAATGGTCATCGGGGTATCTCCTTTATGAGTCATGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATGaataat                                                                                                                                 <  5:41930/148‑1 (MQ=255)
                    ccAGCCCACGGCAGCACAGAAAGTTGTTCACCCGTGGCAGGATTTATCGAAATTGCATGAGTTGCCGGAGTAATGGTCATCGGGGTATCTCCTTTATGAGTCATGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATGaataa                                                                                                                                  >  8:159500/1‑147 (MQ=255)
                      cgCCCACGGCAGCACAGAAAGTTGTTCACCCGTGGCAGGATTTATCGAAATTGCATGAGTGGCCGGAGTAATGGTCATCGGGGTATCTCCTTTATGAGTCAGGGTATGAAGATACGCAGATTTACTCTTGCTTCAAAATGAATActat                                                                                                                               >  4:326213/2‑148 (MQ=255)
                        cccACGGCAGCACAGAAAGTTGTTCACCCGTGGCAGGATTTATCGAAATTGCATGAGTTGCCGGAGTAATGGTCATCGGGGTATCTCCTTTATGAGTCATGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATAtt                                                                                                                              >  7:184924/1‑147 (MQ=255)
                        cccACGGCAGCACAGAAAGTTGTTCACCCGTGGCAGGATTTATCGAAATTGCATGAGTTGCCGGAGTAATGGTCATCGGGGTATCTCCTTTAAGAGTCATGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTa                                                                                                                             <  2:276678/148‑1 (MQ=255)
                            cGGCAGCACAGAAAGTTGTTCACCCGTGGCAGGATTTATCGAAATTGCATGAGTTGCCGGAGTAATGGTCATCGGGGTATCTCCTTTATGAGTCATGGTATGAAGATACGCAGCTTTACTCTTGCTTTAAAATGAATAATATTAAGc                                                                                                                          >  6:191524/1‑147 (MQ=255)
                             ggCAGCACAGAAAGTTGTTCACCCGTGGCAGGATTTATCGAAATTGCATGAGTTGCCGGAGTAATGGTCATCGGGGTATCTCCTTTATGAGTCATGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCa                                                                                                                        <  6:356734/148‑1 (MQ=255)
                                  cacaGAAAGTTGTTCACCCGTGGCAGGATTTATCGAAATTGCATGAGTTGCCGGAGTAATGGTCATCGGGGTATCTCCTTTATGAGTCATGGTATGAAGATACGCAGATGTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTAt                                                                                                                   >  2:324676/1‑148 (MQ=255)
                                    caGAAAGTTGTTCACCCGTGGCAGGATTTATCGAAATTGCATGAGTTGCAGGAGTAATTGTCATCGGTGTATCTCCTTTATGAGTCATGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATAAGCCACTTATTCa                                                                                                                <  2:236564/148‑1 (MQ=255)
                                    caGAAAGTTGTTCACACGTGGCAGGATTTATCGAAATTGCATGAGTTGCCGGAGTAATGGTCATCGGGGTATCTCCTTTATGAGTCATGGTATGAAGATACGCAGATTTACCCTTGCTTTAAAATGAATAATATTAAGCCACTTAtt                                                                                                                  >  7:159205/1‑147 (MQ=255)
                                     aGAAAGTTGTTCACCCGTGGCAGGATTTATCGAAATTGGATGAGTTGCCGGAGTAATGGTCATCGGGGTATCTCCTTTATGAGTCAGGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCa                                                                                                                >  1:210290/1‑148 (MQ=255)
                                     aGAAAGTTGTTCACCCGTGGCAGGATTTATCGAAATTGCATGAGTTGCCGGAGTAATGGTCATCGGGGTATCTCCTTTATGAGTCATGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCa                                                                                                                <  7:320008/148‑1 (MQ=255)
                                              ttCACCCGTGTCAGGATTTATCGAAATTGCATGAGTTGCCGGAGTAATGGTCATCGGGGTATCTCCTTTATGAGTCATGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCACGAATCga                                                                                                        <  4:83266/147‑1 (MQ=255)
                                                cACCCGTGGCAGGATTTATCGAAATTGCATGAGTTGCCGGAGTAATGGTCATCGGGGTATCTCCTTTATGAGTCATGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCACGAATCgaga                                                                                                      <  2:49477/147‑1 (MQ=255)
                                                cACCCGTGGCAGGATTTATCGAAATTGCATGAGTTGCCGGAGTAATGGTCATCGGGGTATCTCCTTTATGAGTCATGGTATGAAGATACGCAGATTTACCCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCACGAATCgaga                                                                                                      <  8:159205/147‑1 (MQ=255)
                                                    cGTGGCAGGATTTATCGAAATTGCATGAGTTGCCGGAGTAATGGTCATCGGGGTATCTCCTTTATGAGTCATGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCACGAATCGAGAATGc                                                                                                  >  6:375766/1‑147 (MQ=255)
                                                    cGTGGCAGGATTTATCGAAATTGCATGAGTTGCCGGAGTAATGGTCATCGGGGTATCTCCTTTATGAGTCATGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCACGAATCGAGAATGCt                                                                                                 <  7:357983/148‑1 (MQ=255)
                                                     gTGGCAGTATTTATCGAAATTGCATGAGTTGCCGGAGTAATGGTCATCGGGGTATCTCCTTTATGAGTCATGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCACGAATCGAGAATGCTa                                                                                                <  4:97307/148‑1 (MQ=255)
                                                     gTGGCAGGATTTATCGAAATTGCATGAGTTGCCGGAGTAATGGTCATCGGGGTATCTCCTTTATGAGTCATGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCACGAATCGAGAATGCTa                                                                                                <  8:184924/148‑1 (MQ=255)
                                                       ggCAGGATTTATCGAAATTGCATGAGTTGCCGGAGTAATGGTCATCGGGGTATCTCCTTTATGAGTCATGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCACGAATCGAGAATGCTATg                                                                                              >  8:322978/1‑148 (MQ=255)
                                                        gCAGGATTTATCGAAATTGCATGAGTTGCCGGAGTAATGGTCATCGGGGTATCTCCTTTATGAGTCATGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCACGAATCGAGAATGCTATgg                                                                                             <  5:275422/148‑1 (MQ=255)
                                                           ggATTTATCGAAATTGCATGAGTTGCCGGAGTAATGGTCATCGGGGTATCTCCTTTATGAGTCATGGTATGAAGATACGCAGATTTACTCTTGCGTTAAAATGAATAATATTAAGCCACTTATTCACGAATCGAGAATGCTATGGATc                                                                                          <  2:233950/148‑1 (MQ=255)
                                                             aTTTATCGAAATTGCATGAGTTGCCGGAGTAATGGTCATCGGGGTATCTCCTTTATGAGTCATTGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCACGAATCGAGAATGCTATGGATCTg                                                                                        <  4:149682/148‑1 (MQ=255)
                                                                          gCATGAGTTGCCGGAGTAATGGTCATCGGGGTATCTCCTTTATGAGTCATGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCACGAATCGAGAATGCTATGGATCTGACCCAACTGgag                                                                            <  8:368916/147‑1 (MQ=255)
                                                                           cATGAGTTGCCGGAGTAATGGTCATCGGGGTATCTCCTTTATGAGTCATGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCACGAATCGAGAATGCTATGGATCTGACCCAACTGGAGAt                                                                          >  1:335165/1‑148 (MQ=255)
                                                                            aTGAGTTGCCGGAGTAATGGTCATCGGGGTATCTCCTTTATGAGTCATGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCACGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATg                                                                         >  1:224764/1‑148 (MQ=255)
                                                                             tGAGTTGCCGGAGTAATGGTCATCGGGGTATCTCCTTTATGAGTCATGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCACGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGt                                                                        <  4:360418/148‑1 (MQ=255)
                                                                             tGAGTTGCCGGAGTAATGGTCATCGGGGTATCTCCTTTATGAGTCATGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCACGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGt                                                                        >  7:107973/1‑148 (MQ=255)
                                                                             tGAGTTGCCGGAGTAATGGTCATCGGGGTATCTCCTTTATGAGTCATGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCACGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGt                                                                        >  4:103036/1‑148 (MQ=255)
                                                                              gAGTTGCCGGAGTAATGGTCATCGGGGTATCTCCTTTATGAGTCATGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCACGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGtt                                                                       >  5:239580/1‑148 (MQ=255)
                                                                               aGTTGCCGGAGTAATGGTCATCGGGGTTTCTCCTTTATGAGTCATGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCACGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGtt                                                                       <  3:184219/147‑1 (MQ=255)
                                                                                 ttGCCGGAGTAATGGTCATCGGGGTATCTCCTTTATGAGTCATGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCACGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCaa                                                                    <  5:299490/148‑1 (MQ=255)
                                                                                 ttGCCGGAGTAATGGTCATCGGGGTATCTCCTTTATGAGTCATGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCACGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCa                                                                     <  5:251271/147‑1 (MQ=255)
                                                                                  tGCCGGAGTAATGGTCATCGGGGTATCTCCTTTATGAGTCATGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCACGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAAc                                                                   >  4:486/1‑148 (MQ=255)
                                                                                      ggAGTAATGGTCATCGGGGTATCTCCTTTATGAGTCATGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCACGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCgg                                                               <  8:145203/148‑1 (MQ=255)
                                                                                      ggAGTAATGGTCATCGGGGTATCTCCTTTATGAGTCATGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCACGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCgg                                                               <  6:110999/148‑1 (MQ=255)
                                                                                           aaTGGTCATCGGGGTATCTCCTTTATGAGTCATGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCACGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGcc                                                          >  1:336506/1‑148 (MQ=255)
                                                                                                     ggggTATCTCCTTTATGAGTCATGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCACGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAAGCTGgc                                                 <  1:221272/147‑1 (MQ=255)
                                                                                                          aTCTCCTTTATGAGTCATGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCACGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAAGCTGGCAGCAt                                            >  6:191612/1‑147 (MQ=255)
                                                                                                            ctcCTTTATGAGTCATGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATg                                                                                                                                       >  2:252421/1‑54 (MQ=255)
                                                                                                            ctcCTTTATGAGTCATGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATg                                                                                                                                       <  1:252421/54‑1 (MQ=255)
                                                                                                            ctcCTTTATGAGTCATGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCACGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAAGCTGGCAGCATAAc                                         >  8:173441/1‑148 (MQ=255)
                                                                                                                tttATGAGTCATGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCACGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAAGCTGGCAGCATAAc                                         <  3:296834/144‑1 (MQ=255)
                                                                                                                tttATGAGTCATGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCACGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAAGCTGGCAGCATAAc                                         >  4:296834/1‑144 (MQ=255)
                                                                                                                  tATGAGTCATGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCACGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAAGCTGGCAGCATAACCCAGGc                                   >  1:60639/1‑148 (MQ=255)
                                                                                                                      aGTCATGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCACGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAAGCTGGCAGCATAACCCAGGCTgc                                <  6:239580/147‑1 (MQ=255)
                                                                                                                        tCATGGTATGAAGATACGCAGATTTACCCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCACGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAAGCTGGCAGCATAACCCAGGCTgc                                >  6:133311/1‑145 (MQ=255)
                                                                                                                        tCATGGTATGAAGATACGCAGATTTACCCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCACGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAAGCTGGCAGCATAACCCAGGCTgc                                <  5:133311/145‑1 (MQ=255)
                                                                                                                         cATGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCACGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAAGCTGGCAGCATAACCCAGGCTgcagc                             >  7:312129/1‑147 (MQ=255)
                                                                                                                         cATGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCACGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAAGCTGGAAGCATAACCCAGGCTgcagc                             >  3:188067/1‑147 (MQ=255)
                                                                                                                           tGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCACGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAAGCTGGCAGCATAACCCAGGCTGCAGCaaa                          <  2:60639/148‑1 (MQ=255)
                                                                                                                            tgtttGAAGATACGCAGATTTACTCTTGCTTTAAAATTAATAATATTAATCCACTTATTCACGAATCGAGAATGCTATGGATCTGACCCAACTGTAGATGTTCAACGCGGTTGCCGAAGCTGGCAGCATAACCCTTGCTGCAGCaaaa                         <  8:338812/144‑1 (MQ=255)
                                                                                                                            ggTATGAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCACGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAAGCTGGCAGCATAACCCAGGCTGCAGCaa                           >  8:65711/1‑146 (MQ=255)
                                                                                                                              tATGAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCACGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAAGCTGGCAGCATAACCCAGGCTGCAGCAAAAGt                       <  1:207468/148‑1 (MQ=255)
                                                                                                                              tATGAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCACGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAAGCTGGCAGCATAACCCAGGCTGCAGCAAAAGt                       <  5:375766/148‑1 (MQ=255)
                                                                                                                               aTGAAGATACTCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCACGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAATCTGTCAGCATAACCCAGGCTGCAGCAAAAg                        >  7:338812/1‑146 (MQ=255)
                                                                                                                               aTGAAGATACGCAGATTTACTCTTGCTGTAAAATGAATAATATTAAGCCACTTATTCACGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCACCGCGGTTGCCGCAGCTGGCAGCATAACCCAGGCTGCAGCAAACGTg                      >  2:375700/1‑148 (MQ=255)
                                                                                                                                 gAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTAACGAATCGAGAATGCTATTGATCTGACCCAACTGGAGATGTTAAACGCGGTTGCCGAAGCTGGCAGCATACCCCAGGCTGCAGCAAAAGTGc                     >  5:101858/1‑147 (MQ=255)
                                                                                                                                   aGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCACGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAAGCTGGCAGCATAACCCAGGCTGCAGCAAAAGTGCATc                  >  7:293770/1‑148 (MQ=255)
                                                                                                                                    gATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCACGAATCGAGAATGCt                                                                                                 >  6:252294/1‑68 (MQ=255)
                                                                                                                                    gATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCACGAATCGAGAATGCt                                                                                                 <  5:252294/68‑1 (MQ=255)
                                                                                                                                    gATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCACGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCATCGCGGTTGCCGAAGCTGGCAGCATAACCCAGGCTGCAGCAAAAGTGCATcg                 >  6:326013/1‑148 (MQ=255)
                                                                                                                                       aCGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCACGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAAGCTGGCAGCATAACCCAGGCTGCAGCAAAAGTGCATCGCGt              >  6:193229/1‑148 (MQ=255)
                                                                                                                                          cAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCCCTTATTCACGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAAGCTGGCAGCATAACCCAGGCTGCAGCAAAAGTGCATCGCGTGcc           >  5:190652/1‑148 (MQ=255)
                                                                                                                                           aGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCACGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAAGCTGGCAGCATAACCCAGGCTGCAGCAAAAGTGCATCGCGTGCCg          >  7:237196/1‑148 (MQ=255)
                                                                                                                                              tttAGTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCACGAATCgaga                                                                                                      >  5:264604/1‑53 (MQ=255)
                                                                                                                                              tttACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCACGAATCgaga                                                                                                      <  6:264604/53‑1 (MQ=255)
                                                                                                                                              tttACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCACGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAAGCTGGCAGCATAACCCAGGCTGCAGCAAAAGTGCATCGCGTGCCGTc        >  4:172974/1‑147 (MQ=255)
                                                                                                                                              tttACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCACGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAAGCTGGCAGCATAACCCAGGCTGCAGCAAAAGTGCATCGCGTGCCGTCg       <  7:322978/148‑1 (MQ=255)
                                                                                                                                               ttACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCACGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAAGCTGGCAGCATAACCCAGGCTGCAGCAAAAGTGCATCGCGTGCCGTCg       >  4:110250/1‑147 (MQ=255)
                                                                                                                                               ttACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCACGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAAGCTGGCAGCATAACCCAGGCTGCAGCAAAAGTGCATCGCGTGCCGTCg       >  2:46325/1‑147 (MQ=255)
                                                                                                                                                  ctctTGCTTTAAAATGAATAATATTAAGCCACTTATTCACGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAAGCTGGCAGCATAACCCAGGCTGCAGCAAAAGTGCATCGCGTGCCGTCGAAtc   <  5:176083/148‑1 (MQ=255)
                                                                                                                                                  ctctTGCTTTAAAATGAATAATATTAAGCCACTTATTCACGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAAGCTGGCAGCATAACCCAGGCTGCAGCAAAAGTGCATCGCGTGCCGTCGAAtc   <  1:237642/148‑1 (MQ=255)
                                                                                                                                                  ctctTGCTTTAAAATGAATAATATTAAGCCACTTATTCACGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAAGCTGGCAGCATAACCCAGGCTGCAGCAAAAGTGCATCGCGTGCCGTCGAAt    <  2:335165/147‑1 (MQ=255)
                                                                                                                                                   tctTGCTTTAAAATGAATAATATTAAGCCACTTATTCACGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAAGCTGGCAGCATAACCCAGGCTGCAGCAAAAGTGCATCGCGTGCCGTCGAAtct  >  3:192702/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CGTTTTCGATATCGTCAGCGCCAGCCCACGGCAGCACAGAAAGTTGTTCACCCGTGGCAGGATTTATCGAAATTGCATGAGTTGCCGGAGTAATGGTCATCGGGGTATCTCCTTTATGAGTCATGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCACGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAAGCTGGCAGCATAACCCAGGCTGCAGCAAAAGTGCATCGCGTGCCGTCGAATCT  >  NC_000913/1614604‑1614898

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