Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,579,570 A→C 6.4% intergenic (‑228/+63) 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,5700AC 6.4% 332.1 / 11.0 110intergenic (‑228/+63)yddA/ydeMputative multidrug ABC transporter permease/ATPase/putative YdeN‑specific sulfatase‑maturating enzyme
Reads supporting (aligned to +/- strand):  ref base A (50/53);  new base C (4/3);  total (54/56)
Fisher's exact test for biased strand distribution p-value = 7.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.34e-01

Only 100 of 110 total aligned reads displayed.

ATATGAAATATAAAGGAATATTTCATTTTCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCGGTACTCTGTTCTTAGCTAACTCTACCATGGCGTTCATA  >  NC_000913/1579423‑1579717
                                                                                                                                                   |                                                                                                                                                   
atatGAAATATAAATGAAAATTTCATTTTCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCAAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATg                                                                                                                                                   <  1:345494/150‑1 (MQ=255)
atatGAAATATAAAGGAATATTTCATTTTCAGTGGTTATGGTTATTCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATg                                                                                                                                                   <  2:1260321/150‑1 (MQ=255)
     aaaTATAAAGGAATATTTCATTTTCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCggg                                                                                                                                              <  1:1122984/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)
                          tttCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATaaga                                                                                                                         >  2:1205358/1‑150 (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)
                                    tatggttatgCTTTTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACTAAACACCCGACATGAATATATGTCGGTTGCCTTGTTAAGGTCTTAAGAATGAGGCTaa                                                                                                               >  2:549506/1‑150 (MQ=255)
                                     atggttatgCTATTGATAATAAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGAAGAAACACCCGACAGGAATTTATGTCGGGGGCCTTGTTAAGGTCATAAGAAGGAGGCTAAg                                                                                                              >  2:154766/1‑150 (MQ=255)
                                       ggttatgCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGaa                                                                                                            <  2:1143574/150‑1 (MQ=255)
                                       ggttatgCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGaa                                                                                                            >  2:269387/1‑150 (MQ=255)
                                        gttatgCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAAt                                                                                                           >  1:762498/1‑150 (MQ=255)
                                            tgCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAg                                                                                                       <  1:930681/150‑1 (MQ=255)
                                            tgCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAg                                                                                                       <  1:1004091/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:1256783/1‑150 (MQ=255)
                                                aTTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTaa                                                                                                   <  1:1680559/150‑1 (MQ=255)
                                                 ttGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTaaa                                                                                                  >  2:1831240/1‑150 (MQ=255)
                                                  tGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTaaaa                                                                                                 >  1:1064148/1‑150 (MQ=255)
                                                  tGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTaaaa                                                                                                 >  1:1340621/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:1020252/1‑150 (MQ=255)
                                                                     tAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAattatt                                                                              >  1:1800239/1‑150 (MQ=255)
                                                                          cccTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCa                                                                         <  2:494229/150‑1 (MQ=255)
                                                                          cccTGATGCATACCACGGAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGCGCGTTATTAATTATTCTACa                                                                         >  1:238396/1‑150 (MQ=255)
                                                                            cTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCAtt                                                                       >  1:878225/1‑150 (MQ=255)
                                                                            cTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCAtt                                                                       >  2:7589/1‑150 (MQ=255)
                                                                            cTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCAAtt                                                                       >  2:726699/1‑150 (MQ=255)
                                                                              gATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTg                                                                     >  1:1701521/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)
                                                                                            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)
                                                                                                  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                        <  2:1704389/150‑1 (MQ=255)
                                                                                                                           tgaCGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCGGTACTCTGTTCTTAGc                        <  1:881640/150‑1 (MQ=255)
                                                                                                                           tgaCGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCGGTACTCTGTTCTTAGc                        >  1:264460/1‑150 (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)
                                                                                                                                   cacCCGACATAAATTCCTGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAATAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCGGTACTCTGTTCTTAGCTAACTCTa                >  1:154017/1‑150 (MQ=255)
                                                                                                                                   cacCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCGGTACTCTGTTCTTAGCTAACTCTa                <  2:1788211/150‑1 (MQ=255)
                                                                                                                                   cacCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCGGTACTCTGTTCTTAGCTAACTCTa                <  2:1762333/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)
                                                                                                                                              aaTTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCGGTACTCTGTTCTTAGCTAACTCTACCATGGCGTTc     <  2:1800239/150‑1 (MQ=255)
                                                                                                                                                 ttATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCGGTACTCTGTTCTTAGCTAACTCTACCATGGCGTTCata  <  2:554674/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
ATATGAAATATAAAGGAATATTTCATTTTCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCGGTACTCTGTTCTTAGCTAACTCTACCATGGCGTTCATA  >  NC_000913/1579423‑1579717

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: