Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,705,9780AC32.5% 353.7 / 238.3 238intergenic (‑188/+120)dppB/dppAdipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base A (68/71);  new base C (63/4);  total (163/75)
Fisher's exact test for biased strand distribution p-value = 1.83e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.92e-01

Only 100 of 243 total aligned reads displayed.

GAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCAT  >  NC_000913/3705879‑3706077
                                                                                                   |                                                                                                         
gAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTa                                                                                                           >  2:1010725/1‑100 (MQ=255)
 naGTATTGCTGCCGGATGGCGCTGCACGCAGAGGNACANCGCCNCCNGGTANACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTaa                                                                                                          >  1:1402848/2‑100 (MQ=255)
  aGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAAc                                                                                                         >  2:1916493/1‑100 (MQ=255)
    tATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGc                                                                                                       >  2:2523086/1‑100 (MQ=255)
    tATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGc                                                                                                       >  1:2233660/1‑100 (MQ=255)
     aTTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGcc                                                                                                      <  1:1602278/100‑1 (MQ=255)
      ttGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCt                                                                                                     <  1:2826148/100‑1 (MQ=255)
        gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                                   >  1:2180198/1‑100 (MQ=255)
         ctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTAt                                                                                                  <  2:944258/100‑1 (MQ=255)
          tgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATc                                                                                                 <  2:1570983/100‑1 (MQ=255)
           gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                                >  1:2669894/1‑100 (MQ=255)
           gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                                >  1:1812156/1‑100 (MQ=255)
            cctgATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCg                                                                                               <  1:57057/97‑1 (MQ=255)
            ccGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCg                                                                                               >  1:377671/1‑100 (MQ=255)
              ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGc                                                                                             <  2:2049852/100‑1 (MQ=255)
               gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                            <  2:689082/100‑1 (MQ=255)
               gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                            <  1:202900/100‑1 (MQ=255)
               gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                            <  2:1095002/100‑1 (MQ=255)
                aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCt                                                                                           <  1:2860538/100‑1 (MQ=255)
                 tGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTa                                                                                          >  2:1322059/1‑100 (MQ=255)
                   gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACa                                                                                        <  1:1420361/100‑1 (MQ=255)
                    cgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACaa                                                                                       <  1:2191731/100‑1 (MQ=255)
                    cgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACaa                                                                                       >  1:2638991/1‑100 (MQ=255)
                    cgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACaa                                                                                       <  2:797524/100‑1 (MQ=255)
                       tGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAAt                                                                                    >  1:1007048/1‑100 (MQ=255)
                            gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGc                                                                               >  1:271997/1‑100 (MQ=255)
                             cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCa                                                                              >  2:1065124/1‑100 (MQ=255)
                             cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCa                                                                              <  2:522884/100‑1 (MQ=255)
                                   nacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCa                                                                        >  1:397039/2‑100 (MQ=255)
                                      aGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGTCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATa                                                                     <  2:2704773/100‑1 (MQ=255)
                                         gcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACAt                                                                  <  1:791655/100‑1 (MQ=255)
                                           cACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTg                                                                >  1:1263404/1‑100 (MQ=255)
                                            aCCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGc                                                               <  2:711688/100‑1 (MQ=255)
                                              cGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑Ag                                                          >  2:1185729/1‑100 (MQ=255)
                                                ggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGaa                                                        >  2:2846732/1‑100 (MQ=255)
                                                 gTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAAt                                                       <  2:233905/100‑1 (MQ=255)
                                                   ataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAAttt                                                    >  2:2602900/1‑100 (MQ=255)
                                                     aCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTc                                                  <  2:1263636/100‑1 (MQ=255)
                                                          gcataaaCAAAGCGCACTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAgg                                           >  2:97915/7‑100 (MQ=18)
                                                               gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGa                                      >  1:1797744/1‑100 (MQ=255)
                                                               gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGa                                      >  1:1091185/1‑100 (MQ=255)
                                                                     gTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGcn                                <  2:397038/100‑2 (MQ=255)
                                                                       gcgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGa                              <  2:1482778/100‑1 (MQ=255)
                                                                         gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAgc                            <  1:216799/100‑1 (MQ=255)
                                                                          gcttttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGC‑‑AGAGTTCAC‑‑GTAGGCCTGATAAGCGTAgcg                           >  2:2197671/11‑100 (MQ=17)
                                                                          cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAgcg                           >  1:1876506/1‑100 (MQ=255)
                                                                           gttttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTGTGCAAATTCAATATATTGC‑‑AGGAAACAC‑‑GTAGGCCTGATAAGCGAAgcgc                          >  1:2581205/10‑100 (MQ=18)
                                                                             tttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTTTGCAAATTCAATATATTGC‑‑AATCTCCGTGTAGGCCTGATAAGCGTAGCGCAt                        >  1:1181570/8‑100 (MQ=18)
                                                                             tGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCACTACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCAt                        >  2:2594305/1‑100 (MQ=255)
                                                                             tGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCAt                        >  2:757035/1‑100 (MQ=255)
                                                                              ttgtaggCCGGATAAGGCGTTCGCGCCGCATCCGGCATACAAAAGTTTGCAAATTCAATAAATTGC‑‑AGAATGCAT‑‑GTAGGCCTGATAAGCGAAGCACATc                       >  2:2452989/7‑100 (MQ=12)
                                                                               tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCa                      >  1:2695557/1‑100 (MQ=255)
                                                                               tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCa                      >  1:2072212/1‑100 (MQ=255)
                                                                               tgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGC‑‑AGAAACCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCa                      >  1:944187/6‑100 (MQ=18)
                                                                                ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAg                     <  2:421558/100‑1 (MQ=255)
                                                                                ttATGCCGGATGCGGCGTTAACGCCGTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAg                     <  2:500621/100‑1 (MQ=255)
                                                                                gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATCGCTTAAATTCAATATATTGC‑‑ACCATCCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                     >  2:2422358/5‑100 (MQ=21)
                                                                                 taggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAACAGTGCAAATTCAATATATTGC‑‑AATCTTCTG‑‑GTAGGCCTGATAAGCATCGCGCATCAgg                    >  1:1280926/4‑100 (MQ=18)
                                                                                 taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAACCTGCAAATTCAATATATTGC‑‑AGGATACGT‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                    >  1:2459294/4‑100 (MQ=18)
                                                                                 taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAACCTGCAAATTCAATATATTGC‑‑AGGATACGT‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                    >  2:2018793/4‑100 (MQ=17)
                                                                                 taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAATCATGCAAATCCAAAATATTGC‑‑AGGAGTTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                   >  1:356265/4‑100 (MQ=21)
                                                                                  aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAATTCTGCAAATTCAACAAATTGC‑‑GAATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                  >  2:1655796/3‑100 (MQ=18)
                                                                                  aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGC‑‑AGAGTTCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                   >  1:680191/3‑100 (MQ=17)
                                                                                  aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGC‑‑AGTTTATTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGc                   >  1:2344368/3‑100 (MQ=25)
                                                                                  aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATCGCTTAAATTCAATATATTGC‑‑ACCATCCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                   >  1:2434943/3‑100 (MQ=21)
                                                                                  aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACCTTGCTAATTCAATATATTGC‑‑AGGGACTAT‑‑GTAGGCCTGATAAGCATAGCGCATCAGGc                   >  2:1784239/3‑100 (MQ=18)
                                                                                  aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAATCATGCAAATCCAAAATATTGC‑‑AGGAGTTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                  >  1:2127215/3‑100 (MQ=21)
                                                                                  aggCCGGATAAGGCGTTAACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAACTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGc                   >  2:1243322/3‑100 (MQ=17)
                                                                                   ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGC‑‑AGAGTTCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                  >  1:1685046/2‑100 (MQ=17)
                                                                                   ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGC‑‑AGAGTTCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                  >  2:2830079/2‑100 (MQ=17)
                                                                                   ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGC‑‑AGAGTTCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                  >  2:315588/2‑100 (MQ=17)
                                                                                   ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAGATCGTGCAAATTCAATATATTGC‑‑ATGAGCCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                  >  1:1055559/2‑100 (MQ=17)
                                                                                   ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGC‑‑AGTTTATTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCa                  >  1:1927959/2‑100 (MQ=21)
                                                                                   ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACCTTGCTAATTCAATATATTGC‑‑AGGGACTAT‑‑GTAGGCCTGATAAGCATAGCGCATCAGGCa                  >  2:2094932/2‑100 (MQ=18)
                                                                                   ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAACCTGCAAATTCAATATATTGC‑‑AGGATACGT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                  >  2:600718/2‑100 (MQ=18)
                                                                                   ggCCGGATAAGGCGTTAACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAATACTTTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGAGCATCAGGc                   >  2:2493784/2‑99 (MQ=17)
                                                                                    gCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt                 >  1:310663/1‑100 (MQ=255)
                                                                                    gCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGC‑‑AGAGTTCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                 >  2:1347844/1‑99 (MQ=17)
                                                                                    gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATCGCTTAAATTCAATATATTGC‑‑ACCATCCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa                 >  1:1572285/1‑99 (MQ=21)
                                                                                    gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAACCGCGCAAATTCAATATATTGCGGAGAAAAT‑‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGctgt                >  2:2659404/1‑97 (MQ=21)
                                                                                     ccggataaggcgTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGC‑‑AGTTTATTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaat                >  1:2310095/9‑98 (MQ=25)
                                                                                     ccggataaggcgTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGC‑‑AGTTTATTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaat                >  1:1773462/9‑98 (MQ=21)
                                                                                     ccggataaggcgTTCACGCCGCATCCGGCCTACAAAAACCTGCAAATTCAATATATTGC‑‑AGGATACGT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaat                >  1:2784144/9‑98 (MQ=18)
                                                                                     ccGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtt                >  1:1874648/1‑100 (MQ=255)
                                                                                     ccGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGC‑‑AGAAACCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                 >  1:2224785/1‑98 (MQ=17)
                                                                                     ccGGATAAGGCGTTCACGCCGCATCCGGCCTACAAACCGCGCAAATTCAATATATTGCGGAGAAAAT‑‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGctgtt               >  1:477658/1‑96 (MQ=21)
                                                                                     caggataaggcgTTCACGCCGCATCCTGCCTACAAAAGCACGCAAATTCAATATATTGC‑‑AGAGATGGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtt                >  1:819304/9‑100 (MQ=14)
                                                                                      cggataaggcgTTTACGCCGCATCCGGCATACAAAAGCATGCAAATTCAATATATTGC‑‑AGAGTTCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagtt               >  2:1429748/8‑97 (MQ=17)
                                                                                      cgatctaggcgTTCACGCCGCATCCGACCCACAAAGTCTTGCAAAATCAATAAATTGC‑‑AGGAACGAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaact               <  1:1146744/93‑4 (MQ=21)
                                                                                       ggataaggcgTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGC‑‑AGAAACCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt              >  1:2001385/7‑100 (MQ=18)
                                                                                       ggataaggcgTTCACGCCGCATCCGGCCTACAAATCATGCAAATCCAAAATATTGC‑‑AGGAGTTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtttat             >  2:2271105/7‑98 (MQ=21)
                                                                                       ggataaggcgTTCACGCCGCATCCGGCCTACAAATCATGCAAATCCAAAATATTGC‑‑AGGAGTCAT‑‑GTACGCCTGATACGCGTAGCGCATCAGGCAtttat             >  2:1952617/7‑98 (MQ=21)
                                                                                        gataaggcgTTCACGCCGCATCCGGCCTACAAAATCGCTTAAATTCAATATATTGC‑‑ACCATCCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaattct             >  1:1640765/6‑95 (MQ=25)
                                                                                             ggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctt        >  1:828002/1‑100 (MQ=255)
                                                                                               cgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg      <  2:1630585/100‑1 (MQ=255)
                                                                                               cgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg      >  1:2159031/1‑100 (MQ=255)
                                                                                                gTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgc     <  2:1876493/100‑1 (MQ=255)
                                                                                                gTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgc     <  2:1221602/100‑1 (MQ=255)
                                                                                                 ttAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgcc    >  1:134289/1‑100 (MQ=255)
                                                                                                   aaCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCAt  >  1:1967762/1‑100 (MQ=255)
                                                                                                   |                                                                                                         
GAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCAT  >  NC_000913/3705879‑3706077

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