Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,238,464 A→C 5.9% intergenic (+167/‑116) 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,4640AC 5.9% 228.8 / 4.4 85intergenic (+167/‑116)ygjR/alxputative NAD(P)‑dependent dehydrogenase/putative membrane‑bound redox modulator
Reads supporting (aligned to +/- strand):  ref base A (45/32);  new base C (2/3);  total (47/38)
Fisher's exact test for biased strand distribution p-value = 6.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.65e-01

GTAAAAGGTATGTAACAGACCATTGACTGGCTGAATGGTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTCTATGAATACTGTCGGCACGCCGTTGCTATGGGGC  >  NC_000913/3238316‑3238612
                                                                                                                                                    |                                                                                                                                                    
gtAAAAGGTATGTAACAGACCATTGACTGGCTGAATGGTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAg                                                                                                                                                     >  2:88419/1‑150 (MQ=255)
      ggTATGTAACAGACCATTGACTGGCTGAATGGTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAACGGAATCCGCAGACGTAt                                                                                                                                               >  1:1010379/1‑150 (MQ=255)
      ggTATGTAACAGACCATTGACTGGCTGAATGGTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTAt                                                                                                                                               >  2:214046/1‑150 (MQ=255)
       gTATGTAACAGACCATTGACTGGCTGAATGGTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTAtt                                                                                                                                              <  2:1139924/150‑1 (MQ=255)
          tGTAACAGACCATTGACTGGCTGAATGGTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCt                                                                                                                                           <  2:1294147/150‑1 (MQ=255)
            tAACAGACCATTGACTGGCTGAATGGTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCttt                                                                                                                                         >  1:1098215/1‑150 (MQ=255)
             aaCAGACCATTGACTGGCTGAATGGTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCtttt                                                                                                                                        <  2:553221/150‑1 (MQ=255)
             aaCAGACCATTGACTGGCTGAATGGTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGAAACCCGAAAGGAATACGCAGACGTATTCCtttt                                                                                                                                        >  1:747806/1‑150 (MQ=255)
              aCAGACCATTGACTGGCTGAATGGTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGTCTGCGTATTCCtttcg                                                                                                                                       <  2:748385/150‑3 (MQ=255)
              aCAGACCATTGACTGGCTGAATGGTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGTCTGCGTATTCCtttcg                                                                                                                                       <  1:941206/150‑3 (MQ=255)
              aCAGACCATTGACTGGCTGAATGGTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCttttt                                                                                                                                       <  2:1054526/150‑1 (MQ=255)
                  gaccttGACTGGCTGAATGGTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATTTGAGAAAAAACACATCAGGTCACAGTGCAACCCGAAACGAAAACGCAGACGAAATCCTTTTttgtt                                                                                                                                   >  2:263827/5‑150 (MQ=255)
                  aCCATTGACTGGCTGAATGGTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTttgtt                                                                                                                                   <  2:564044/150‑1 (MQ=255)
                   ccATTGACTGGCTGAATGGTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTttgttg                                                                                                                                  >  2:85264/1‑150 (MQ=255)
                    cATTGACTGGCTGAATGGTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTATTtgttgt                                                                                                                                 >  2:856966/1‑150 (MQ=255)
                     aTTGACTGGCTGAATGGTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGGCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTa                                                                                                                                <  2:677267/150‑1 (MQ=255)
                      ttGACTGGCTGAATGGTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTaa                                                                                                                               <  2:392188/150‑1 (MQ=255)
                       tGACTGGCTGAATGGTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAg                                                                                                                              <  2:841291/150‑1 (MQ=255)
                            ggCTGAATGGTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTgaga                                                                                                                         >  2:1286698/1‑150 (MQ=255)
                              cTGAATGGTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGAcc                                                                                                                       >  1:1036738/1‑150 (MQ=255)
                               tGAATGGTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCt                                                                                                                      <  1:654511/150‑1 (MQ=255)
                                     gtatGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATAGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGCGCAACCCTAAAGGAATACGCAGACGTAATCCTGTTATGTTCTAAGTGAGACATTGCCgg                                                                                                                >  1:1139693/4‑150 (MQ=255)
                                     gTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCgg                                                                                                                <  1:317805/150‑1 (MQ=255)
                                     gTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCAGAAAGGAATACGCAGACGTATTCCTTTTTTCTTGTAAGTGAGACATTGACgt                                                                                                                >  1:515923/1‑149 (MQ=255)
                                      tCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTTTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGa                                                                                                               <  1:966738/150‑1 (MQ=255)
                                        tGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCACCCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGCCCTTGCCGGAAg                                                                                                             >  2:129785/1‑150 (MQ=255)
                                         gTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAgg                                                                                                            >  1:74510/1‑150 (MQ=255)
                                           cATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCg                                                                                                          <  1:732895/150‑1 (MQ=255)
                                               cTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGCGACCTTGCCGGAAGGCGAGGt                                                                                                      >  2:797250/1‑150 (MQ=255)
                                                tttGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAATGCGAGGTc                                                                                                     >  2:843012/1‑150 (MQ=255)
                                                tttGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTc                                                                                                     <  2:713410/150‑1 (MQ=255)
                                                  tGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTa                                                                                                   <  2:571390/150‑1 (MQ=255)
                                                  tGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTa                                                                                                   <  1:755896/150‑1 (MQ=255)
                                                     tACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGGCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTCTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGc                                                                                                >  1:75865/1‑150 (MQ=255)
                                                            aaaGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGACAAAACACATCCGGTCACCGGGCAACCCGCAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGCAGGCGAGGTCTATGCATaaaaa                                                                                         <  1:908213/150‑1 (MQ=255)
                                                             aaGGGGAGTCACTTCATTGCCGCTCGATCCTCATTACGATGTTTGAAAAAACACATGCGGTCACCGGGCAACCCGTAAGGAATACGCTGACGTATTCCTTTTTTGTTGTAAGTGACACCTTGCCGGATGGCGAGGTCTATGCATAAAAag                                                                                        <  1:211531/150‑1 (MQ=255)
                                                             aaGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAag                                                                                        <  2:281098/150‑1 (MQ=255)
                                                              aGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAagc                                                                                       <  1:223055/150‑1 (MQ=255)
                                                              aGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAagc                                                                                       >  1:302729/1‑150 (MQ=255)
                                                                 ggAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAagcagc                                                                                    >  1:986615/1‑150 (MQ=255)
                                                                 ggAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAagcagc                                                                                    >  1:522868/1‑150 (MQ=255)
                                                                 ggAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAagcagc                                                                                    <  1:531891/150‑1 (MQ=255)
                                                                 ggAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAagcagc                                                                                    <  2:766095/150‑1 (MQ=255)
                                                                 ggAGTAACTTCAGTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAagcagc                                                                                    <  1:582464/150‑1 (MQ=255)
                                                                   aGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCgg                                                                                  <  2:1054012/150‑1 (MQ=255)
                                                                       aCTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGa                                                                              <  2:873601/150‑1 (MQ=255)
                                                                          tCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGt                                                                           >  1:1111953/1‑150 (MQ=255)
                                                                           cATTGCCTGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTc                                                                          <  2:1154088/150‑1 (MQ=255)
                                                                            aTTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACTTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCt                                                                         >  1:180720/1‑150 (MQ=255)
                                                                               gCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTcc                                                                      <  1:846101/150‑1 (MQ=255)
                                                                               gCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTcc                                                                      >  1:1155716/1‑150 (MQ=255)
                                                                                  ggTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCTAc                                                                   >  1:898743/1‑150 (MQ=255)
                                                                                  ggTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGAc                                                                   >  2:218107/1‑150 (MQ=255)
                                                                                     cGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGtt                                                                >  1:143661/1‑150 (MQ=255)
                                                                                        tCGTCATTACGATGCGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGc                                                             <  2:963497/150‑1 (MQ=255)
                                                                                         cGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGcc                                                            <  1:684241/150‑1 (MQ=255)
                                                                                          gTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCg                                                           >  2:1112682/1‑150 (MQ=255)
                                                                                          gTCATTACGATGTGGGAAAAAACCCATCCGGGCAACGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCg                                                           <  1:204741/150‑1 (MQ=255)
                                                                                             aTTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGttt                                                        >  2:318889/1‑150 (MQ=255)
                                                                                                aCGATGTGCGAAAAAACACATCCGGACACCGGGCAACCCGAAAGGAATACTCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGTAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGtttttt                                                     <  1:472490/150‑1 (MQ=255)
                                                                                                      tgtgAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTAtgtg                                               >  2:727643/1‑150 (MQ=255)
                                                                                                       gtgAAACAACACATCCGGGCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATgtgt                                              <  1:851799/150‑1 (MQ=255)
                                                                                                       gtgAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGttttttt                                                    >  1:457167/1‑144 (MQ=255)
                                                                                                       gtgAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGttttttt                                                    <  2:457167/144‑1 (MQ=255)
                                                                                                       gtgAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATgtgt                                              >  1:567330/1‑150 (MQ=255)
                                                                                                       gtgAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGACGGAAGGCGACGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATgtgt                                              >  1:425297/1‑150 (MQ=255)
                                                                                                        tgAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTa                                             >  2:1079350/1‑150 (MQ=255)
                                                                                                          aaaaaaCACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGGCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAg                                           >  1:748934/1‑150 (MQ=255)
                                                                                                          aaaaaaCACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAACGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCAGACGTTGGCCGTTTTATTATGTGTatg                                           >  1:219698/1‑148 (MQ=255)
                                                                                                                 acaTCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTc                                    >  1:117609/1‑150 (MQ=255)
                                                                                                                    tCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTCTAt                                 <  1:924553/150‑1 (MQ=255)
                                                                                                                    tCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTCTAt                                 >  1:988345/1‑150 (MQ=255)
                                                                                                                       ggTCACCGGGCAACCCGAACGGAATACGCTGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTCTATGaa                              <  2:1065394/150‑1 (MQ=255)
                                                                                                                         tCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTCTATGAATa                            >  1:900071/1‑150 (MQ=255)
                                                                                                                                gCAACCCGAAAGGAATACGACGACGTATTACTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTCTATGAATACTGTCgg                     <  1:1287529/150‑1 (MQ=255)
                                                                                                                                  aaCCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTCTATGAATACTGTCGGCa                   >  2:768254/1‑150 (MQ=255)
                                                                                                                                   aCCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTCTATGAATACTGTCGGCAc                  >  2:333659/1‑150 (MQ=255)
                                                                                                                                     aaaaaaGGAATACGTCTGCGTATTCCTTTTTTGTTGTAAGTGAGACATTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGTCTGACGTCTTCCGCCGTTGGCCGTTTTTTTATGTGTAAGGAACTTCTATGAATACTGTCGGCACGc                >  2:1113591/4‑150 (MQ=255)
                                                                                                                                        aaaGGAATACGTCTGCGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTCTATGAATACTGTCGGCACGCCGt             >  1:1299716/1‑150 (MQ=255)
                                                                                                                                         aaGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTCTATGAATACTGTCGGCACGCCGtt            >  1:883646/1‑150 (MQ=255)
                                                                                                                                         aaGGAATACGCAGACGTAGTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTCTATGAATACTGTCGGCACGCCGtt            >  1:79659/1‑150 (MQ=255)
                                                                                                                                           tgAATACCCTTACGTATTCCTTTATTGTTCTAAGTGAGACCTTGCAGGAAGGCGAGGTCTATGCATAAAAAGCAGCCGCTGACGACTACCGTCGTTGGCCGTTTTTTTATGTGTAAGGAACTGCTATGATTACTCTCGGCACGCCGTTGc          <  2:658345/149‑1 (MQ=255)
                                                                                                                                             aaTACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCAGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTCTATGAATACTGTCGGCACGCCGTTGCTa        >  2:1181741/1‑150 (MQ=255)
                                                                                                                                                   cAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTCTATGAATACTGTCGGCACGCCGTTGCTATGGGGc  >  1:1120921/1‑150 (MQ=255)
                                                                                                                                                   cAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTCTATGAATACTGTCGGCACGCCGTTGCTATGGGGc  >  1:1115735/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GTAAAAGGTATGTAACAGACCATTGACTGGCTGAATGGTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTCTATGAATACTGTCGGCACGCCGTTGCTATGGGGC  >  NC_000913/3238316‑3238612

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