Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,705,984 T→G 30.1% intergenic (‑194/+114) dppB ← / ← dppA dipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,705,9840TG30.1% 496.9 / 277.0 287intergenic (‑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 (93/107);  new base G (76/10);  total (169/118)
Fisher's exact test for biased strand distribution p-value = 6.48e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.29e-05

Only 100 of 289 total aligned reads displayed.

TTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCT  >  NC_000913/3705885‑3706083
                                                                                                   |                                                                                                        
  gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCATTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                                        >  1:1993627/1‑100 (MQ=255)
    tgcCGGATGGCGCTGCACGCAGAGGCATAGCGTCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATc                                                                                                      <  1:522855/100‑1 (MQ=255)
    tgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATc                                                                                                      <  1:2067364/100‑1 (MQ=255)
     gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCATTGTTTATACCGGATGCGACGTTAACGCCTTATcc                                                                                                     >  1:1943126/1‑100 (MQ=255)
       cGGATGGCGCTGTACGCAGAGTCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCgg                                                                                                   <  2:3680894/100‑1 (MQ=255)
       cGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCACATTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCgg                                                                                                   >  1:2094199/1‑100 (MQ=255)
        ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGc                                                                                                  >  1:1585419/1‑100 (MQ=255)
        ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTGATCCGGc                                                                                                  <  2:1118786/100‑1 (MQ=255)
        ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGc                                                                                                  <  2:3618483/100‑1 (MQ=37)
         gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                                 <  2:251542/100‑1 (MQ=255)
          aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                                 <  1:3780093/99‑1 (MQ=255)
             gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGCATaaa                                                                                             <  1:2142847/100‑3 (MQ=25)
              cgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑aa                                                                                           >  2:779110/1‑100 (MQ=255)
               gcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑aaa                                                                                          >  2:1200600/1‑100 (MQ=255)
                cTGCACGCAGAGGCACAGCGCCACAGGATATACAGATTGTTGACAAAGGCCGCGTCGTTTACGCCGCATGCGGCGTTAACGCCTCATCCGGACTAC‑aaa                                                                                          >  2:1713954/1‑99 (MQ=18)
                  gCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATc                                                                                       >  1:908798/1‑100 (MQ=255)
                   cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCg                                                                                      <  1:415612/100‑1 (MQ=255)
                     cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTg                                                                                    >  2:2074319/1‑100 (MQ=255)
                      gCAGAGGCACAGCTCTACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGc                                                                                   <  1:3488582/100‑1 (MQ=255)
                      gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGc                                                                                   <  2:2838934/100‑1 (MQ=255)
                      gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGc                                                                                   <  2:3339669/100‑1 (MQ=255)
                           ggCACAGCGCCACCGGGTATACAGATTGTTTACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAAtt                                                                              <  1:1107384/100‑1 (MQ=255)
                           ggCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAAtt                                                                              >  2:3782376/1‑100 (MQ=255)
                            gCACAGCGCAACCGGGTATACAGATTGTTTACAAAGTGTGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAAtcc                                                                             <  1:2655158/100‑3 (MQ=18)
                             cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCa                                                                            <  2:414433/100‑1 (MQ=255)
                             cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCa                                                                            >  2:1227680/1‑100 (MQ=255)
                             cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCa                                                                            <  1:1743513/100‑1 (MQ=255)
                                 gcgcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATAc                                                                        >  2:1802213/1‑100 (MQ=255)
                                 gcgcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATAc                                                                        >  2:571346/1‑100 (MQ=255)
                                   gcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACAt                                                                      >  1:781340/1‑100 (MQ=255)
                                    ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACAtt                                                                     <  1:3307646/100‑1 (MQ=255)
                                     cACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTg                                                                    >  1:2690150/1‑100 (MQ=255)
                                       ccGGTTATACAGATTGTTGACAAATTTCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGc                                                                   <  1:3398284/99‑1 (MQ=255)
                                        cGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAg                                                                 <  1:2261436/100‑1 (MQ=255)
                                          ggcataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGTAGa‑a                                                              >  2:2942291/4‑100 (MQ=255)
                                               aCAGATTGTTGACAAAGTGCGCTATGTTTATGCCGGATGCAGAGTTAACGCCTAATCCGGCCTAA‑AAAATCGTGCAAATTAAATATATTGCAGA‑AT‑TTTc                                                        >  2:860138/1‑100 (MQ=17)
                                                cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑g                                                     <  1:2027931/100‑1 (MQ=255)
                                                  gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTa                                                   <  1:1683849/100‑1 (MQ=255)
                                                     tgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGc                                                <  1:1906276/100‑1 (MQ=255)
                                                      gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGcc                                               <  1:2914105/100‑1 (MQ=255)
                                                       ttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATTGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCt                                              >  2:48179/1‑100 (MQ=255)
                                                       ttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCt                                              >  2:1525194/1‑100 (MQ=255)
                                                       atctcacaGTCCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTCTCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCt                                              <  2:3271554/93‑1 (MQ=18)
                                                         gACAAAGTGCGCTTTGTTTATGCCGGATGCGGTGTTAACACCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGa                                            <  2:2752734/100‑1 (MQ=255)
                                                         gACAAAGTGCGCTGTGTTTATGCCGGATGCGGAGTAAACGCCTTATCGGGCGTAG‑AAAATCGTGGAAATTCAACATATGGCACA‑GT‑TATT‑‑GTAgg                                                 >  2:2468206/1‑95 (MQ=14)
                                                          aCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGAt                                           <  1:3460860/100‑1 (MQ=255)
                                                           cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCtagat                                          >  1:3271586/1‑96 (MQ=255)
                                                           cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATa                                          <  2:2889410/100‑1 (MQ=255)
                                                               gTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCg                                      <  1:246085/100‑1 (MQ=255)
                                                                tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAAACGAAAAAATTAAATACATTACAGA‑AA‑TTTC‑‑GTAGGCCTGATAAGCg                                      >  1:230252/1‑99 (MQ=17)
                                                                   gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAgc                                  >  2:2302199/1‑100 (MQ=255)
                                                                   gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAgc                                  <  2:3330813/100‑1 (MQ=255)
                                                                    cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAgcg                                 >  1:1172899/1‑100 (MQ=255)
                                                                    cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAgcg                                 >  1:1099675/1‑100 (MQ=255)
                                                                    cTTTGTTTATGACGGATGCGGCGTTAACGCCATATCCAGTCTAA‑AAAATCGTGCAAATTGAATACATTGCAGA‑AT‑TTTT‑‑GTAGGCGTGATAAGCGGAgt                                  <  2:587324/99‑2 (MQ=25)
                                                                     tttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAgcgc                                >  1:2432653/1‑100 (MQ=255)
                                                                     tttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAgcgc                                <  1:1822863/100‑1 (MQ=255)
                                                                     cttttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAGCATGAAAATTCAATATATTGCAGG‑AG‑CTGC‑‑GTAGGCCTGATAAGCGTAgcgc                                >  1:2784324/10‑100 (MQ=17)
                                                                      ttttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAGG‑AC‑CAT‑‑‑GTAGGCCTGATAAGCGTAGCGCa                               >  1:3615220/9‑100 (MQ=18)
                                                                       tGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCAt                              >  2:1263790/1‑100 (MQ=255)
                                                                        ttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAACATGCAAATTCAATAGATTGCAGAGAT‑TAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATc                             <  1:2862793/94‑1 (MQ=17)
                                                                        gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑ATTC‑‑GTAGGCCTGATAAGCGGAGCGCATc                             >  1:2634574/1‑100 (MQ=255)
                                                                         tgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGCAGA‑GT‑TCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCa                            >  1:691795/6‑100 (MQ=14)
                                                                          gtaggCCGGATAAGGCGTTAACGCCGCATCCGGCATAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAg                           >  1:3028124/5‑100 (MQ=18)
                                                                           taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑CAAACTTTGCAAATTCAATATATTGC‑‑A‑AT‑CTCCGTGTAGGCCTGATAAGCGTAGCGCATCAgg                          >  2:2095046/4‑100 (MQ=18)
                                                                            aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCATGAAAATTCAATACATTGCAAG‑AT‑TTTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                         >  1:1522145/3‑100 (MQ=14)
                                                                            aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACCTTGCTAATTCAATATATTGCAGG‑GA‑CTAT‑‑GTAGGCCTGATAAGCATAGCGCATCAGGc                         >  1:3165231/3‑100 (MQ=18)
                                                                             tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCa                        >  2:3561352/1‑100 (MQ=255)
                                                                             ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATGTTGCAGA‑AA‑CCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                        >  2:2504201/2‑100 (MQ=18)
                                                                             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAGG‑AC‑CAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                        >  2:2806148/2‑100 (MQ=18)
                                                                             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAGA‑A‑‑CCAA‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                        >  1:858175/2‑100 (MQ=17)
                                                                             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAGTTCGTGCAAATTCAATAAATTGCA‑A‑TA‑TGAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGcga                       >  1:2364922/2‑98 (MQ=18)
                                                                             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATTGTGTAAATTCAATAAATTGATCG‑AT‑CCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                        >  1:2602120/2‑100 (MQ=21)
                                                                             ggCCGGATAAGGCGTTAACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCACATCAGGCa                        >  2:2172666/2‑100 (MQ=18)
                                                                               ccggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGCACA‑AT‑CCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaat                      >  1:1647455/9‑98 (MQ=18)
                                                                               ccGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtt                      >  1:3102731/1‑100 (MQ=255)
                                                                                cggataaggcgTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑GA‑ATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtttg                    >  1:2950735/8‑99 (MQ=18)
                                                                                cGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAttt                     >  1:60220/1‑100 (MQ=255)
                                                                                 ggataaggcgTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑GA‑ATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtttgc                   >  2:1933976/7‑98 (MQ=18)
                                                                                 ggataaggcgTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAGG‑AC‑CAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                    >  1:1605903/7‑100 (MQ=18)
                                                                                 ggataaggcgTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAGG‑AC‑CAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                    >  2:314289/7‑100 (MQ=18)
                                                                                 ggataaggcgTTCACGCCGCATACGGCCTAC‑AAAATTGTGTAAATTGAATAAATTGATCG‑AT‑CCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                    >  2:2683186/7‑100 (MQ=18)
                                                                                 agataaggcgTTTCCGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTACATATATTGCAGA‑AA‑CCTT‑‑GTCGGCATTATAAGCGTAGCACATCAGGCa                        >  2:2206371/7‑96 (MQ=12)
                                                                                  gataaggcgTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGCACC‑AT‑CCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaattct                   >  2:1985494/6‑95 (MQ=25)
                                                                                  gataaggcgTTCACGCCGCATCCGGCCTAC‑‑AAATCATGCAAATCCAAAATATTGCAGG‑AG‑TTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtttatg                  >  1:1355078/6‑97 (MQ=21)
                                                                                   aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtttt                    >  2:1622565/1‑98 (MQ=255)
                                                                                   aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                  <  1:2571630/100‑1 (MQ=255)
                                                                                   aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                  <  1:1197727/100‑1 (MQ=255)
                                                                                      cggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGTCATTtttgcct               >  2:3175897/1‑100 (MQ=255)
                                                                                        gcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgccctt             <  2:1021377/100‑1 (MQ=255)
                                                                                         cgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg            <  1:2405534/100‑1 (MQ=255)
                                                                                          gTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAGTTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgc           <  1:1265379/100‑1 (MQ=255)
                                                                                            tAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTCCCa         >  2:3023718/1‑100 (MQ=255)
                                                                                            tAACGCCTAATCCGGCCTAC‑AAAATCGTACAAATTCAACACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGCGCACATCAGGCATCTTTTCCTTTGCCa         >  2:3916067/1‑100 (MQ=18)
                                                                                             aaCGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCAt        >  2:1431861/1‑100 (MQ=255)
                                                                                             aaCGCCTTATCCGACCTAC‑AAAATCGTGCAAATTCAATTCATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCAt        >  2:2408127/1‑100 (MQ=255)
                                                                                              aCGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc       <  1:2239539/100‑1 (MQ=255)
                                                                                               cGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGTAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCa      <  1:46533/100‑1 (MQ=255)
                                                                                               cGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATAAGGCATTTTTGCCTTTGCCATCa      >  1:1015696/1‑100 (MQ=255)
                                                                                                gCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAg     <  1:1463593/100‑1 (MQ=255)
                                                                                                   |                                                                                                        
TTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCT  >  NC_000913/3705885‑3706083

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