Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,238,465 G→T 7.2% intergenic (+168/‑115) ygjR → / → alx putative NAD(P)‑dependent dehydrogenase/putative membrane‑bound redox modulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,238,4650GT 7.2% 340.6 / 13.8 97intergenic (+168/‑115)ygjR/alxputative NAD(P)‑dependent dehydrogenase/putative membrane‑bound redox modulator
Reads supporting (aligned to +/- strand):  ref base G (47/43);  new base T (3/4);  total (50/47)
Fisher's exact test for biased strand distribution p-value = 7.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.61e-01

AAAGGTATGTAACAGACCATTGACTGGCTGAATGGTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTCTATGAATACTGTCGGCACGCCGTTGCTAT  >  NC_000913/3238319‑3238607
                                                                                                                                                  |                                                                                                                                              
aaaGGTATGTAACAGACCATTGACTGGCTGAATGGTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGAc                                                                                                                                              <  2:221331/149‑1 (MQ=255)
       tGTAACAGACCATTGACTGGCTGAATGGTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTcc                                                                                                                                       >  4:667648/1‑149 (MQ=255)
       tGTAACAGACCATTGACTGGCTGAATGGTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTcc                                                                                                                                       <  1:290844/149‑1 (MQ=255)
        gTAACAGACCATTGACTGGCTGAATGGTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCt                                                                                                                                      >  3:500402/1‑149 (MQ=255)
        gTAACAGACCATTGACTGGCTGAATGGTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCt                                                                                                                                      <  1:426593/149‑1 (MQ=255)
         tAACAGACCATTGACTGGCTGAATGGTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTc                                                                                                                                        >  3:23383/1‑146 (MQ=255)
         tAACAGACCATTGACTGGCTGAATGGTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTc                                                                                                                                        <  4:23383/146‑1 (MQ=255)
          aaCAGACCATTGACTGGCTGAATGGTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGTCTGCGTATTCCttt                                                                                                                                    <  1:238874/149‑1 (MQ=255)
          aaCAGACCATTGACTGGCTGAATGGTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGTCTGCGTATTCCttt                                                                                                                                    <  3:5281/149‑1 (MQ=255)
           aCAGACCATTGACTGGCTGAATGGTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCCAAAGGAATACGCAGACGTATTCCtttt                                                                                                                                   >  2:676225/1‑149 (MQ=255)
            cAGACCATTGACTGGCTGAATGGTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCttttt                                                                                                                                  >  1:582331/1‑149 (MQ=255)
             aGACCATTGACTGGCTGAATGGTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCtttttt                                                                                                                                 <  1:201002/149‑1 (MQ=255)
             aGACCATTGACTGGCTGAATGGTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCCAAAGGAATACGCAGACGTATTCCtttttt                                                                                                                                 <  1:676225/149‑1 (MQ=255)
                 cATTGACTGGCTGAATGGTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTttgttg                                                                                                                             >  4:671049/1‑149 (MQ=255)
                  aTTGACTGGCTGAATGGTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTcc                                                                                                                                       >  3:299569/1‑138 (MQ=255)
                  aTTGACTGGCTGAATGGTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTcc                                                                                                                                       <  4:299569/138‑1 (MQ=255)
                  aTTGACTGGCTGAATGGTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTtgttgt                                                                                                                            <  1:310064/149‑1 (MQ=255)
                    tGACTGGCTGAATGGTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTaa                                                                                                                          <  2:351649/149‑1 (MQ=255)
                    tGACTGGCTGAATGGTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTaa                                                                                                                          <  3:667648/149‑1 (MQ=255)
                    tGACTGGCTGAATGGTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTaa                                                                                                                          <  1:384847/149‑1 (MQ=255)
                    tGACTGGCTGAATGGTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTaa                                                                                                                          <  3:383809/149‑1 (MQ=255)
                       cTGGCTGAATGGTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTg                                                                                                                       >  1:285259/1‑149 (MQ=255)
                       cTGGCTGAATGGTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTg                                                                                                                       >  4:625881/1‑149 (MQ=255)
                        tGGCTGAATGGTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTga                                                                                                                      <  2:27867/149‑1 (MQ=255)
                            tGAATGGTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGAcc                                                                                                                  >  2:611481/1‑149 (MQ=255)
                             gAATGGTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTttgtt                                                                                                                              >  4:685241/1‑136 (MQ=255)
                             gAATGGTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTttgtt                                                                                                                              <  3:685241/136‑1 (MQ=255)
                                  gTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCttt                                                                                                                                    >  1:662300/1‑125 (MQ=255)
                                  gTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCttt                                                                                                                                    <  2:662300/125‑1 (MQ=255)
                                      gTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAg                                                                                                        >  3:466288/1‑149 (MQ=255)
                                      gTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAg                                                                                                        <  2:285259/149‑1 (MQ=255)
                                      gTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAg                                                                                                        <  1:120117/149‑1 (MQ=255)
                                           aCTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAg                                                                                                   <  4:500402/149‑1 (MQ=255)
                                           aCTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAg                                                                                                   <  1:313367/149‑1 (MQ=255)
                                             tttGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGt                                                                                                 >  2:405424/1‑149 (MQ=255)
                                             tttGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGt                                                                                                 >  1:658788/1‑149 (MQ=255)
                                                 ttACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTAt                                                                                             <  3:625881/149‑1 (MQ=255)
                                                 ttACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTAt                                                                                             >  4:546844/1‑149 (MQ=255)
                                                 ttACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTAt                                                                                             >  1:442314/1‑149 (MQ=255)
                                                 ttACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGAGCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTAt                                                                                             <  2:446501/149‑1 (MQ=255)
                                                   aCCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCGGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGc                                                                                           <  2:120861/149‑1 (MQ=255)
                                                     cTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGa                                                                                                    >  1:273426/1‑138 (MQ=255)
                                                     cTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGa                                                                                                    <  2:273426/138‑1 (MQ=255)
                                                     cTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCAt                                                                                         >  2:205538/1‑149 (MQ=255)
                                                     cTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCAt                                                                                         >  4:124495/1‑149 (MQ=255)
                                                          aaGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATaaaaa                                                                                    >  2:486964/1‑149 (MQ=255)
                                                              ggAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAg                                                                                                        >  3:67904/1‑125 (MQ=255)
                                                              ggAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAg                                                                                                        <  4:67904/125‑1 (MQ=255)
                                                               gAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAagcagc                                                                               >  1:524033/1‑149 (MQ=255)
                                                                 gTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCgg                                                                             >  3:151765/1‑149 (MQ=255)
                                                                  tAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGc                                                                            >  3:310224/1‑149 (MQ=255)
                                                                   aaCTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCt                                                                           >  4:599344/1‑149 (MQ=255)
                                                                    aCTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTg                                                                          >  2:384453/1‑149 (MQ=255)
                                                                      ttCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGAc                                                                        >  2:317590/1‑149 (MQ=255)
                                                                      ttCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGAc                                                                        >  3:572387/1‑149 (MQ=255)
                                                                      ttCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGAc                                                                        >  1:399386/1‑149 (MQ=255)
                                                                      ttCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCCATGCATAAAAAGCAGCGGCTGAc                                                                        >  4:212340/1‑149 (MQ=255)
                                                                             ccGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTcc                                                                 <  4:310224/149‑1 (MQ=255)
                                                                                     tCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGTCTGCGTATTCCtttcgg                                                                                                                                 >  2:595114/1‑74 (MQ=25)
                                                                                     tCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGTCTGCGTATTCCtttcgg                                                                                                                                 <  1:595114/77‑4 (MQ=25)
                                                                                     tCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTgg                                                         >  1:243928/1‑149 (MQ=255)
                                                                                        tCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAagcagc                                                                               <  3:441880/124‑1 (MQ=255)
                                                                                        tCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAagcagc                                                                               >  4:441880/1‑124 (MQ=255)
                                                                                        tCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCg                                                      >  4:300143/1‑149 (MQ=255)
                                                                                              cGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGtttttt                                                <  1:405424/149‑1 (MQ=255)
                                                                                               gATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGCCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGttttttt                                               >  2:570009/1‑149 (MQ=255)
                                                                                                 tgtgtgAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTAt                                             <  4:466288/149‑1 (MQ=255)
                                                                                                   tgtgAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTAtgt                                           <  1:486964/149‑1 (MQ=255)
                                                                                                    gtgAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTAtgag                                          >  3:606343/1‑147 (MQ=255)
                                                                                                     tgAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATgtgt                                         >  1:7352/1‑149 (MQ=255)
                                                                                                           aaCACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGaa                                   <  3:212340/149‑1 (MQ=255)
                                                                                                           aaCACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGaa                                   <  3:673136/149‑1 (MQ=255)
                                                                                                            acacATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAAc                                  >  2:496180/1‑149 (MQ=255)
                                                                                                               caTCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGAGTAAGGAACTTc                               <  4:606343/149‑1 (MQ=255)
                                                                                                                aTCCGGTCACCGGGCAACCCGAAAGGAATACGCCGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGACGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAAGAACTTCt                              >  2:244217/1‑149 (MQ=255)
                                                                                                                aTCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTCt                              >  2:665779/1‑149 (MQ=255)
                                                                                                                aTCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTCt                              >  1:663195/1‑149 (MQ=255)
                                                                                                                 tCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGtttttt                                                <  4:466824/130‑1 (MQ=255)
                                                                                                                 tCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGtttttt                                                >  3:466824/1‑130 (MQ=255)
                                                                                                                  ccGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTCTAt                            <  1:570009/149‑1 (MQ=255)
                                                                                                                  ccGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTCTAt                            <  1:205538/149‑1 (MQ=255)
                                                                                                                  ccGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTCTAt                            <  1:189746/149‑1 (MQ=255)
                                                                                                                    ggTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTCTATGa                          <  1:611481/149‑1 (MQ=255)
                                                                                                                     gTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTCTATGaa                         <  3:300143/149‑1 (MQ=255)
                                                                                                                        aCCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTCTATGAATAc                      >  2:220234/1‑149 (MQ=255)
                                                                                                                               aaCCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTa                                                                                              <  4:124885/70‑1 (MQ=255)
                                                                                                                               aaCCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTa                                                                                              >  3:124885/1‑70 (MQ=255)
                                                                                                                               aaCCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTCTATGAATACTGTCGGc               >  1:672496/1‑149 (MQ=255)
                                                                                                                                  aaaaaaGGAATACGTCTGCGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTAt                                             <  2:231204/113‑1 (MQ=39)
                                                                                                                                  aaaaaaGGAATACGTCTGCGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTAt                                             >  1:231204/4‑116 (MQ=39)
                                                                                                                                   cGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTCTATGAATACTGTCGGCACGc           >  1:31114/1‑149 (MQ=255)
                                                                                                                                   cGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTCTATGAATACTGTCGGCACGc           <  2:7352/149‑1 (MQ=255)
                                                                                                                                      aaGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTCTATGAATACTGTCGGCACGCCGt        >  4:255878/1‑149 (MQ=255)
                                                                                                                                       aGGAATACGTCTGCGTATTCCTTTCGGGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTCTATGAATACTGTCGGCACGCCGtt       >  3:266164/1‑149 (MQ=255)
                                                                                                                                       aGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTCTATGAATACTGTCGGCACGCCGtt       <  3:599344/149‑1 (MQ=255)
                                                                                                                                            tACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTCTATGAATACTGTCGGCACGCCGTTGCTAt  <  1:496180/149‑1 (MQ=255)
                                                                                                                                            tACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTCTATGAATACTGTCGGCACGCCGTTGCTAt  <  3:671049/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                              
AAAGGTATGTAACAGACCATTGACTGGCTGAATGGTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTCTATGAATACTGTCGGCACGCCGTTGCTAT  >  NC_000913/3238319‑3238607

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