Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,579,564 G→A 16.9% 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,5640GA16.9% 372.9 / 57.8 159intergenic (‑222/+69)yddA/ydeMputative multidrug ABC transporter permease/ATPase/putative YdeN‑specific sulfatase‑maturating enzyme
Reads supporting (aligned to +/- strand):  ref base G (59/73);  new base A (14/13);  total (73/86)
Fisher's exact test for biased strand distribution p-value = 5.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.33e-01

Only 100 of 165 total aligned reads displayed.

TATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCC  >  NC_000913/1579465‑1579663
                                                                                                   |                                                                                                   
tatgCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAgg                                                                                                     >  2:4053300/1‑100 (MQ=255)
    cTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAAtt                                                                                                 <  2:3505869/100‑1 (MQ=255)
    cTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAAtt                                                                                                 <  1:339794/100‑1 (MQ=255)
     tATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAAttt                                                                                                >  1:1019793/1‑100 (MQ=255)
     tATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAAttt                                                                                                >  2:3008030/1‑100 (MQ=255)
        tGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATg                                                                                             >  2:1978418/1‑100 (MQ=255)
        tGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATg                                                                                             <  1:1802103/100‑1 (MQ=255)
         gATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGGCGAAACACCCGACAGGAATTTATGt                                                                                            >  1:1592649/1‑100 (MQ=255)
          aTAACCAAGGATAGTGCTAATGCGCTGATTCATACCACATAAGCCAGGATTTTAGCAAGGGTAGGGGATGATGACGAAACACCCGACAGGAATTTAt                                                                                              <  1:3167691/97‑1 (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:1645295/100‑1 (MQ=255)
              tCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACATAAATTCCTGTCGGGt                                                                                       <  2:947489/100‑1 (MQ=255)
              tCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACATAAATTCCTGTCGGGt                                                                                       <  1:1224811/100‑1 (MQ=255)
              tCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACATAAATTCCTGTCGGGt                                                                                       <  1:1754840/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)
                    atatTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACATAAATTCCTGTCGGGTgtttcg                                                                                 <  2:2311406/100‑6 (MQ=37)
                    atatTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTACCTTg                                                                                 <  1:3587024/100‑1 (MQ=255)
                     tatTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGt                                                                                >  2:1373168/1‑100 (MQ=255)
                     tatTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGt                                                                                <  2:4029375/100‑1 (MQ=255)
                     tatTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGt                                                                                >  2:356637/1‑100 (MQ=255)
                       ttGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAGCACCCGACATAAATTCCTGTCGGGTgtttcgtca                                                                              <  2:2705889/100‑9 (MQ=25)
                          cTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAgg                                                                           <  1:1523218/100‑1 (MQ=255)
                          cTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAgg                                                                           <  1:1741560/100‑1 (MQ=255)
                           tAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAgg                                                                           <  2:1801697/99‑1 (MQ=255)
                           tAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGt                                                                          <  1:3293018/100‑1 (MQ=255)
                           tAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGt                                                                          >  2:3251694/1‑100 (MQ=255)
                            aaTGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTc                                                                         >  2:2903923/1‑100 (MQ=255)
                             aTGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCa                                                                        <  2:2140526/100‑1 (MQ=255)
                               gCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATa                                                                      <  1:941133/100‑1 (MQ=255)
                                   tGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGCAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATaagaa                                                                  >  1:3711361/1‑100 (MQ=255)
                                     aTGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAagg                                                                <  2:426098/100‑1 (MQ=255)
                                     aTGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAagg                                                                <  2:3605552/100‑1 (MQ=255)
                                     aTGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAagg                                                                <  2:4131649/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)
                                          tACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCt                                                           >  1:1932547/1‑100 (MQ=255)
                                                  aaGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATgg                                                   >  2:2756333/1‑100 (MQ=255)
                                                   aGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGa                                                  <  2:3741955/100‑1 (MQ=255)
                                                    gCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAg                                                 <  2:2702964/100‑1 (MQ=255)
                                                    gCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAg                                                 <  1:995505/100‑1 (MQ=255)
                                                     ccAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTCAGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGCAGt                                                >  1:4171020/1‑100 (MQ=255)
                                                     ccAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGt                                                >  1:4019508/1‑100 (MQ=255)
                                                     ccAGGATCTTCGCAAGGAAAGATGATGATGACGAAACACCCGACCGGAATTTATGTCGGGTGCCTTGTCAAGGTCATAAGAAGGAGGCTAAGAATGGAg                                                 <  1:1930972/99‑1 (MQ=255)
                                                       aGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTa                                              <  1:707821/100‑1 (MQ=255)
                                                       aGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTa                                              >  2:3581585/1‑100 (MQ=255)
                                                       aGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTa                                              <  1:2903907/100‑1 (MQ=255)
                                                         gATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTaaa                                            <  2:3291529/100‑1 (MQ=255)
                                                             ttCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAagag                                        <  2:139371/100‑1 (MQ=255)
                                                             ttCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAagag                                        >  2:1604933/1‑100 (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                                      >  2:4320062/1‑100 (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)
                                                                      gAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGttatt                               <  2:634120/100‑1 (MQ=255)
                                                                      gAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGttatt                               <  2:2940097/100‑1 (MQ=255)
                                                                         gatgatgatGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTattaat                            <  1:4055092/100‑1 (MQ=255)
                                                                         gatgatgatGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTattaat                            <  1:845101/100‑1 (MQ=255)
                                                                         gatgatgatGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTattaan                            <  2:3922355/100‑2 (MQ=255)
                                                                          atgatgatgACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTattaatt                           <  1:4031858/100‑1 (MQ=255)
                                                                           tgatgatgaCGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTattaatta                          >  2:1124875/1‑100 (MQ=255)
                                                                            gatgatgaCGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAattat                         >  2:2837990/1‑100 (MQ=255)
                                                                              tgatgaCGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATtc                       <  2:3770741/100‑1 (MQ=255)
                                                                                  gaCGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCa                   >  2:2754285/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)
                                                                                     aaggcacCCGACATAAATTCCTGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGACTGGAGTTAAAAGAGAGCGTTATTAATTATTCTACAttt                >  1:351365/5‑100 (MQ=37)
                                                                                      aaaCACCCTACAGGAATTTATGTCGGGTGCCGTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCAtt                 <  1:1422281/98‑1 (MQ=255)
                                                                                      aaaCACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTg               >  2:1285213/1‑100 (MQ=255)
                                                                                      aaaCACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTg               <  2:3151625/100‑1 (MQ=255)
                                                                                        gcacCCGACATAAATTCCTGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGtt             >  1:966751/2‑100 (MQ=38)
                                                                                        gcacCCGACATAAATTCCTGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGtt             >  1:2647824/2‑100 (MQ=38)
                                                                                        acacCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGtt             <  2:3119241/100‑1 (MQ=255)
                                                                                         cacCCGACATAAATTCCTGTCGGGTGCCTTGTTACGGCCATAAGAAGGAGGATAAGACTGGAGTTAAAAGCGAGCGTTATTAATTATTCTCCATTTGt              >  1:2409201/1‑98 (MQ=25)
                                                                                         cacCCGACATAAATTCCTGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGttt            >  2:3968126/1‑100 (MQ=38)
                                                                                         cacCCGACATAAATTCCTGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGttt            >  1:1134056/1‑100 (MQ=38)
                                                                                         cacCCGACATAAATTCCTGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGttt            >  1:2392758/1‑100 (MQ=38)
                                                                                         cacCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGttt            <  2:1034190/100‑1 (MQ=255)
                                                                                         cacCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGttt            >  2:3215008/1‑100 (MQ=255)
                                                                                          acCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGtttt           <  2:509540/100‑1 (MQ=255)
                                                                                          acCCGACACGAATTTATGTCGGGCGCCTTGTTAAGGCCCTAAGAAGGAGGGTAAGAATGGATTTAAAAGAGAGCGTTATTAATTCTTCTCCATTTGtttt           >  1:4118852/1‑100 (MQ=255)
                                                                                            ccGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGc         <  1:1627682/100‑1 (MQ=255)
                                                                                            ccGACAGGAATTTATGTCGGGTGACTTGTTACGGTCATAAAAAGGACTCTAAGAATACAGTTAAAAGAGAGCGTCATTAATTATTCTCGATTTGtttt           >  1:4226133/1‑98 (MQ=255)
                                                                                                cAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACa     >  2:918218/1‑100 (MQ=255)
                                                                                                cAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACa     <  1:1646810/100‑1 (MQ=255)
                                                                                                  ggAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATc   >  2:1115760/1‑100 (MQ=255)
                                                                                                  ggAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATc   >  2:837152/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
TATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCC  >  NC_000913/1579465‑1579663

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: