Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,356,445 A→T 15.4% intergenic (‑34/+25) dtpC ← / ← cadA dipeptide and tripeptide permease/lysine decarboxylase, acid‑inducible

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,356,4450AT15.4% 579.2 / 88.1 234intergenic (‑34/+25)dtpC/cadAdipeptide and tripeptide permease/lysine decarboxylase, acid‑inducible
Reads supporting (aligned to +/- strand):  ref base A (114/83);  new base T (16/20);  total (131/103)
Fisher's exact test for biased strand distribution p-value = 1.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.27e-02

Only 100 of 239 total aligned reads displayed.

GAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATAT  >  NC_000913/4356346‑4356543
                                                                                                   |                                                                                                  
gAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCa                                                                                                    >  1:3154127/1‑100 (MQ=255)
    tACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGcc                                                                                                <  2:1575165/100‑1 (MQ=255)
    tACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGcc                                                                                                >  2:3287535/1‑100 (MQ=255)
    tACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTAAAAAAAAAAAAAAAAGGGAAGTGGCAAGcc                                                                                                <  1:1580005/100‑1 (MQ=255)
      cTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGcc                                                                                                <  1:2725514/98‑1 (MQ=255)
      cTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCAc                                                                                              >  1:1474716/1‑100 (MQ=255)
        cgcAGATTTGGATGGCCACGATATAGTATATCGCGCGCGGCTATGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCTTGCCACtt                                                                                            >  2:2974765/3‑100 (MQ=255)
        cccAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACtt                                                                                            >  2:3138192/1‑100 (MQ=255)
        cccAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACtt                                                                                            >  1:939214/1‑100 (MQ=255)
          cAGATTTGGATCGCCACGATATAGTCTATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGTGCAAAAAAGGGACGTGGCAAGCCCCTTcc                                                                                          >  1:2481498/1‑100 (MQ=255)
           aGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTTCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTccc                                                                                         <  1:2354458/100‑1 (MQ=255)
           aGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCTTGCCACTTccc                                                                                         <  2:3564938/100‑1 (MQ=255)
           aGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCTTGCCACTTccc                                                                                         <  2:3965818/100‑1 (MQ=255)
           aGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTACTTATGAGCAAAAAAGGGAAGTGGCTTGCCACTTccc                                                                                         <  2:4066702/100‑1 (MQ=255)
           aGATTGGGATCGCCACGATATAGAATATCGCGCGCGGATGTGAGGGTGTTTTCATGTGTTCTCATAATAAAAAAAAAAGGGAAGTGGCTTGCCACTTccc                                                                                         <  1:2960374/100‑1 (MQ=255)
            gATTGGGAGCGTCGCGATATATTATATCGCGCTCGGCTGTGGGTTTGTTTTCATGTTTTCTCCTTATGATCAAAAAAGGGAAGTGGCTTGCCACTTccc                                                                                         <  1:3640011/99‑1 (MQ=255)
              tttGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTg                                                                                      <  1:4288216/100‑1 (MQ=255)
              tttGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTAAGAGCAAAAAAGGGAAGTGGCTTGCCACTTCCCtt                                                                                       <  1:3471370/99‑1 (MQ=39)
               ttGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGt                                                                                     <  2:459425/100‑1 (MQ=255)
                tGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTa                                                                                    >  1:3972193/1‑100 (MQ=255)
                tGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTa                                                                                    <  2:450046/100‑1 (MQ=255)
                 ggATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCTTGCCACTTCCCtttt                                                                                     <  1:3906974/98‑3 (MQ=255)
                  gATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACg                                                                                  >  2:3870703/1‑100 (MQ=255)
                   aTCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCTTGCCACTTCCCttttttgc                                                                                 <  1:2022852/100‑7 (MQ=38)
                     tgCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGc                                                                               <  2:2238008/99‑1 (MQ=255)
                      gCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCt                                                                              >  1:577755/1‑100 (MQ=255)
                       ccACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTa                                                                             >  2:724117/1‑100 (MQ=255)
                           gATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAatta                                                                         <  1:631310/100‑1 (MQ=255)
                            atatAGTATATCGCGCGCGGCTGTGAGGGTGTTCTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAattat                                                                        >  2:3056467/1‑100 (MQ=255)
                               taGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTAtttt                                                                     >  2:1678307/1‑100 (MQ=255)
                                aGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTAttttt                                                                    >  1:1870801/1‑100 (MQ=255)
                                aGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTAttttt                                                                    <  1:46267/100‑1 (MQ=255)
                                aGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTAttt                                                                      >  1:1693525/1‑98 (MQ=255)
                                   atatCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGc                                                                 <  2:492910/100‑1 (MQ=255)
                                   atatCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGc                                                                 <  1:1807397/100‑1 (MQ=255)
                                    tatCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCt                                                                >  1:3304452/1‑100 (MQ=255)
                                     atCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCtt                                                               >  2:1838653/1‑100 (MQ=255)
                                     atCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCtt                                                               >  1:1207541/1‑100 (MQ=255)
                                      tCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCttt                                                              >  2:338045/1‑100 (MQ=255)
                                       cgcgcgcgGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCttt                                                              <  1:2399051/99‑1 (MQ=255)
                                       cgcgcgcgGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTttc                                                             <  1:1761037/100‑1 (MQ=255)
                                          gcgcgGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTttcttc                                                          <  1:1526570/100‑1 (MQ=255)
                                             cgGCTGTGAGGGTGGTTTCATGTGGTCTCCTTATGAGCAAAAAAGGGAAGTGGTAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCttt                                                       <  1:1370560/100‑1 (MQ=255)
                                              ggCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTc                                                      <  2:1270934/100‑1 (MQ=255)
                                               gCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCa                                                     >  2:2140389/1‑100 (MQ=255)
                                               gCTGTGAGATTGTTTTCATGTGCTCTCCTTATGAGCAAAAAATTGAAGTGTCAAGCCACTTCCCTTCTACGAGCTAATTATTCTTTGCTTTCTTCTTTCa                                                     <  1:3737769/100‑1 (MQ=255)
                                                cTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCaa                                                    <  1:2876203/100‑1 (MQ=255)
                                                 tgtgAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAAt                                                   <  2:2726867/100‑1 (MQ=255)
                                                   tgAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCACGCCACTTCCCTTGTCCGAGCTAATTATTTTTTGCTTTCTTCTTTCAATAc                                                 >  1:4202349/1‑100 (MQ=255)
                                                    gAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATAcc                                                >  1:3943521/1‑100 (MQ=255)
                                                      gggTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCtt                                              <  2:754244/100‑1 (MQ=255)
                                                        gtgtTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTaa                                            <  2:3243737/100‑1 (MQ=255)
                                                         tgtTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAAc                                           >  2:4369704/1‑100 (MQ=255)
                                                          gtTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAAc                                           >  1:2595833/1‑99 (MQ=255)
                                                          gtTTGCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTTGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCGTTCAATACCTTAACg                                          <  1:1002663/100‑1 (MQ=255)
                                                           ttttCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACgg                                         <  2:630603/100‑1 (MQ=255)
                                                             ttCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGta                                       <  1:3243887/100‑1 (MQ=255)
                                                              tCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGtat                                      <  2:1482119/100‑1 (MQ=255)
                                                              tCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGtat                                      >  2:2530999/1‑100 (MQ=255)
                                                               cATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGtata                                     >  2:1302289/1‑100 (MQ=255)
                                                                 tgtgTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGc                                   <  2:4153590/100‑1 (MQ=255)
                                                                  gtgtTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCg                                  >  2:1747900/1‑100 (MQ=255)
                                                                   tgtTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCgg                                 >  1:3251078/1‑100 (MQ=255)
                                                                    gtTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATCATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGc                                >  2:250646/1‑100 (MQ=255)
                                                                       ctcCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCAt                             >  2:2579679/1‑100 (MQ=255)
                                                                         ccTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCa                           <  1:3063240/100‑1 (MQ=255)
                                                                         ccTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCa                           >  1:3443919/1‑100 (MQ=255)
                                                                         ccTTATGAGCAAAAAAGGGAAGTGGCAAGCAACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCa                           <  2:194135/100‑1 (MQ=255)
                                                                          cTTATGAGCAAAAAAAAAAAAAAGAAAACCACTTCACAAGTACAAGCTAATTATTTTTTACTTACTTATTTCAATACATTAACGGtataac                                   >  1:1318491/1‑89 (MQ=255)
                                                                           ttATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGc                         >  1:2045906/1‑100 (MQ=255)
                                                                            tATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGcc                        <  2:4300585/100‑1 (MQ=255)
                                                                            tATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGcc                        <  2:2667884/100‑1 (MQ=255)
                                                                            tATGAAAAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGcc                        <  1:696580/100‑1 (MQ=255)
                                                                              tGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCt                       >  1:3871314/1‑99 (MQ=255)
                                                                               gAGCAAAAAAGGGAAGAGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGGTTTCTTCTTTCAATACCTTAACGGTATAGTGGCCATCAGCGTg                      >  1:4301500/1‑99 (MQ=255)
                                                                                 tcgtgcaaGGGAAGTGGCTTTCAACTTCCCTTGGAGGGGCTAATTATTTTTTGCTTTATACTTTCAATACATTACCGGTATAGCGGCCATCAGCCTGACg                   <  1:2515086/94‑1 (MQ=255)
                                                                                 gCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACg                   >  1:3258697/1‑100 (MQ=255)
                                                                                 gCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACg                   >  2:1695281/1‑100 (MQ=255)
                                                                                 gCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACg                   <  2:1565309/100‑1 (MQ=255)
                                                                                  cAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACgg                  >  2:2804042/1‑100 (MQ=255)
                                                                                   aaaaaaGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGt                 <  1:893658/100‑1 (MQ=255)
                                                                                   aaaaaaGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGt                 >  2:2884837/1‑100 (MQ=255)
                                                                                    aaaaaGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTTTTTTTTTTTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTa                >  1:1056883/1‑100 (MQ=255)
                                                                                     aaaaGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTAt               <  2:3553660/100‑1 (MQ=255)
                                                                                       aaGGGAAGTGGCTTGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGc             >  2:233215/1‑100 (MQ=255)
                                                                                       aaGGGAAGTGGCTTGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGc             >  2:1229668/1‑100 (MQ=255)
                                                                                         gggAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCAc           >  2:4025352/1‑100 (MQ=255)
                                                                                         gggAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTATTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCAc           >  2:4113780/1‑100 (MQ=255)
                                                                                          ggAAGTGGCTTGCCACTTCCCTTGTACGAGCTAATTATTTTTTTTTTTTTTCTTTCAATACACTAACGGTATAGAATCCATTAGCCTGACGGTAGGcccc          >  1:2183505/1‑100 (MQ=255)
                                                                                          ggAAGTGGCTTGCCACTTCCCTTGTACGAGCTAATTATTTTTTGTTTTTTTTTTTCAATACATTAACGGTATAGCGGCCATCAGCTTGACGGTATGCAcc          >  1:4397362/1‑100 (MQ=255)
                                                                                          ggAAGTGGCTTGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCAcc          >  2:1153067/1‑100 (MQ=255)
                                                                                          ggAAGTGGCTTGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCAc           >  2:3495201/1‑99 (MQ=255)
                                                                                          ggAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCAcc          <  1:1711710/100‑1 (MQ=255)
                                                                                           gAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCg         >  2:3510146/1‑100 (MQ=255)
                                                                                            aaGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGt        <  2:2000792/100‑1 (MQ=255)
                                                                                            aaGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGt        <  1:714109/100‑1 (MQ=255)
                                                                                             aGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTg       <  2:1351704/100‑1 (MQ=255)
                                                                                                ggCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAat    >  2:2097805/1‑100 (MQ=255)
                                                                                                ggCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAat    <  2:3444404/100‑1 (MQ=255)
                                                                                                ggCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAat    >  2:3608171/1‑100 (MQ=255)
                                                                                                   |                                                                                                  
GAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATAT  >  NC_000913/4356346‑4356543

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: