Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,705,9840TG26.2% 231.5 / 94.5 107intergenic (‑194/+114)dppB/dppAdipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base T (41/38);  new base G (26/2);  total (67/40)
Fisher's exact test for biased strand distribution p-value = 7.35e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.35e-02

Only 100 of 107 total aligned reads displayed.

GCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGT  >  NC_000913/3705887‑3706081
                                                                                                 |                                                                                                       
gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                                       >  2:297678/1‑100 (MQ=255)
gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                                       <  2:1681662/100‑1 (MQ=255)
    ccGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCg                                                                                                   >  1:821922/1‑100 (MQ=255)
     cGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCgg                                                                                                  >  1:1236090/1‑100 (MQ=255)
      ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGc                                                                                                 <  1:264075/100‑1 (MQ=255)
      ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGc                                                                                                 >  1:1702187/1‑100 (MQ=255)
       gATGGCGCTGCACGCAGAGGCACAGTGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                                <  2:677213/100‑1 (MQ=255)
             gcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACaaa                                                                                          >  2:1256370/1‑100 (MQ=255)
             gcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACaaa                                                                                          >  2:1068255/1‑100 (MQ=255)
              cTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGcataaacaa                                                                                         <  1:1763549/100‑9 (MQ=25)
              cTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCCTGGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACaaa                                                                                          >  1:58529/1‑99 (MQ=255)
                 cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCg                                                                                      <  2:1149932/100‑1 (MQ=255)
                 cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCg                                                                                      <  1:1100084/100‑1 (MQ=255)
                 cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGAATCCGGCGTTAACACCTTATCCGGCCTACAAAATCg                                                                                      >  2:592340/1‑100 (MQ=255)
                   cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTg                                                                                    <  2:734117/100‑1 (MQ=255)
                    gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGc                                                                                   <  2:706941/100‑1 (MQ=255)
                        agGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAAt                                                                               <  1:1606729/100‑1 (MQ=255)
                         ggCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAAtt                                                                              >  2:1532130/1‑100 (MQ=255)
                         ggCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAAtt                                                                              >  2:1523415/1‑100 (MQ=255)
                           cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCa                                                                            >  2:1680175/1‑100 (MQ=255)
                           cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCa                                                                            >  2:488539/1‑100 (MQ=255)
                            acaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCa                                                                            <  2:1476944/99‑1 (MQ=255)
                             caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAAt                                                                          >  2:67366/1‑100 (MQ=255)
                             caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAAt                                                                          <  2:338379/100‑1 (MQ=255)
                               gcgcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATAc                                                                        <  1:1033347/100‑1 (MQ=255)
                               gcgcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATAc                                                                        <  1:1032586/100‑1 (MQ=255)
                                cgcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATAc                                                                        <  2:1047255/99‑1 (MQ=255)
                                cgcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACa                                                                       >  2:911581/1‑100 (MQ=255)
                                    aCCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACAtt                                                                     <  2:1293695/98‑1 (MQ=255)
                                    aCCGGGGATAGAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGc                                                                   <  1:143326/100‑1 (MQ=255)
                                      cGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑Ag                                                              >  1:1762033/1‑100 (MQ=255)
                                      cGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑Ag                                                              <  2:541449/100‑1 (MQ=255)
                                       gggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGa                                                             >  1:503843/1‑100 (MQ=255)
                                        ggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGaa                                                            >  2:1191166/1‑100 (MQ=255)
                                         gTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAAt                                                           <  2:170582/100‑1 (MQ=255)
                                          tataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCGGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAAtt                                                         <  1:620334/100‑1 (MQ=255)
                                          tataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAAtt                                                         <  1:595451/100‑1 (MQ=255)
                                              cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑TTCg                                                   <  1:1744610/100‑1 (MQ=255)
                                                   tgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑TTCGTAGGc                                              <  2:1740260/100‑1 (MQ=255)
                                                    gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑TTCGTAGGcc                                             >  1:103363/1‑100 (MQ=255)
                                                    gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑TTCGTAGGcc                                             >  1:1169348/1‑100 (MQ=255)
                                                       gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑TTCGTAGGCCTGa                                          >  2:1340399/1‑100 (MQ=255)
                                                       gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑TTCGTAGGCCTGa                                          >  1:1683413/1‑100 (MQ=255)
                                                       gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑TTCGTAGGCCTGa                                          <  2:1702200/100‑1 (MQ=255)
                                                         cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑TTCGTAGGCCTGATa                                        >  2:901313/1‑100 (MQ=255)
                                                          aaaGTGCTCTTTGTTTATGCCGGATGCGGCGTTAAAGGCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑TTCGTAGGCCTGATa                                        >  1:1055188/1‑99 (MQ=255)
                                                          aaaGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑TTCGTAGGCCTGATaa                                       >  1:1055189/1‑100 (MQ=255)
                                                            aGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑TTCGTAGGCCTGATAAGc                                     <  1:520659/100‑1 (MQ=255)
                                                            aGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑TTCGTAGGCCTGATAAGc                                     >  1:12707/1‑100 (MQ=255)
                                                                cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑TTCGTAGGCCTGATAAGCGGAg                                 <  1:679691/100‑1 (MQ=255)
                                                                cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑TTCGTAGGCCTGATAAGCGGAg                                 <  1:1536712/100‑1 (MQ=255)
                                                                cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑TTCGTAGGCCTGATAAGCGGAg                                 >  1:1010374/1‑100 (MQ=255)
                                                                 gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑TTCGTAGGCCTGATAAGCGGAgc                                <  1:363970/100‑1 (MQ=255)
                                                                  tgatttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAATCATGCAAATCCAAAATATTGC‑‑AGGAGT‑‑TATGTAGGCCTGATAAGCGTAgcgc                              >  1:551126/11‑100 (MQ=18)
                                                                  cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑TTCGTAGGCCTGATAAGCGGAgcg                               >  2:1548186/1‑100 (MQ=255)
                                                                  cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑TTCGTAGGCCTGATAAGCGGAgcg                               >  2:1438777/1‑100 (MQ=255)
                                                                    ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑TTCGTAGGCCTGATAAGCGGAGCGCa                             >  1:10159/1‑100 (MQ=255)
                                                                    ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑TTCGTAGGCCTGATAAGCGGAGCGCa                             <  2:1037093/100‑1 (MQ=255)
                                                                      gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCTTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑TTCGTAGGCCTGATAAGCGGAGCGCATc                           <  2:1502798/100‑1 (MQ=255)
                                                                       tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCa                          <  2:1144574/100‑1 (MQ=255)
                                                                       tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCa                          <  1:59539/100‑1 (MQ=255)
                                                                       tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCa                          <  1:874024/100‑1 (MQ=255)
                                                                        gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAACATGCAAATTCAATAGATTGC‑‑AGAGAT‑‑TATGTAGGCCTGATAAGCGTAGCGCATCAg                         >  1:311414/5‑100 (MQ=17)
                                                                         taggCCGGATAAGGCGTTAACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAgg                        >  1:424834/4‑100 (MQ=18)
                                                                          aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGC‑‑ATGGGCC‑‑ATGTAGGCCTGATAAGCGTAGCGCATCAGGc                       >  1:576835/3‑100 (MQ=17)
                                                                          aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGC‑‑AGAGTT‑‑CACGTAGGCCTGATAAGCGTAGCGCATCAGGc                       >  1:1304631/3‑100 (MQ=17)
                                                                          aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGC‑‑AGAGTT‑‑CACGTAGGCCTGATAAGCGTAGCGCATCAAGc                       >  2:1127876/3‑100 (MQ=14)
                                                                          aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAACAGTGCAAATTCAATATATTGC‑‑AATCTT‑‑CTGGTAGGCCTGATAAGCATCGCGCATCAGGc                       >  2:1387766/3‑100 (MQ=18)
                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATTGTGTAAATTCAATAAATTGA‑‑TCGATC‑‑CACGTAGGCCTGATAAGCGTAGCGCATCAgg                        >  2:80290/3‑99 (MQ=18)
                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAACCTGCAAATTCAATATATTGC‑‑AGGATA‑‑CGTGTAGGCCTGATAAGCGTAGCGCATCAGGc                       >  2:863322/3‑100 (MQ=17)
                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAATCATGCAAATCCAAAATATTGC‑‑AGGAGT‑‑TATGTAGGCCTGATAAGCGTAGCGCATCAGGCa                      >  1:1016666/3‑100 (MQ=21)
                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAATCATGCAAATCCAAAATATTGC‑‑AGGAGT‑‑TATGTAGGCCTGATAAGCGTAGCGCATCAGGCa                      >  1:1634581/3‑100 (MQ=21)
                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAATCATGCAAATCCAAAATATTGC‑‑AGGAGT‑‑TATGTAGGCCTGATAAGCGTAGCGCATCAGGCa                      >  1:371322/3‑100 (MQ=21)
                                                                           tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCa                      >  1:1423901/1‑100 (MQ=255)
                                                                           tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCa                      >  2:353382/1‑100 (MQ=255)
                                                                           ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAA‑‑CTTGTAGGCCTGATAAGCGCAGCGCATCAGGc                       >  2:1356032/2‑100 (MQ=18)
                                                                           ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGC‑‑AATTCAACTAATAGGCCTGATAAGCGCAGCGCATCAGGc                       >  2:1382502/2‑100 (MQ=18)
                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATCATGAAAATTCAATACATTGC‑‑AAGATT‑‑TTCGTAGGCCTGATAAGAGTAGCGCATCAGGCa                      >  2:1792775/2‑100 (MQ=17)
                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTTTGCAAATTCAATATATTGC‑‑AATCTC‑‑CGTGTAGGCCTGATAAGCGTAGCGCATCAGGCa                      >  2:972699/2‑100 (MQ=18)
                                                                            gtcggataaggcgTTCACGCCGCATCCGACCCACAAAGTCTTGCAAAATCAATAAATTGC‑‑AGGAAC‑‑GATGTAGGCCTGATAAGCGTAGCGCATCAGGCaa                     >  2:1383905/10‑99 (MQ=21)
                                                                            gCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGC‑‑AGAAAC‑‑CTTGTAGGACTGATAAGAGTAGCGCATAAggaa                      >  2:785541/1‑97 (MQ=17)
                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACGAAATCGTGCTAATTCAATATATTGT‑‑GGGCTC‑‑CACGTAGGCCTGATAAGCGTAGCGCATCAGGCag                     >  2:1642074/1‑99 (MQ=18)
                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGC‑‑AGTTTA‑‑TTGGTAGGCCTGATAAGCGCAGCGCATCAGGCaa                     >  1:466815/1‑99 (MQ=21)
                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAACCGCGCAAATTCAATATATTGCGGAGAAAA‑‑T‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGctgt                    >  2:1815709/1‑97 (MQ=21)
                                                                            aaagcacctggcgGTCACGCCGCATCCGGCCTACGAAACCGCAAGAATTCAATATATTGC‑‑AGGAGC‑‑GGTGTAGGCCTGATAAGCGTAGCGCATCAAGCag                     <  1:799084/91‑2 (MQ=18)
                                                                             ccggataaggcgTTTACGCCGCATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGC‑‑ATGGGCC‑‑ATGTAGGCCTGATAAGCGTAGCGCATCAGGCaat                    >  1:1461589/9‑98 (MQ=17)
                                                                              cGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAtt                    >  2:729353/1‑99 (MQ=255)
                                                                               ggataaggcgTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGC‑‑AGAGTT‑‑CACGTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                  >  1:1562793/7‑100 (MQ=17)
                                                                               ggataaggcgTTTACGCCGCAACCGGCCTACAAATTCTGCAAATTCAACAAATTGC‑‑GAATCC‑‑CTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAtttg                  >  1:1809809/7‑98 (MQ=18)
                                                                               ggataaggcgTTCACGCCTTATCCGGCCTACAAAATCGCTTAAATTCAATATATTGC‑‑ACAATA‑‑CTCGTAGGCCTAATAAGCGTAGAGCATCAGGCaat                    >  2:228972/7‑96 (MQ=17)
                                                                               ggataaggcgTTCACGCCGCATCCGGCCTACAAAATCATGAAAATTCAATACATTGC‑‑AAGATT‑‑TTCGTAGGCCTGATAAGCGTAGCGCATCAGGCacgct                  >  2:1450541/7‑96 (MQ=17)
                                                                               ggataaggcgTTCACGCCGCATCCGGCCTACAAAACTGTGCAAATTCAATATATTGC‑‑AGGAAA‑‑CACGTAGGCCTGATAAGCGAAGCGCATCAGGCAgttt                  >  2:128892/7‑100 (MQ=18)
                                                                                     ggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctt            >  2:142002/1‑100 (MQ=255)
                                                                                       cgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg          <  1:1523404/100‑1 (MQ=255)
                                                                                       cgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg          <  1:138339/100‑1 (MQ=255)
                                                                                         ttAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgcc        <  1:436529/100‑1 (MQ=255)
                                                                                          tAACGCCTTATCCGACCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa       >  1:1143213/1‑100 (MQ=255)
                                                                                            aCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc     <  2:1449394/100‑1 (MQ=255)
                                                                                             cGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCa    <  2:37410/100‑1 (MQ=255)
                                                                                                cTTATCCGGCCTACAAAATCGTGCAACTTCAATACATTGC‑‑AGAATT‑‑TTCGTAGGCCTGATAAGCGGAACGCATCAAACATATTTGCCTTTGCCATAAGt  >  2:179884/1‑99 (MQ=18)
                                                                                                 |                                                                                                       
GCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGT  >  NC_000913/3705887‑3706081

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: