Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,238,463 C→T 7.1% intergenic (+166/‑117) 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,4630CT 7.1% 310.4 / 11.1 99intergenic (+166/‑117)ygjR/alxputative NAD(P)‑dependent dehydrogenase/putative membrane‑bound redox modulator
Reads supporting (aligned to +/- strand):  ref base C (47/45);  new base T (3/4);  total (50/49)
Fisher's exact test for biased strand distribution p-value = 7.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.08e-01

GTAAAAGGTATGTAACAGACCATTGACTGGCTGAATGGTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTCTATGAATACTGTCGGCACGCCGTTGCTAT  >  NC_000913/3238316‑3238607
                                                                                                                                                   |                                                                                                                                                
gtAAAAGGTATGTAACAGACCATTGACTGGCTGAATGGTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCa                                                                                                                                                 <  4:322549/149‑1 (MQ=255)
gtAAAAGGTATGTAACAGACCATTGACTGGCTGAATGGTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCa                                                                                                                                                 <  1:297972/149‑1 (MQ=255)
   aaaGGTATGTAACAGACCATTGACTGGCTGAATGGTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGAc                                                                                                                                              <  2:221331/149‑1 (MQ=255)
          tGTAACAGACCATTGACTGGCTGAATGGTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTcc                                                                                                                                       <  1:290844/149‑1 (MQ=255)
          tGTAACAGACCATTGACTGGCTGAATGGTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTcc                                                                                                                                       >  4:667648/1‑149 (MQ=255)
           gTAACAGACCATTGACTGGCTGAATGGTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCt                                                                                                                                      <  1:426593/149‑1 (MQ=255)
           gTAACAGACCATTGACTGGCTGAATGGTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCt                                                                                                                                      >  3:500402/1‑149 (MQ=255)
            tAACAGACCATTGACTGGCTGAATGGTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTc                                                                                                                                        <  4:23383/146‑1 (MQ=255)
            tAACAGACCATTGACTGGCTGAATGGTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTc                                                                                                                                        >  3:23383/1‑146 (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                                                                                                                                       <  4:299569/138‑1 (MQ=255)
                     aTTGACTGGCTGAATGGTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTcc                                                                                                                                       >  3:299569/1‑138 (MQ=255)
                     aTTGACTGGCTGAATGGTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTtgttgt                                                                                                                            <  1:310064/149‑1 (MQ=255)
                       tGACTGGCTGAATGGTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTaa                                                                                                                          <  3:667648/149‑1 (MQ=255)
                       tGACTGGCTGAATGGTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTaa                                                                                                                          <  3:383809/149‑1 (MQ=255)
                       tGACTGGCTGAATGGTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTaa                                                                                                                          <  2:351649/149‑1 (MQ=255)
                       tGACTGGCTGAATGGTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTaa                                                                                                                          <  1:384847/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                                                                                                                              <  3:685241/136‑1 (MQ=255)
                                gAATGGTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTttgtt                                                                                                                              >  4:685241/1‑136 (MQ=255)
                                     gTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCttt                                                                                                                                    >  1:662300/1‑125 (MQ=255)
                                     gTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCttt                                                                                                                                    <  2:662300/125‑1 (MQ=255)
                                         gTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAg                                                                                                        <  1:120117/149‑1 (MQ=255)
                                         gTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAg                                                                                                        >  3:466288/1‑149 (MQ=255)
                                         gTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAg                                                                                                        <  2:285259/149‑1 (MQ=255)
                                              aCTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAg                                                                                                   <  1:313367/149‑1 (MQ=255)
                                              aCTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAg                                                                                                   <  4:500402/149‑1 (MQ=255)
                                                tttGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGt                                                                                                 >  1:658788/1‑149 (MQ=255)
                                                tttGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGt                                                                                                 >  2:405424/1‑149 (MQ=255)
                                                    ttACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTAt                                                                                             >  4:546844/1‑149 (MQ=255)
                                                    ttACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTAt                                                                                             <  3:625881/149‑1 (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                                                                                         >  4:124495/1‑149 (MQ=255)
                                                        cTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCAt                                                                                         >  2:205538/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                                                                        >  3:572387/1‑149 (MQ=255)
                                                                         ttCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGAc                                                                        >  1:399386/1‑149 (MQ=255)
                                                                         ttCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGAc                                                                        >  2:317590/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:205538/149‑1 (MQ=255)
                                                                                                                     ccGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTCTAt                            <  1:570009/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)
                                                                                                                                                   |                                                                                                                                                
GTAAAAGGTATGTAACAGACCATTGACTGGCTGAATGGTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTCTATGAATACTGTCGGCACGCCGTTGCTAT  >  NC_000913/3238316‑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: