Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,579,570 A→C 16.6% 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,5700AC16.6% 389.1 / 66.1 160intergenic (‑228/+63)yddA/ydeMputative multidrug ABC transporter permease/ATPase/putative YdeN‑specific sulfatase‑maturating enzyme
Reads supporting (aligned to +/- strand):  ref base A (59/74);  new base C (16/11);  total (75/85)
Fisher's exact test for biased strand distribution p-value = 2.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 162 total aligned reads displayed.

TGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATG  >  NC_000913/1579473‑1579669
                                                                                                 |                                                                                                   
tGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATg                                                                                                   <  1:1720687/100‑1 (MQ=255)
 gagaaTCAAGGATATTGCTAATGCCCGGATGCATACCACGTAATCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGt                                                                                                  >  1:4179283/4‑100 (MQ=255)
 gATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGGCGAAACACCCGACAGGAATTTATGt                                                                                                  >  1:1592649/1‑100 (MQ=255)
 gATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGt                                                                                                  >  2:4273760/1‑100 (MQ=255)
 gATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGt                                                                                                  >  1:1889616/1‑100 (MQ=255)
     aTCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAGGATGATGATGACGAAACACCCGACAGGAATTTATGTCggg                                                                                              >  2:1172701/1‑100 (MQ=255)
     aTCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACATAAATTCCTGTCggg                                                                                              <  1:4335713/100‑1 (MQ=255)
      tCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGTATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACATAAATTCCTGTCGGGt                                                                                             <  1:3344257/100‑1 (MQ=255)
      tCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACATAAATTCCTGTCGGGt                                                                                             <  2:947489/100‑1 (MQ=255)
      tCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACATAAATTCCTGTCGGGt                                                                                             <  2:1645295/100‑1 (MQ=255)
       cAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTg                                                                                            >  2:3501661/1‑100 (MQ=255)
        aaGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGTATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACATAAATTCCTGTCGGGTg                                                                                            <  1:916729/99‑1 (MQ=38)
         aGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGcc                                                                                          <  2:491393/100‑1 (MQ=255)
          ggATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCt                                                                                         >  1:563317/1‑100 (MQ=255)
            atatTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTACCTTg                                                                                       <  1:3587024/100‑1 (MQ=255)
             tatTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGt                                                                                      >  2:1373168/1‑100 (MQ=255)
             tatTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGt                                                                                      >  2:356637/1‑100 (MQ=255)
              atTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAATGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGtt                                                                                     <  1:3187328/100‑1 (MQ=255)
               ttGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAGCACCCGACATAAATTCCTGTCGGGTgtttcgtca                                                                                    <  2:2705889/100‑9 (MQ=25)
                  cTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAgg                                                                                 <  1:1741560/100‑1 (MQ=255)
                   tAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGt                                                                                >  2:3251694/1‑100 (MQ=255)
                   tAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGt                                                                                <  2:4065253/100‑1 (MQ=255)
                    aaTGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTc                                                                               >  2:2903923/1‑100 (MQ=255)
                     aTGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTc                                                                               <  2:1391972/99‑1 (MQ=255)
                     aTGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCa                                                                              <  2:2140526/100‑1 (MQ=255)
                      tGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCAt                                                                             <  2:2794719/100‑1 (MQ=255)
                      tGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCAt                                                                             >  1:1522626/1‑100 (MQ=255)
                       gCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATa                                                                            <  1:941133/100‑1 (MQ=255)
                        cccTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATaa                                                                           >  2:3690681/1‑100 (MQ=255)
                           tGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGCAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATaagaa                                                                        >  1:3711361/1‑100 (MQ=255)
                            gATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATaagaag                                                                       <  2:270405/100‑1 (MQ=255)
                             aTGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAagg                                                                      <  2:426098/100‑1 (MQ=255)
                                cATACCACGTAAGCCAGGATTTTCGCAAGGGAGGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAaggagg                                                                   <  1:1172696/100‑1 (MQ=255)
                                cATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAaggagg                                                                   <  2:4369770/100‑1 (MQ=255)
                                 aTACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGc                                                                  <  2:1846505/100‑1 (MQ=255)
                                   aCCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTa                                                                >  2:1700002/1‑100 (MQ=255)
                                        gTAAGCCAGGAGTTTCGCAAGAGATGAGGCTGATGACTAAATACCCGACAAAAAGTTATGTCGGGTGCCTTGTTAAGGTCATAATAAGGAGGCTAAGaa                                                            <  1:3527428/99‑1 (MQ=255)
                                         tAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATg                                                          <  1:1558896/100‑1 (MQ=255)
                                           agtcAGTATTTTCACAAGTGTAGATGATGACGACGAAAAACCCCACAAAAATTCATGTCAAGTGCCTTGTTAATGTCATAAGAAGGAGGCTAAGaagcga                                                        <  1:1760320/97‑5 (MQ=255)
                                           aGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGa                                                        <  2:3741955/100‑1 (MQ=255)
                                            gCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAGCACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAg                                                       <  2:3097009/100‑1 (MQ=255)
                                            gCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGan                                                       <  1:123247/100‑2 (MQ=255)
                                            gCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAg                                                       <  1:995505/100‑1 (MQ=255)
                                            gCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCGTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAg                                                       <  2:1081998/100‑1 (MQ=255)
                                             ccAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTCAGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGCAGt                                                      >  1:4171020/1‑100 (MQ=255)
                                             ccAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGt                                                      >  2:1535269/1‑100 (MQ=255)
                                             ccAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGt                                                      >  1:4019508/1‑100 (MQ=255)
                                             ccAGGATCTTCGCAAGGAAAGATGATGATGACGAAACACCCGACCGGAATTTATGTCGGGTGCCTTGTCAAGGTCATAAGAAGGAGGCTAAGAATGGAg                                                       <  1:1930972/99‑1 (MQ=255)
                                               aGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTa                                                    >  2:3581585/1‑100 (MQ=255)
                                               aGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTa                                                    <  1:2903907/100‑1 (MQ=255)
                                                    tttCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAaga                                               <  2:2242376/100‑1 (MQ=255)
                                                     ttCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAagag                                              >  2:2530345/1‑100 (MQ=255)
                                                      tCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAagaga                                             <  2:487224/100‑1 (MQ=255)
                                                      tCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAagag                                              <  1:1439910/99‑1 (MQ=255)
                                                       cGCAAGGGAAGATTGTGGTGACGAAACACCCTACAGGAATTTATGTCGGGTGGCTTGTTAAGCTCATAAGAAGTAGGTTAAGAAGTGAGGAAAAagc                                               <  1:2160870/97‑2 (MQ=255)
                                                       cGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAagagag                                            >  1:4258217/1‑100 (MQ=255)
                                                       cGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAagagag                                            <  2:2697586/100‑1 (MQ=255)
                                                        gCAAGGGGAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGc                                           >  1:4124586/1‑100 (MQ=255)
                                                        gCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAagagaac                                           >  1:707987/1‑98 (MQ=255)
                                                        gCAAGGGAAGATGATGATGACGAAAAACCCGACAGGAATTTATGTCGGGTGCCTAGTTAAGGTCCTAAGAAGTAGGCTAAGAATGGAGTTAAAAGAGAGc                                           >  1:1015915/1‑100 (MQ=255)
                                                           aGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGtt                                        >  1:3549143/1‑100 (MQ=255)
                                                            gggAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGtta                                       <  1:3573417/100‑1 (MQ=255)
                                                             ggAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGttat                                      >  2:11821/1‑100 (MQ=255)
                                                              gAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGttatt                                     <  2:2000596/100‑1 (MQ=255)
                                                              gAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGttatt                                     <  2:634120/100‑1 (MQ=255)
                                                              gAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGttatt                                     <  1:1646293/100‑1 (MQ=255)
                                                                 gatgatgatGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTattaat                                  <  1:845101/100‑1 (MQ=255)
                                                                 gatgatgatGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTattaat                                  <  1:4055092/100‑1 (MQ=255)
                                                                  atgatgatgACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTattaatt                                 <  1:4031858/100‑1 (MQ=255)
                                                                   tgatgatgaCGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTattaatta                                >  2:1124875/1‑100 (MQ=255)
                                                                    gatgatgaCGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAattat                               >  2:2837990/1‑100 (MQ=255)
                                                                     atgatgaCGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAattatt                              <  2:3959257/100‑1 (MQ=255)
                                                                     atgatgaCGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAattatt                              >  1:1200003/1‑100 (MQ=255)
                                                                      tgatgaCGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATtc                             >  2:331404/1‑100 (MQ=255)
                                                                          gaCGAAACACCCGACAGGAATTTCTGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCa                         >  1:2040981/1‑100 (MQ=255)
                                                                          gaCGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCa                         >  2:2754285/1‑100 (MQ=255)
                                                                            cGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCAtt                       >  2:4034694/1‑100 (MQ=255)
                                                                            cGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCAtt                       >  1:1128709/1‑100 (MQ=255)
                                                                             gAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGTAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCAttt                      <  1:3244475/100‑1 (MQ=255)
                                                                             gAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCAttt                      >  2:3356905/1‑100 (MQ=255)
                                                                             aaggcacCCGACATAAATTCCTGTCGGGTGTCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCAttt                      >  2:4148844/5‑100 (MQ=37)
                                                                              aaaCACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTg                     <  2:3151625/100‑1 (MQ=255)
                                                                               aaCACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGt                    <  1:331403/100‑1 (MQ=255)
                                                                                gcacCCGACATAAATTCCTGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGtt                   >  1:2647824/2‑100 (MQ=38)
                                                                                 cacCCGACATAAATTCCTGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAACAGAGAGCGTTATTAATTATTCTCCATTTGttt                  >  2:3830662/1‑100 (MQ=38)
                                                                                 cacCCGACATAAATTCCTGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGttt                  >  2:3968126/1‑100 (MQ=38)
                                                                                 cacCCGACATAAATTCCTGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGttt                  >  1:2428963/1‑100 (MQ=38)
                                                                                 cacCCGACATAAATTCCTGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGttt                  >  1:2392758/1‑100 (MQ=38)
                                                                                 cacCCGACATAAATTCCTGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGttt                  >  1:952075/1‑100 (MQ=38)
                                                                                 cacCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGttt                  >  2:3193968/1‑100 (MQ=255)
                                                                                 cacCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGttt                  <  2:1034190/100‑1 (MQ=255)
                                                                                  acCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGtttt                 <  2:509540/100‑1 (MQ=255)
                                                                                  acCCGACACGAATTTATGTCGGGCGCCTTGTTAAGGCCCTAAGAAGGAGGGTAAGAATGGATTTAAAAGAGAGCGTTATTAATTCTTCTCCATTTGtttt                 >  1:4118852/1‑100 (MQ=255)
                                                                                     cGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCa              <  2:1840604/100‑1 (MQ=255)
                                                                                        cAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACa           >  2:918218/1‑100 (MQ=255)
                                                                                          ggAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATc         >  2:837152/1‑100 (MQ=255)
                                                                                          ggAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATc         >  2:3004671/1‑100 (MQ=255)
                                                                                            aaTTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCa       >  2:2702546/1‑100 (MQ=255)
                                                                                                 aTGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAAtg  >  2:4398596/1‑100 (MQ=255)
                                                                                                 aTGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAAtg  <  1:3359896/100‑1 (MQ=255)
                                                                                                 |                                                                                                   
TGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATG  >  NC_000913/1579473‑1579669

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