Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,579,564 G→A 5.3% 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 5.3% 577.3 / 14.6 169intergenic (‑222/+69)yddA/ydeMputative multidrug ABC transporter permease/ATPase/putative YdeN‑specific sulfatase‑maturating enzyme
Reads supporting (aligned to +/- strand):  ref base G (78/82);  new base A (3/6);  total (81/88)
Fisher's exact test for biased strand distribution p-value = 4.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.03e-01

Only 100 of 175 total aligned reads displayed.

TTTAACAATATGAAATATAAAGGAATATTTCATTTTCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCGGTACTCTGTTCTTAGCTAACTCTACCATGGCG  >  NC_000913/1579416‑1579711
                                                                                                                                                    |                                                                                                                                                   
tttAACAATATGAAATATAAAGGAATATTTCATTTTCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAgg                                                                                                                                                     >  2:235241/1‑149 (MQ=255)
  tAACAATATGAAATATAAAGGAATATCTCATTTTCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGaa                                                                                                                                                   <  1:96456/149‑1 (MQ=255)
    aCAATATGAAATATAAAGGAATATTTCATTTTCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAAtt                                                                                                                                                 <  1:986706/149‑1 (MQ=255)
      aaTATGAAATATAAAGGAATATTTCATTTTCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGaa                                                                                                                                                   <  2:839189/145‑1 (MQ=255)
      aaTATGAAATATAAAGGAATATTTCATTTTCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGaa                                                                                                                                                   >  1:839189/1‑145 (MQ=255)
            aaaTATAAAGGAATATTTCATTTTCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACATAAATTCCTGTCgg                                                                                                                                         >  3:1269961/1‑149 (MQ=255)
            aaaTATAAAGGAATATTTCATTTTCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCgg                                                                                                                                         >  4:1458641/1‑149 (MQ=255)
             aaTATAAAGGAATATTTCATTTTCAGTGGTTATGGTTATGCTATTGACAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACATAAATTCCTGTCggg                                                                                                                                        <  4:1269961/149‑1 (MQ=255)
              atatAAAGGAATATTTCATTTTCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACATAAATTCCTGTCGGGt                                                                                                                                       <  1:1850281/149‑1 (MQ=255)
              atatAAAGGAATATTTCATTTTCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACATAAATTCCTGTCGGGt                                                                                                                                       <  3:129891/149‑1 (MQ=255)
              atatAAAGGAATATTTCATTTTCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGt                                                                                                                                       <  3:1108378/149‑1 (MQ=255)
                   aaGGAATATTTCATTTTCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCtt                                                                                                                                  >  3:465103/1‑149 (MQ=255)
                   aaGGAATATTTCATTTTCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCtt                                                                                                                                  >  2:1422203/1‑149 (MQ=255)
                        atatTTCATTTTCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTaa                                                                                                                             <  4:465103/149‑1 (MQ=255)
                        atatTTCATTTTCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTaa                                                                                                                             <  1:1422203/149‑1 (MQ=255)
                        atatTTCATTTTCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTaa                                                                                                                             <  2:1264076/149‑1 (MQ=255)
                           tttCATTTTCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGt                                                                                                                          >  2:323116/1‑149 (MQ=255)
                             tCATTTTCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCtt                                                                                                                                  <  4:1316985/139‑1 (MQ=255)
                             tCATTTTCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCtt                                                                                                                                  >  3:1645644/1‑139 (MQ=255)
                             tCATTTTCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCtt                                                                                                                                  <  4:1645644/139‑1 (MQ=255)
                             tCATTTTCAGTGGTTATGGTTATGCTATTGATAATCAAGGACATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCtt                                                                                                                                  >  3:1316985/1‑139 (MQ=255)
                             tCATTTTCACTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCtt                                                                                                                                  <  4:1287567/139‑1 (MQ=255)
                                ttttCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATaa                                                                                                                     <  1:235241/149‑1 (MQ=255)
                                   tCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATaagaa                                                                                                                  >  2:299916/1‑149 (MQ=255)
                                     aGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCACTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAagg                                                                                                                <  1:777835/149‑1 (MQ=255)
                                         gttatggttatgCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACATAAATTCCTGTCggg                                                                                                                                        >  3:680522/1‑121 (MQ=255)
                                         gttatggttatgCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACATAAATTCCTGTCggg                                                                                                                                        <  4:680522/121‑1 (MQ=255)
                                          ttatggttatgCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGt                                                                                                                          <  3:355092/134‑1 (MQ=255)
                                          ttatggttatgCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCt                                                                                                           >  4:748674/1‑149 (MQ=255)
                                          ttatggttatgCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCt                                                                                                           >  1:1455409/1‑149 (MQ=255)
                                                     cTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATgg                                                                                                   >  3:649175/1‑146 (MQ=255)
                                                      tATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAaggagg                                                                                                             >  2:826100/1‑135 (MQ=255)
                                                       aTTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAaggag                                                                                                              >  2:832634/1‑133 (MQ=255)
                                                       aTTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAaggag                                                                                                              <  1:832634/133‑1 (MQ=255)
                                                        ttGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAg                                                                                                                            >  1:252226/1‑118 (MQ=255)
                                                        ttGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTaa                                                                                                         <  1:343581/137‑1 (MQ=255)
                                                         tGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTaaa                                                                                            <  3:1458641/149‑1 (MQ=255)
                                                          gATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGc                                                                                                            <  3:985914/132‑1 (MQ=255)
                                                          gATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGc                                                                                                            >  4:985914/1‑132 (MQ=255)
                                                          gATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTaaaa                                                                                           <  1:452696/149‑1 (MQ=255)
                                                           aTAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCgg                                                                                                                                         <  3:1467396/102‑1 (MQ=255)
                                                           aTAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCgg                                                                                                                                         >  4:1467396/1‑102 (MQ=255)
                                                           aTAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAag                                                                                          >  3:414551/1‑149 (MQ=255)
                                                               tCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAagagag                                                                                      >  1:764849/1‑149 (MQ=255)
                                                                 aaGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGa                                                                                                  >  2:38762/1‑135 (MQ=255)
                                                                 aaGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTaaaa                                                                                           >  1:123983/1‑142 (MQ=255)
                                                                 aaGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTaaaa                                                                                           <  2:123983/142‑1 (MQ=255)
                                                                 aaGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCg                                                                                    >  4:232471/1‑149 (MQ=255)
                                                                    gATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGa                                                                                                  <  1:683674/132‑1 (MQ=255)
                                                                      tatTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAagaga                                                                                       <  2:196180/141‑1 (MQ=255)
                                                                      tatTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAagaga                                                                                       >  1:196180/1‑141 (MQ=255)
                                                                            tAATGCCCTGATGCATACCGCGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAattat                                                                         >  3:1118069/1‑149 (MQ=255)
                                                                                gCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATaagaag                                                                                                                 <  1:1711866/105‑1 (MQ=255)
                                                                                gCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAg                                                                                                 <  3:849661/121‑1 (MQ=255)
                                                                                gCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATtctc                                                                     <  3:232471/149‑1 (MQ=255)
                                                                                 cccTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACATAAAttc                                                                                                                                                <  4:1822288/73‑2 (MQ=255)
                                                                                  ccTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATaa                                                                                                                     <  4:844959/99‑1 (MQ=255)
                                                                                  ccTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGtt                                                                                               <  2:1577893/121‑1 (MQ=255)
                                                                                  ccTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCa                                                                   <  4:414551/149‑1 (MQ=255)
                                                                                   cTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCAt                                                                  >  1:1197266/1‑149 (MQ=255)
                                                                                       tGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTaaaa                                                                                           >  3:101716/1‑120 (MQ=255)
                                                                                       tGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGt                                                              <  1:11657/149‑1 (MQ=255)
                                                                                           tACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATg                                                                                                    >  4:529575/1‑107 (MQ=255)
                                                                                           tACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATg                                                                                                    <  3:529575/107‑1 (MQ=255)
                                                                                              cACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAagaga                                                                                       >  3:344059/1‑117 (MQ=255)
                                                                                              cACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAagaga                                                                                       <  4:345896/117‑1 (MQ=255)
                                                                                              cACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAagaga                                                                                       <  4:344059/117‑1 (MQ=255)
                                                                                              cACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCAtt                                                                 <  3:1841928/139‑1 (MQ=255)
                                                                                              cACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCAtt                                                                 >  4:1841928/1‑139 (MQ=255)
                                                                                                cGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAg                                                                                                                            <  1:1364033/78‑1 (MQ=255)
                                                                                                cGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTattaat                                                                            <  3:1688221/126‑1 (MQ=255)
                                                                                                   aaGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTaa                                                                                             >  3:252708/1‑106 (MQ=255)
                                                                                                      ccAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGa                                                                                                       <  2:1388496/93‑1 (MQ=255)
                                                                                                        aGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTaa                                                                                                         >  3:1133856/1‑89 (MQ=255)
                                                                                                        aGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTaa                                                                                                         <  4:1133856/89‑1 (MQ=255)
                                                                                                                        aaGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGa                                                                                                  >  2:1827824/1‑80 (MQ=255)
                                                                                                                           atgatgatgACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGtttt                                                           <  1:544788/116‑1 (MQ=255)
                                                                                                                           atgatgatgACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGtttt                                                           >  2:544788/1‑116 (MQ=255)
                                                                                                                            tgatgatgaCGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCGGTACTCTGt                         >  4:1501534/1‑149 (MQ=255)
                                                                                                                            tgatgatgaCGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCGGTACTCTGt                         >  4:729075/1‑149 (MQ=255)
                                                                                                                             gatgatgaCGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTattaa                                                                             <  4:542006/96‑1 (MQ=255)
                                                                                                                             gatgatgaCGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCGGTACTCTGtt                        >  4:835567/1‑149 (MQ=255)
                                                                                                                                gatgaCGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAagaga                                                                                       <  3:286860/83‑1 (MQ=255)
                                                                                                                                gatgaCGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAagaga                                                                                       >  4:286860/1‑83 (MQ=255)
                                                                                                                                gatgaCGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCGGTACTCTGTTCtt                     <  3:835567/149‑1 (MQ=255)
                                                                                                                                      gAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCGGTACTCTGTTCTTAGCTaa               >  4:1791228/1‑149 (MQ=255)
                                                                                                                                      gAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCGGTACTCTGTTCTTAGCTaa               <  3:1791228/149‑1 (MQ=255)
                                                                                                                                      gAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCGGTACTCTGTTCTTAGCTaa               <  1:299916/149‑1 (MQ=255)
                                                                                                                                           acCCGACATAAATTCCTGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATtc                                                                       <  3:736877/88‑1 (MQ=38)
                                                                                                                                            cccGGCAGGAATTCATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCGGTACTCTGTTCTTAGCTAACTCTAc         >  2:1484899/1‑149 (MQ=255)
                                                                                                                                             ccGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTattaatt                                                                           <  4:1076082/82‑1 (MQ=255)
                                                                                                                                             ccGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGtttt                                                           <  2:640785/98‑1 (MQ=255)
                                                                                                                                             ccGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGtttt                                                           >  1:640785/1‑98 (MQ=255)
                                                                                                                                              cGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCGGTACTCTGTTCTTAGCTAACTCTACCa       <  1:589953/149‑1 (MQ=255)
                                                                                                                                                aCAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTaaa                                     >  1:156947/1‑117 (MQ=255)
                                                                                                                                                  aGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCGGTACTCTGtt                        >  2:622011/1‑128 (MQ=255)
                                                                                                                                                   ggAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCGGTACTCTGTTCTTAGCTAACTCTACCATGGCg  >  2:661667/1‑149 (MQ=255)
                                                                                                                                                    gAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATc                                                   >  2:698394/1‑99 (MQ=255)
                                                                                                                                                    gAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATc                                                   <  1:1868192/99‑1 (MQ=255)
                                                                                                                                                    gAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATc                                                   >  2:1868192/1‑99 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TTTAACAATATGAAATATAAAGGAATATTTCATTTTCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCGGTACTCTGTTCTTAGCTAACTCTACCATGGCG  >  NC_000913/1579416‑1579711

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