Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,705,985 T→C 41.0% intergenic (‑195/+113) 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,9850TC41.0% 323.1 / inf 322intergenic (‑195/+113)dppB/dppAdipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base T (97/93);  new base C (122/10);  total (219/103)
Fisher's exact test for biased strand distribution p-value = 1.63e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.96e-03

Only 100 of 323 total aligned reads displayed.

GCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CAGA‑AT‑T‑‑‑T‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCT  >  NC_000913/3705887‑3706083
                                                                                                  |                                                                                                          
gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGCGCGCTGTGTTTATGACGGATGCGGCGTTAACGCCTTa                                                                                                           >  2:2825249/1‑100 (MQ=255)
 ctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTAt                                                                                                          >  2:3816725/1‑100 (MQ=255)
   gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                                        <  2:2707694/100‑1 (MQ=255)
     cGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCgg                                                                                                      >  1:146899/1‑100 (MQ=255)
      ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGc                                                                                                     <  1:451843/100‑1 (MQ=37)
       gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                                    <  2:1047547/100‑1 (MQ=255)
         tGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGCATa                                                                                                  <  2:266267/100‑1 (MQ=37)
          ggCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAc                                                                                                 >  2:2122988/1‑100 (MQ=255)
          ggCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGCATaa                                                                                                 <  2:1746596/100‑2 (MQ=25)
           gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑a                                                                                               >  2:305772/1‑100 (MQ=255)
             gcTGCACGCAGAGGCACAGCGCCACCGGGTATCCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑aaa                                                                                             >  1:3975185/1‑100 (MQ=255)
             gcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑aaa                                                                                             <  2:3473602/100‑1 (MQ=255)
               tGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAAt                                                                                           <  1:712390/100‑1 (MQ=255)
                   cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTg                                                                                       >  2:3014917/1‑100 (MQ=255)
                   cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTg                                                                                       <  1:2092396/100‑1 (MQ=255)
                    gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGc                                                                                      <  2:3953558/100‑1 (MQ=255)
                    gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGc                                                                                      >  1:1012826/1‑100 (MQ=255)
                    gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGc                                                                                      <  1:948507/100‑1 (MQ=255)
                    gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGc                                                                                      <  1:2334001/100‑1 (MQ=255)
                      agagGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCaa                                                                                    <  1:1961493/100‑1 (MQ=255)
                             caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAAt                                                                             <  1:1799076/100‑1 (MQ=255)
                              aGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATa                                                                            <  2:1945320/100‑1 (MQ=255)
                                     acGGGTATACAAATAGTTCACAGAGTGCGCTTAGTTTCTGCAGGATGCGGCGTTAAAGCCTTATCCGACCTAC‑AAAGTCGTTCAAGTTCAATACATTG‑Ca                                                                    >  2:1161765/2‑100 (MQ=21)
                                         gTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑c                                                                     <  2:3732767/95‑1 (MQ=255)
                                              cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CAGA‑AT‑T‑‑‑T‑TCg                                                     <  2:2642480/100‑1 (MQ=255)
                                                gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CAGA‑AT‑T‑‑‑T‑TCGTa                                                   >  1:3090983/1‑100 (MQ=255)
                                                   tgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CAGA‑AT‑T‑‑‑T‑TCGTAGGc                                                <  1:860543/100‑1 (MQ=255)
                                                    gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CAGA‑AT‑T‑‑‑T‑TCGTAGGcc                                               >  1:1093793/1‑100 (MQ=255)
                                                    gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTACTCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CAGA‑AT‑T‑‑‑T‑TCGCAGGcc                                               <  2:1686147/100‑1 (MQ=255)
                                                      tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CAGA‑AT‑T‑‑‑T‑TCGTAGGCCTg                                             >  2:4065163/1‑100 (MQ=255)
                                                       gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CAGA‑AT‑T‑‑‑T‑TCGTAGGCCTGa                                            >  1:1554201/1‑100 (MQ=255)
                                                       gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CAGA‑AT‑T‑‑‑T‑TCGTAGGCCTGa                                            >  1:3346544/1‑100 (MQ=255)
                                                         cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CAGA‑AT‑T‑‑‑T‑TCGTAGGCCTGATa                                          <  2:2174638/100‑1 (MQ=255)
                                                            aGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CAGA‑AT‑T‑‑‑T‑TCGTAGGCCTGATAAGc                                       >  2:3349954/1‑100 (MQ=255)
                                                              tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CAGA‑AT‑T‑‑‑T‑TCGTAGGCCTGATAAGCgg                                     <  1:2174010/100‑1 (MQ=255)
                                                                cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CAGA‑AT‑T‑‑‑T‑TCGTAGGCCTGATAAGCGGAg                                   <  2:512306/100‑1 (MQ=255)
                                                                cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACACTG‑CAGA‑AT‑T‑‑‑T‑TCGTAGGTCTGATAAGCgggg                                   <  1:2825256/100‑3 (MQ=255)
                                                                 gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CAGA‑AT‑T‑‑‑T‑TCGTAGGCCTGATAAGCGGAgc                                  <  1:2467864/100‑1 (MQ=255)
                                                                    ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CAGA‑AT‑T‑‑‑T‑TCGTAGGCCTGATAAGCGGAGCGCa                               <  2:36542/100‑1 (MQ=255)
                                                                     tGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CAGA‑AT‑T‑‑‑T‑TCGTAGGCCTGATAAGCGGAGCGCAt                              >  1:3657878/1‑100 (MQ=255)
                                                                     tGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CAGA‑AT‑T‑‑‑T‑TCGTAGGCCTGATAAGCGGAGCGCAt                              >  2:1606354/1‑100 (MQ=255)
                                                                     tGTCTATGCCAGAGGGGGCCTTAACGGCGTATCCGGCCTAC‑AAGAGCGTGCAAATTCAATTCGTTG‑CAGA‑AG‑T‑‑‑T‑TCGTAGGCCTGATAAGCGGTgc                                  <  2:202789/96‑1 (MQ=12)
                                                                      gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CAGA‑AT‑T‑‑‑T‑TCGTAGGCCTGATAAGCGGAGCGCATc                             >  1:484109/1‑100 (MQ=255)
                                                                       tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CAGA‑AT‑T‑‑‑T‑TCGTAGGCCTGATAAGCGGAGCGCATCa                            <  2:3456550/100‑1 (MQ=255)
                                                                       tgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTG‑CAGA‑GT‑T‑‑‑C‑ACGTAGGCCTGATAAGCGTAGCGCATCa                            >  2:2963406/6‑100 (MQ=17)
                                                                        ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CAGA‑AT‑T‑‑‑T‑TCGTAGGCCTGATAAGCGGAGCGCATCAg                           <  1:3412513/100‑1 (MQ=255)
                                                                        gtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAACATATTG‑CA‑A‑TTCT‑‑‑C‑TTGTAGGCCTGATAAGCGTAGCGCATCAg                           >  1:252764/5‑100 (MQ=18)
                                                                          aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTG‑CAGA‑AA‑C‑‑‑C‑TTGTAGGCCTGATAAGCGTAGCGCATCAGGc                         >  1:2978924/3‑100 (MQ=18)
                                                                          aTGACGGATAAGGCGTTCACGCCGCATCCGTCCTAC‑AAGTTCATGCAATTTTAATAAGTTGCCAGA‑TA‑T‑‑‑C‑CTGTAGGCCTGATAAGCGTAGCGCATCAgg                          >  2:1527434/1‑100 (MQ=25)
                                                                           tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CAGA‑AT‑T‑‑‑T‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGCa                        >  2:1385350/1‑100 (MQ=255)
                                                                           ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTG‑C‑GA‑ATCC‑‑‑C‑TTGTAGGCCTGATAAGCGTAGCGCATCAGGCAt                       >  1:3628725/2‑100 (MQ=18)
                                                                           ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTG‑C‑GA‑ATCC‑‑‑C‑TTGTAGGCCTGATAAGCGTAGCGCATCAGGCAt                       >  2:1352664/2‑100 (MQ=18)
                                                                           ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTG‑C‑GA‑AT‑C‑‑‑CCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAt                       >  1:864509/2‑100 (MQ=18)
                                                                           ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTG‑CAGA‑AA‑C‑‑‑C‑TTGTAGGCCTGATAAGCGTAGCGCATCAGGCa                        >  1:3237597/2‑100 (MQ=17)
                                                                           ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTG‑CAGA‑GT‑T‑‑‑C‑ACGTAGGCCTGATAAGCGTAGCGCATCAGGCa                        >  1:3904255/2‑100 (MQ=17)
                                                                           ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTG‑CAGA‑GT‑T‑‑‑C‑ACGTAGGCCTGATAAGCGTAGCGCATCAGGCa                        >  2:1436629/2‑100 (MQ=17)
                                                                           ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAATGTGCAAATTCAATAAATTG‑C‑‑A‑AT‑TCAAC‑TTGTAGGCCTGATAAGCGCAGCGCATCAGGc                         >  1:3585531/2‑100 (MQ=21)
                                                                           ggCCGGATAAGGCGTTTACGCCGCATCCGGCATAC‑AAAATCGTGCAAATTCAACATATTG‑CAAT‑TC‑T‑‑‑C‑TTGTAGGCCTGATAAGCGTAGCGCATCAgg                          >  2:776051/2‑98 (MQ=18)
                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑GAAATCGTGCTAATTCAATATATTG‑TGGG‑CT‑C‑‑‑C‑ACGTAGGCCTGATAAGCGTAGCGCATCAGGCa                        >  1:4050688/2‑100 (MQ=21)
                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑GAAATCGTGCTAATTCAATATATTG‑TGGG‑CT‑C‑‑‑C‑ACGTAGGCCTGATAAGCGTAGCGCATCAGGCa                        >  1:1502194/2‑100 (MQ=18)
                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCATGAAAATTCAATACATTG‑CAAG‑AT‑T‑‑‑T‑TCGTAGGCCTGATAAGCGTAGCGCATCAGGCa                        >  1:1881861/2‑100 (MQ=17)
                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCATAC‑AAAACTGTGCAAATTCAATATATTG‑CAGG‑AA‑A‑‑‑C‑ACGTAGGCCTGATAAGCGAAGCGCATCAGGCa                        >  1:1060097/2‑100 (MQ=18)
                                                                            gCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CAGA‑AT‑T‑‑‑T‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGCAt                       >  1:2650065/1‑100 (MQ=255)
                                                                            gCCGGATAAGGCGTTCGCGCCGCATCCGGCCTAC‑AAAAGTTTGCAAATTCAATAAATTG‑CAGA‑AT‑T‑‑‑C‑ATGTAGGCCTGATAAGCGAAGCGCATCAGGCAt                       >  2:4191360/1‑100 (MQ=17)
                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAGTTCGTGCAAATTCAATAAATTG‑CA‑A‑TA‑T‑‑‑G‑ACGTAGGCCTGATAAGCGTAGCGCATCAGGcgat                      >  1:1712766/1‑97 (MQ=21)
                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAACATGCAAATTCAATAGATTG‑CAGAGAT‑T‑‑‑A‑T‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa                       >  1:1199639/1‑99 (MQ=17)
                                                                            gCCGGATAAGGCGTTAACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CAGA‑AT‑T‑‑‑T‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGc                         >  2:3690750/1‑98 (MQ=17)
                                                                             ccggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTG‑CAGA‑AA‑C‑‑‑C‑TTGTAGGCCTGATAAGCGTAGCGCATCAGGCagt                      >  1:2446499/9‑98 (MQ=18)
                                                                             ccggataaggcgTTCACGCCGCATCCGGCCTAC‑GAAATCGTGCTAATTCAATATATTG‑TGGG‑CT‑C‑‑‑C‑ACGTAGGCCTGATAAGCGTAGCGCATCAGGCagt                      >  2:254173/9‑98 (MQ=18)
                                                                             ccggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAATTCAATGTATTG‑CAGT‑TT‑A‑‑‑T‑TGGTAGGCCTGATAAGCGCAGCGCATCAGGCaat                      >  1:2244663/9‑98 (MQ=25)
                                                                             ccggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTG‑CACC‑AT‑C‑‑‑C‑TCGTAGGCCTGATAAGCGTAGCGCATCAGGCaat                      >  2:914270/9‑98 (MQ=21)
                                                                             ccGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CAGA‑AT‑T‑‑‑T‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGCAtt                      <  2:2502378/100‑1 (MQ=255)
                                                                             ccGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CAGA‑AT‑T‑‑‑T‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGCAtt                      >  1:1329031/1‑100 (MQ=255)
                                                                             ccGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CAGA‑AT‑T‑‑‑T‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGCAtt                      >  1:638783/1‑100 (MQ=255)
                                                                             ccGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CAGA‑AT‑T‑‑‑T‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGCAtt                      >  1:4144737/1‑100 (MQ=255)
                                                                              cGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CAGA‑AT‑T‑‑‑T‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGCAttt                     <  2:3259963/100‑1 (MQ=255)
                                                                               ggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTG‑CAGA‑AA‑C‑‑‑C‑TTGTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                    >  2:2346684/7‑100 (MQ=18)
                                                                               ggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTG‑CAGA‑AA‑C‑‑‑C‑TTGTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                    >  2:1209954/7‑100 (MQ=18)
                                                                               ggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTG‑CAGA‑GT‑T‑‑‑C‑ACGTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                    >  2:3983671/7‑100 (MQ=17)
                                                                               ggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTG‑CAGA‑GT‑T‑‑‑C‑ACGTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                    >  1:3278712/7‑100 (MQ=17)
                                                                               ggataaggcgTTCACGCCTTATCCTGCCTAC‑AAAAGCACGCAAATTCAATATATTG‑CAGA‑GA‑T‑‑‑G‑GCGTAGGCCTGATAAGCGTAGCGCATCAGGCAtttt                    >  1:3540081/7‑100 (MQ=14)
                                                                               ggataaggcgTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTG‑CAGG‑AC‑C‑‑‑A‑T‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                    >  1:3879331/7‑100 (MQ=18)
                                                                               ggataaggcgTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTG‑CAGG‑AC‑C‑‑‑A‑T‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                    >  1:2400924/7‑100 (MQ=18)
                                                                               ggataaggcgTTCACGCCGCATCCGGCCTAC‑GAAATCGTGCTAATTCAATATATTG‑TGGG‑CT‑C‑‑‑C‑ACGTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                    >  1:2149184/7‑100 (MQ=21)
                                                                               ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAGATCGTGCAAATTCAATATATTG‑CATG‑AG‑C‑‑‑C‑ACGTAGGCCTGATAAGCGTAGCACATAAGGCAgttt                    >  1:745872/7‑100 (MQ=17)
                                                                               ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAATTCAATGTATTG‑CAGT‑TT‑A‑‑‑T‑TGGTAGGCCTGATAAGCGCAGCGCATCAGGCaatct                    >  1:2568786/7‑96 (MQ=25)
                                                                               ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTG‑CACC‑AT‑C‑‑‑C‑TCGTAGGCCTGATAAGCGTAGCGCATCAGGCaattc                    >  1:549514/7‑96 (MQ=21)
                                                                               ggataaggcgTTCACGCCGCATCCGACCCAC‑AAAGTCTTGCAAAATCAATAAATTG‑CAGG‑AA‑C‑‑‑G‑ATGTAGGCCTGATAAGCGTAGCGCATCAGGCaacta                    >  1:1435705/7‑96 (MQ=25)
                                                                                gctaaggcgTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTG‑C‑GA‑ATCC‑‑‑C‑TTGTAGGCCTGATAAGCGTAGCGCATCAGGCAtttgcg                  >  1:2848181/6‑97 (MQ=21)
                                                                                 aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CAGA‑AT‑T‑‑‑T‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                  <  1:116989/100‑1 (MQ=255)
                                                                                  tGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CAGA‑AT‑T‑‑‑T‑TCGTTGGCCTGATAAGCGGAGCGCATCAGGCATTtttgc                 >  1:3771128/1‑100 (MQ=255)
                                                                                       cgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CAGA‑AT‑T‑‑‑T‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg            <  1:2731698/100‑1 (MQ=255)
                                                                                       cgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CAGA‑AT‑T‑‑‑T‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg            <  2:3091136/100‑1 (MQ=255)
                                                                                        gTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CAGA‑AT‑T‑‑‑T‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgc           >  2:810154/1‑100 (MQ=255)
                                                                                         ttAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CAGA‑AT‑T‑‑‑T‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgcc          <  1:1620320/100‑1 (MQ=255)
                                                                                          tAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CAGA‑AT‑T‑‑‑T‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa         >  2:1805779/1‑100 (MQ=255)
                                                                                           aaCGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CAGA‑AT‑T‑‑‑T‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCAt        >  2:3610401/1‑100 (MQ=255)
                                                                                            aCGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CAGA‑AT‑T‑‑‑T‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc       <  2:1332259/100‑1 (MQ=255)
                                                                                              gCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CAGA‑AT‑T‑‑‑T‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAg     <  1:3713940/100‑1 (MQ=255)
                                                                                                cTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CAGA‑AT‑T‑‑‑T‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTc   >  1:131155/1‑100 (MQ=255)
                                                                                                  |                                                                                                          
GCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CAGA‑AT‑T‑‑‑T‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCT  >  NC_000913/3705887‑3706083

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