Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,579,564 G→A 6.6% intergenic (‑222/+69) yddA ← / ← ydeM putative multidrug ABC transporter permease/ATPase/putative YdeN‑specific sulfatase‑maturating enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,579,5640GA 6.6% 313.3 / 8.6 107intergenic (‑222/+69)yddA/ydeMputative multidrug ABC transporter permease/ATPase/putative YdeN‑specific sulfatase‑maturating enzyme
Reads supporting (aligned to +/- strand):  ref base G (51/48);  new base A (4/3);  total (55/52)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.18e-01

Only 100 of 108 total aligned reads displayed.

TAACAATATGAAATATAAAGGAATATTTCATTTTCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCGGTACTCTGTTCTTAGCTAACTCTACCATGGCGTT  >  NC_000913/1579418‑1579713
                                                                                                                                                  |                                                                                                                                                     
tAACAATATGAAATATAAAGGAATATTTCATTTTCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAAt                                                                                                                                                    >  1:785481/1‑150 (MQ=255)
     atatGAAATATAAATGAAAATTTCATTTTCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCAAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATg                                                                                                                                               <  1:345494/150‑1 (MQ=255)
     atatGAAATATAAAGGAATATTTCATTTTCAGTGGTTATGGTTATTCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATg                                                                                                                                               <  2:1260321/150‑1 (MQ=255)
          aaaTATAAAGGAATATTTCATTTTCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCggg                                                                                                                                          <  1:1122984/150‑1 (MQ=255)
            atatAAAGGAATATTTCATTTTCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACATAAATTCCTGTCGGGTg                                                                                                                                        <  1:1581928/150‑1 (MQ=255)
            atatAAAGGAATATTTCATTTTCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACATAAATTCCTGTCGGGTg                                                                                                                                        <  1:539077/150‑1 (MQ=255)
             tataAAGGAATATTTCATTTTCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTGATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACATAAATTCCTGTCGGGTgt                                                                                                                                       <  2:1766400/150‑2 (MQ=255)
               taAAGGAATATTTCATTTTCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCATGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCt                                                                                                                                     >  1:1862472/1‑150 (MQ=255)
               taAAGGAATATTTCATTTTCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCt                                                                                                                                     >  2:1281694/1‑150 (MQ=255)
                   ggAATATTTCATTTTCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGtt                                                                                                                                 >  2:1680698/1‑150 (MQ=255)
                     aaTATTTCATTTTCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTaa                                                                                                                               <  1:631497/150‑1 (MQ=255)
                      atctttCATTTTCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACTTAATCCAGGATTTTCGCAAGGGAAGATGTTGATGACGAAACACCCGACAGGAATTTATGTCGGGAGCCTTGTTAAg                                                                                                                              <  2:227983/147‑1 (MQ=255)
                        atTTCATTTTCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGt                                                                                                                            <  1:1704460/150‑1 (MQ=255)
                        atTTCATTTTCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACAACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTTCCTTGTTAAGGt                                                                                                                            >  1:1762333/1‑150 (MQ=255)
                         tttCATTTTCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTc                                                                                                                           <  2:1787895/150‑1 (MQ=255)
                         tttCATTTTCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTc                                                                                                                           >  2:820247/1‑150 (MQ=255)
                         tttCATTTTCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTACGTCGGGTGCCTTGTTAAGGTc                                                                                                                           <  1:542730/150‑1 (MQ=255)
                            cATTTTCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATa                                                                                                                        >  2:60708/1‑150 (MQ=255)
                            cATTTTCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATa                                                                                                                        >  2:881640/1‑150 (MQ=255)
                            aaTTTTCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATa                                                                                                                        <  2:931011/149‑1 (MQ=255)
                                  cAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAagg                                                                                                                  <  2:536225/150‑1 (MQ=255)
                                  cAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAagg                                                                                                                  <  2:831144/150‑1 (MQ=255)
                                      ggttatggttatGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGc                                                                                                              <  1:1785706/150‑1 (MQ=255)
                                      ggttatggttatGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGc                                                                                                              >  1:1758554/1‑150 (MQ=255)
                                      ggttatggttatGCTATTGATAATCAAGGATATTGCTAAGCCCCTGATGCATACCACGTTAGCCAGGATTTTCGCAAGGGAAGATGATTATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGGCATAAGAAGGAGGc                                                                                                              <  2:93917/150‑1 (MQ=255)
                                          atggttatgCTATTGATAATAAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGAAGAAACACCCGACAGGAATTTATGTCGGGGGCCTTGTTAAGGTCATAAGAAGGAGGCTAAg                                                                                                          >  2:154766/1‑150 (MQ=255)
                                            ggttatgCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGaa                                                                                                        >  2:269387/1‑150 (MQ=255)
                                            ggttatgCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGaa                                                                                                        <  2:1143574/150‑1 (MQ=255)
                                             gttatgCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAAt                                                                                                       >  1:762498/1‑150 (MQ=255)
                                              ttatgCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATg                                                                                                      <  2:60507/150‑1 (MQ=255)
                                                 tgCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAg                                                                                                   <  1:1004091/150‑1 (MQ=255)
                                                 tgCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAg                                                                                                   <  1:930681/150‑1 (MQ=255)
                                                 tgCTATTGATAATCAAGGATATTGCTAATGCCCTGATGAATACCACGTAAGACAGGATTTTCGCTATGGAAGATGATGATGACGATACCCCCTACAGGATTTTATTTCGGGTGCCTTGTTACGGTCATAAGATGGAGGCTAAGAATGGAg                                                                                                   <  1:1508725/150‑1 (MQ=255)
                                                  gCTATTGATAATCAAGGATATTGCTAATGACCTGATGCATACCCCGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGACTTGTTAAGGTCATAAGAAGGACGCTAAGAATGGAGt                                                                                                  >  1:1116062/1‑150 (MQ=255)
                                                   cTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGtt                                                                                                 >  1:1600108/1‑150 (MQ=255)
                                                    tATTGTTAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTa                                                                                                <  1:157733/150‑1 (MQ=255)
                                                     aTTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTaa                                                                                               <  1:1680559/150‑1 (MQ=255)
                                                     aTTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTaa                                                                                               >  1:1256783/1‑150 (MQ=255)
                                                      ttGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTaaa                                                                                              >  2:1831240/1‑150 (MQ=255)
                                                       tGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTaaaa                                                                                             >  1:1064148/1‑150 (MQ=255)
                                                        gATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAag                                                                                            <  1:367587/150‑1 (MQ=255)
                                                          tAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAagag                                                                                          >  1:1452510/1‑150 (MQ=255)
                                                            aTCAAGGATATTGCTAATTCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAagagag                                                                                        <  1:762724/150‑1 (MQ=255)
                                                            aTCAAGGATATTGCTAAAGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCAGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAagagag                                                                                        >  1:201116/1‑150 (MQ=255)
                                                              cAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCg                                                                                      >  1:1294375/1‑150 (MQ=255)
                                                                aGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGtt                                                                                    >  1:1788211/1‑150 (MQ=255)
                                                                     atTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTattaa                                                                               >  1:1803650/1‑150 (MQ=255)
                                                                       tGCTAATGACCTTATGAATAACACGTTAGACAGTATTTTTGCAAGTGAAGATTATGATGAAGAAACACCCGACATGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGATGCTAAGAATGGAGTTAAAAGAGAGCGTTattaatt                                                                             <  1:1004595/150‑1 (MQ=255)
                                                                         cTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAattat                                                                           <  1:322366/150‑1 (MQ=255)
                                                                          tAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAattatt                                                                          >  1:1800239/1‑150 (MQ=255)
                                                                          tAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAattatt                                                                          >  1:1020252/1‑150 (MQ=255)
                                                                               cccTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCa                                                                     <  2:494229/150‑1 (MQ=255)
                                                                               cccTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCa                                                                     >  1:1264238/1‑150 (MQ=255)
                                                                               cccTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGAAGAAACACCCGACATGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCa                                                                     <  1:320793/150‑1 (MQ=255)
                                                                                 cTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCAtt                                                                   >  1:878225/1‑150 (MQ=255)
                                                                                 cTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCAtt                                                                   >  2:7589/1‑150 (MQ=255)
                                                                                 cTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCAAtt                                                                   >  2:726699/1‑150 (MQ=255)
                                                                                    aTGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGt                                                                <  1:208267/150‑1 (MQ=255)
                                                                                     tGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGtt                                                               >  1:1395025/1‑150 (MQ=255)
                                                                                       cATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGtttt                                                             <  1:1321092/150‑1 (MQ=255)
                                                                                           ccACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCaa                                                         >  1:1385419/1‑150 (MQ=255)
                                                                                            cACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTAATCTCCATTTGTTTTGCAAc                                                        <  2:394812/150‑1 (MQ=255)
                                                                                                tAAGCCAGGATTTTCGCAAGGGAAGATGATGATGATGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATcc                                                    >  2:1048979/1‑150 (MQ=255)
                                                                                                 aaGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGATTTAAAAGAGAGCGTTTTTAATTATTCTCCATTTGTTTTGCAACATCCa                                                   >  1:941473/1‑150 (MQ=255)
                                                                                                    ccAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATaa                                                >  1:482462/1‑150 (MQ=255)
                                                                                                     cAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTTTTCTCCATTTGTTTTGCAACATCCATAAt                                               >  2:1486405/1‑150 (MQ=255)
                                                                                                      aGGATTTTCGCAAGGGAAGATTATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCATCATCCATAAtg                                              <  1:34221/150‑1 (MQ=255)
                                                                                                      aGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTTGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGCGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAAtg                                              >  2:240000/1‑150 (MQ=255)
                                                                                                       ggATTTTCGCAAGTGAAGATGATGTTTACGAAACACCCGACATGAATTTATGTCGTGTGCCTTGTTAAGGTCATAAGAATGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTCTTAATTATTCTCCTTTTGTTTTGCAACATCCATAAtgt                                             <  1:97463/150‑1 (MQ=255)
                                                                                                       ggATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAAtgt                                             <  1:1680698/150‑1 (MQ=255)
                                                                                                        gATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAAtgtg                                            <  1:1655859/150‑1 (MQ=255)
                                                                                                        gATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAAtgtg                                            >  1:175614/1‑150 (MQ=255)
                                                                                                         aTTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTgg                                           >  2:49717/1‑150 (MQ=255)
                                                                                                                 aaGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCgg                                   <  2:563353/150‑1 (MQ=255)
                                                                                                                  aGGGAATATGATGATTACGAAACACCCGACAGTAATTTATGTCGGGTGCCTTGTTAATGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTACATTTGTTTTGCAACATCCATAATGTGGTAAAGCGGt                                  <  1:549506/150‑1 (MQ=255)
                                                                                                                    ggAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCGGTAc                                >  2:622826/1‑150 (MQ=255)
                                                                                                                     taaGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCGGTAct                               <  2:542582/149‑1 (MQ=255)
                                                                                                                        catgatgatgACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCGGTACTCTg                            <  1:1708115/149‑1 (MQ=255)
                                                                                                                           gatgatgaCGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCGGTACTCTGTTc                         <  1:1205358/150‑1 (MQ=255)
                                                                                                                             tgatgaCGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCGGTACTCTGTTCtt                       <  1:1098021/150‑1 (MQ=255)
                                                                                                                              gatgaCGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCGGTACTCTGTTCTTa                      <  1:680650/150‑1 (MQ=255)
                                                                                                                               atgaCGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCGGTACTCTGTTCTTAg                     <  2:762498/150‑1 (MQ=255)
                                                                                                                                tgaCGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCGGTACTCTGTTCTTAGc                    <  1:881640/150‑1 (MQ=255)
                                                                                                                                tgaCGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCGGTACTCTGTTCTTAGc                    >  1:264460/1‑150 (MQ=255)
                                                                                                                                tgaCGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCGGTACTCTGTTCTTAGc                    <  2:1704389/150‑1 (MQ=255)
                                                                                                                                tgaCGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGAGGTACTCTGTTCTTatc                    >  2:1737560/1‑148 (MQ=255)
                                                                                                                                 gaCGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCGGTACTCTGTTCTTAGCt                   <  2:1587849/150‑1 (MQ=255)
                                                                                                                                 gaCGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGAGTTATTATTTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCGGTACTCTGTTCTTAGCt                   >  1:1175677/1‑150 (MQ=255)
                                                                                                                                  acaaggcacCCGACATAAATTCCTGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCGGTACTCTGTTCTTAGCTa                  >  1:1833641/7‑150 (MQ=255)
                                                                                                                                      ggcacCCGACATAAATTCCTGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCGGTACTCTGTTCTTAGCTAActc              >  2:1470049/3‑150 (MQ=255)
                                                                                                                                       gcacCCGACATAAATTCCTGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATATTGTGGTAAAGCGGTACTCTGTTCTTAGCTAActct             >  1:1459588/2‑150 (MQ=255)
                                                                                                                                       acacCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCGGTACTCTGTTCTTAGCTAActct             >  2:1608189/1‑150 (MQ=255)
                                                                                                                                        cacCCGACATAAATTCCTGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAATAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCGGTACTCTGTTCTTAGCTAACTCTa            >  1:154017/1‑150 (MQ=255)
                                                                                                                                        cacCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCGGTACTCTGTTCTTAGCTAACTCTa            <  2:1788211/150‑1 (MQ=255)
                                                                                                                                           ccGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCGGTACTCTGTTCTTAGCTAACTCTACCa         >  1:472253/1‑150 (MQ=255)
                                                                                                                                             gACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCGGTACTCTGTTCTTAGCTAACTCTACCATg       <  2:1516723/150‑1 (MQ=255)
                                                                                                                                              aCAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAACAGAGCGGAATTAATAATTCTCCATTTGTTTTGCTACATCCATTATGTGTTAAAGCGGTACTCTGTTCTTAGCTAACTCTACCATgg      <  2:25891/150‑1 (MQ=255)
                                                                                                                                               cAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCGGTACTCTGTTCTTAGCTAACTCTACCATGGc     <  2:1498822/150‑1 (MQ=255)
                                                                                                                                                aGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACCTCCATAATGTGGTAAAGCGGTACTCTGTTCTTAGCTAACTCTACCATGGCg    >  1:920704/1‑150 (MQ=255)
                                                                                                                                                  gAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCGGTACTCTGTTCTTAGCTAACTCTACCATGGCGtt  <  1:49717/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
TAACAATATGAAATATAAAGGAATATTTCATTTTCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCGGTACTCTGTTCTTAGCTAACTCTACCATGGCGTT  >  NC_000913/1579418‑1579713

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: